Pith. sign in

REVIEW 3 major objections 3 minor 1 cited by

Once a protostar forms, the unaccreted gas lifetime scales with final stellar mass: low-mass stars accrete briefly from local gas, high-mass stars for millions of years from a larger volume.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-13 09:01 UTC pith:JFBI2GOX

load-bearing objection Useful STARFORGE Lagrangian tracking of mass-dependent accretion, but the post-selected cell definition makes the lifetime–mass trend partly tautological. the 3 major comments →

arxiv 2604.06471 v2 pith:JFBI2GOX submitted 2026-04-07 astro-ph.GA

The Evolution of Star-Forming Gas in STARFORGE: From Clouds, to Cores, to Stars

classification astro-ph.GA
keywords star formationgiant molecular cloudsdense coresprotostellar accretionturbulencemagnetic fieldsSTARFORGELagrangian tracking
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Star formation is known to happen in dense parts of giant molecular clouds, but it has been unclear how gas gathers into individual stars and what role dense cores play. This paper tracks the gas that ends up in stars by following the fluid elements in three cloud simulations that differ mainly in magnetic field strength. It finds that after a protostar appears, the remaining unaccreted gas lives for a time that rises with the star's final mass: stars below half a solar mass finish accreting in roughly half a million years from nearby gas, while stars above two solar masses keep accreting for three to five million years from a much larger volume. At the moment of formation that gas already obeys the same linewidth-size and mass-size relations seen in real dense cores, and its radii, velocities, and energy ratios barely depend on the cloud's global field. Low- and intermediate-mass stars accrete fairly steadily and can be matched by standard theoretical models, but many high-mass stars accrete in bursts that none of those models describe. The result matters because it shows that the gas reservoir is more extended than the cores observers usually catalogue, yet is still shaped by turbulence and feedback in ways that look like real cores.

Core claim

Once a protostar forms, the lifetime of its unaccreted gas correlates with final stellar mass: low-mass stars (M_* < 0.5 M_☉) accrete for 0.5–0.6 Myr from a relatively local reservoir, while high-mass stars (M_* > 2 M_☉) accrete over 3.3–4.7 Myr from a much larger volume. At formation the unaccreted gas follows turbulently regulated core scalings σ_v ∝ R^{0.47–0.55} and M ∝ R^{1.0–1.1}, and high-mass accretion histories are not well-fit by isothermal-sphere, turbulent-core, or competitive-accretion models.

What carries the argument

Lagrangian cell tracking of unaccreted gas in the STARFORGE radiation-magnetohydrodynamic cloud simulations: every fluid element that will later join a sink particle is followed from the moment of protostar formation, yielding the mass-dependent lifetime, radius, and accretion history of the true reservoir.

Load-bearing premise

That the fluid cells labelled as unaccreted gas in these three simulated clouds are a faithful stand-in for the real physical reservoir that feeds a protostar.

What would settle it

Direct comparison of measured accretion durations and reservoir sizes for a statistical sample of low-mass versus high-mass protostars in nearby star-forming regions; if the durations do not rise from ~0.5 Myr to several Myr with final mass, or if the gas does not obey the reported linewidth-size and mass-size relations at formation, the claim fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 3 minor

Summary. The manuscript uses Lagrangian cell tracking in three STARFORGE GMC simulations (varying magnetic field strength) to follow gas that eventually accretes onto sink particles. It reports that post-formation accretion duration and reservoir size scale with final stellar mass (low-mass stars accrete ~0.5–0.6 Myr from local gas; high-mass stars ~3.3–4.7 Myr from larger volumes), that at protostar formation the unaccreted gas obeys turbulently regulated core scalings (σ_v ∝ R^{0.47–0.55}, M ∝ R^{1.0–1.1}), that low/intermediate-mass accretion histories are reasonably fit by isothermal-sphere, turbulent-core, or competitive-accretion models while many high-mass histories are intermittent and poorly fit by any of them, and that radii, velocity dispersions, virial parameters, and magnetic energy ratios of the accreting gas are largely insensitive to global cloud B. The authors conclude that star-forming gas is more extended than classical dense cores yet has core-like physical properties regulated by turbulence and feedback.

