Population III Stars: Early Cosmic Pioneers

• Population III stars are the first generation of metal-free stars formed in primordial minihalos, initiating reionization and chemical enrichment.
• Simulations show that their formation is clustered with a top-heavy initial mass function and significant disk fragmentation leading to multiple protostars.
• Observational constraints using EMP stars, 21-cm cosmology, and binary evolution models provide insights into their feedback, mass limits, and surviving low-mass remnants.

Population III (Pop III) stars comprise the first stellar generation to form from metal-free primordial gas, initiating the cosmic transition out of the Dark Ages. Theoretical and numerical studies demonstrate that Pop III stars are critical agents of early cosmic feedback, responsible for the onset of cosmological reionization and chemical enrichment. Despite their pivotal role, their properties—especially their mass function (IMF), multiplicity, and feedback—remain imprecisely constrained due to intrinsic observational inaccessibility, rapid self-enrichment of their host environments, and limitations in simulation resolution. Advances in 3D radiation-hydrodynamical simulations and “stellar archaeology” via extremely metal-poor (EMP) stars are providing increasingly stringent constraints.

1. Formation Physics and Early Environments

Pop III stars formed at redshifts $z\sim 30$–$20$ within minihalos of total mass $\sim 2\times10^5$–$1\times10^6\ M_\odot$ collapsing from primordial cosmological fluctuations (Glover et al., 18 Sep 2025). The dominant cooling pathway in these environments is via molecular hydrogen, as primordial gas lacks metals and thus cannot utilize fine-structure or dust cooling channels (Glover et al., 18 Sep 2025, Latif et al., 2021). The critical halo mass for efficient H$_2$ cooling is set by the requirement that the associated cooling time is less than the Hubble time,

$$M_{\rm crit}(z) \approx 3\times10^5\ M_\odot\ \left(\frac{1+z}{30}\right)^{-3/2},$$

with the actual threshold increased by baryon–dark-matter streaming and Lyman–Werner UV backgrounds (Glover et al., 18 Sep 2025, Latif et al., 2021). Gas collapse proceeds via isothermal contraction down to $T\sim200$ K and $n\sim10^4$ cm$^{-3}$, yielding a characteristic fragmentation mass scale $M_J\sim10^3\ M_\odot$ (Klessen et al., 2023).

Disk formation occurs from the conservation of angular momentum during collapse, and the characteristic infall rates are high,

$\dot{M} \sim c_s^3 / G \sim 10^{-3}$–$10^{-2}\ M_\odot\ {\rm yr}^{-1},$

where $c_s$ is the sound speed in 200–1000 K gas (Glover et al., 18 Sep 2025).

2. Multiplicity, Fragmentation and the Pop III Initial Mass Function

3D cosmological RHD simulations, as in (Latif et al., 2021), consistently show that Pop III star formation is inherently multiple, with each minihalo producing $\sim$10–20 protostellar fragments over $10^5$ years. Disk fragmentation occurs due to gravitational instability (Toomre $Q < 1$), resulting in the formation of small stellar clusters rather than isolated monolithic stars (Glover et al., 18 Sep 2025, Latif et al., 2021, Klessen et al., 2023). The birth-mass function derived from such simulations is well described by a power law,

$\frac{dN}{dM_*} \propto M_*^{-1.18},$

over the range 0.1–37 $M_\odot$. The mass distribution is top-heavy compared to that of present-day stars (Salpeter slope $-2.35$), but still numerically dominated by stars of characteristic mass $1$–$10\ M_\odot$ (median $\sim5\ M_\odot$) (Latif et al., 2021).

Dynamical effects are central: up to 70% of protostellar fragments are ejected from the disk via N-body or three-body interactions. Ejection terminates accretion for low-mass fragments, leading to a low-mass tail in the IMF. Ionizing feedback from massive protostars deforms the disk, induces outflows, and sets an upper-mass ceiling of $\sim30$–$40\ M_\odot$ in small primordial halos (Latif et al., 2021).

