Search for: All records

Atomic hydrogen (H i) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of H i near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts $z \lesssim 3$, the low-mass end remains as-of-yet unexplored in observational and computational surveys. We employ a cosmological, hydrodynamical simulation (FIREbox) supplemented with zoom-in simulations (MassiveFIRE) from the Feedback In Realistic Environments (FIRE) project to investigate the H i covering fraction of Lyman Limit Systems ($N_{{\text{H}}\, \rm{{\small I}}} \gtrsim 10^{17.2}$ cm$^{-2}$) across a wide range of redshifts ($z=0-6$) and halo masses ($10^8-10^{13} \, \,\mathrm{ M}_{\odot }$ at $z=0$, $10^8-10^{11}\, \,\mathrm{ M}_{\odot }$ at $z=6$) in the absence of feedback from active galactic nuclei. We find that the covering fraction inside haloes exhibits a strong increase with redshift, with only a weak dependence on halo mass for higher mass haloes. For massive haloes ($M_{\mathrm{vir}} \sim 10^{11}-10^{12} \,\mathrm{ M}_{\odot }$), the radial profiles showcase scale-invariance and remain independent of mass. The radial dependence is well captured by a fitting function. The covering fractions in our simulations are in good agreement with measurements of the covering fraction in LBGs. Our comprehensive analysis unveils a complex dependence with redshift and halo mass for haloes with $M_{\mathrm{vir}} \lesssim 10^{10} \,\mathrm{ M}_{\odot }$ that future observations aim to constrain, providing key insights into the physics of structure formation and gas assembly.

 

Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$_2$ gas mass of the region using the CO-to-H$_2$ conversion factor ($X_{\rm{{\small CO}}}$). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the fire-2 subgrid models coupled to the chimes non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H$_2$ column densities, and the relation between $X_{\rm{{\small CO}}}$ and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H$_2$ mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5 per cent) for our high-mass galaxies, though still affects the H$_2$ mass and $L_{\rm{{\small CO}}}$ by $\approx$25  per cent.

 

The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium metallicity measurements and enabled characterization of the gas-phase mass–metallicity relation (MZR) for galaxies atz≳ 5 for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations atz= 5–12 and for stellar massesM_*∼ 10⁶–10¹⁰M_⊙. These simulations implement a multichannel stellar feedback model and produce broadly realistic galaxy properties, including when evolved toz= 0. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than 0.01 dex as redshift decreases fromz= 12 toz= 5. The median MZR in the simulations is well approximated as a constant power-law relation across this redshift range given by$\log (Z/Z_{\odot })=0.37\log (M_{*}/{\mathrm{M}}_{\odot })-4.3$. We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR atz> 5 contrasts with the evolution atz≲ 3, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations.

 

Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyse how the MZR evolves over z = 0–3, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas-phase metallicity, and star formation rate.

 

Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the H i properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and Firebox). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the H i content of galaxies. We find that the galaxy-wide H i properties agree with each other and with observations. However, differences appear for small-scale properties. The thin H i discs observed in the local universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of H i discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the H i discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies.

 

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.

 

An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21 cm hyperfine transition of diffuse neutral hydrogen. Under the basic assumption that the 21 cm signal is activated by the Lyαphoton field produced by metal-poor stellar systems, we show that a detection at the low frequencies of the EDGES and SARAS3 experiments may be expected from a simple extrapolation of the declining UV luminosity density inferred atz≲ 14 from JWST early galaxy data. Accounting for an early radiation excess above the cosmic microwave background suggests a shallower or flat evolution to simultaneously reproduce low- and high-zcurrent UV luminosity density constraints, which cannot be entirely ruled out, given the large uncertainties from cosmic variance and the faint-end slope of the galaxy luminosity function at cosmic dawn. Our findings raise the intriguing possibility that a high star formation efficiency at early times may trigger the onset of intense Lyαemission at redshiftz≲ 20 and produce a cosmic 21 cm absorption signal 200 Myr after the Big Bang.

 

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. 

Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C\, \small {II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C\, \small {II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C\, \small {II}]}{}/{}\rm SFR$ than local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origins of this ‘$\rm [C\, \small {II}]$ deficit’ is unclear. In this work, we study the $L_{\rm [C\, \small {II}]}$-SFR relation of galaxies using a sample of z = 0 − 8 galaxies with M* ≈ 107 − 5 × 1011 M⊙ extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (fire) project. We find a simple analytic expression for $L_{\rm [C\, \small {II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C\, \small {II}]$-emitting gas ($f_{\rm [C\, \small {II}]}$), gas metallicity (Zgas), gas density (ngas) and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes: $\rm H_2$-rich galaxies where tdep is the main driver of the $\rm [C\, \small {II}]$ deficit and $\rm H_2$-poor galaxies where Zgas is the main driver. The observed $\rm [C\, \small {II}]$ deficit of IR-luminous galaxies and early EoR galaxies, corresponding to the two different regimes, is due to short gas depletion time and low gas metallicity, respectively. Our result indicates that the $\rm [C\, \small {II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C\, \small {II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C\, \small {II}]$ line intensity mapping experiments.

 

Accurately reproducing the thin cold gas discs observed in nearby spiral galaxies has been a long standing issue in cosmological simulations. Here, we present measurements of the radially resolved H i scale height in 22 non-interacting Milky Way-mass galaxies from the FIREbox cosmological volume. We measure the H i scale heights using five different approaches commonly used in the literature: fitting the vertical volume density distribution with a Gaussian, the distance between maximum and half-maximum of the vertical volume density distribution, a semi-empirical description using the velocity dispersion and the galactic gravitational potential, the analytic assumption of hydrostatic equilibrium, and the distance from the midplane which encloses ≳60 per cent of the H i mass. We find median H i scale heights, measured using the vertical volume distribution, that range from ∼100 pc in the galactic centres to ∼800 pc in the outskirts and are in excellent agreement with recent observational results. We speculate that the presence of a realistic multiphase interstellar medium, including cold gas, and realistic stellar feedback are the drivers behind the realistic H i scale heights.