scholarly journals EDGE: a new approach to suppressing numerical diffusion in adaptive mesh simulations of galaxy formation

2020 ◽  
Vol 501 (2) ◽  
pp. 1755-1765
Author(s):  
Andrew Pontzen ◽  
Martin P Rey ◽  
Corentin Cadiou ◽  
Oscar Agertz ◽  
Romain Teyssier ◽  
...  

ABSTRACT We introduce a new method to mitigate numerical diffusion in adaptive mesh refinement (AMR) simulations of cosmological galaxy formation, and study its impact on a simulated dwarf galaxy as part of the ‘EDGE’ project. The target galaxy has a maximum circular velocity of $21\, \mathrm{km}\, \mathrm{s}^{-1}$ but evolves in a region that is moving at up to $90\, \mathrm{km}\, \mathrm{s}^{-1}$ relative to the hydrodynamic grid. In the absence of any mitigation, diffusion softens the filaments feeding our galaxy. As a result, gas is unphysically held in the circumgalactic medium around the galaxy for $320\, \mathrm{Myr}$, delaying the onset of star formation until cooling and collapse eventually triggers an initial starburst at z = 9. Using genetic modification, we produce ‘velocity-zeroed’ initial conditions in which the grid-relative streaming is strongly suppressed; by design, the change does not significantly modify the large-scale structure or dark matter accretion history. The resulting simulation recovers a more physical, gradual onset of star formation starting at z = 17. While the final stellar masses are nearly consistent ($4.8 \times 10^6\, \mathrm{M}_{\odot }$ and $4.4\times 10^6\, \mathrm{M}_{\odot }$ for unmodified and velocity-zeroed, respectively), the dynamical and morphological structure of the z = 0 dwarf galaxies are markedly different due to the contrasting histories. Our approach to diffusion suppression is suitable for any AMR zoom cosmological galaxy formation simulations, and is especially recommended for those of small galaxies at high redshift.

2010 ◽  
Vol 6 (S277) ◽  
pp. 263-266
Author(s):  
Bruno Thooris ◽  
Daniel Pomarède

AbstractOur understanding of the structuring of the Universe from large-scale cosmological structures down to the formation of galaxies now largely benefits from numerical simulations. The RAMSES code, relying on the Adaptive Mesh Refinement technique, is used to perform massively parallel simulations at multiple scales. The interactive, immersive, three-dimensional visualization of such complex simulations is a challenge that is addressed using the SDvision software package. Several rendering techniques are available, including ray-casting and isosurface reconstruction, to explore the simulated volumes at various resolution levels and construct temporal sequences. These techniques are illustrated in the context of different classes of simulations. We first report on the visualization of the HORIZON Galaxy Formation Simulation at MareNostrum, a cosmological simulation with detailed physics at work in the galaxy formation process. We then carry on in the context of an intermediate zoom simulation leading to the formation of a Milky-Way like galaxy. Finally, we present a variety of simulations of interacting galaxies, including a case-study of the Antennae Galaxies interaction.


2004 ◽  
Vol 21 (4) ◽  
pp. 360-365 ◽  
Author(s):  
J. I. Davies ◽  
S. Sabatini ◽  
S. Roberts

AbstractLow luminosity (dwarf) galaxies play a crucial role in our current theories of galaxy and large scale structure formation. In the hierarchical picture they are the building blocks from which other structures form. These theories in their basic form overpredict the numbers of small dark matter halos (dwarf galaxies?) unless some form of star formation supression is invoked. In this paper we describe observations of dwarf galaxies in a range of different environments. We find that there are far too few dwarf galaxies in low density environments to be compatible with the theories. These observations are not consistent with an environment-independent mechanism suppressing dwarf galaxy formation. It is also not clear how these mechanisms can supress star formation if dwarf galaxies have large mass-to-light ratios (≈100). Either the whole idea of hierarchical galaxy formation has to be rejected or other environmentally dependent physical processes have to be invoked. We suggest that small, gas-rich dI galaxies have their evolution rapidly advanced as they move into the dense cluster environment.


