ABSTRACT We use a particle tracking analysis to study the origins of the circumgalactic medium (CGM), separating it into (1) accretion from the intergalactic medium (IGM), (2) wind from the central galaxy, and (3) gas ejected from other galaxies. Our sample consists of 21 FIRE-2 simulations, spanning the halo mass range Mh ∼ 1010–1012 M⊙, and we focus on z = 0.25 and z = 2. Owing to strong stellar feedback, only ∼L⋆ haloes retain a baryon mass $\gtrsim\! 50\hbox{ per cent}$ of their cosmic budget. Metals are more efficiently retained by haloes, with a retention fraction $\gtrsim\! 50\hbox{ per cent}$. Across all masses and redshifts analysed $\gtrsim \!60\hbox{ per cent}$ of the CGM mass originates as IGM accretion (some of which is associated with infalling haloes). Overall, the second most important contribution is wind from the central galaxy, though gas ejected or stripped from satellites can contribute a comparable mass in ∼L⋆ haloes. Gas can persist in the CGM for billions of years, resulting in well mixed-halo gas. Sightlines through the CGM are therefore likely to intersect gas of multiple origins. For low-redshift ∼L⋆ haloes, cool gas (T < 104.7 K) is distributed on average preferentially along the galaxy plane, however with strong halo-to-halo variability. The metallicity of IGM accretion is systematically lower than the metallicity of winds (typically by ≳1 dex), although CGM and IGM metallicities depend significantly on the treatment of subgrid metal diffusion. Our results highlight the multiple physical mechanisms that contribute to the CGM and will inform observational efforts to develop a cohesive picture.
more »
« less
The fates of the circumgalactic medium in the FIRE simulations
ABSTRACT We analyse the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z = 0.25 and 2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range $M_{\rm h}(z=0) \sim 10^{10} \!-\! 10^{12} \rm {\,M}_\odot$. We analyse both where the gas ends up after first leaving the CGM (its ‘proximate’ fate) and its location at z = 0 (its ‘ultimate’ fate). Of the CGM at z = 2, about half is found in the ISM or stars of the central galaxy by z = 0 in Mh(z = 2) ∼ 5 × 1011 M⊙ haloes, but most of the CGM in lower mass haloes ends up in the intergalactic medium (IGM). This is so even though most of the CGM in Mh(z = 2) ∼ 5 × 1010 M⊙ haloes first accretes on to the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z = 0.25 remains in the CGM by z = 0 at all halo masses analysed. Of the CGM gas that subsequently accretes on to the central galaxy in the progenitors of Mh(z = 0) ∼ 1012 M⊙ haloes, most of it is cool (T ∼ 104 K) at z = 2 but hot (∼Tvir) at z ∼ 0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z = 0.25 and $2 \, {\gtrsim} 80{{\ \rm per\ cent}}$ of the cool gas in $M_{\rm h} \gtrsim 10^{11} \rm {M}_\odot$ haloes will accrete on to a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM.
more »
« less
- NSF-PAR ID:
- 10184318
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 494
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 3581 to 3595
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $\sim 10^{12}\, {\rm M_{\odot }}$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of ${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $\gtrsim {\dot{M}}_{\rm crit}$ the volume-filling gas phase can in principle be ‘transonic’ – virialized in the outer halo but cool and free-falling near the galaxy. We compare ${\dot{M}}_{\rm crit}$ to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that ${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$ at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of $\sim 10^{12}\, {\rm M_{\odot }}$.more » « less
-
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
-
ABSTRACT Recent searches for the hosts of z ∼ 4 damped Ly α absorbers (DLAs) have detected bright galaxies at distances of tens of kpc from the DLA. Using the FIRE-2 cosmological zoom simulations, we argue that these relatively large distances are due to a predominantly cool and neutral inner circumgalactic medium (CGM) surrounding high-redshift galaxies. The inner CGM is cool because of the short cooling time of hot gas in ${\lesssim}10^{12}\, {\rm M_{\odot }}$ haloes, which implies that accretion and feedback energy are radiated quickly, while it is neutral due to high volume densities and column densities at high redshift that shield cool gas from photoionization. Our analysis predicts large DLA covering factors (${\gtrsim}50{{\ \rm per\ cent}}$) out to impact parameters ∼0.3[(1 + z)/5]3/2Rvir from the central galaxies at z ≳ 1, equivalent to a proper distance of ${\sim}21\, M_{12}^{1/3} \left(\left(1+z\right)/5\right)^{1/2}\, {\rm kpc}$ (Rvir and M12 are the halo virial radius and mass in units of $10^{12}\, {\rm M_{\odot }}$, respectively). This implies that DLA covering factors at z ∼ 4 may be comparable to unity out to a distance ∼10 times larger than stellar half-mass radii. A predominantly neutral inner CGM in the early universe suggests that its mass and metallicity can be directly constrained by absorption surveys, without resorting to the large ionization corrections as required for ionized CGM.more » « less
-
null (Ed.)ABSTRACT This paper presents a survey of Mg ii absorbing gas in the vicinity of 380 random galaxies, using 156 background quasi-stellar objects (QSOs) as absorption-line probes. The sample comprises 211 isolated (73 quiescent and 138 star-forming galaxies) and 43 non-isolated galaxies with sensitive constraints for both Mg ii absorption and H α emission. The projected distances span a range from d = 9 to 497 kpc, redshifts of the galaxies range from z = 0.10 to 0.48, and rest-frame absolute B-band magnitudes range from MB = −16.7 to −22.8. Our analysis shows that the rest-frame equivalent width of Mg ii, Wr(2796), depends on halo radius (Rh), B-band luminosity(LB), and stellar mass (Mstar) of the host galaxies, and declines steeply with increasing d for isolated, star-forming galaxies. At the same time, Wr(2796) exhibits no clear trend for either isolated, quiescent galaxies or non-isolated galaxies. In addition, the covering fraction of Mg ii absorbing gas 〈κ〉 is high with 〈κ〉 ≳ 60 per cent at <40 kpc for isolated galaxies and declines rapidly to 〈κ〉 ≈ 0 at d ≳ 100 kpc. Within the gaseous radius, the incidence of Mg ii gas depends sensitively on both Mstar and the specific star formation rate inferred from H α. Different from what is known for massive quiescent haloes, the observed velocity dispersion of Mg ii absorbing gas around star-forming galaxies is consistent with expectations from virial motion, which constrains individual clump mass to $m_{\rm cl} \gtrsim 10^5 \, \rm M_\odot$ and cool gas accretion rate of $\sim 0.7\!-\!2 \, \mathrm{ M}_\odot \, \rm yr^{-1}$. Finally, we find no strong azimuthal dependence of Mg ii absorption for either star-forming or quiescent galaxies. Our results demonstrate that multiple parameters affect the properties of gaseous haloes around galaxies and highlight the need of a homogeneous, absorption-blind sample for establishing a holistic description of chemically enriched gas in the circumgalactic space.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/staa902
