arxiv: 2604.26006 · v1 · submitted 2026-04-28 · 🌌 astro-ph.GA · astro-ph.CO

Recognition: unknown

Towards A Universal Analytical Model of Population III Star Formation: A Bridge Between Cosmological Scales and Protostars

James Gurian , Boyuan Liu , Donghui Jeong , Takashi Hosokawa , Shingo Hirano , Volker Bromm , Naoki Yoshida

Authors on Pith no claims yet

Pith reviewed 2026-05-07 15:22 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO

keywords Population III starsstar formation efficiencyanalytical modelLyman-Werner radiationprotostellar diskscosmological scaleshalo-scale star formation

0 comments

The pith

An analytical model connects cosmological radiation backgrounds to sub-AU protostellar disk fragmentation for Population III star formation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper builds a chain of analytical models that together span from the large-scale cosmological environment down through dark-matter halos and gas clouds to the small-scale fragmentation inside accreting protostellar disks. Each segment is calibrated against existing simulations and then linked to the next, producing a single, fast-to-evaluate description of how pristine gas turns into the first stars. The resulting framework shows that the fraction of halo gas converted into Population III stars can change by more than two orders of magnitude depending on halo mass and the strength of the Lyman-Werner radiation field, with abrupt jumps when cooling switches from HD to H2 to atomic lines.

Core claim

By stitching together scale-specific analytical prescriptions, the authors obtain a computationally inexpensive model that reproduces simulation results across the full dynamic range from cosmological backgrounds to sub-AU disks. The model predicts halo-scale star-formation efficiencies ranging from roughly 0.001 to 0.5, set by the transitions among HD cooling, molecular-hydrogen cooling, and atomic cooling, while the efficiency inside individual unstable clouds remains high, above 0.2.

What carries the argument

A linked chain of three analytical modules: (1) cosmological radiation background to halo properties, (2) halo to star-forming cloud, and (3) cloud to fragmenting protostellar disk, each validated separately against simulations before being joined.

If this is right

Halo-scale Pop III star-formation efficiency changes by more than two orders of magnitude with Lyman-Werner intensity and virial temperature.
Sharp efficiency jumps appear exactly where cooling mechanisms switch from HD to H2 to atomic lines.
Cloud-scale conversion of unstable gas into stars remains efficient (greater than or equal to 0.2) across the explored conditions.
The full model remains inexpensive enough to scan wide ranges of radiation background and halo properties.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The framework could be used to generate large statistical samples of Pop III formation sites for input into reionization or chemical-enrichment calculations without running new simulations for every parameter combination.
If the model holds, the sharp efficiency transitions may imprint observable features on the spatial distribution or metallicity patterns of the earliest stars.
Extending the same modular approach to include later metal-line cooling or different initial mass functions would test whether the same linking logic remains valid once the first supernovae enrich the gas.

Load-bearing premise

That interactions between the different scales can be ignored and that agreement of each piece with its own simulations is enough to trust the combined prediction under all conditions.

What would settle it

A single high-resolution simulation that evolves one halo from cosmological scales through disk fragmentation and reports a star-formation efficiency differing by more than a factor of a few from the model's output for the same Lyman-Werner flux and halo mass.

Figures

Figures reproduced from arXiv: 2604.26006 by Boyuan Liu, Donghui Jeong, James Gurian, Naoki Yoshida, Shingo Hirano, Takashi Hosokawa, Volker Bromm.

**Figure 1.** Figure 1: Connections between the cosmological context (red) and the three models of relevant scales in this work: the halo scale (orange), the cloud scale (green), and the disk scale (blue). The solid arrows show the workflow of our combined model: The output of a larger-scale model serves as key input for the subsequent smaller-scale model, and a minimal set of additional free parameters are considered to capture … view at source ↗

