An official website of the United States government
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
A Dusty Locale: evolution of galactic dust populations from Milky Way to dwarf-mass galaxies
ABSTRACT Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments project, spanning $$M_{\rm vir}=10^{9-12}{M}_{\odot };\, M_{*}=10^{6-11}\, {M}_{\odot }$$. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species. All galactic dust populations in our suite exhibit similar evolutionary histories, with gas–dust accretion being the dominant producer of dust mass for all but the most metal-poor galaxies. Similar to previous works, we find the onset of efficient gas–dust accretion occurs above a ‘critical’ metallicity threshold (Zcrit). Due to this threshold, our simulations reproduce observed trends between galactic D/Z and metallicity and element depletion trends in the interstellar medium. However, we find Zcrit varies between dust species due to differences in key element abundances, dust physical properties, and life cycle processes resulting in $$Z_{\rm crit}\sim 0.05{\rm Z}_{\odot },\, 0.2{\rm Z}_{\odot },\, 0.5{\rm Z}_{\odot }$$ for metallic iron, silicates, and carbonaceous dust, respectively. These variations could explain the lack of small carbonaceous grains observed in the Magellanic Clouds. We also find a delay between the onset of gas–dust accretion and when a dust population reaches equilibrium, which we call the equilibrium time-scale (τequil). The relation between τequil and the metal enrichment time-scale of a galaxy, determined by its recent evolutionary history, can contribute to the scatter in the observed relation between galactic D/Z and metallicity.
more »
« less
- PAR ID:
- 10496062
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 529
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2356-2378
- Size(s):
- p. 2356-2378
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (∼1000 galaxies with $$M_{\rm star}\gt 10^8\, M_\odot$$ at z = 0) at a resolution ($$\Delta {}x\sim {}20\, {\rm pc}$$ , $$m_{\rm b}\sim {}6\times {}10^4\, M_\odot$$ ) comparable to state-of-the-art galaxy zoom-in simulations. FIREbox captures the multiphase nature of the interstellar medium in a fully cosmological setting (L = 22.1 Mpc) thanks to its exceptionally high dynamic range (≳106) and the inclusion of multichannel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that star formation rates, gas masses, and metallicities of simulated galaxies with $$M_{\rm star}\lt 10^{10.5-11}\, M_\odot$$ broadly agree with observations. These galaxy scaling relations extend to low masses ($$M_{\rm star}\sim {}10^7\, M_\odot$$ ) and follow a (broken) power-law relationship. Also reproduced are the evolution of the cosmic HI density and the HI column density distribution at z ∼ 0–5. At low z , FIREbox predicts a peak in the stellar-mass–halo-mass relation but also a higher abundance of massive galaxies and a higher cosmic star formation rate density than observed, showing that stellar feedback alone is insufficient to reproduce the properties of massive galaxies at late times. Given its high resolution and sample size, FIREbox offers a baseline prediction of galaxy formation theory in a ΛCDM Universe while also highlighting modelling challenges to be addressed in next-generation galaxy simulations.more » « less
-
ABSTRACT Recent observations indicate that galactic outflows are ubiquitous in high-redshift (high-z) galaxies, including normal star-forming galaxies, quasar hosts, and dusty star-forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star-formation histories and galactic outflow properties of galaxies in massive haloes ($$10^{12}\, {\rm M}_{\odot }\ \lt\ M_{\rm vir}\ \lt\ 5\times 10^{12}\, {\rm M}_{\odot }$$) at z ≳ 5.5 in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei. The simulated galaxies resemble z > 4 DSFGs, with star-formation rates of $$\sim\!{1000}\ {\rm M}_{\odot }\, \rm yr^{-1}$$ and molecular gas masses of Mmol ∼ 1010 M⊙. However, the simulated galaxies are characterized by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of ∼0.1, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium. We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of 30 per cent. Finally, compared with FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high-z sample. We argue that both effects originate from the higher total and gas surface densities that characterize high-z massive systems.more » « less
-
ABSTRACT We analyse the rest-optical emission-line ratios of z ∼ 1.5 galaxies drawn from the Multi-Object Spectrometer for Infra-Red Exploration Deep Evolution Field (MOSDEF) survey. Using composite spectra, we investigate the mass–metallicity relation (MZR) at z ∼ 1.5 and measure its evolution to z = 0. When using gas-phase metallicities based on the N2 line ratio, we find that the MZR evolution from z ∼ 1.5 to z = 0 depends on stellar mass, evolving by $$\Delta \rm log(\rm O/H) \sim 0.25$$ dex at M*< $$10^{9.75}\, \mathrm{M}_{\odot }$$ down to $$\Delta \rm log(\rm O/H) \sim 0.05$$ at M* ≳ $$10^{10.5}\, \mathrm{M}_{\odot }$$. In contrast, the O3N2-based MZR shows a constant offset of $$\Delta \rm log(\rm O/H) \sim 0.30$$ across all masses, consistent with previous MOSDEF results based on independent metallicity indicators, and suggesting that O3N2 provides a more robust metallicity calibration for our z ∼ 1.5 sample. We investigated the secondary dependence of the MZR on star formation rate (SFR) by measuring correlated scatter about the mean M*-specific SFR and M*−$$\log (\rm O3N2)$$ relations. We find an anticorrelation between $$\log (\rm O/H)$$ and sSFR offsets, indicating the presence of a M*−SFR−Z relation, though with limited significance. Additionally, we find that our z ∼ 1.5 stacks lie along the z = 0 metallicity sequence at fixed μ = log (M*/M⊙) − 0.6 × $$\log (\rm SFR / M_{\odot } \, yr^{-1})$$ suggesting that the z ∼ 1.5 stacks can be described by the z = 0 fundamental metallicity relation (FMR). However, using different calibrations can shift the calculated metallicities off of the local FMR, indicating that appropriate calibrations are essential for understanding metallicity evolution with redshift. Finally, understanding how [N ii]/H α scales with galaxy properties is crucial to accurately describe the effects of blended [N ii] and H α on redshift and H α fiux measurements in future large surveys utilizing low-resolution spectra such as with Euclid and the Roman Space Telescope.more » « less
-
ABSTRACT We introduce the thesan project, a suite of large volume ($$L_\mathrm{box} = 95.5 \, \mathrm{cMpc}$$) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the intergalactic medium during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter and baryonic mass resolutions of $$3.1 \times 10^6\, {\rm M_\odot }$$ and $$5.8 \times 10^5\, {\rm M_\odot }$$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z = 5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (arepo-rt) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying dark matter models, and numerical convergence. The fiducial simulation and model variations are calibrated to produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including the stellar-to-halo-mass relation, galaxy stellar mass function, star formation rate density, and the mass–metallicity relation, despite the galaxy formation model being mainly calibrated at z = 0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21 cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21 cm experiments to accurately characterize the sources that dominate the ionizing photon budget.more » « less
