{"id":"c195af6b-a9b2-4a14-9797-a46ddc9a11dc","arxiv_id":"2507.16944","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A framework linking protostellar evolution tracks to radiative-transfer SED grids predicts that isothermal-sphere, turbulent-core, and competitive accretion histories produce distinct 100-micron and 3-millimeter flux signatures.","lead":"This paper builds a way to predict the infrared and millimeter light that forming stars should emit, by connecting two existing families of models: one for how the star grows and one for the surrounding dust. The authors show that three competing theories of how stars gain mass predict different apparent brightness over time, which could let astronomers distinguish the theories with telescopes.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Distinguishability of accretion histories is asserted without a quantitative test against the framework's own flux uncertainties.","rationale":"The reader's weakest assumption focuses on the fidelity of the M10 accretion prescriptions (Section 2.2), which is indeed a serious limitation. My concern is related but distinct: even granting those prescriptions, the paper does not demonstrate that the resulting flux tracks are separable given the framework's own quantified uncertainties. Section 3.2 reports ~20–30% accuracy and ~50% neighbor scatter at 1 mm, yet Section 3.1 claims a 'clear distinction' based on visual inspection of Figure 3 without any statistical separation metric or classification test. This is the load-bearing gap because the abstract and Section 3.1 explicitly promise observational distinguishability, and that promise is the paper's central contribution over prior work. The proposed Monte Carlo classifier test would directly settle whether the predicted differences are larger than the noise. Since this is a call for additional analysis rather than a demonstration of an error, the existing CONDITIONAL verdict remains appropriate; I therefore recommend UNCHANGED.","tokens_in":37190,"tokens_out":5715,"duration_ms":69934,"concrete_test":"For a grid of final masses and ages, draw synthetic fluxes from the framework's predicted S100µm and S3mm tracks for IS, TC, and CA, adding noise drawn from the empirically measured recovery-error distribution (Section 3.2, Figure 5) at each wavelength. Then run a simple classifier (e.g., quadratic discriminant or nearest-centroid) on the noised fluxes to see whether accretion history can be recovered at greater than chance. If classification accuracy is near chance for realistic noise levels, the central distinguishability claim fails; if it is high, the claim is supported. A complementary check: compute the separation (e.g., Bhattacharyya distance) between the three history distributions in the S100–S3mm plane and compare it to the within-history scatter.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that IS, TC, and CA accretion histories can be observationally distinguished via 100-µm and 3-mm fluxes—is not supported by a quantitative separation test. Figure 3 shows tracks that appear distinct, but the framework's own error budget (Section 3.2) gives a 1-mm flux recovery accurate to only ~20–30% (16th to 84th percentiles) with a neighbor-scatter σ_MAD averaging ~50% of the flux. The paper never computes whether the track-to-track separations in the S100µm–S3mm plane exceed these uncertainties for realistic observed populations, nor how inclination, disk fraction, and distance affect the separation. Because the three histories also differ primarily in accretion timescale (Section 3.1), the claimed distinguishability hinges on the assumption that these timescales and the resulting tracks are robust; the authors acknowledge (Section 2.2) that the M10 prescriptions are steady, nonepisodic, and not fully consistent with modern star-formation understanding. Thus the headline result is an untested prediction rather than a demonstrated observable distinction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a framework that links the R24 radiative-transfer YSO model grid to protostellar evolutionary tracks generated with a modified Klassen et al. (2012) code, for isothermal-sphere (IS), turbulent-core (TC), and competitive accretion (CA) histories. The framework predicts flux evolution at 100 um and 3 mm for stellar masses from 0.2 to 50 solar masses, claims that these accretion histories may be observationally distinguished, characterizes prediction uncertainty through leave-one-out recovery, compares the framework against the Zhang & Tan (2018) grid to attribute a systematic offset to dust-model differences, and constructs Class/Stage confusion matrices for interpreting observed YSO populations. The paper also critically evaluates alternative Class and Stage definitions.","tokens_in":37476,"tokens_out":3135,"duration_ms":37831,"significance":"If the central distinguishability claim is established, the framework would be a valuable theory-agnostic bridge between protostellar evolution models and direct observables over a much wider mass range than most prior work. The paper has notable strengths: it ships public confusion-matrix data, it performs a careful leave-one-out validation against its own grid, it quantifies wavelength-dependent accuracy and precision honestly, and it tests the impact of dust opacities by rerunning a subset of models with the ZT18 dust configuration. These elements make the methodology reusable and the caveats transparent. However, the headline claim that IS, TC, and CA histories can be observationally distinguished is not supported by a quantitative test against the framework's own flux uncertainties, and the underlying accretion models are acknowledged by the authors to be simplified in ways that directly affect the claimed distinguishing timescales.","major_comments":[{"comment":"The central claim that IS, TC, and CA accretion histories can be observationally distinguished via 100-um and 3-mm fluxes is asserted from visual separation of tracks in Figure 3, but the paper never performs a quantitative separation test against its own uncertainty budget. Section 3.2 and Figure 5 show that the 16th-84th percentile recovery ratio at 1 mm is 0.83-1.28 (about 20-30% spread), and Figure 4 reports a mean fractional sigma_MAD of approximately 0.5. The tracks in Figure 3 are also smoothed with a rolling median, which can reduce apparent scatter. The paper should compute, for representative masses and ages, whether the inter-history separations in the S100um-S3mm plane exceed the propagated neighbor-scatter uncertainties, and how inclination, distance, and population mixture affect the distinguishability. Without such a test, the statement in Section 3.1 that 'it is theoretically possible to distinguish the mechanism of accretion at play through observation' remains an untested prediction rather than a demonstrated result.","section":"Section 3.1 and Figure 3"},{"comment":"The distinguishability claim depends on the accretion timescales predicted by the three M10 prescriptions, but the authors themselves state that 'the IS/TC/CA models utilized here are therefore not, as implemented, fully consistent with the modern understanding of star formation.' In particular, the implementation assumes steady, nonepisodic accretion, monolithic collapse, and a finite mass reservoir, and the CA model is implemented under monolithic-collapse assumptions that are at odds with the hierarchical nature of competitive accretion. Because the flux tracks in Figure 3 differ primarily in the timescale of mass assembly, these acknowledged simplifications are load-bearing for the distinguishability claim. The paper should either test the sensitivity of the predicted tracks to episodic or tapered accretion and to reservoir replenishment, or explicitly reframe the results as a proof-of-concept under the stated simplified assumptions rather than as a robust observational discriminant.","section":"Section 2.2"}],"minor_comments":[{"comment":"In the first paragraph, 'following the the Offner et al. (2009) implementation' contains a duplicated article and should read 'following the Offner et al. (2009) implementation.'","section":"Section 2.2"},{"comment":"The paragraph beginning 'Overestimates in the flux recovery distribution are thick disks —' is a sentence fragment and reads as an incomplete heading; it should be rewritten as a complete sentence or integrated into the surrounding text.","section":"Section 3.2"},{"comment":"The caption contains the typo 'efficiency' (non-standard ligature and spelling); it should be 'efficiency'.","section":"Figure 13 caption"},{"comment":"The abstract states that the histories 'may be observationally distinguished,' while Section 3.1 states that it is 'theoretically possible to distinguish' them; the conclusion then states that the framework 'found' observable differences. These phrasings should be made consistent with the level of quantitative support actually provided, especially given the lack of a separation test.","section":"Abstract and Section 3.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid methodology contribution with honest uncertainty quantification and useful public products. My main concern is that the abstract and conclusion overstate the distinguishability result relative to what is demonstrated; the quantitative separation test suggested in the major comments should be feasible with the existing machinery and would either strengthen or appropriately temper the headline claim. The acknowledged simplifications in the accretion prescriptions also warrant a sensitivity test or a clear reframing. These issues are addressable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the paper does something genuinely new: it links a theory-agnostic radiative-transfer grid (R24) to protostellar evolutionary tracks (modified K12) via nearest-neighbor matching in T*, L*, Mcore space, so you can predict fluxes for arbitrary accretion histories without building a bespoke grid. Second, the headline claim that IS/TC/CA histories can be distinguished with 100 µm and 3 mm fluxes is not actually demonstrated; it's an interesting prediction, but the paper never checks whether the track separations exceed its own substantial flux uncertainties.\n\nWhat's good: The internal leave-one-out validation is honest and useful—long-wavelength recovery within ~20-30%, with a clear caveat that IR fluxes are much worse. The comparison to ZT18 is careful, and the dust-model rerun experiment that explains the 30% offset is a nice piece of detective work. The confusion matrices between Class and Stage are a practical tool, and the release of those matrices is a plus. The authors also state their limitations plainly: the M10 prescriptions are steady, nonepisodic, and not fully consistent with modern star formation, so the tracks are simplified toys.\n\nThe soft spots: The central distinguishability claim is the load-bearing result, and it's untested. The paper's own σ_MAD at 1 mm averages about 50% of the flux, and the recovery scatter is 20-30%; the tracks in Figure 3 look separated, but no test shows that a realistic population with inclinations, distances, and disk fractions would land in distinct regions. That's an addressable gap, not a fatal one—the framework can produce the tracks; the authors just need to run the separation test. Also, the framework code isn't released, only the confusion matrices, which limits reproducibility. And the Stage II phase is underrepresented because accretion stops at envelope depletion, which the authors note.\n\nWho this is for: YSO observers who want to interpret photometry without assuming an accretion history, and modelers who want to compare theories directly to fluxes. It deserves a serious referee—the method is novel and the execution is mostly careful—but it needs a revision that quantifies whether the claimed distinguishability survives the error budget.\n\nBottom line: send it out, but ask for the separation test and the code release before publication.","headline":"A solid, genuinely new framework that overreaches slightly in its headline claim about distinguishing accretion histories—the separation isn't tested against the paper's own error bars.","tokens_in":37957,"tokens_out":2677,"would_cite":true,"duration_ms":29510,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that the way a young star accretes its mass leaves a distinguishable imprint in its 100-micron and 3-millimeter fluxes, so surveys can test competing star-formation theories directly.","keywords":["young stellar objects","protostellar accretion","radiative transfer","spectral energy distributions","star formation histories","Class and Stage classification","flux predictions"],"falsifier":"Take a nearby embedded cluster with a well-determined age, distance, and membership, measure each YSO at 100 $\\mu$m and 3 mm, and plot the sources against the predicted flux tracks and isochrones for the appropriate age and mass. If the observed population does not fall into the predicted IS/TC/CA patterns, or if a single cluster's sources scatter across the whole plane rather than following one history, the claim that the accretion mechanism can be distinguished this way would be falsified.","tokens_in":36990,"feed_emoji":"🌟","tokens_out":8480,"duration_ms":85441,"temperature":0.7,"pith_summary":"This paper develops a framework for predicting how a young stellar object's observable flux evolves as it accretes mass, without committing in advance to any one theory of star formation. The framework joins a large, theory-agnostic set of radiative-transfer spectral energy distribution (SED) models to protostellar evolutionary tracks, so that any assumed accretion history can be turned into predicted fluxes. Applying it to isothermal-sphere, turbulent-core, and competitive accretion histories for birth masses from 0.2 to 50 solar masses, the authors argue that the three histories leave distinguishable tracks in the 100-micron versus 3-millimeter flux plane, primarily because they assemble mass on different timescales. If correct, this gives observers a direct, multi-wavelength test of which accretion mechanism actually operates, and it provides physically grounded confusion matrices for translating observed YSO classes into evolutionary stages.","feed_headline":"Three star-formation histories leave distinct far-IR and mm tracks","feed_subtitle":"Predicting 100-micron and 3-mm fluxes from theory lets surveys test how young stars actually accrete.","key_machinery":"The carrying mechanism is a three-parameter nearest-neighbor matching procedure. Each protostellar evolutionary track is converted into a sequence of points in the space of source temperature $T_\\star$, total luminosity $L_\\star$, and circumstellar mass $M_{\\rm core}$ within roughly 10,000 au; a quantile transform maps each coordinate to a uniform distribution so that distances in different physical units are comparable, and the ten nearest radiative-transfer models are averaged as the median SED separately by inclination to produce a predicted SED at each time step. The accretion-rate parameterization behind the tracks has the common form $\\dot{m} = \\dot{m}_1 (m/m_f)^j m_f^{j_f}$, with history-specific scaling parameters, so the same machinery can host any accretion prescription. The 100-micron versus 3-mm plane is the diagnostic output because the former traces the luminosity that heats the dust while the latter traces the mass of optically thin dust.","core_discovery":"The central claim is that the mechanism by which a protostar gains mass can be read off from long-wavelength observations of the forming star. By matching evolutionary tracks that follow the standard isothermal-sphere, turbulent-core, and competitive accretion rate laws to a broad grid of radiative-transfer SED templates in the space of stellar temperature, total luminosity, and circumstellar mass, the framework produces, for each history and for final stellar masses between 0.2 and 50 $M_\\odot$, a time-ordered sequence of 100-micron and 3-mm fluxes. The three histories trace similar shapes in this plane but advance along it at very different speeds: low-mass stars deplete their dust reservoirs on timescales that differ by factors of two to three, while high-mass stars take several million years to accrete under isothermal-sphere conditions but only a few hundred thousand years under turbulent-core or competitive accretion. Because 100-micron flux tracks luminosity and 3-mm flux tracks dust mass, the authors argue that these timescale differences translate into observationally distinguishable flux tracks, and that a mixed population using different accretion modes in different mass regimes would appear distinct as well.","pith_inferences":["My inference: because the distinguishing power comes from accretion timescales, adding episodic or tapered accretion would blur or shift the predicted tracks; a synthetic-population test mixing histories would show how much separation survives real variability.","My inference: the same quantile-transform matching could attach other evolutionary tracks, such as externally fed or episodic accretion models, to the existing SED grid, an extension the authors note is deferred to future work.","My inference: the confusion matrices imply that published Class 0/I lifetimes inferred from Class counts may be systematically misestimated if the true population is not isothermal-sphere-like; the paper's c2d example points in this direction but does not make it a final conclusion.","My inference: the strong inclination dependence seen in reproducing an independent grid suggests that edge-on, disk-dominated YSOs may be systematically overbright at millimeter wavelengths, so flux-plane diagnostics should be applied to populations or with inclination control rather than to individual sources."],"forward_implications":["Survey data at 100 $\\mu$m and 3 mm can, in principle, distinguish whether a YSO population accretes via isothermal-sphere, turbulent-core, or competitive accretion, with the largest separation for massive stars whose accretion timescales differ by roughly an order of magnitude.","Low-mass stars of 0.2 $M_\\odot$ are predicted to be visible at long wavelengths for differing durations of about 0.13 Myr (IS), 0.25 Myr (TC), and 0.36 Myr (CA), so the same flux measurement implies different ages under different histories.","The Stage 0/I boundary defined by half-mass assembly corresponds to a knee in the predicted flux tracks, offering a photometric way to separate the earliest evolutionary phases rather than relying solely on Class.","Class-to-Stage confusion matrices let observers convert observed Class counts into Stage counts under each accretion history; applying them to a c2d-like sample shifts the inferred population strongly toward Stage 0/I.","The framework's own 1-mm flux recovery is accurate to roughly 20-30% on average, and the systematic offset found when reproducing an independent turbulent-core grid can be explained by differences in dust opacity and disk structure."],"supporting_citations":[{"why":"It supplies the theory-agnostic radiative-transfer YSO model grid and the SEDs the framework selects from.","marker":"R24"},{"why":"It provides the protostellar evolution code whose tracks are modified to generate isothermal-sphere, turbulent-core, and competitive accretion histories.","marker":"K12"},{"why":"It prescribes the accretion-rate parameterization used for all three accretion histories in the generated tracks.","marker":"M10"},{"why":"It defines the isothermal-sphere accretion model that serves as one of the three baseline scenarios.","marker":"Shu 1977"},{"why":"It defines the turbulent-core accretion model used as the second baseline scenario.","marker":"McKee & Tan 2002"},{"why":"It introduces the competitive accretion scenario used as the third baseline history.","marker":"Bonnell et al. 1997"},{"why":"It provides the original model grid and Hyperion-based radiative transfer methods that R24 extends.","marker":"R17"},{"why":"It is the independent turbulent-core model grid used to test whether the framework reproduces externally built predictions.","marker":"ZT18"}],"fun_headline_variants":["Accretion history leaves distinct 100-μm and 3-mm tracks","Far-IR and mm fluxes betray how young stars accrete","Three star-formation histories show separate far-IR and mm tracks","Long-wavelength light reveals the accretion recipe of forming stars","How a protostar gains mass is readable in its 100-μm and 3-mm glow"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the standard steady, non-episodic accretion prescriptions used here, combined with the assumption that each star forms from a single collapsing core with a finite mass reservoir, describe the observable behavior of real protostars well enough for the predicted flux tracks to be meaningful; if real accretion is episodic or slows over time, the tracks and their distinguishability could change.","fun_headline_variants_meta":{"raw":{"variants":["Accretion history leaves distinct 100-μm and 3-mm tracks","Far-IR and mm fluxes betray how young stars accrete","Three star-formation histories show separate far-IR and mm tracks","Long-wavelength light reveals the accretion recipe of forming stars","How a protostar gains mass is readable in its 100-μm and 3-mm glow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000726,"raw_usage":{"total_tokens":3320,"prompt_tokens":1077,"completion_tokens":2243,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":693,"completion_tokens_details":{"reasoning_tokens":2147}},"tokens_in":693,"tokens_out":2243,"duration_ms":17785,"temperature":1.0,"reasoning_tokens":2147,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:00:32.470839+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a nearby embedded cluster with a well-determined age, distance, and membership, measure each YSO at 100 $\\mu$m and 3 mm, and plot the sources against the predicted flux tracks and isochrones for the appropriate age and mass. If the observed population does not fall into the predicted IS/TC/CA patterns, or if a single cluster's sources scatter across the whole plane rather than following one history, the claim that the accretion mechanism can be distinguished this way would be falsified.","supporting_citations":[],"review_version":1}