Significance. If the mass-dependent lifetimes, core-like scalings at formation, and model-fit dichotomy survive scrutiny, the work would clarify how much of a star’s mass is set by a local core versus larger-scale, longer-lived accretion, and would provide a concrete simulation-based challenge to applying isothermal-sphere / turbulent-core / competitive-accretion prescriptions uniformly across the IMF. The use of the STARFORGE suite (full feedback, MHD, sink particles) and explicit comparison to analytic accretion models are genuine strengths. The result that many high-mass accretion histories are intermittent and poorly fit by standard models is a falsifiable, observationally relevant claim.

major comments (3)
  1. [Methods / Lagrangian tracking definition] The central mass–lifetime and mass–radius trends rest on defining the reservoir as the set of Lagrangian cells that are eventually accreted by each sink. By construction, the “lifetime of the unaccreted gas” is then essentially the accretion duration, and high-mass sinks that continue accreting for several Myr automatically sample larger, later-arriving volumes. The abstract itself notes that the accreting gas is “more extended than typically-defined dense cores.” Without controls that (i) identify a coherent, bound or contiguous reservoir at the instant of sink formation independent of final membership, (ii) compare to randomly selected or non-accreting cell sets of equal mass, or (iii) show that the same scalings appear when the reservoir is defined by density/velocity thresholds at formation time alone, the reported correlations risk being tautological rather than physical. This defin
  2. [Results / sample of high-mass stars] Only three GMC realizations (varying B) are used. High-mass stars (M_* > 2 M_⊙) are rare; the reported 3.3–4.7 Myr range and the statement that “many” high-mass histories are intermittent and poorly fit therefore rest on a small sample. The paper needs to report the number of sinks in each mass bin per run, the distribution of accretion durations (not only the quoted ranges), and whether the high-mass conclusions survive leave-one-out or bootstrap tests across the three clouds.
  3. [Methods / numerical resolution and sink rules] Sink accretion rules and finite resolution couple cell membership to numerical parameters. The abstract asserts that radii, σ_v, virial parameters, and magnetic energy ratios are “largely insensitive” to global B, yet provides no resolution or accretion-threshold convergence tests for the Lagrangian-tracked quantities. Without those tests, it remains possible that the mass-dependent lifetimes and the failure of analytic models for high-mass stars are simulation-specific rather than physical.
minor comments (3)
  1. [Abstract / Results] The abstract quotes power-law indices (0.47–0.55, 1.0–1.1) without uncertainties or the mass/radius range over which they are measured; these should be stated with fit errors and the fitting procedure.
  2. [Methods] Clarify whether “unaccreted gas” at formation includes only cells that will later be accreted by that specific sink, or also cells that remain unbound / are accreted by other sinks; the distinction matters for comparison to observed cores.
  3. [Results / accretion history fits] The claim that no single analytic model fits all masses is useful; a quantitative goodness-of-fit metric (e.g., reduced χ² or residual time series) per mass bin would make the “well-fit” / “not well-fit” statements reproducible.

Circularity Check

0 steps flagged

No significant circularity: empirical Lagrangian tracking results are measured outputs, not forced by construction or self-citation.

full rationale

This is a simulation-analysis paper that reports measured properties of gas tracked via Lagrangian cells in STARFORGE runs, not a first-principles derivation whose conclusions reduce to its inputs. The reservoir is defined as cells that are eventually accreted; the reported lifetimes, enclosing radii, velocity dispersions, virial parameters, and power-law indices (σ_v ∝ R^{0.47–0.55}, M ∝ R^{1.0–1.1}) are then measured on that set at and after protostar formation. Those quantities are free to take any values consistent with the hydrodynamics, magnetic fields, and feedback; nothing in the abstract or method forces the mass-dependent timescales (0.5–0.6 Myr vs 3.3–4.7 Myr), the larger volumes for high-mass stars, or the specific exponents. Comparison of accretion histories to external analytic models (isothermal sphere, turbulent core, competitive accretion) yields both successes (low/intermediate-mass) and failures (many high-mass stars), which is independent grounding rather than a fitted-input-called-prediction. Use of the authors’ own STARFORGE suite is ordinary self-citation of a numerical experiment and does not invoke uniqueness theorems or smuggled ansatzes. Post-selection of eventually-accreted cells raises legitimate interpretive caveats about whether the set constitutes a coherent physical core at formation, but that is a methodological limitation, not circular reduction of a claimed derivation. Score 0 is therefore appropriate.

