We present a semi-analytic, computationally inexpensive model to identify halos capable of forming a Population III star in cosmological simulations across a wide range of times and environments. This allows for a much more complete and representative set of Population III star forming halos to be constructed, which will lead to Population III star formation simulations that more accurately reflect the diversity of Population III stars, both in time and halo mass. This model shows that Population III and chemically enriched stars coexist beyond the formation of the first generation of stars in a cosmological simulation until at least z ∼ 10, and likely beyond, though Population III stars form at rates that are 4-6 orders of magnitude lower than chemically enriched stars by z = 10. A catalog of more than 40,000 candidate Population III forming halos were identified, with formation times temporally ranging from z = 30 to z = 10, and ranging in mass from 2.3 × 10 5 M to 1.2 × 10 10 M . At early times, the environment that Population III stars form in is very similar to that of halos hosting chemically enriched star formation. At later times Population III stars are found to form in low-density regions that are not yet chemically polluted due to a lack of previous star formation in the area. Population III star forming halos become increasingly spatially isolated from one another at later times, and are generally closer to halos hosting chemically enriched star formation than to another halo hosting Population III star formation by z ∼ 10.

This paper presents the first results from a model for chemical evolution that can be applied to N-body cosmological simulations and quantitatively compared to measured stellar abundances from large astronomical surveys. This model convolves the chemical yield sets from a range of stellar nucleosynthesis calculations (including asymptotic giant branchstars, Type Ia and II supernovae, and stellar wind models) with a user-specified stellar initial mass function (IMF) and metallicity to calculate the time-dependent chemical evolution model for a "simple stellar population" (SSP) of uniform metallicity and formation time. These SSP models are combined with a semianalytic model for galaxy formation and evolution that uses merger trees from N-body cosmological simulations to track several α-and iron-peak elements for the stellar and multiphase interstellar medium components of several thousand galaxies in the early (z 6) universe. The simulated galaxy population is then quantitatively compared to two complementary datasets of abundances in the Milky Way stellar haloand is capable of reproducing many of the observed abundance trends. The observed abundance ratio distributions arebest reproduced with a Chabrier IMF, a chemically enriched star formation efficiency of 0.2, and a redshift of reionization of 7. Many abundances are qualitatively well matched by our model, but our model consistently overpredicts the carbon-enhanced fraction of stars at low metallicities, likely owingto incomplete coverage of Population III stellar yields and supernova modelsand the lack of dust as a component of our model.

No abstract