-
LightAMR format standard and lossless compression algorithms for adaptive mesh refinement grids: RAMSES use case
Authors:
Loïc Strafella,
Damien Chapon
Abstract:
The evolution of parallel I/O library as well as new concepts such as 'in transit' and 'in situ' visualization and analysis have been identified as key technologies to circumvent I/O bottleneck in pre-exascale applications. Nevertheless, data structure and data format can also be improved for both reducing I/O volume and improving data interoperability between data producer and data consumer. In t…
▽ More
The evolution of parallel I/O library as well as new concepts such as 'in transit' and 'in situ' visualization and analysis have been identified as key technologies to circumvent I/O bottleneck in pre-exascale applications. Nevertheless, data structure and data format can also be improved for both reducing I/O volume and improving data interoperability between data producer and data consumer. In this paper, we propose a very lightweight and purpose-specific post-processing data model for AMR meshes, called lightAMR. Based on this data model, we introduce a tree pruning algorithm that removes data redundancy from a fully threaded AMR octree. In addition, we present two lossless compression algorithms, one for the AMR grid structure description and one for AMR double/single precision physical quantity scalar fields. Then we present performance benchmarks on RAMSES simulation datasets of this new lightAMR data model and the pruning and compression algorithms. We show that our pruning algorithm can reduce the total number of cells from RAMSES AMR datasets by 10-40% without loss of information. Finally, we show that the RAMSES AMR grid structure can be compacted by ~ 3 orders of magnitude and the float scalar fields can be compressed by a factor ~ 1.2 for double precision and ~ 1.3 - 1.5 in single precision with a compression speed of ~ 1 GB/s.
△ Less
Submitted 25 August, 2022;
originally announced August 2022.
-
Formation of compact galaxies in the Extreme-Horizon simulation
Authors:
Solène Chabanier,
Frédéric Bournaud,
Yohan Dubois,
Sandrine Codis,
Damien Chapon,
David Elbaz,
Christophe Pichon,
Olivier Bressand,
Julien Devriendt,
Raphael Gavazzi,
Katarina Kraljic,
Taysun Kimm,
Clotilde Laigle,
Jean-Baptiste Lekien,
Garreth Martin,
Nathalie Palanque-Delabrouille,
Sébastien Peirani,
Pierre-Franck Piserchia,
Adrianne Slyz,
Maxime Trebitsch,
Christophe Yèche
Abstract:
We present the Extreme-Horizon (EH) cosmological simulation: EH models galaxy formation with stellar and AGN feedback and uses a very high resolution in the intergalactic and circumgalactic medium. The high resolution in low-density regions results in smaller-size massive galaxies at redshift $z=2$, in better agreement with observations compared to other simulations. This results from the improved…
▽ More
We present the Extreme-Horizon (EH) cosmological simulation: EH models galaxy formation with stellar and AGN feedback and uses a very high resolution in the intergalactic and circumgalactic medium. The high resolution in low-density regions results in smaller-size massive galaxies at redshift $z=2$, in better agreement with observations compared to other simulations. This results from the improved modeling of cold gas flows accreting onto galaxies. Besides, the EH simulation forms a population of particularly compact galaxies with stellar masses of $10^{10-11}$\,M$_\sun$ that are reminiscent of observed ultracompact galaxies at $z\simeq2$. These objects form mainly through repeated major mergers of low-mass progenitors, independently of baryonic feedback mechanisms. This formation process can be missed in simulations using a too low resolution in low-density intergalactic regions.
△ Less
Submitted 9 July, 2020;
originally announced July 2020.
-
Boosting I/O and visualization for exascale era using Hercule: test case on RAMSES
Authors:
Loic Strafella,
Damien Chapon
Abstract:
It has been clearly identified that I/O is one of the bottleneck to extend application for the exascale era. New concepts such as 'in transit' and 'in situ' visualization and analysis have been identified as key technologies to circumvent this particular issue. A new parallel I/O and data management library called Hercule, developed at CEA-DAM, has been integrated to Ramses, an AMR simulation code…
▽ More
It has been clearly identified that I/O is one of the bottleneck to extend application for the exascale era. New concepts such as 'in transit' and 'in situ' visualization and analysis have been identified as key technologies to circumvent this particular issue. A new parallel I/O and data management library called Hercule, developed at CEA-DAM, has been integrated to Ramses, an AMR simulation code for self-gravitating fluids. Splitting the original Ramses output format in Hercule database formats dedicated to either checkpoints/restarts (HProt format) or post-processing (HDep format) not only improved I/O performance and scalability of the Ramses code but also introduced much more flexibility in the simulation outputs to help astrophysicists prepare their DMP (Data Management Plan). Furthermore, the very lightweight and purpose-specific post-processing format (HDep) will significantly improve the overall performance of analysis and visualization tools such as PyMSES 5. An introduction to the Hercule parallel I/O library as well as I/O benchmark results will be discussed.
△ Less
Submitted 4 June, 2020;
originally announced June 2020.
-
A sub-parsec resolution simulation of the Milky Way: Global structure of the ISM and properties of molecular clouds
Authors:
F. Renaud,
F. Bournaud,
E. Emsellem,
B. Elmegreen,
R. Teyssier,
J. Alves,
D. Chapon,
F. Combes,
A. Dekel,
J. Gabor,
P. Hennebelle,
K. Kraljic
Abstract:
We present a self-consistent hydrodynamical simulation of a Milky Way-like galaxy, at the resolution of 0.05 pc. The model includes star formation and a new implementation of stellar feedback through photo-ionization, radiative pressure and supernovae. The simulation resolves the structure of the interstellar medium at subparsec resolution for a few cloud lifetimes, and at 0.05 pc for about a clou…
▽ More
We present a self-consistent hydrodynamical simulation of a Milky Way-like galaxy, at the resolution of 0.05 pc. The model includes star formation and a new implementation of stellar feedback through photo-ionization, radiative pressure and supernovae. The simulation resolves the structure of the interstellar medium at subparsec resolution for a few cloud lifetimes, and at 0.05 pc for about a cloud crossing time. Turbulence cascade and gravitation from the kpc scales are de facto included in smaller structures like molecular clouds. We show that the formation of a bar influences the dynamics of the central ~100 pc by creating resonances. At larger radii, the spiral arms host the formation of regularly spaced clouds: beads on a string and spurs. These instabilities pump turbulent energy into the gas, generally in the supersonic regime. Because of asymmetric drift, the supernovae explode outside of their gaseous nursery, which diminishes the effect of feedback on the structure of clouds. The evolution of clouds is thus mostly due to fragmentation and gas consumption, regulated mainly by supersonic turbulence. The transition from turbulence supported to self-gravitating gas is detected in the gas density probability distribution function at ~2000 cm^-3. The power spectrum density suggests that gravitation governs the hierarchical organisation of structures from the galactic scale down to a few parsecs.
△ Less
Submitted 6 September, 2013; v1 submitted 22 July, 2013;
originally announced July 2013.
-
Beyond the nuclear starburst? Clustered star formation in major mergers
Authors:
Leila C. Powell,
Frederic Bournaud,
Damien Chapon,
Romain Teyssier
Abstract:
Recent simulation work has successfully captured the formation of the star clusters that have been observed in merging galaxies. These studies, however, tend to focus on studying extreme starbursts, such as the Antennae galaxies. We aim to establish whether there is something special occurring in these extreme systems or whether the mechanism for cluster formation is present in all mergers to a gr…
▽ More
Recent simulation work has successfully captured the formation of the star clusters that have been observed in merging galaxies. These studies, however, tend to focus on studying extreme starbursts, such as the Antennae galaxies. We aim to establish whether there is something special occurring in these extreme systems or whether the mechanism for cluster formation is present in all mergers to a greater or lesser degree. We undertake a general study of merger-induced star formation in a sample of 5 pc resolution adaptive mesh refinement simulations of low redshift equal-mass mergers with randomly-chosen orbital parameters. We find that there is an enhanced mass fraction of very dense gas that appears as the gas density probability density function evolves during the merger. This finding has implications for the interpretation of some observations; a larger mass fraction of dense gas could account for the enhanced HCN/CO ratios seen in ULIRGs and predicts that alpha_CO is lower in mergers, as for a given mass of H_2, CO emission will increase in a denser environment. We also find that as the star formation rate increases, there is a correlated peak in the velocity dispersion of the gas, which we attribute to increasing turbulence driven by the interaction itself. Star formation tends to be clumpy: in some cases there is extended clumpy star formation, but even when star formation is concentrated within the inner kpc (i.e. what may be considered a nuclear starburst) it still often has a clumpy, rather than a smooth, distribution. We find no strong evidence for a clear bimodality in the Kennicutt-Schmidt relation for the average mergers simulated here. Instead, they are typically somewhat offset above the predicted quiescent relation during their starbursts.
△ Less
Submitted 11 June, 2013;
originally announced June 2013.
-
Hydrodynamics of galaxy mergers with supermassive black holes: is there a last parsec problem ?
Authors:
Damien Chapon,
Lucio Mayer,
Romain Teyssier
Abstract:
We study the formation of a supermassive black hole (SMBH) binary and the shrinking of the separation of the two holes to sub-pc scales starting from a realistic major merger between two gas-rich spiral galaxies with mass comparable to our Milky Way. The simulations, carried out with the Adaptive Mesh Refinement (AMR) code RAMSES, are capable of resolving separations as small as 0.1 pc. The collis…
▽ More
We study the formation of a supermassive black hole (SMBH) binary and the shrinking of the separation of the two holes to sub-pc scales starting from a realistic major merger between two gas-rich spiral galaxies with mass comparable to our Milky Way. The simulations, carried out with the Adaptive Mesh Refinement (AMR) code RAMSES, are capable of resolving separations as small as 0.1 pc. The collision of the two galaxies produces a gravo-turbulent rotating nuclear disk with mass (10^9 Msun) and size (60 pc) in excellent agreement with previous SPH simulations with particle splitting that used a similar setup (Mayer et al. 2007) but were limited to separations of a few parsecs. The AMR results confirm that the two black holes sink rapidly as a result of dynamical friction onto the gaseous background, reaching a separation of 1 pc in less than 10^7 yr. We show that the dynamical friction wake is well resolved by our model and we find good agreement with analytical predictions of the drag force as a function of the Mach number. Below 1 pc, black hole pairing slows down significantly, as the relative velocity between the sinking SMBH becomes highly subsonic and the mass contained within their orbit falls below the mass of the binary itself, rendering dynamical friction ineffective. In this final stage, the black holes have not opened a gap as the gaseous background is highly pressurized in the center. Non-axisymmetric gas torques do not arise to restart sinking in absence of efficient dynamical friction, at variance with previous calculations using idealized equilibrium nuclear disk models. (abridged)
△ Less
Submitted 24 December, 2012; v1 submitted 27 October, 2011;
originally announced October 2011.
-
The impact of ISM turbulence, clustered star formation and feedback on galaxy mass assembly through cold flows and mergers
Authors:
Leila C. Powell,
Frederic Bournaud,
Damien Chapon,
Julien Devriendt,
Adrianne Slyz,
Romain Teyssier
Abstract:
Two of the dominant channels for galaxy mass assembly are cold flows (cold gas supplied via the filaments of the cosmic web) and mergers. How these processes combine in a cosmological setting, at both low and high redshift, to produce the whole zoo of galaxies we observe is largely unknown. Indeed there is still much to understand about the detailed physics of each process in isolation. While thes…
▽ More
Two of the dominant channels for galaxy mass assembly are cold flows (cold gas supplied via the filaments of the cosmic web) and mergers. How these processes combine in a cosmological setting, at both low and high redshift, to produce the whole zoo of galaxies we observe is largely unknown. Indeed there is still much to understand about the detailed physics of each process in isolation. While these formation channels have been studied using hydrodynamical simulations, here we study their impact on gas properties and star formation (SF) with some of the first simulations that capture the multiphase, cloudy nature of the interstellar medium (ISM), by virtue of their high spatial resolution (and corresponding low temperature threshold). In this regime, we examine the competition between cold flows and a supernovae(SNe)-driven outflow in a very high-redshift galaxy (z {\approx} 9) and study the evolution of equal-mass galaxy mergers at low and high redshift, focusing on the induced SF. We find that SNe-driven outflows cannot reduce the cold accretion at z {\approx} 9 and that SF is actually enhanced due to the ensuing metal enrichment. We demonstrate how several recent observational results on galaxy populations (e.g. enhanced HCN/CO ratios in ULIRGs, a separate Kennicutt Schmidt (KS) sequence for starbursts and the population of compact early type galaxies (ETGs) at high redshift) can be explained with mechanisms captured in galaxy merger simulations, provided that the multiphase nature of the ISM is resolved.
△ Less
Submitted 21 February, 2011;
originally announced February 2011.
-
Star formation in galaxy mergers: ISM turbulence, dense gas excess, and scaling relations for disks and starbusts
Authors:
Frederic Bournaud,
Leila C. Powell,
Damien Chapon,
Romain Teyssier
Abstract:
Galaxy interactions and mergers play a significant, but still debated and poorly understood role in the star formation history of galaxies. Numerical and theoretical models cannot yet explain the main properties of merger-induced starbursts, including their intensity and their spatial extent. Usually, the mechanism invoked in merger-induced starbursts is a global inflow of gas towards the central…
▽ More
Galaxy interactions and mergers play a significant, but still debated and poorly understood role in the star formation history of galaxies. Numerical and theoretical models cannot yet explain the main properties of merger-induced starbursts, including their intensity and their spatial extent. Usually, the mechanism invoked in merger-induced starbursts is a global inflow of gas towards the central kpc, resulting in a nuclear starburst. We show here, using high-resolution AMR simulations and comparing to observations of the gas component in mergers, that the triggering of starbursts also results from increased ISM turbulence and velocity dispersions in interacting systems. This forms cold gas that are denser and more massive than in quiescent disk galaxies. The fraction of dense cold gas largely increases, modifying the global density distribution of these systems, and efficient star formation results. Because the starbursting activity is not just from a global compacting of the gas to higher average surface densities, but also from higher turbulence and fragmentation into massive and dense clouds, merging systems can enter a different regime of star formation compared to quiescent disk galaxies. This is in quantitative agreement with recent observations suggesting that disk galaxies and starbursting systems are not the low-activity end and high-activity end of a single regime, but actually follow different scaling relations for their star formation.
△ Less
Submitted 23 December, 2010;
originally announced December 2010.
-
Galactic star formation in parsec-scale resolution simulations
Authors:
Leila C. Powell,
Frederic Bournaud,
Damien Chapon,
Julien Devriendt,
Adrianne Slyz,
Romain Teyssier
Abstract:
The interstellar medium (ISM) in galaxies is multiphase and cloudy, with stars forming in the very dense, cold gas found in Giant Molecular Clouds (GMCs). Simulating the evolution of an entire galaxy, however, is a computational problem which covers many orders of magnitude, so many simulations cannot reach densities high enough or temperatures low enough to resolve this multiphase nature. Therefo…
▽ More
The interstellar medium (ISM) in galaxies is multiphase and cloudy, with stars forming in the very dense, cold gas found in Giant Molecular Clouds (GMCs). Simulating the evolution of an entire galaxy, however, is a computational problem which covers many orders of magnitude, so many simulations cannot reach densities high enough or temperatures low enough to resolve this multiphase nature. Therefore, the formation of GMCs is not captured and the resulting gas distribution is smooth, contrary to observations. We investigate how star formation (SF) proceeds in simulated galaxies when we obtain parsec-scale resolution and more successfully capture the multiphase ISM. Both major mergers and the accretion of cold gas via filaments are dominant contributors to a galaxy's total stellar budget and we examine SF at high resolution in both of these contexts.
△ Less
Submitted 24 September, 2010;
originally announced September 2010.
-
Hydrodynamics of high-redshift galaxy collisions: From gas-rich disks to dispersion-dominated mergers and compact spheroids
Authors:
Frederic Bournaud,
Damien Chapon,
Romain Teyssier,
Leila C. Powell,
Bruce G. Elmegreen,
Debra Meloy Elmegreen,
Pierre-Alain Duc,
Thierry Contini,
Benoit Epinat,
Kristen L. Shapiro
Abstract:
Disk galaxies at high redshift (z~2) are characterized by high fractions of cold gas, strong turbulence, and giant star-forming clumps. Major mergers of disk galaxies at high redshift should then generally involve such turbulent clumpy disks. Merger simulations, however, model the ISM as a stable, homogeneous, and thermally pressurized medium. We present the first merger simulations with high frac…
▽ More
Disk galaxies at high redshift (z~2) are characterized by high fractions of cold gas, strong turbulence, and giant star-forming clumps. Major mergers of disk galaxies at high redshift should then generally involve such turbulent clumpy disks. Merger simulations, however, model the ISM as a stable, homogeneous, and thermally pressurized medium. We present the first merger simulations with high fractions of cold, turbulent, and clumpy gas. We discuss the major new features of these models compared to models where the gas is artificially stabilized and warmed. Gas turbulence, which is already strong in high-redshift disks, is further enhanced in mergers. Some phases are dispersion-dominated, with most of the gas kinetic energy in the form of velocity dispersion and very chaotic velocity fields, unlike merger models using a thermally stabilized gas. These mergers can reach very high star formation rates, and have multi-component gas spectra consistent with SubMillimeter Galaxies. Major mergers with high fractions of cold turbulent gas are also characterized by highly dissipative gas collapse to the center of mass, with the stellar component following in a global contraction. The final galaxies are early-type with relatively small radii and high Sersic indices, like high-redshift compact spheroids. The mass fraction in a disk component that survives or re-forms after a merger is severely reduced compared to models with stabilized gas, and the formation of a massive disk component would require significant accretion of external baryons afterwards. Mergers thus appear to destroy extended disks even when the gas fraction is high, and this lends further support to smooth infall as the main formation mechanism for massive disk galaxies.
△ Less
Submitted 13 January, 2011; v1 submitted 24 June, 2010;
originally announced June 2010.
-
The driving mechanism of starbursts in galaxy mergers
Authors:
Romain Teyssier,
Damien Chapon,
Frederic Bournaud
Abstract:
We present hydrodynamic simulations of a major merger of disk galaxies, and study the ISM dynamics and star formation properties. High spatial and mass resolutions of 12pc and 4x10^4 M_sol allow to resolve cold and turbulent gas clouds embedded in a warmer diffuse phase. We compare to lower resolution models, where the multiphase ISM is not resolved and is modeled as a relatively homogeneous and s…
▽ More
We present hydrodynamic simulations of a major merger of disk galaxies, and study the ISM dynamics and star formation properties. High spatial and mass resolutions of 12pc and 4x10^4 M_sol allow to resolve cold and turbulent gas clouds embedded in a warmer diffuse phase. We compare to lower resolution models, where the multiphase ISM is not resolved and is modeled as a relatively homogeneous and stable medium. While merger-driven bursts of star formation are generally attributed to large-scale gas inflows towards the nuclear regions, we show that once a realistic ISM is resolved, the dominant process is actually gas fragmentation into massive and dense clouds and rapid star formation therein. As a consequence, star formation is more efficient by a factor of up to 10 and is also somewhat more extended, while the gas density probability distribution function (PDF) rapidly evolves towards very high densities. We thus propose that the actual mechanism of starburst triggering in galaxy collisions can only be captured at high spatial resolution and when the cooling of gas is modeled down to less than 10^3 K. Not only does our model reproduce the properties of the Antennae system, but it also explains the ``starburst mode'' revealed recently in high-redshift mergers compared to quiescent disks.
△ Less
Submitted 24 June, 2010;
originally announced June 2010.