2019 ◽  
Vol 15 (S341) ◽  
pp. 253-256 ◽  
Author(s):  
Li-Hsin Chen ◽  
Ke-Jung Chen

AbstractModern cosmological simulations suggest that the hierarchical assembly of dark matter halos provided the gravitational wells that allowed the primordial gases to form stars and galaxies inside them. The first galaxies comprised of the first systems of stars gravitationally bound in dark matter halos are naturally recognized as the building blocks of early Universe. To understand the formation of the first galaxies, we use an adaptive mesh refinement (AMR) cosmological code, Enzo to simulate the formation and evolution of the first galaxies. We first model an isolated galaxy by considering much microphysics such as star formation, stellar feedback, and primordial gas cooling. To examine the effect of Pop III stellar feedback to the first galaxy formation, we adjust the initial temperature, density distribution and metallicity distributions by assuming different IMFs of the first stars. Our results suggest that star formation in the first galaxies is sensitive to the initial conditions of Pop III supernovae and their remnants. Our study can help to correlate the populations of the first stars and supernovae to star formation inside these first galaxies which may be soon observed by the (James Webb Space Telescope JWST).


2008 ◽  
Vol 4 (S255) ◽  
pp. 111-115
Author(s):  
Britton D. Smith ◽  
Matthew J. Turk ◽  
Steinn Sigurdsson ◽  
Brian W. O'Shea ◽  
Michael L. Norman

AbstractIt is generally accepted that the very first stars in the universe were significantly more massive and formed much more in isolation than stars observed today. This suggests that there was a transition in star formation modes that was most likely related to the metallicity of the star-forming environment. We study how the addition of heavy elements alters the dynamics of collapsing gas by performing a series of numerical simulations of primordial star formation with various levels of pre-enrichment, using the adaptive mesh refinement, hydrodynamic + N-body code, Enzo. At high redshifts, the process of star formation is heavily influenced by the cosmic microwave background (CMB), which creates a temperature floor for the gas. Our results show that cloud-collapse can follow three distinct paths, depending on the metallicity. For very low metallicities (log10(Z/Z⊙) < −3.5), star formation proceeds in the primordial mode, producing only massive, singular objects. For high metallicities (log10(Z/Z⊙) > −3), efficient cooling from the metals cools the gas to the CMB temperature when the core density is still very low. When the gas temperature reaches the CMB temperature, the core becomes very thermally stable, and further fragmentation is heavily suppressed. In our simulations with log10(Z/Z⊙) > −3, only a single object forms with a mass-scale of a few hundred M⊙. We refer to this as the CMB-regulated star formation mode. For metallicities between these two limits (−3.5 < log10(Z/Z⊙) < −3), the gas cools efficiently, but never reaches the CMB temperature. In this mode, termed the metallicity-regulated star formation mode, the minimum gas temperature is reached at much higher densities, allowing the core to fragment and form multiple objects with mass-scales of only a few M⊙. Our results imply that the stellar initial mass function was top-heavy at very high redshift due to stars forming in the CMB-regulated mode. As the CMB temperature lowers with time, the metallicity-regulated star formation mode (producing multiple low-mass stars) operates at higher metallicities and eventually becomes the sole mode of star formation.


2007 ◽  
Vol 3 (S245) ◽  
pp. 471-476
Author(s):  
Aprajita Verma ◽  
Matthew Lehnert ◽  
Natascha Förster Schreiber ◽  
Malcolm Bremer ◽  
Laura Douglas