**Figure 2.** Figure 2: The density-temperature evolution for four cases which illustrate the possible cooling channels, labeled by the parameters which differ from their fiducial values ( view at source ↗

**Figure 3.** Figure 3: The minimum temperature of the collapsing gas for various halo masses, redshifts, and radiation field intensities. Note that white corresponds approximately to the minimum temperature attainable by H2 cooling, so that dark blue indicates HD cooling while light red indicates weak molecular cooling and dark red, atomic cooling. In the bottom right corner of the top panels, the virial temperature is so low th… view at source ↗

**Figure 4.** Figure 4: Left: the ratio of the mass enclosed to the modified Bonnor-Ebert mass corresponding to the thermal evolution tracks shown in view at source ↗

**Figure 5.** Figure 5: The simulations of Sugimura et al. (2023) compared with the model results for three clouds labeled by the turbulent Mach number and rotation parameter. Left, center, and right correspond respectively to the high, intermediate, and low M˙ cases of Sugimura et al. (2023). Top row: the density-temperature relationships from the simulations of Sugimura et al. (2023) (thick, solid lines), extrapolated to lower … view at source ↗

**Figure 6.** Figure 6: The stellar mass in the single star case at the time of disk photoevaporation, the end of the stellar lifetime, or the onset of general-relativistic instability, for the rescaled infall rates for the molecular cooling case (solid) and the atomic cooling case (dashed). Also shown are the simulation samples of Hirano et al. (2014) and Toyouchi et al. (2023). with Liu et al. (2024), we label the fiducial inf… view at source ↗

**Figure 7.** Figure 7: The total stellar mass (left), most massive star (center), and mean stellar mass (right) from the simulations of Hirano et al. (2014); Toyouchi et al. (2023); Sugimura et al. (2023) and the model of this work in the case with a single star, and with power law mass function α = −0.9, and α = −2.3. For the case with α = −0.9, we show both m1 = 1 M⊙ and m1 = 15 M⊙. The discontinuities in the model results are… view at source ↗

**Figure 9.** Figure 9: The stellar masses predicted by the combined model of this work (background color) assuming no fragmentation as in Section 2.3.1, compared with the simulations of Hirano et al. (2014) (colored points). The transition from orange to purple is due to the delayed collapse and concomitant turnon of HD cooling in the slowly cooling halos, as determined by Eq. (2.1) (compare with view at source ↗

**Figure 10.** Figure 10: The predicted total stellar mass for different values of J21 under the fiducial fragmentation model, as well as the bounds on the region where star formation is most likely to occur as in view at source ↗

**Figure 11.** Figure 11: The halo-scale SFE for the same cases as view at source ↗

**Figure 12.** Figure 12: The cloud-scale SFE for the same cases as view at source ↗

read the original abstract

We construct an analytical model of Population III star formation that connects the cosmological radiation background to sub-AU protostellar disk fragmentation, a dynamic range inaccessible to any single simulation. Our approach is based on combining separate models of the disparate relevant scales: from the cosmological environment to the host-halo scale, from the halo scale to the scale of the star-forming cloud, and from the cloud scale to the fragmenting, accreting protostellar disk. Individually and collectively, the models agree well with the predictions of state of the art simulations, while remaining computationally inexpensive and physically transparent. As an example of the applicability of the model, we study the effects of varying the Lyman-Werner flux on the Pop. III star formation efficiency. We show that depending on the halo properties and the strength of the dissociating radiation field, the halo-scale Pop. III star formation efficiency varies by more than two orders of magnitude from $\varepsilon_{\rm SFE,H} \approx 10^{-3}$ to $\varepsilon_{\rm SFE, H} \approx 0.5$. The abrupt transitions between hydrogen-deuteride cooling (in low virial temperature mini-halos subjected to low radiation backgrounds), molecular hydrogen cooling (at intermediate temperatures and radiation intensities), and atomic cooling (in higher temperature halos exposed to strong radiation fields) produces sharp features in the halo-scale star formation efficiency as a function of the halo properties. Meanwhile, at the scale of individual star-forming clouds, the star formation efficiency is $\varepsilon_{\rm SFE,c} \gtrsim 0.2$. That is, pristine gas in a halo is converted into unstable clouds at a wide range of efficiencies, and these unstable clouds are efficiently converted into Pop. III stars.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper chains separate analytical models to estimate Pop III efficiencies from cosmology down to protostellar disks, which is practical for parameter sweeps even if the module handoffs are not fully re-validated.

read the letter

The main takeaway is that this paper builds a chained analytical model for Pop III star formation efficiencies that spans from the Lyman-Werner background all the way down to protostellar disks. It reproduces trends from simulations at each step and runs cheaply enough to scan parameters. They do a good job showing how efficiency can swing by orders of magnitude with halo properties and radiation strength, with clear jumps at the HD, H2, and atomic cooling regimes. The cloud-scale efficiency stays high at over 20 percent, which is a clean result. The potential issue is whether the modules connect smoothly. Each part is tuned to its own simulations, but the paper may not have cross-checked the full pipeline against integrated runs that cover multiple scales at once. That leaves room for the reported sharp transitions to shift if the assumed distributions at the interfaces don't hold. This kind of tool is aimed at people building semi-analytic models of the early universe who need fast estimates for star formation rates in pristine halos. It could feed into reionization calculations without much cost. I would send it to referees. The synthesis is new enough and the comparisons are solid enough to warrant a full review, even with some work needed on validating the links.

Referee Report

2 major / 1 minor

Summary. The manuscript constructs an analytical model of Population III star formation by combining separate scale-specific models (cosmological radiation background to halo virial temperature and cooling, halo to star-forming cloud, and cloud to sub-AU protostellar disk fragmentation). It reports good agreement with state-of-the-art simulations individually and collectively, and as an application shows that halo-scale star formation efficiency ε_SFE,H varies by more than two orders of magnitude (from ~10^{-3} to ~0.5) with Lyman-Werner flux and halo properties, with sharp transitions at HD/H2/atomic cooling boundaries, while cloud-scale efficiency ε_SFE,c ≳ 0.2.

Significance. If the chained model holds, it supplies a computationally inexpensive and physically transparent framework for predicting Pop III efficiencies over a dynamic range inaccessible to single simulations, which could be incorporated into semi-analytic models or large-volume cosmological runs. The reported efficiency variations and sharp cooling-regime features are concrete, falsifiable outputs that highlight the role of the dissociating background. The approach of linking calibrated modules is a strength when interfaces are shown to be robust.

major comments (2)

[model construction and example application sections] The central claim that the linked model bridges scales without missing interactions rests on the untested premise that outputs from one module (e.g., halo-scale SFE and cloud mass function) are statistically compatible with inputs to the next (e.g., disk initial conditions). No re-validation against any integrated multi-scale simulation is reported, so mismatches in angular momentum or density PDF could propagate and alter the final efficiencies or stellar masses. This directly affects the 'universal' and 'bridge' assertions.
[results on efficiency variation] The abrupt jumps in ε_SFE,H at the HD/H2/atomic cooling boundaries are presented as physical features, but they inherit directly from the separate cooling-regime models; without explicit propagation of uncertainties or interface consistency checks, it is unclear whether these discontinuities survive in the full chain.

minor comments (1)

Notation for efficiencies (ε_SFE,H vs ε_SFE,c) is clear in the abstract but should be defined explicitly at first use in the main text with reference to the relevant equations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and insightful comments, which have helped us improve the clarity and robustness of the manuscript. We address each major comment point by point below, with revisions made where appropriate to strengthen the presentation of the model interfaces and the physical nature of the efficiency features.

read point-by-point responses

Referee: [model construction and example application sections] The central claim that the linked model bridges scales without missing interactions rests on the untested premise that outputs from one module (e.g., halo-scale SFE and cloud mass function) are statistically compatible with inputs to the next (e.g., disk initial conditions). No re-validation against any integrated multi-scale simulation is reported, so mismatches in angular momentum or density PDF could propagate and alter the final efficiencies or stellar masses. This directly affects the 'universal' and 'bridge' assertions.

Authors: We agree that a direct end-to-end validation against a single integrated simulation spanning all scales would provide the strongest test of interface compatibility. However, no such simulation exists or is feasible at present due to the prohibitive computational cost of resolving the full dynamic range from cosmological to sub-AU scales. Our model instead validates each module individually against dedicated simulations at its respective scale and ensures that the statistical outputs (e.g., cloud mass functions, densities, and angular momenta) from one module lie within the validated input ranges of the next. In the revised manuscript, we have added a new subsection explicitly discussing the interface assumptions, including quantitative checks on angular momentum and density PDF consistency between modules. We have also moderated the language around the 'universal' and 'bridge' claims to emphasize that the model is a chained framework with validated components rather than a fully self-consistent multi-scale simulation. These changes constitute a partial revision. revision: partial
Referee: [results on efficiency variation] The abrupt jumps in ε_SFE,H at the HD/H2/atomic cooling boundaries are presented as physical features, but they inherit directly from the separate cooling-regime models; without explicit propagation of uncertainties or interface consistency checks, it is unclear whether these discontinuities survive in the full chain.

Authors: The sharp transitions originate from the distinct and well-validated cooling physics (HD, H2, and atomic) that define each regime, as confirmed by the individual module comparisons to simulations. To demonstrate that these features persist in the chained model, the revised manuscript now includes an explicit uncertainty propagation analysis: we vary key inputs (Lyman-Werner flux, halo virial temperature, and radiation background) within their plausible ranges and show that the discontinuities remain robust features rather than artifacts. This analysis is supported by the interface consistency checks added in response to the first comment. We have updated the results section and added a supplementary figure illustrating the sensitivity of ε_SFE,H to input variations. revision: yes

Circularity Check

0 steps flagged

No circularity: composite model chains externally validated scale-specific components

full rationale

The derivation chains independent analytical models across cosmological, halo, cloud, and disk scales, each reported to match separate simulation suites. No equation or claim reduces an output (such as ε_SFE,H or ε_SFE,c) to a fitted input by construction, nor does any load-bearing step rely on a self-citation whose content is itself unverified or defined circularly. The reported variations with Lyman-Werner flux and cooling regime boundaries follow from the explicit combination of the modules rather than from renaming or self-definition. The interfaces are presented as direct mappings whose validity is asserted via per-module simulation agreement, leaving the overall construction self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The model rests on the assumption that component models at each scale are sufficiently accurate and combinable; no free parameters or invented entities are mentioned in the abstract.

axioms (1)

domain assumption Separate analytical models for cosmological, halo, cloud, and disk scales can be linked without significant missing cross-scale physics.
Invoked by the claim that the combined model connects all scales and matches simulations.

pith-pipeline@v0.9.0 · 5654 in / 1137 out tokens · 47748 ms · 2026-05-07T15:22:08.269721+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Metal Enrichment by the First Stars Exploding at the Lower Energy Limit of Pair-Instability Supernovae
astro-ph.GA 2026-05 unverdicted novelty 5.0

Low-energy PISNe from 140 solar-mass Pop III stars produce second-generation stars at median [Fe/H] ~ -5.5 with odd-even patterns, but their absence from EMP observations disfavors PISNe as the main early enrichment channel.

Reference graph

Works this paper leans on

3 extracted references · 2 canonical work pages · cited by 1 Pith paper

[1]

E., 2018, ApJ, 857, 138 Beers T

Becerra F., Marinacci F., Inayoshi K., Bromm V., Hernquist L. E., 2018, ApJ, 857, 138 Beers T. C., Christlieb N., 2005, ARA&A, 43, 531 Bonnor W. B., 1956, MNRAS, 116, 351 Brauer K., et al., 2025a, The Astrophysical Journal, 980, 41 Brauer K., et al., 2025b, The Astrophysical Journal, 993, 2 Bromm V., 2013, Reports on Progress in Physics, 76, 112901 Bromm ...

work page arXiv 2018
[2]

L., Bromm V., 2009, MNRAS, 392, L50 Fujimoto S., et al., 2025a, arXiv e-prints, p

pp 533–557, doi:10.1016/B978-0-443-21439-4.00052-3 Frebel A., Johnson J. L., Bromm V., 2009, MNRAS, 392, L50 Fujimoto S., et al., 2025a, arXiv e-prints, p. arXiv:2512.11790 Fujimoto S., et al., 2025b, ApJ, 989, 46 Gessey-Jones T., et al., 2022, Monthly Notices of the Royal Astronomical Society, 516, 841–860 Grassi T., Bovino S., Schleicher D. R. G., Priet...

work page doi:10.1016/b978-0-443-21439-4.00052-3 2009
[3]

(A13) in the left hand side, and substituting Eq

By writing ˙m2mαc(t)−˙m1mαc(t) = ˙m2mαc(t)−˙N,(A17) expanding the accretion law Eq. (A13) in the left hand side, and substituting Eq. (2.20), one can show that ˙m2 dN dm ⏐⏐⏐⏐ m− 2 = ˙M∗−¯m˙N m2−¯m,(A18) where ¯m= ∫m2 m1 dmm dN dm∫m2 m1 dm dN dm .(A19) This paper has been typeset from a TEX/LATEX file prepared by the author. MNRAS000, 1–21 (2026)

2026