ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($$M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$$, $$M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($$\gtrsim 1\,$$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($$M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$$, $$M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($$M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $$\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $$\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant.
more »
« less
Be it therefore resolved: cosmological simulations of dwarf galaxies with 30 solar mass resolution
ABSTRACT We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies ($$M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$$), run to z = 0 with $$30\, \mathrm{M}_{\odot }$$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $$M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$$ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $$M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$$) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs ($$M_{\ast }\gt 10^{5}\, \mathrm{M}_{\odot }$$) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M*/Mhalo > 10−4 to form a dark matter core $${\gt}200\rm \, pc$$, while lower mass UFDs exhibit cusps down to $${\lesssim}100\rm \, pc$$, as expected from energetic arguments. Our dwarfs with $$M_{\ast }\gt 10^{4}\, \mathrm{M}_{\odot }$$ have half-mass radii (R1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).
more »
« less
- PAR ID:
- 10123453
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 490
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 4447-4463
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT Carbon enhanced metal poor (CEMP)-no stars, a subset of CEMP stars ($$\rm [C/Fe]\ge 0.7$$ and $$\rm [Fe/H]\lesssim -1$$) have been discovered in ultra-faint dwarf (UFD) galaxies, with $$M_{\rm vir}\approx 10^8{\, \mathrm{ M}_\odot }$$ and $$M_{\ast }\approx 10^3-10^4{\, \mathrm{ M}_\odot }$$ at z = 0, as well as in the halo of the Milky Way (MW). These CEMP-no stars are local fossils that may reflect the properties of the first (Pop III) and second (Pop II) generation of stars. However, cosmological simulations have struggled to reproduce the observed level of carbon enhancement of the known CEMP-no stars. Here, we present new cosmological hydrodynamic zoom-in simulations of isolated UFDs that achieve a gas mass resolution of $$m_{\rm gas}\approx 60{\, \mathrm{ M}_\odot }$$. We include enrichment from Pop III faint supernovae (SNe), with ESN = 0.6 × 1051 erg, to understand the origin of CEMP-no stars. We confirm that Pop III and Pop II stars are mainly responsible for the formation of CEMP and C-normal stars, respectively. New to this study, we find that a majority of CEMP-no stars in the observed UFDs and the MW halo can be explained by Pop III SNe with normal explosion energy (ESN = 1.2 × 1051 erg) and Pop II enrichment, but faint SNe might also be needed to produce CEMP-no stars with $$\rm [C/Fe]\gtrsim 2$$, corresponding to the absolute carbon abundance of $$\rm A(C)\gtrsim 6.0$$. Furthermore, we find that while we create CEMP-no stars with high carbon ratio $$\rm [C/Fe]\approx 3-4$$, by adopting faint SNe, it is still challenging to reproduce CEMP-no stars with extreme level of carbon abundance of $$\rm A(C)\approx 7.0-7.5$$, observed both in the MW halo and UFDs.more » « less
-
ABSTRACT We study star formation histories (SFHs) of 500 dwarf galaxies (stellar mass $$M_\ast =10^5\!-\!10^9\, \rm {M}_\odot$$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of ‘satellite versus central’ and ‘LG versus individual MW versus isolated dwarf central.’ Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the ‘near field’ have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated (true field) dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for M* = 106–$$10^7\, \rm {M}_\odot$$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG.more » « less
-
ABSTRACT The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in z = 0 galaxies at $$M_{*}=10^{5\!-\!10}\, \rm {M}_{\odot }$$ around 14 MW-mass hosts from the Feedback In Realistic Environments 2 (FIRE-2) simulations. Lower mass satellites ($$M_{*}\lesssim 10^7\, \rm {M}_{\odot }$$) are mostly quiescent and higher mass satellites ($$M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$$) are mostly star forming, with intermediate-mass satellites ($$M_{*}\approx 10^{7\!-\!8}\, \rm {M}_{\odot }$$) split roughly equally between quiescent and star forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($$M_{*}=10^{8\!-\!9}\, \rm {M}_{\odot }$$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and Large Magellanic Cloud (LMC)-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as pre-processing in low-mass groups prior to MW infall. Lower mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced ≥1–2 pericentre passages (≈2.5–5 Gyr) within a MW-mass halo. Most galaxies with $$M_{*}\gtrsim 10^{6.5}\, \rm {M}_{\odot }$$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the Satellites Around Galactic Analogs (SAGA) survey’s quiescent fraction at $$M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$$, they do not generally reproduce SAGA’s turnover at lower masses.more » « less
-
ABSTRACT We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $$(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$$ at z = 0) as part of the ‘Feedback in Realistic Environment’ (FIRE) project. Compact (r1/2 < 30 pc), relatively massive (0.5 × 105 ≲ M⋆/M⊙ ≲ 5 × 105), self-bound stellar clusters form at 11 ≳ z ≳ 5 in progenitors with $$M_{\rm vir} \approx 10^9\, \mathrm{M}_{\odot }$$. Cluster formation is triggered when at least $$10^7\, \mathrm{M}_{\odot }$$ of dense, turbulent gas reaches $$\Sigma _{\rm gas} \approx 10^4\, {\rm M}_\odot \, {\rm pc}^{-2}$$ as a result of the compressive effects of supernova feedback or from cloud–cloud collisions. The clusters can survive for $$2-3\, {\rm Gyr}$$; absent numerical effects, they could possibly survive substantially longer, perhaps to z = 0. The longest lived clusters are those that form at significant distance – several hundreds of pc – from their host galaxy. We therefore predict that globular clusters forming in progenitors of present-day dwarf galaxies will be offset from any pre-existing stars within their host dark matter haloes as opposed to deeply embedded within a well-defined galaxy. Properties of the nascent clusters are consistent with observations of some of the faintest and most compact high-redshift sources in Hubble Space Telescope lensing fields and are at the edge of what will be detectable as point sources in deep imaging of non-lensed fields with JWST. By contrast, the star clusters’ host galaxies will remain undetectable.more » « less