Axiom & Free-Parameter Ledger

3 free parameters · 3 axioms · 0 invented entities

The central claims rest on the STARFORGE simulation physics and on the operational definition of unaccreted Lagrangian gas as the accretion reservoir. Free parameters include the three cloud magnetic-field setups and the mass bins used to quote lifetimes. Domain assumptions include that the simulated turbulence, feedback, and sink accretion capture the relevant real-cloud physics. No new physical entities are invented; the work analyzes existing simulation machinery.

free parameters (3)
  • GMC magnetic field strengths (three runs)
    The study design varies B across three GMCs; the specific field values and how they are normalized are free simulation choices that define the sample over which ‘insensitive to global cloud properties’ is claimed.
  • Stellar mass bins (M_* < 0.5 M_⊙ and M_* > 2 M_⊙)
    Lifetime ranges are quoted for these cuts; the bin edges are analysis choices that structure the main correlation claim.
  • Definition of unaccreted / accreting gas reservoir (Lagrangian cell set)
    Which cells count as the reservoir after protostar formation (density threshold, binding, membership over time) is an operational free choice that directly sets measured lifetimes, radii, and virial parameters.
axioms (3)
  • domain assumption STARFORGE GMC simulations with sink particles and included feedback adequately represent the mass assembly of real protostars.
    All lifetime and model-fit conclusions are drawn from these runs; the abstract treats them as the physical system under study.
  • domain assumption At protostar formation, unaccreted gas should be compared to turbulently regulated isothermal dense-core scalings (σ_v–R, M–R).
    The abstract uses those characteristic relations as the benchmark for ‘core-like’ structure.
  • domain assumption Isothermal-sphere, turbulent-core, and competitive-accretion models are the appropriate analytic baselines for continuous accretion histories.
    Model-fit success/failure statements depend on this choice of comparison set.

pith-pipeline@v1.1.0-grok45 · 6780 in / 3081 out tokens · 38172 ms · 2026-07-13T09:01:01.435239+00:00 · methodology

0 comments
read the original abstract

Star formation occurs within dense regions of giant molecular clouds (GMCs), however, exactly how gas collects and evolves to form individual stars and what role dense cores play remains unclear. We use the Lagrangian cell information in the STARFORGE simulation suite to track star-forming gas in three GMCs with varying magnetic field strengths. We find that, once a protostar forms, the lifetime of the unaccreted gas correlates with the final stellar mass, where low-mass stars ($M_*$ < 0.5 M$_\odot$) accrete for 0.5-0.6 Myr from a relatively local reservoir of gas, and high-mass stars ($M_*$ > 2 M$_\odot$) accrete over 3.3-4.7 Myr from a much larger volume. Although the protostellar accretion time increases weakly with magnetic field strength, the accreting gas radii, velocity dispersions, virial parameters, and magnetic energy ratios are largely insensitive to the global cloud properties. At the time of protostar formation, the unaccreted gas exhibits linewidth-size and mass-size relations characteristic of turbulently regulated, isothermal dense cores, following $\sigma_v \propto R^{0.47-0.55}$ and $M \propto R^{1.0-1.1}$, respectively. Low- and intermediate-mass stars undergo relatively continuous accretion and their accretion histories are well-fit by either isothermal sphere, turbulent core, or competitive accretion models, where no one model fits all masses. However, many high-mass stars experience intermittent accretion and their accretion histories are not well-fit by any of these models. While the distribution of accreting gas is more extended than typically-defined dense cores, the physical properties and structure of the star-forming gas resemble those of observed cores and are largely regulated by turbulence and feedback.

Figures

