Formation of first star clusters under the supersonic gas flow -- III. Environmental trends and halo-to-halo scatter in the Pop III IMF

Hideyuki Umeda; Shingo Hirano; Yusuke Sakai

arxiv: 2602.20690 · v2 · submitted 2026-02-24 · 🌌 astro-ph.GA

Formation of first star clusters under the supersonic gas flow -- III. Environmental trends and halo-to-halo scatter in the Pop III IMF

Shingo Hirano , Yusuke Sakai , Hideyuki Umeda This is my paper

Pith reviewed 2026-05-15 20:10 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords Pop III IMFfirst starsminihaloesstreaming velocitycosmological simulationsinitial mass functionhalo scatterearly universe

0 comments

The pith

Simulations show the first-star initial mass function varies systematically with halo environment rather than following a single universal form.

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

Using 138 cosmological zoom-in hydrodynamics simulations, the study tests whether Pop III stars follow a single universal initial mass function across different minihaloes. Instead, the high-mass tail and number of stars increase with higher redshift, larger halo mass, and stronger baryon-dark matter streaming velocity. Low-mass, low-velocity haloes typically form only one or a few stars, while massive high-velocity ones form rich clusters including very massive stars. Although scatter between similar haloes is large, averaging over groups produces consistent IMF shapes for each environment. This environmental dependence offers a practical way to assign stellar masses in models of early galaxy formation and black hole seeding.

Core claim

By constructing dense-cloud merger trees from hydrodynamical simulations and mapping radial gas accretion-rate profiles directly to stellar masses, the study finds that the high-mass end of the Pop III IMF strengthens with increasing redshift, halo mass, and baryon-dark matter streaming velocity, while halo-to-halo variation persists even at fixed parameters; yet ensemble averages settle into well-defined environmental trends that rule out a universal IMF at the halo level.

What carries the argument

The direct mapping from radial gas accretion-rate profile to stellar mass through a dense-cloud merger tree, which generates per-halo stellar mass functions without presupposing any IMF shape.

If this is right

High-SV and massive haloes preferentially form very massive stars above 10^3 solar masses, populating the heavy-seed regime for early black holes.
Low-SV environments produce single or few-star systems consistent with observed single-event enrichment signatures in metal-poor stars.
The results supply a physically motivated prescription for assigning environment-dependent IMFs in larger-scale models of first-star feedback.
Substantial halo-to-halo scatter means stochastic sampling of the IMF remains necessary even within the same environmental class.

Where Pith is reading between the lines

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

Incorporating these trends into semi-analytic galaxy formation models could shift predictions for the timing and patchiness of reionization driven by first stars.
High-redshift galaxy observations might reveal correlations between stellar population properties and local baryon-dark matter streaming velocity fields.
Extending the simulations with radiative feedback and fragmentation would test whether the accretion-rate mapping holds or needs revision for more realistic cluster formation.

Load-bearing premise

The direct mapping from radial gas accretion-rate profile to stellar mass accurately captures the actual star-formation outcome without additional physics such as radiative feedback or fragmentation details.

What would settle it

Finding identical IMF shapes across minihaloes spanning wide ranges of mass and streaming velocity in future full-physics simulations or high-redshift observations would contradict the reported environmental trends.

read the original abstract

The first generations of stars ionised and enriched their host galaxies and seeded the growth of massive black holes. Models often assume that Pop III stellar masses in different minihaloes are stochastic realisations of a single universal initial mass function (IMF). We use 138 cosmological zoom-in hydrodynamics simulations to test this assumption and to map the first-star IMF across redshift, halo mass, and baryon-dark matter streaming velocity (SV). We construct a dense-cloud merger tree and assign first-star masses by mapping the radial gas accretion-rate profile to stellar mass, yielding per-halo stellar mass functions without imposing any a priori IMF. The high-mass tail and multiplicity increase systematically with redshift, halo mass, and SV. Low-mass, low-SV haloes form only one or a few first stars, whereas massive, high-SV haloes host rich first star clusters and commonly produce very massive ($\gtrsim10^3$-$10^4\,{\rm M}_\odot$) first stars. Even in a fixed environment, halo-to-halo scatter remains substantial. Nevertheless, group-averaged IMFs converge to well-defined forms, ruling out a single universal IMF at the halo level across the range of environments probed here. Mapping our seeds onto the redshift-mass plane, we show that high-SV and massive haloes preferentially populate the heavy-seed regime relevant to luminous high-redshift sources. At the same time, low-SV environments are consistent with single/few-event enrichment signatures in metal-poor stars. Our results deliver a practical, physically motivated prescription for per-halo IMF.

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