Table 1: Pop III Mass Spectrum Characteristics

Characteristic	Value(s)	Reference
Mass range	0.1–37 $M_\odot$ (simulated)	(Latif et al., 2021)
Power-law IMF slope	$-1.18$ (birth function, simulation)	(Latif et al., 2021)
Massive star fraction	$\sim$20% with $M>10\ M_\odot$	(Latif et al., 2021)
Ejection fraction	50–70% of fragments	(Latif et al., 2021)

3. Nucleosynthetic Feedback and Chemical Enrichment

Pop III stars end the cosmic Dark Ages by injecting metals and energetic feedback into their surroundings, enabling rapid self-enrichment and facilitating the transition to Population II (Pop II) star formation (Wise et al., 2010). Stars in the 8.5–90 $M_\odot$ range yield core-collapse supernovae, with no robust evidence for pair-instability SN (PISN) yields seen in the abundance patterns of EMP stars (Fraser et al., 2015). Simulations and chemical archaeology establish a rapid self-enrichment timescale; a single core-collapse SN can raise a minihalo’s metallicity to $Z_{\rm crit}\sim10^{-4}–10^{-3.5}\ Z_\odot$ in $10^6$–$10^7$ yr. This critical metallicity triggers a transition to Pop II formation via more efficient cooling and low-mass fragmentation (Wise et al., 2010, Kulkarni et al., 2013).

Owing to the rapid enrichment and inefficient mixing, Pop III stars dominate feedback at $z\gtrsim30$ but become negligible contributors to ionizing photon budgets by $z=10$, where their fractional ionizing output is $<1\%$ (Kulkarni et al., 2013).

4. Observational Constraints and Indirect Signatures

Direct detection of Pop III stars is precluded by their high redshift and faintness, but indirect constraints are robust and multi-faceted:

a. EMP and CEMP Stars

The abundance patterns of EMP stars ([Fe/H]$<-3.5$) in the Milky Way allow for a reconstruction of the Pop III progenitor mass function (Fraser et al., 2015). The best-fit IMF is Salpeter-like, with

$$\frac{dN}{dM} \propto M^{-2.35},\quad M_{\rm min} = 8.5^{+0.2}_{-0.4}\ M_\odot,\quad M_{\rm max} = 87^{+13}_{-33}\ M_\odot$$

(Fraser et al., 2015). No chemical or statistical evidence exists for significant numbers of progenitors above $\sim$120 $M_\odot$; PISN signatures are absent. The high [C/Fe] ratios observed in CEMP-no stars match yields from $20$–$120\ M_\odot$ Pop III progenitors (Sarmento et al., 2019); the CEMP fraction in EMPs constrains the minimum Pop III mass to 15–27 $M_\odot$ if the entire EMP population is attributed to Pop III pollution (Tanikawa et al., 2022).

b. Reionization and 21-cm Cosmology

Pop III stars generate copious H-ionizing photons, but self-limited SFRs and rapid enrichment restrict their role in cosmic reionization to $z>15$–$20$ (Tanaka et al., 2020, Kulkarni et al., 2013). Feedback processes—especially the escape fractions of ionizing and Lyman–Werner photons—are strongly mass- and environment-dependent. Accurate modelling of their impact is essential for interpreting 21-cm signals from cosmic dawn and reionization (Tanaka et al., 2020).

c. Gravitational Wave and High-z Transients

Massive Pop III binaries may produce long gamma-ray bursts and binary black hole mergers visible in gravitational waves, especially at $z>10$ where intrinsic event rates could be as high as 40% of all LGRBs (Campisi et al., 2011). Extreme-value statistics applied to high-resolution simulations indicate that the most massive Pop III stars at $z=10$–$20$ can reach $10^3$–$10^4\ M_\odot$, potentially acting as black hole seeds for supermassive quasars (Chantavat et al., 2023). Prospects for direct transient detection (e.g. PISNe) are restricted but potentially within JWST or large synoptic surveys’ capabilities for rare, high-mass progenitors (Glover et al., 18 Sep 2025).

