skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, December 13 until 2:00 AM ET on Saturday, December 14 due to maintenance. We apologize for the inconvenience.


Title: FIREbox: simulating galaxies at high dynamic range in a cosmological volume
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
Award ID(s):
2108230 1652522 1715216 2108314 2108318 1713353 2107872
PAR ID:
10428376
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
522
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3831 to 3860
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We investigate galaxy sizes at redshift $z\gtrsim 6$ with the cosmological radiation-magnetohydrodynamic simulation suite thesan(-hr). These simulations simultaneously capture reionization of the large-scale intergalactic medium and resolved galaxy properties. The intrinsic sizes ($r^{\ast }_{1/2}$) of simulated galaxies increase moderately with stellar mass at $M_{\ast } \lesssim 10^{8}{\, \rm M_\odot}$ and decrease fast at larger masses, resulting in a hump feature at $M_{\ast }\sim 10^{8}{\, \rm M_\odot}$ that is insensitive to redshift. Low-mass galaxies are in the initial phase of size growth and are better described by a spherical shell model with feedback-driven outflows competing with the cold inflowing gas streams. In contrast, massive galaxies fit better with the disc formation model. They generally experience a phase of rapid compaction and gas depletion, likely driven by internal disc instability rather than external processes. We identify four compact quenched galaxies in the $(95.5\, {\rm cMpc})^{3}$ volume of thesan-1 at $z\simeq 6$ and their quenching follows reaching a characteristic stellar surface density akin to the massive compact galaxies at cosmic noon. Compared to observations, we find that the median ultraviolet effective radius ($R^{\rm UV}_{\rm eff}$) of simulated galaxies is at least three times larger than the observed ones at $M_{\ast }\lesssim 10^{9}{\, \rm M_\odot}$ or $M_{\rm UV}\gtrsim -20$ at $6 \lesssim z \lesssim 10$. The population of compact galaxies ($R^{\rm UV}_{\rm eff}\lesssim 300\, {\rm pc}$) galaxies at $M_{\ast }\sim 10^{8}{\, \rm M_\odot}$ is missing in our simulations. This inconsistency persists across many other cosmological simulations with different galaxy formation models and demonstrates the potential of using galaxy morphology to constrain physics of galaxy formation at high redshifts.

     
    more » « less
  2. ABSTRACT

    The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at $z \gtrsim 6$ is an open question for understanding which galaxies primarily drove cosmic reionization. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful ‘near-far’ technique for studying the reionization-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass ($M_{\rm star}\lesssim 10^{9}\, \rm {M_\odot }$) galaxies across 11 Milky Way (MW)- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterize their progenitors at $z=6\!-\!9$, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with $M_{\rm star}\sim 10^5\, \rm {M_\odot }$ ($\sim 10^9\, \rm {M_\odot }$) had $\approx 1$ ($\approx 30$) progenitors at $z\approx 7$, and its main progenitor comprised $\approx 100~{{\ \rm per\ cent}}$ ($\approx 10~{{\ \rm per\ cent}}$) of the total stellar mass of all its progenitors at $z\approx 7$. We show that although only $\sim 15~{{\ \rm per\ cent}}$ of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at $z \sim 6 \!-\! 9$. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just low-mass galaxies at $z = 0$ and not the MW/M31 hosts correctly recovers the slope of the SMF down to $M_{\rm star} \sim 10^{4.5}\, \rm {{\rm M}_{\odot }}$ at $z \gtrsim 6$. Thus, we validate the ‘near-far’ approach as an unbiased tool for probing low-mass reionization-era galaxies.

     
    more » « less
  3. 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
  4. 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
  5. null (Ed.)
    ABSTRACT The galaxy size–stellar mass and central surface density–stellar mass relationships are fundamental observational constraints on galaxy formation models. However, inferring the physical size of a galaxy from observed stellar emission is non-trivial due to various observational effects, such as the mass-to-light ratio variations that can be caused by non-uniform stellar ages, metallicities, and dust attenuation. Consequently, forward-modelling light-based sizes from simulations is desirable. In this work, we use the skirt  dust radiative transfer code to generate synthetic observations of massive galaxies ($M_{*}\sim 10^{11}\, \rm {M_{\odot }}$ at z = 2, hosted by haloes of mass $M_{\rm {halo}}\sim 10^{12.5}\, \rm {M_{\odot }}$) from high-resolution cosmological zoom-in simulations that form part of the Feedback In Realistic Environments project. The simulations used in this paper include explicit stellar feedback but no active galactic nucleus (AGN) feedback. From each mock observation, we infer the effective radius (Re), as well as the stellar mass surface density within this radius and within $1\, \rm {kpc}$ (Σe and Σ1, respectively). We first investigate how well the intrinsic half-mass radius and stellar mass surface density can be inferred from observables. The majority of predicted sizes and surface densities are within a factor of 2 of the intrinsic values. We then compare our predictions to the observed size–mass relationship and the Σ1−M⋆ and Σe−M⋆ relationships. At z ≳ 2, the simulated massive galaxies are in general agreement with observational scaling relations. At z ≲ 2, they evolve to become too compact but still star forming, in the stellar mass and redshift regime where many of them should be quenched. Our results suggest that some additional source of feedback, such as AGN-driven outflows, is necessary in order to decrease the central densities of the simulated massive galaxies to bring them into agreement with observations at z ≲ 2. 
    more » « less