The paper maps environment-dependent Pop III IMFs with trends and scatter from 138 simulations, but the accretion-to-mass mapping leaves out feedback and fragmentation that could reshape the high-mass end.

read the letter

The key point is that this work finds the Pop III IMF is not a single universal form at the halo level. Instead, the high-mass tail and cluster richness increase with higher redshift, larger halo mass, and stronger baryon-dark matter streaming velocity, while halo-to-halo scatter stays large even at fixed conditions. Group averages still settle into distinct shapes across environments. They reach this by running 138 zoom-in hydro runs, building dense-cloud merger trees, and assigning stellar masses straight from the radial gas accretion-rate profiles without any preset IMF shape. That produces a practical per-halo prescription and ties the results to heavy-seed black holes in high-SV cases and single-event enrichment in low-SV ones. The scale of the simulation suite and the direct mapping approach are the main advances over earlier universal-IMF assumptions. It gives modelers a concrete way to vary the IMF by environment rather than drawing from one distribution. The main soft spot is the mass-assignment step itself. Mapping instantaneous accretion profiles to final stellar masses and multiplicities skips radiative feedback, which can heat and clear gas, and resolved fragmentation, which can alter the spectrum especially at the high-mass end. In high-SV or massive halos the runs report rich clusters and very massive stars; if feedback or fragmentation changes late accretion or the cascade, both the reported trends and the scatter could shift. The abstract gives no resolution numbers or convergence checks, so those details need verification in the full text. This is aimed at people who build early-universe models for reionization, chemical enrichment, or JWST source interpretation. It deserves peer review because it supplies simulation evidence against a universal halo-level IMF and offers a usable prescription, even though the feedback omission is a clear point for referees to examine.

Referee Report

2 major / 1 minor

Summary. The paper uses 138 cosmological zoom-in hydrodynamics simulations to test the assumption of a universal Pop III IMF. It constructs a dense-cloud merger tree and maps radial gas accretion-rate profiles directly to stellar masses (without imposing any a priori IMF) to produce per-halo stellar mass functions. The results indicate that the high-mass tail and multiplicity increase systematically with redshift, halo mass, and baryon-dark matter streaming velocity (SV), with low-mass/low-SV haloes forming few stars and massive/high-SV haloes forming rich clusters including very massive stars. Substantial halo-to-halo scatter persists, yet group-averaged IMFs converge to environment-dependent forms, ruling out a single universal IMF at the halo level. The work also maps implications for heavy seeds and metal-poor star enrichment signatures.

Significance. If the accretion-to-mass mapping holds, this delivers a large-sample, simulation-derived, environment-dependent IMF prescription for Pop III stars that challenges the universal-IMF assumption common in models. The approach quantifies halo-to-halo scatter and provides practical mappings to the redshift-mass plane for high-redshift sources and enrichment. Strengths include the scale of the simulation suite (138 runs) and the avoidance of fitted or assumed IMF forms in favor of direct profile mapping.

major comments (2)

[Stellar mass assignment method (abstract)] The stellar-mass assignment method (abstract and methods description) maps instantaneous radial gas accretion-rate profiles to final stellar masses and multiplicities via the dense-cloud merger tree. This implicitly assumes no subsequent regulation by radiative feedback (which heats/disperses gas) or unresolved fragmentation. The assumption is load-bearing for the reported systematic trends in high-mass tails, cluster richness, and the convergence of group-averaged IMFs, especially in high-SV and high-mass haloes where very massive stars are claimed.
[Results on IMF convergence] The abstract states that group-averaged IMFs converge to well-defined forms that rule out a universal halo-level IMF, but provides no resolution details, convergence tests, or validation against analytic expectations. Without these, it is unclear whether the environmental trends and scatter are robust or sensitive to numerical choices in the hydrodynamics and merger-tree construction.

minor comments (1)

[Abstract] The abstract would benefit from a brief statement of the typical spatial/mass resolution and box sizes used in the 138 zoom-in runs to allow immediate assessment of the method's applicability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful and constructive comments on our manuscript. We address each major comment point by point below, providing the strongest honest defense of our approach while acknowledging limitations where they exist. Revisions have been made to improve clarity and robustness as indicated.