AbstractHigh redshift galaxies play a key role in our developing understanding of galaxy formation and evolution. Since such galaxies are being studied within a Gyr of the big bang, they provide a unique probe of the physics of one of the first generations of large-scale star-formation. We have performed a complete statistical study of the physical properties of a robust sample of z~5 UV luminous galaxies selected using the Lyman-break technique. The characteristic properties of this sample differ from LBGs at z~3 of comparable luminosity in that they are a factor of ten less massive (~few×109 M⊙) and the majority (~70%) are considerably younger (<100Myr). Our results support no more than a modest decline in the global star formation rate density at high redshifts and suggest that ~1% of the stellar mass density of the universe had already assembled at z~5. The constraint derived for the latter is affected by their young ages and short duty cycles which imply existing z~5 LBG samples may be highly incomplete. These intense starbursts have high unobscured star formation rate surface densities (~100s M⊙ yr−1 kpc−2), suggesting they drive outflows and winds that enrich the intra- and inter-galactic media with metals. These properties imply that the majority of z~5 LBGs are in formation meaning that most of their star-formation has likely occurred during the last few crossing times. They are experiencing their first (few) generations of large-scale star formation and are accumulating their first significant stellar mass. As such, z~5 LBGs are the likely progenitors of the spheroidal components of present-day massive galaxies (supported by their high stellar mass surface densities and their core phase-space densities).


2020 ◽  
Vol 638 ◽  
pp. A123 ◽  
Author(s):  
Gohar Dashyan ◽  
Yohan Dubois

The regulation of the baryonic content in dwarf galaxies is a long-standing problem. Supernovae (SNe) are supposed to play a key role in forming large-scale galactic winds by removing important amounts of gas from galaxies. SNe are efficient accelerators of non-thermal particles, so-called cosmic rays (CRs), which can substantially modify the dynamics of the gas and conditions to form large-scale galactic winds. We investigate how CR injection by SNe impacts the star formation and the formation of large-scale winds in dwarf galaxies, and whether it can produce galaxy star-formation rates (SFR) and wind properties closer to observations. We ran CR magneto-hydrodynamical simulations of dwarf galaxies at high resolution (9 pc) with the adaptive mesh refinement code RAMSES. Those disc galaxies are embedded in isolated halos of mass of 1010 and 1011 M⊙, and CRs are injected by SNe. We included CR isotropic and anisotropic diffusion with various diffusion coefficients, CR radiative losses, and CR streaming. The injection of CR energy into the interstellar medium smooths out the highest gas densities, which reduces the SFR by a factor of 2–3. Mass outflow rates are significantly greater with CR diffusion, by 2 orders of magnitudes for the higher diffusion coefficients. Without diffusion and streaming, CRs are inefficient at generating winds. CR streaming alone allows for the formation of winds but which are too weak to match observations. The formation of galactic winds strongly depends on the diffusion coefficient: for low coefficients, CR energy stays confined in high density regions where CR energy losses are highest, and higher coefficients, which allow for a more efficient leaking of CRs out of dense gas, produce stronger winds. CR diffusion leads to colder and denser winds than without CRs, and brings outflow rates and mass loading factors much closer to observations.


2020 ◽  
Vol 499 (2) ◽  
pp. 2648-2661
Author(s):  
Aaron A Dutton ◽  
Tobias Buck ◽  
Andrea V Macciò ◽  
Keri L Dixon ◽  
Marvin Blank ◽  
...  

ABSTRACT We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the response of cold dark matter (CDM) haloes to baryonic processes. Previous work has shown that the halo response is primarily a function of the ratio between galaxy stellar mass and total virial mass, and the density threshold above which gas is eligible to form stars, n[cm−3]. At low n all simulations in the literature agree that dwarf galaxy haloes are cuspy, but at high n ≳ 100 there is no consensus. We trace halo contraction in dwarf galaxies with n ≳ 100 reported in some previous simulations to insufficient spatial resolution. Provided the adopted star formation threshold is appropriate for the resolution of the simulation, we show that the halo response is remarkably stable for n ≳ 5, up to the highest star formation threshold that we test, n = 500. This free parameter can be calibrated using the observed clustering of young stars. Simulations with low thresholds n ≤ 1 predict clustering that is too weak, while simulations with high star formation thresholds n ≳ 5, are consistent with the observed clustering. Finally, we test the CDM predictions against the circular velocities of nearby dwarf galaxies. Low thresholds predict velocities that are too high, while simulations with n ∼ 10 provide a good match to the observations. We thus conclude that the CDM model provides a good description of the structure of galaxies on kpc scales provided the effects of baryons are properly captured.


