skip to main content

Explore Research Products in the NSF-PAR

Title: Quasi-equilibrium models of high-redshift disc galaxy evolution
ABSTRACT In recent years, simple models of galaxy formation have been shown to provide reasonably good matches to available data on high-redshift luminosity functions. However, these prescriptions are primarily phenomenological, with only crude connections to the physics of galaxy evolution. Here, we introduce a set of galaxy models that are based on a simple physical framework but incorporate more sophisticated models of feedback, star formation, and other processes. We apply these models to the high-redshift regime, showing that most of the generic predictions of the simplest models remain valid. In particular, the stellar mass–halo mass relation depends almost entirely on the physics of feedback (and is thus independent of the details of small-scale star formation) and the specific star formation rate is a simple multiple of the cosmological accretion rate. We also show that, in contrast, the galaxy’s gas mass is sensitive to the physics of star formation, although the inclusion of feedback-driven star formation laws significantly changes the naive expectations. While these models are far from detailed enough to describe every aspect of galaxy formation, they inform our understanding of galaxy formation by illustrating several generic aspects of that process, and they provide a physically grounded basis for extrapolating predictions to faint galaxies and high redshifts currently out of reach of observations. If observations show violations from these simple trends, they would indicate new physics occurring inside the earliest generations of galaxies.  more » « less
Award ID(s):
1812458
NSF-PAR ID:
10293904
Author(s) / Creator(s):
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
500
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3394 to 3412
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ( , The Astrophysical Journal)
    Abstract

    Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. The complexity of these physics can prove difficult to disentangle and obfuscate how mass-dependent trends in the galaxy population originate. Here, we train a machine-learning method called Explainable Boosting Machines (EBMs) to infer how the stellar mass and star formation rate of nearly 6 million galaxies simulated by the Cosmic Reionization on Computers project depend on the physical properties of halo mass, the peak circular velocity of the galaxy during its formation historyvpeak, cosmic environment, and redshift. The resulting EBM models reveal the relative importance of these properties in setting galaxy stellar mass and star formation rate, withvpeakproviding the most dominant contribution. Environmental properties provide substantial improvements for modeling the stellar mass and star formation rate in only ≲10% of the simulated galaxies. We also provide alternative formulations of EBM models that enable low-resolution simulations, which cannot track the interior structure of dark matter halos, to predict the stellar mass and star formation rate of galaxies computed by high-resolution simulations with detailed baryonic physics.

     
    more » « less
  2. ; ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    We post-process galaxies in the IllustrisTNG simulations with skirt radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at z ≥ 4. The rest-frame K- and z-band galaxy luminosity functions from TNG are overall consistent with observations, despite ${\sim}0.5\, \mathrm{dex}$ underprediction at z = 4 for MK ≲ −25 and Mz ≲ −24. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. Based on theoretical estimations, we show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of $500\, {\rm arcmin}^{2}$ at z = 6 (z = 8). As opposed to the consistency in the UV, optical, and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshift are qualitatively consistent with observations, the peak dust temperature of z ≥ 6 galaxies are overestimated by about $20\, {\rm K}$. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts.

     
    more » « less
  3. ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    The connection between the escape fraction of ionizing radiation (fesc) and the properties of galaxies, such as stellar mass ($\rm M_{\rm *}$), age, star-formation rate (SFR), and dust content, are key inputs for reionization models, but many of these relationships remain untested at high redshift. We present an analysis of a sample of 96 $z$ ∼ 3 galaxies from the Keck Lyman Continuum Spectroscopic Survey (KLCS). These galaxies have both sensitive Keck/LRIS spectroscopic measurements of the Lyman continuum (LyC) region, and multiband photometry that places constraints on stellar population parameters. We construct composite spectra from subsamples binned as a function of galaxy property and quantify the ionizing-photon escape for each composite. We find a significant anti-correlation between fesc and $\rm M_{\rm *}$, consistent with predictions from cosmological zoom-in simulations. We also find significant anti-correlation between fesc and E(B−V), encoding the underlying physics of LyC escape in our sample. We also find no significant correlation between fesc and either stellar age or specific SFR (= SFR/$\rm M_{\rm *}$), challenging interpretations that synchronize recent star formation and favorable conditions for ionizing escape. The galaxy properties now shown to correlate with fesc in the KLCS are Lyα equivalent width, UV Luminosity, $\rm M_{\rm *}$, SFR, and E(B−V), but not age or sSFR. This comprehensive analysis of galaxy properties and LyC escape at high redshift will be used to guide future models and observations of the reionization epoch.

     
    more » « less
  4. ; ; ; ; ; ; ; ( , The Astrophysical Journal)
    Abstract

    The epoch of reionization (EoR) offers a unique window into the dawn of galaxy formation, through which high-redshift galaxies can be studied by observations of both themselves and their impact on the intergalactic medium. Line intensity mapping (LIM) promises to explore cosmic reionization and its driving sources by measuring intensity fluctuations of emission lines tracing the cosmic gas in varying phases. Using LIMFAST, a novel seminumerical tool designed to self-consistently simulate LIM signals of multiple EoR probes, we investigate how building blocks of galaxy formation and evolution theory, such as feedback-regulated star formation and chemical enrichment, might be studied with multitracer LIM during the EoR. On galaxy scales, we show that the star formation law and the feedback associated with star formation can be indicated by both the shape and redshift evolution of LIM power spectra. For a baseline model of metal production that traces star formation, we find that lines highly sensitive to metallicity are generally better probes of galaxy formation models. On larger scales, we demonstrate that inferring ionized bubble sizes from cross-correlations between tracers of ionized and neutral gas requires a detailed understanding of the astrophysics that shape the line luminosity–halo mass relation. Despite various modeling and observational challenges, wide-area, multitracer LIM surveys will provide important high-redshift tests for the fundamentals of galaxy formation theory, especially the interplay between star formation and feedback by accessing statistically the entire low-mass population of galaxies as ideal laboratories, complementary to upcoming surveys of individual sources by new-generation telescopes.

     
    more » « less
  5. ; ; ; ; ( , The Astrophysical Journal)
    Abstract

    Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy’s star-forming efficiency. Mass, energy, and metal outflows (ηM,ηE, andηZ, here normalized by the star formation rate, the SNe energy, and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (CGM), preventing future accretion. Traditionally, models have assumed that galaxies self-regulate by ejecting a large fraction of the gas, which enters the interstellar medium (ISM), although whether such high mass loadings agree with observations is still unclear. To better understand how the relative importance of ejective (i.e., high mass loading) versus preventative (i.e., high energy loading) feedback affects the present-day properties of galaxies, we develop a simple gas-regulator model of galaxy evolution, where the stellar mass, ISM, and CGM are modeled as distinct reservoirs which exchange mass, metals, and energy at different rates within a growing halo. Focusing on the halo mass range from 1010to 1012M⊙, we demonstrate that, with reasonable parameter choices, we can reproduce the stellar-to-halo mass relation and the ISM-to-stellar mass relation with low-mass-loaded (ηM∼ 0.1–10) but high-energy-loaded (ηE∼ 0.1–1) winds, with self-regulation occurring primarily through heating and cooling of the CGM. We show that the model predictions are robust against changes to the mass loading of outflows but are quite sensitive to our choice of the energy loading, preferringηE∼ 1 for the lowest-mass halos and ∼0.1 for Milky Way–like halos.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email