Abstract:ABSTRACT JWST has brought us new insights into Cosmic Dawn with tentative detection of the unique signatures of metal-free Population III (Pop III) stars, such as strong He II emission, extremely blue ultraviolet spectrum, and enhanced nitrogen abundance. Self-consistent theoretical predictions of the formation rates, sites, and masses of Pop III stars are crucial for interpreting the observations, but are challenging due to complex physical processes operating over the large range of length-scales involved. One solution is to combine analytical models for the small-scale star formation process with cosmological simulations that capture the large-scale physics such as structure formation, radiation backgrounds, and baryon-dark matter streaming motion that regulate the conditions of Pop III star formation. We build an analytical model to predict the final masses of Pop III stars/clusters from the properties of star-forming clouds, based on the key results of small-scale star formation simulations and stellar evolution models. Our model for the first time considers the interplay between feedback and fragmentation and covers different modes of Pop III star formation ranging from ordinary small ($\sim\!{10{-}2000}\ \rm M_\odot$) clusters in molecular-cooling clouds to massive ($\gtrsim\!{10^{4}}\ \rm M_\odot$) clusters containing supermassive ($\sim\!{10^{4}{-}3}\times 10^{5}\ \rm M_\odot$) stars under violent collapse of atomic-cooling clouds with large gas accretion rates of $\gtrsim\!{0.1}\ \rm M_\odot \ yr^{-1}$. As an example, the model is applied to the Pop III star-forming clouds in the progenitors of typical haloes hosting high-z luminous quasars ($M_{\rm h}\sim 10^{12}\ \rm M_\odot$ at $z\sim 6$), which shows that formation of Pop III massive clusters is common ($\sim\!{20{-}70}{{\ \rm per\ cent}}$) in such biased ($\sim\!{4}\sigma$) regions, and the resulting heavy black hole seeds from supermassive stars can account for a significant fraction of observed luminous ($\gtrsim\!{10^{46}}\ \rm erg\ s^{-1}$) quasars at $z\sim 6$.

Building up the Population III initial mass function from cosmological initial conditions

How the Population III Initial Mass Function Governs the Properties of the First Galaxies

Population III star formation in the presence of turbulence, magnetic fields and ionizing radiation feedback

A self-consistent semi-analytic model for Population III star formation in minihalos

The Birth of a Galaxy: Primordial Metal Enrichment and Stellar Populations

Towards a universal analytical model for Population III star formation: interplay between feedback and fragmentation

Following the Cosmic Evolution of Pristine Gas III: The Observational Consequences of the Unknown Properties of Population III Stars

Population III Star Formation In Large Cosmological Simulations I. Halo Temporal and Physical Environment

Explaining the high nitrogen abundances observed in high-z galaxies via population III stars of a few thousand solar masses

Simulations of early structure formation: Properties of halos that host primordial star formation

Self-consistent semi-analytic models of the first stars

The evolution of accreting population III stars at 10^{-6} - 10^3 Ms yr^{-1}

Semi-analytic modelling of Pop. III star formation and metallicity evolution - I. Impact on the UV luminosity functions at z = 9 − 16

Kindling the First Stars. I. Dependence of Detectability of the First Stars with JWST on the Population III Stellar Masses

Rise of the First Super-Massive Stars

Galaxy formation with cold gas accretion and evolving stellar initial mass function

Semi-analytic modelling of Pop. III star formation and metallicity evolution -- I. Impact on the UV luminosity functions at z = 9-16

The Growth of Black Holes from Population III Remnants in the Renaissance Simulations

The first stars: formation, properties, and impact