51ÁÔÆæÈë¿Ú

Abstract:

We use intergalactic medium (IGM) metallicity distributions from several state-of-the-art cosmological simulations of Milky Way analogs and a semi-analytic model of ultra-faint dwarf galaxy (UFD) formation to model the stellar metallicities of UFDs in MW-like environments. We study simulations with different treatments of star formation, stellar feedback, and Population III enrichment, and in all cases, we find that only a few percent of the IGM accretable by UFD progenitors is enriched to metallicities [Fe/H] ≥ −4. When the metallicity of accreted IGM in the semi-analytic galaxy formation model is set using these IGM metallicity distributions, the model underpredicts UFD metallicities and their scatter compared to the observed luminosity–metallicity relation. Our results indicate that IGM enrichment is not the dominant mechanism setting metallicities of UFD stars. Instead, UFD stellar metallicity is determined primarily by the interplay between internal enrichment and metal loss through feedback-driven outflows. We examine models with different values of the maximum outflow mass loading factor nmax and show that the full range of average stellar metallicities of UFDs at Mv < −7 can be reproduced if the maximum mass loading factor varies in the range 200 ≲ nmax ≲ 2000. We also consider stellar metallicity distribution functions (MDFs) within individual model galaxies with different assumptions about IGM enrichment and nmax. We find that all considered models are in reasonable agreement with observed UFD MDFs, with model differences less than the uncertainties of current metallicity measurements.

Abstract:

<jats:p>Star formation is a key process that governs the baryon cycle within galaxies, however, the question of how it controls their growth remains elusive due to modeling uncertainties. To understand the impact of star formation models on galaxy evolution, we performed cosmological zoom-in radiation-hydrodynamic simulations of a dwarf dark matter halo, with a virial mass of <jats:italic>M</jats:italic><jats:sub>vir</jats:sub> ∼ 10<jats:sup>9</jats:sup> <jats:italic>M</jats:italic><jats:sub>⊙</jats:sub> at <jats:italic>z</jats:italic> = 6. We compared two different star formation models: a multi-freefall model combined with a local gravo-thermo-turbulent condition and a more self-consistent model based on a sink particle algorithm, where gas accretion and star formation are directly controlled by the gas kinematics. As the first study in this series, we used cosmological zoom-in simulations with different spatial resolutions and found that star formation is more bursty in the runs with the sink algorithm, generating stronger outflows than in the runs with the gravo-thermo-turbulent model. The main reason for the increased burstiness is that the gas accretion rates on the sinks are high enough to form stars on very short timescales, leading to more clustered star formation. As a result, the star-forming clumps are disrupted more quickly in the sink run due to more coherent radiation and supernova feedback. The difference in burstiness between the two star formation models becomes even more pronounced when the supernova explosion energy is artificially increased. Our results suggest that improving the modeling of star formation on small, sub-molecular cloud scales can significantly impact the global properties of simulated galaxies.</jats:p>

Abstract:

Constraining the escape fraction of Lyman Continuum (LyC) photons from high-redshift galaxies is crucial to understanding reionization. Recent observations have demonstrated that various characteristics of the Ly α emission line correlate with the inferred LyC escape fraction (f LyC esc ) of low-redshift galaxies. Using a data set of 9600 mock Ly α spectra of star-forming galaxies at 4.64 ≤ z ≤ 6 from the SPHINX20 cosmological radiation hydrodynamical simulation, we study the physics controlling the escape of Ly α and LyC photons. We find that our mock Ly α observations are representative of high-redshift observations and that typical observational methods tend to overpredict the Ly α escape fraction (f Ly α esc ) by as much as 2 dex. We investigate the correlations between f LyC esc and f Ly α esc , Ly α equivalent width (Wλ(Ly α)), peak separation (vsep), central escape fraction (fcen), and red peak asymmetry (Ared f ). We find that f Ly α esc and fcen are good diagnostics for LyC leakage, selecting for galaxies with lower neutral gas densities and less UV attenuation that have recently experienced supernova feedback. In contrast, Wλ(Ly α) and vsep are found to be necessary but insufficient diagnostics, while Ared f carries little information. Finally, we use stacks of Ly α, H α, and F150W mock surface brightness profiles to find that galaxies with high f LyC esc tend to have less extended Ly α and F150W haloes but larger H α haloes than their non-leaking counterparts. This confirms that Ly α spectral profiles and surface brightness morphology can be used to better understand the escape of LyC photons from galaxies during the epoch of reionization.