Author(s):  
Marta B. Silva ◽  
Ely D. Kovetz ◽  
Garrett K. Keating ◽  
Azadeh Moradinezhad Dizgah ◽  
Matthieu Bethermin ◽  
...  

AbstractThis paper outlines the science case for line-intensity mapping with a space-borne instrument targeting the sub-millimeter (microwaves) to the far-infrared (FIR) wavelength range. Our goal is to observe and characterize the large-scale structure in the Universe from present times to the high redshift Epoch of Reionization. This is essential to constrain the cosmology of our Universe and form a better understanding of various mechanisms that drive galaxy formation and evolution. The proposed frequency range would make it possible to probe important metal cooling lines such as [CII] up to very high redshift as well as a large number of rotational lines of the CO molecule. These can be used to trace molecular gas and dust evolution and constrain the buildup in both the cosmic star formation rate density and the cosmic infrared background (CIB). Moreover, surveys at the highest frequencies will detect FIR lines which are used as diagnostics of galaxies and AGN. Tomography of these lines over a wide redshift range will enable invaluable measurements of the cosmic expansion history at epochs inaccessible to other methods, competitive constraints on the parameters of the standard model of cosmology, and numerous tests of dark matter, dark energy, modified gravity and inflation. To reach these goals, large-scale structure must be mapped over a wide range in frequency to trace its time evolution and the surveyed area needs to be very large to beat cosmic variance. Only a space-borne mission can properly meet these requirements.


2020 ◽  
Vol 500 (3) ◽  
pp. 3394-3412
Author(s):  
Steven R Furlanetto

ABSTRACT In recent years, simple models of galaxy formation have been shown to provide reasonably good matches to available data on high-redshift luminosity functions. However, these prescriptions are primarily phenomenological, with only crude connections to the physics of galaxy evolution. Here, we introduce a set of galaxy models that are based on a simple physical framework but incorporate more sophisticated models of feedback, star formation, and other processes. We apply these models to the high-redshift regime, showing that most of the generic predictions of the simplest models remain valid. In particular, the stellar mass–halo mass relation depends almost entirely on the physics of feedback (and is thus independent of the details of small-scale star formation) and the specific star formation rate is a simple multiple of the cosmological accretion rate. We also show that, in contrast, the galaxy’s gas mass is sensitive to the physics of star formation, although the inclusion of feedback-driven star formation laws significantly changes the naive expectations. While these models are far from detailed enough to describe every aspect of galaxy formation, they inform our understanding of galaxy formation by illustrating several generic aspects of that process, and they provide a physically grounded basis for extrapolating predictions to faint galaxies and high redshifts currently out of reach of observations. If observations show violations from these simple trends, they would indicate new physics occurring inside the earliest generations of galaxies.


2020 ◽  
Vol 498 (1) ◽  
pp. 164-180 ◽  
Author(s):  
Harley Katz ◽  
Dominika Ďurovčíková ◽  
Taysun Kimm ◽  
Joki Rosdahl ◽  
Jeremy Blaizot ◽  
...  

ABSTRACT Identifying low-redshift galaxies that emit Lyman continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly ‘analogues’ of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at z ∼ 3 and z ∼ 0. We find that the simulated galaxies with high LyC escape fractions (fesc) often exhibit high O32 and populate the same regions of the R23–O32 plane as z ∼ 3 LyC leakers. However, we show that viewing angle, metallicity, and ionization parameter can all impact where a galaxy resides on the O32–fesc plane. Based on emission line diagnostics and how they correlate with fesc, lower metallicity LyC leakers at z ∼ 3 appear to be good analogues of reionization-era galaxies. In contrast, local [S ii]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the S ii Baldwin–Phillips–Terlevich (BPT) diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using infrared and nebular emission lines. We show that our model using only [C ii]158 μm and [O iii]88 μm can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS 1149_JD1 at z = 9.1 is the most likely galaxy to be actively contributing to the reionization of the Universe.


Sign in / Sign up

Export Citation Format

Share Document