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.
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- PAR ID:
- 10533842
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 533
- Issue:
- 1
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 1059-1074
- Size(s):
- p. 1059-1074
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT JWST has revealed a large population of UV-bright galaxies at $z\gtrsim 10$ and possibly overly massive galaxies at $z\gtrsim 7$, challenging standard galaxy formation models in the ΛCDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tension. Our benchmark model demonstrates excellent agreement with the UV luminosity functions (UVLFs) at $4\lesssim z \lesssim 10$ in both ΛCDM and EDE cosmologies. In the EDE cosmology, the UVLF measurements at $z\simeq 12$ based on spectroscopically confirmed galaxies (eight galaxies at $z\simeq 11\!-\!13.5$) exhibit no tension with the benchmark model. Photometric constraints at $12 \lesssim z\lesssim 16$ can be fully explained within EDE via either moderately increased star-formation efficiencies ($\epsilon _{\ast}\sim 3\!-\!10\ \hbox{per cent}$ at $M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$) or enhanced UV variabilities ($\sigma _{\rm UV}\sim 0.8\!-\!1.3$ mag at $M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$) that are within the scatter of hydrodynamical simulation predictions. A similar agreement is difficult to achieve in $\Lambda$CDM, especially at $z\gtrsim 14$, where the required $\sigma _{\rm UV}$ exceeds the maximum value seen in simulations. Furthermore, the implausibly large cosmic stellar mass densities inferred from some JWST observations are no longer in tension with cosmology when the EDE is considered. Our findings highlight EDE as an intriguing unified solution to a fundamental problem in cosmology and the recent tensions raised by JWST observations. Data at the highest redshifts reached by JWST will be crucial for differentiating modified galaxy formation physics from new cosmological physics.
-
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.
-
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
-
null (Ed.)ABSTRACT We present the escape fraction of hydrogen ionizing photons (fesc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z ≥ 5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 haloes in Mvir ∼ 108–$10^{12}\, M_{\odot }$ (M* ∼ 104–$10^{10}\, M_{\odot }$) at z = 5–12. We find the sample average 〈fesc〉increases with halo mass for Mvir ∼ 108–$10^{9.5}\, M_{\odot }$, becomes nearly constant for 109.5–$10^{11}\, M_{\odot }$, and decreases at ${\gtrsim}10^{11}\, M_{\odot }$. Equivalently, 〈fesc〉 increases with stellar mass up to $M_{\ast }\sim 10^8\, M_{\odot }$ and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate 〈fesc〉 ∼ 0.2 for galaxies at $M_{\ast }\sim 10^8\, M_{\odot }$. Nearly half of the escaped ionizing photons come from stars 1–3 Myr old and the rest from stars 3–10 Myr old. Binaries only have a modest effect, boosting 〈fesc〉 by ∼25–35 per cent and the number of escaped photons by 60–80 per cent. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of 〈fesc〉 at the high-mass end is due to dust attenuation, while at the low-mass end, 〈fesc〉 decreases owing to inefficient star formation and hence feedback. At fixed mass, 〈fesc〉 tends to increase with redshift. Although the absolute 〈fesc〉does not fully converge with resolution in our simulations, the mass- and redshift-dependence of 〈fesc〉 is likely robust. Our simulations produce sufficient ionizing photons for cosmic reionization.more » « less