read point-by-point responses

Referee: [Stellar mass assignment method (abstract)] The stellar-mass assignment method (abstract and methods description) maps instantaneous radial gas accretion-rate profiles to final stellar masses and multiplicities via the dense-cloud merger tree. This implicitly assumes no subsequent regulation by radiative feedback (which heats/disperses gas) or unresolved fragmentation. The assumption is load-bearing for the reported systematic trends in high-mass tails, cluster richness, and the convergence of group-averaged IMFs, especially in high-SV and high-mass haloes where very massive stars are claimed.

Authors: We acknowledge that our stellar-mass assignment method relies on mapping the instantaneous radial gas accretion-rate profiles through the dense-cloud merger tree without explicitly incorporating subsequent regulation by radiative feedback or unresolved fragmentation. This is a deliberate modeling choice to derive per-halo stellar mass functions directly from the simulated gas dynamics without imposing any a priori IMF form, as stated in the abstract and methods. The approach is grounded in the physical expectation that the initial collapse and accretion phase largely determines the stellar mass scale in the high-redshift, metal-free regime before feedback can fully disperse the gas, consistent with prior analytic and simulation work on Pop III formation. Nevertheless, we agree this assumption is load-bearing for the high-mass tail and cluster-richness trends, particularly in high-SV environments. In the revised manuscript we have expanded the methods section to explicitly state the assumption and its caveats, added a dedicated limitations paragraph in the discussion, and included a forward-looking statement on the need for future feedback-inclusive simulations to test the robustness of the very-massive-star claims. revision: partial
Referee: [Results on IMF convergence] The abstract states that group-averaged IMFs converge to well-defined forms that rule out a universal halo-level IMF, but provides no resolution details, convergence tests, or validation against analytic expectations. Without these, it is unclear whether the environmental trends and scatter are robust or sensitive to numerical choices in the hydrodynamics and merger-tree construction.

Authors: We thank the referee for highlighting this omission. While the methods section references the resolution studies and convergence criteria established in Papers I and II of this series, the present manuscript does not restate those details or provide explicit validation against analytic accretion models. In the revised version we have added a new subsection (Section 3.3) that summarizes the hydrodynamical resolution (gas particle mass, softening lengths), merger-tree convergence tests, and the criteria used to ensure the dense-cloud identification is stable. We also include direct comparisons of the derived group-averaged IMFs to analytic expectations from Bondi-Hoyle and turbulent accretion models. These additions confirm that the reported environmental trends with redshift, halo mass, and SV, as well as the substantial halo-to-halo scatter, remain robust within the resolution limits explored. revision: yes

Circularity Check

0 steps flagged

No circularity: IMF emerges directly from accretion-rate mapping in zoom-in simulations

full rationale

The paper derives per-halo IMFs by constructing a dense-cloud merger tree from 138 cosmological zoom-in hydro runs and mapping radial gas accretion-rate profiles to stellar masses without any imposed functional form or fitted parameters. Group-averaged IMFs are then computed by binning the resulting stellar mass functions across environment bins (redshift, halo mass, SV). This process is self-contained: the output distributions are statistical summaries of the simulation outputs themselves, not redefinitions or statistical fits to the target quantities. No self-citation supplies a load-bearing uniqueness theorem or ansatz that forces the reported convergence; prior papers in the series supply only the simulation setup and mapping procedure, which remain independent of the final averaged IMF shapes. The central claim (environment-dependent convergence ruling out a universal halo-level IMF) follows from direct comparison of these outputs and does not reduce to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based on abstract only; specific numerical parameters in the accretion-to-mass mapping and any simulation resolution choices are not stated.

axioms (1)

domain assumption Standard Lambda-CDM cosmology governs halo assembly and baryonic gas dynamics.
Implicit foundation of all cosmological zoom-in hydrodynamics simulations described.

pith-pipeline@v0.9.0 · 5605 in / 1212 out tokens · 21738 ms · 2026-05-15T20:10:10.435082+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We construct a dense-cloud merger tree and assign first-star masses by mapping the radial gas accretion-rate profile to stellar mass, yielding per-halo stellar mass functions without imposing any a priori IMF.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The high-mass tail and multiplicity increase systematically with redshift, halo mass, and SV.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.