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$.

The most massive Population III stars

The first stars: formation, properties, and impact

Characterization of Population III Stars with Stellar Atmosphere and Evolutionary Modeling and Predictions of their Observability with the James Webb Space Telescope

Unveiling the contribution of Pop III stars in primeval galaxies at redshift $\geq 6$

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

Constraining Very High Mass Population III Stars through He II Emission in Galaxy BDF-521 at z = 7.01

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

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

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

The Challenges of Identifying Population III Stars in the Early Universe

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

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

Discrimination of heavy elements originating from Pop III stars in z = 3 intergalactic medium

Detecting Population III Stars through Tidal Disruption Events in the Era of JWST and Roman

The Evolution of Accreting Population III Stars at 10$^{-6}$ - 10$^3$ M$_\odot$/yr

A needle in a haystack? Catching Pop III stars in the Epoch of Reionization: I. Pop III star forming environments

Formation of Massive and Wide First-star Binaries in Radiation Hydrodynamics Simulations

Evolution of low mass population III stars from the pre-main sequence to the white dwarf cooling track

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

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