Title: A Framework to Calibrate a Semianalytic Model of the First Stars and Galaxies to the Renaissance Simulations
Abstract We present a method that calibrates a semianalytic model to the Renaissance Simulations, a suite of cosmological hydrodynamical simulations with high-redshift galaxy formation. This approach combines the strengths of semianalytic techniques and hydrodynamical simulations, enabling the extension to larger volumes and lower redshifts that are inaccessible to simulations due to computational expense. Using a sample of Renaissance star formation histories from an average density region of the Universe, we construct a four-parameter prescription for metal-enriched star formation characterized by an initial bursty stage followed by a steady stage where stars are formed at constant efficiencies. Our model also includes a treatment of Pop III star formation where a minimum halo mass and log-normal distribution of stellar mass are adopted to match the numerical simulations. Star formation is generally well reproduced for halos with masses ≲109M⊙. Between 11 <z< 25 our model produces metal-enriched star formation rate densities (SFRDs) that typically agree with Renaissance within a factor of ∼2 for the average density region. Additionally, the total metal-enriched stellar mass only differs from Renaissance by about 10% atz∼ 11. For regions that are either more overdense or rarefied and not included in the calibration, we produce metal-enriched SFRDs that agree with Renaissance within a factor of ∼2 at high-zbut eventually differ by higher factors for later times. This is likely due to environmental dependencies not included in the model. Our star formation prescriptions can easily be adopted in other analytic or semianalytic works to match our calibration to Renaissance.  more » « less
Award ID(s):
2009309 2108020
PAR ID:
10627195
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Astronomical Society
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
978
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
13
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present a new self-consistent semianalytic model of the first stars and galaxies to explore the high-redshift (z≥ 15) Population III (PopIII) and metal-enriched star formation histories. Our model includes the detailed merger history of dark matter halos generated with Monte Carlo merger trees. We calibrate the minimum halo mass for PopIII star formation from recent hydrodynamical cosmological simulations that simultaneously include the baryon–dark matter streaming velocity, Lyman–Werner (LW) feedback, and molecular hydrogen self-shielding. We find an overall increase in the resulting star formation rate density (SFRD) compared to calibrations based on previous simulations (e.g., the PopIII SFRD is over an order of magnitude higher atz= 35−15). We evaluate the effect of the halo-to-halo scatter in this critical mass and find that it increases the PopIII stellar mass density by a factor ∼1.5 atz≥ 15. Additionally, we assess the impact of various semianalytic/analytic prescriptions for halo assembly and star formation previously adopted in the literature. For example, we find that models assuming smooth halo growth computed via abundance matching predict SFRDs similar to the merger tree model for our fiducial model parameters, but that they may underestimate the PopIII SFRD in cases of strong LW feedback. Finally, we simulate subvolumes of the Universe with our model both to quantify the reduction in total star formation in numerical simulations due to a lack of density fluctuations on spatial scales larger than the simulation box, and to determine spatial fluctuations in SFRD due to the diversity in halo abundances and merger histories. 
    more » « less
  2. Abstract The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star-forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across broad mass (108.0M⊙≤M⋆ ≲ 1012.0M⊙) and redshift (0 ≤z≤ 8) ranges. We find that all simulations predict strong negative (i.e., radially decreasing) metallicity gradients at early cosmic times, likely due to their similar treatments of relatively smooth stellar feedback not providing sufficient mixing to quickly flatten gradients. The strongest redshift evolution occurs in galaxies with stellar masses of 1010.0–1011.0M⊙, while galaxies with stellar mass < 1010M⊙and >1011M⊙exhibit weaker redshift evolution. Our result of negative gradients at high redshift contrast with the many positive and flat gradients in the 1 < z < 4 observational literature. Atz > 6, the negative gradients observed with JWST and the Atacama Large Millimeter/submillimeter Array are flatter than those in simulations, albeit with closer agreement than at lower redshift. Overall, we suggest that these smooth stellar feedback galaxy simulations may not sufficiently mix their metal content radially, and that either stronger stellar feedback or additional subgrid turbulent metal diffusion models may be required to better reproduce observed metallicity gradients. 
    more » « less
  3. Abstract Recent observations suggest a nearly constant gas-phase mass–metallicity relation (MZR) atz ≳ 5, in agreement with many theoretical predictions. This lack of evolution contrasts with observations atz ≲ 3, which find an increasing normalization of the MZR with decreasing redshift. We analyze a high-redshift suite of FIRE-2 cosmological zoom-in simulations to identify the physical drivers of the MZR. Previous studies have explained the weak evolution of the high-redshift MZR in terms of weakly evolving or saturated gas fractions, but we find that this alone does not explain the evolution in FIRE-2. Instead, stellar feedback following intense bursts of star formation drives enriched gas out of galaxies, resetting their interstellar medium and separating their histories into distinct “burst cycles.” We develop the “reduced burst model,” a simplified gas-regulator model that successfully reproduces the simulated MZR and identifies the dominant drivers of its evolution. As redshift decreases, the metallicity of inflows within burst cycles increases at fixed stellar mass due to increased wind recycling of enriched gas. Meanwhile, the metal mass produced by stars per inflowing gas mass within these cycles decreases because of decreased star formation per gas mass inflowing into the galaxy. The effects of these two processes on the median metallicity largely cancel, holding the MZR constant forz = 5–12. At fixed stellar mass, the simulations predict lower gas metallicities at higher Hα-derived star formation rates, in qualitative agreement with the fundamental metallicity relation, but this effect is reduced in rest UV-selected samples. 
    more » « less
  4. We analysed a sample of 552 galaxies from the LEGA-C spectroscopic survey (0.6 < z < 0.77), for which we estimated the stellar population parameters by a Bayesian analysis of the stellar absorption features and photometry. We investigated the effect of the current star formation activity on light-weighted mean stellar ages and metallicities and their median trends with stellar mass or velocity dispersion. The bimodality in the global age–mass relation stems from the different age distributions in the quiescent and star-forming populations. No bimodality is observed in the stellar metallicity-mass relation, although quiescent and star-forming galaxies have different distributions in this parameter space. We identified a high-metallicity sequence populated by quiescent and weakly star-forming galaxies. At masses lower than 1010.8 M⊙, the median stellar metallicity–mass relation of star-forming galaxies steepens as a consequence of the increasing scatter towards lower stellar metallicities for galaxies with an increasing specific star formation rate at fixed mass. Relying on a consistent analysis of SDSS DR7 spectra and accounting for aperture corrections, we quantified the evolution of the volume-weighted stellar age and stellar metallicity scaling relations betweenz = 0.7 and the present. We found negligible evolution in the stellar metallicity–mass relation of quiescent galaxies and forM* > 1011 M⊙galaxies in general. Lower-mass star-forming galaxies instead have typically lower metallicities than their local counterparts, indicating significant enrichment sincez ∼ 0.7 in the low-mass regime. Notably, the median of the stellar ages of the general population and of quiescent galaxies has changed by only 2 Gyr betweenz = 0.7 andz = 0.1, which is less than expected from cosmic ageing. Some quiescent galaxies must evolve passively to reach the old boundary of the local population. In order to explain the evolution of the median trends, however, both individual evolution through rejuvenation and/or minor merging that affects the outer galaxy regions and population evolution through quenching of massive metal-rich star-forming galaxies are required. 
    more » « less
  5. Abstract We investigate the impact of massive primordial black holes (PBHs;mBH ∼ 106M⊙) on the star formation and first galaxy assembly process using high-resolution hydrodynamical simulations fromz= 1100 toz ∼ 9. We find that PBH accretion is self-regulated by feedback, suppressing mass growth unless feedback is weak. PBHs accelerate structure formation by seeding dark matter (DM) halos and gravitationally attracting gas, but strong feedback can delay cooling and suppress star formation. In addition, the presence of baryon-DM streaming creates an offset between the PBH location and the peaks induced in gas density, promoting earlier and more efficient star formation compared to standard ΛCDM. Byz ∼ 10, PBH-seeded galaxies form dense star clusters, with PBH-to-stellar mass ratios comparable to observed high-zactive galactic nuclei like UHZ-1. Our results support PBHs as viable supermassive black hole (SMBH) seeds but do not exclude alternative scenarios. We emphasize that PBH-seeding provides a natural explanation for some of the newly discovered overmassive SMBHs at high redshift, in particular those with extreme ratios of BH-to-dynamical (virial) mass that challenge standard formation channels. Future studies with ultra-deep JWST surveys, the Roman Space Telescope, and radio surveys with facilities such as the Square Kilometre Array and Hydrogen Epoch of Reionization Array will be critical in distinguishing PBH-driven SMBH growth from other pathways. 
    more » « less