5. Binary Evolution and Feedback Enhancement

Modern simulations and population synthesis establish that Pop III star formation yields a non-negligible binary fraction. Binary mass transfer, Roche-lobe overflow, and common envelope evolution modify the radiative and chemical yields of Pop III stars. Binary interactions can extend main-sequence lifetimes, nearly double collective UV photon output, and alter SN type outcomes (e.g., shifting compact remnant masses and suppressing/enabling PISNe) (Tsai et al., 2023). This complexity must be reflected in subgrid models for cosmological reionization and early enrichment.

6. Pop III Star Formation Across Cosmic Time and Spatial Isolation

Pop III star formation persists over an extended redshift window $30>z>6$ in chemically pristine halos, coexisting with enriched star formation but at much lower rates—by $z=10$ Pop III SFRD is $4$–$6$ orders of magnitude lower than that of metal-enriched stars (Crosby et al., 2013). Spatial isolation of Pop III–forming halos increases toward lower $z$, with most late-forming Pop III stars appearing in low-density “pockets” far from bright galaxies (Crosby et al., 2013).

Halo mass thresholds for Pop III formation rise as Lyman–Werner feedback strengthens, from $M_{\rm vir,\ min}\sim6\times10^4\ M_\odot$ at $z=30$ to $\sim10^7\ M_\odot$ at $z=10$ (Crosby et al., 2013). Local feedback, including both UV fields and metal enrichment, tightly regulates the temporal and spatial distribution of Pop III events (Wise et al., 2010, Muratov et al., 2012).

7. Survivors: Prospects for Low-Mass Pop III Detection

Long-lived Pop III stars with $M<0.8\ M_\odot$ remain a major target for discovery in the Milky Way and Local Group (Johnson, 2014, Komiya et al., 2016). Surface pollution models predict that metal-free survivors would generally acquire [Fe/H]$<-5$ due to ISM accretion, but rare ejected “escapees” from mini-halos may retain pristine surfaces, a prospect currently explored via kinematic and chemical tagging of high-latitude, extremely metal-poor stars (Komiya et al., 2016). Characteristic abundance signatures—enhanced [C/Fe], [O/Fe], and [Zn/Fe] due to dust-segregated accretion—could uniquely identify Pop III remnants (Johnson, 2014).

---

References

• (Latif et al., 2021) The Birth Mass Function of Pop III Stars
• (Fraser et al., 2015) The Mass Distribution of Population III Stars
• (Glover et al., 18 Sep 2025) The first stars
• (Klessen et al., 2023) The first stars: formation, properties, and impact
• (Wise et al., 2010) The Birth of a Galaxy: Primordial Metal Enrichment and Stellar Populations
• (Kulkarni et al., 2013) Chemical constraints on the contribution of Population III stars to cosmic reionization
• (Johnson, 2014) The chemical signature of surviving Population III stars in the Milky Way
• (Sarmento et al., 2019) Following the Cosmic Evolution of Pristine Gas III
• (Crosby et al., 2013) Population III Star Formation In Large Cosmological Simulations I
• (Muratov et al., 2012) Revisiting The First Galaxies: The epoch of Population III stars
• (Komiya et al., 2016) Population III stars around the Milky Way
• (Chantavat et al., 2023) The most massive Population III stars
• (Campisi et al., 2011) Population III stars and the Long Gamma Ray Burst rate
• (Tsai et al., 2023) The Evolution of Population III and Extremely Metal-Poor Binary Stars
• (Tanaka et al., 2020) Modelling Population III stars for semi-numerical simulations
• (Tanikawa et al., 2022) Can Population III stars be major origins of both merging binary black holes and extremely metal poor stars?