Figures reproduced from arXiv: 2604.06471 by Ananya Kaalva, Michael Y. Grudic, Nina Filippova, Stella S. R. Offner.

Figure 1
Figure 1. Figure 1: Gas forming typical individual low-mass and high-mass stars during the prestellar phase at three times in an M2e4 cloud. The top row represents the spatial distribution of gas forming a low-mass star with final mass 0.206 M⊙ and the bottom row is the gas for a high-mass star of 9.997 M⊙. Gray areas indicate the spatial density of the cloud gas, colored circles represent the star-forming gas, and the star s… view at source ↗
Figure 2
Figure 2. Figure 2: Gas forming typical individual low-mass and high-mass stars during the prestellar phase at three times in an M2e4 cloud. The top row represents the spatial distribution of gas forming a low-mass star with final mass 0.206 M⊙ and the bottom row is the gas for a high-mass star of 9.997 M⊙. Gray areas indicate the spatial density of the cloud gas, colored circles represent the star-forming gas, and the star s… view at source ↗
Figure 3
Figure 3. Figure 3: Time evolution of the prestellar gas mass function (PGMF) across the three STARFORGE simulations with varying magnetic field strengths (M2e4 mu4.2, M2e4, M2e4 mu0.4). Each panel shows the log-binned distributions of gas masses at five time snapshots, with line shading indicating time. Counts represent the number of prestellar gas subsets present per mass bin at each snapshot [PITH_FULL_IMAGE:figures/full_… view at source ↗
Figure 4
Figure 4. Figure 4: Time evolution of gas properties for the inter￾mediate mass bin (0.5 < M < 2M⊙), averaged across five bins sorted by protostellar duration. The three shades of lines correspond to the three simulation runs, with the darker lines indicating a stronger magnetic field. Panels show (top to bottom): gas mass, velocity dispersion, effective radius, virial parameter, BE/PE ratio, and KE/BE ratio. The ver￾tical do… view at source ↗
Figure 5
Figure 5. Figure 5: Average velocity dispersion and magnetic-to– gravitational energy ratio (BE/PE) across all cores versus time. Solid lines indicate the mean value at each time, and the shaded regions show ±1σ standard deviation. The three magnetic field strengths (M2e4 mu4.2, M2e4, M2e4 mu0.4) are indicated by the different colors. the strong field run, M2e4 mu0.4, as magnetic flux is removed from the gas. The KE/PE ratios… view at source ↗
Figure 6
Figure 6. Figure 6: Log velocity dispersion (top) and log mass (bot￾tom) versus log effective radius at the onset of accretion (start of the protostellar phase) for all stars in each sim￾ulation. In the top panel the points are colored by mass. Solid lines represent the best-fit linear regressions for each simulation, colored by magnetic field strength. The median value is marked in black. The fit parameters are reported in … view at source ↗
Figure 7
Figure 7. Figure 7 [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Distribution of j = 1−1/b parameter values for low- (M < 0.5M⊙), intermediate- (0.5 < M < 2M⊙), and high-mass (M > 2M⊙) stars across the three clouds. A χ 2 cutoff of 3.0 was applied to exclude poorly fit mass accretion histories, and the resulting distributions peak near j ≈ 0.5. be dispersed by stellar feedback (e.g., C. D. Matzner & C. F. McKee 2000; M. N. Machida & T. Hosokawa 2013; S. S. R. Offner & H… view at source ↗
Figure 9
Figure 9. Figure 9: Same as [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Top panels: Same as [PITH_FULL_IMAGE:figures/full_fig_p017_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Top panels: Same as [PITH_FULL_IMAGE:figures/full_fig_p018_11.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

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

  1. Preferential alignment of Class 0, Class I protostellar disks in multiple systems across nine nearby molecular clouds

    astro-ph.SR 2026-07 conditional novelty 6.0

    Disks around young stars in binary and higher-order multiple systems are preferentially aligned out to 6000 AU, implying turbulent fragmentation alone cannot explain how most multiples form.

Reference graph

Works this paper leans on

3 extracted references · 1 linked inside Pith · cited by 1 Pith paper

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

  3. [3]

    Episodic Accretion in Young Stars

    thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...