{"id":"6054a161-3280-465f-aa0a-ad1af45c47a6","arxiv_id":"2505.07764","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A measured correlation between most massive core mass and bound gas mass, combined with an assumed 30% efficiency and a theoretical SFR-most massive star relation, reproduces the approximately linear Gao-Solomon law.","lead":"This paper reports a new correlation between the mass of the most massive dense core in a molecular cloud and the cloud's gravitationally bound gas mass, spanning three orders of magnitude. The authors argue that this correlation, combined with two previously proposed relations, can reproduce the linear star formation law known as the Gao-Solomon relation, suggesting a more deterministic picture of star formation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central comparison assumes Mbound_gas ≈ Mdense_gas (Footnote 3, deferred to a companion paper); if wrong, Figure 3 does not test Gao-Solomon. Eq. 8's support is a self-cited theory plus consistency generated through the same chain.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the paper's central comparison substitutes Mbound_gas for Mdense_gas without quantitative demonstration, deferring to an unavailable companion paper. I do not find an internal inconsistency in the derivation; the Mmax_core–Mbound_gas correlation is a genuine empirical contribution and the random-sampling test adds falsifiable support. However, the headline explanation of the Gao-Solomon relation depends on that substitution, and the associated circularity in treating Eq. 8 as 'empirically constrained by the presented data' reinforces the need for independent validation. This does not justify rejection, because the paper is explicit about its conditional nature and the new correlation stands on its own. It does justify keeping a conditional verdict until the companion comparison and an independent test of Eq. 8 are available. Hence the reader's verdict is unchanged.","tokens_in":24729,"tokens_out":9473,"duration_ms":98593,"concrete_test":"For each of the 23 regions, measure Mdense_gas from an independent dense-gas tracer (e.g., HCN J=1-0 or N2H+ from public surveys, applying Gao-Solomon's definition) and test log(Mbound_gas/Mdense_gas) for zero mean, scatter <0.3 dex, and no trend with Mbound or SFR. Recompute Figure 3 replacing Mbound with Mdense; if the resulting slope or normalization moves by more than ~30% or no longer overlaps the Lada-rescaled Gao-Solomon line, the central claim fails. This directly resolves the dependency on Jiao et al. (submitted).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The Gao-Solomon relation is defined for Mdense_gas. Figure 3 instead plots SFRtheory against Mbound_gas and overlays the Gao-Solomon line, relying on the equivalence asserted in Footnote 3: 'in many but not all star-forming regions, Mbound_gas ~ Mdense_gas' (Jiao et al. submitted). This is an unquantified claim in an unavailable companion paper. If Mbound_gas/Mdense_gas varies with cloud mass or with SFR, the agreement in Figure 3 is not a test of the Gao-Solomon relation; it is an overlay against an arbitrary x-variable. The normalization and slope of the derived relation are sensitive to this substitution. Additionally, Eq. 8 is taken from Yan et al. 2017 by fitting their theoretical model; the paper then suggests in §4.1 that Eq. 8 'can also be regarded as empirically constrained by the presented data,' but the presented data constrain Eq. 8 only through Eq. 7, the assumed 30% SFE, and the same Mbound≈Mdense substitution. That makes part of the consistency check circular. The Mmax_core–Mbound_gas correlation itself appears real and is robust to temperature method, but it is not by itself a measurement of the dense-gas mass used in Gao-Solomon.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an explanation for the linear Gao-Solomon star formation law (SFR proportional to dense gas mass) using a chain of three correlations: SFR vs. the most massive star, the most massive star vs. the most massive core, and the most massive core vs. the gravitationally bound gas mass. Using archival Herschel and ALMA-IMF data for 11 nearby and 12 distant star-forming regions, the authors measure M_max_core and M_bound_gas and report log(M_max_core/M_sun) = 0.506 log(M_bound_gas/M_sun) - 0.32 (Eq. 7). Assuming a 30% efficiency to convert M_max_core to m_max_star and adopting an SFR-m_max_star relation from Yan et al. (2017) (Eq. 8), they derive SFR proportional to M_bound_gas^1.03, matching the Gao-Solomon relation (Fig. 3). They further argue, via Monte Carlo simulations, that the observed M_max_core-M_bound_gas relation is too tight to arise from random sampling of a standard core mass function.","tokens_in":25039,"tokens_out":3459,"duration_ms":34856,"significance":"If the central claim holds, the paper identifies a new empirical scaling relation between the most massive core and the gravitationally bound gas mass, and it offers an alternative to the stochastic-sampling picture for the origin of the Gao-Solomon relation. The analysis has notable strengths: the core extraction is applied consistently at the same physical resolution across both nearby and distant samples, the correlation is robust to the three adopted dust temperature estimators (Appendix C), and the Monte Carlo test in Section 4.2 explicitly quantifies the tension with random CMF sampling. The M_max_core-M_bound_gas correlation itself appears to be a real and useful observational result. However, the conversion of this correlation into a derivation of the Gao-Solomon relation rests on several load-bearing assumptions that are not independently established in this manuscript, so the explanatory claim is not yet secured.","major_comments":[{"comment":"The comparison in Figure 3 tests the Gao-Solomon relation only if M_bound_gas can be substituted for M_dense_gas. This equivalence is asserted in Footnote 3, with the actual validation deferred entirely to Jiao et al. (submitted), a companion paper that is not available for evaluation. The manuscript does not quantify M_bound_gas/M_dense_gas for the sample or show that this ratio is constant over the three orders of magnitude in mass. If the ratio varies with cloud mass or with star formation activity, then both the slope and the normalization of the derived SFR-M_bound_gas relation change, and Figure 3 is not a comparison against the observed Gao-Solomon relation. The authors should either present a direct demonstration of M_bound_gas ~ M_dense_gas for these clouds or reframe the Gao-Solomon comparison as a conditional prediction rather than an empirical test.","section":"§4.1, Footnote 3 and Figure 3"},{"comment":"Equation (8) is taken from a theoretical model (Yan et al. 2017) and then, later in the same section, the paper states that it 'can also be regarded as empirically constrained by the presented data.' But the presented data constrain Eq. (8) only through Eq. (7), the assumed 30% efficiency, and the M_bound_gas ≈ M_dense_gas substitution. Using Eq. (8) both as input to the recipe and as part of the consistency check in Figure 3 makes the agreement partially circular. The authors should separate what is assumed from what is tested; an independent calibration of Eq. (8) from directly measured m_max_star values or from YSO-count SFRs would resolve this concern.","section":"§4.1, Eq. (8)"},{"comment":"The constant 30% star-forming efficiency converting M_max_core to m_max_star is assumed without direct measurement, and it sets the zero-point of the derived SFR. The agreement in Figure 3 therefore does not validate the normalization of the hypothesized relation; a different assumed SFE would shift the derived relation vertically while preserving its slope. The authors should present the sensitivity of the SFR-M_bound_gas comparison to the assumed SFE and provide any empirical constraints that favor ~30% specifically.","section":"§4.1, recipe step 2"},{"comment":"The fiducial slope 0.506 in Eq. (7) is quoted without an uncertainty, and it changes to 0.55 when four ALMA-IMF fields with known missing flux are excluded; the three temperature methods in Appendix C yield slopes of 0.506, 0.554, and 0.575, a spread of about 0.07. Since the exponent chain in Eq. (9) uses 0.506 directly, the claim that the final relation is 'approximately linear' ('slope 1.03') should carry the propagated uncertainty from these choices. Reporting a central slope with a systematic uncertainty and showing how the derived SFR-M_bound_gas slope changes under the alternative assumptions would make the result quantitative rather than suggestive.","section":"§3.3 and Eq. (7)"}],"minor_comments":[{"comment":"The caption states the scaling relation as M_bound_gas ∝ (M_max_core)^0.506, which is inverted relative to Eq. (7); the text and figure should be consistent about which quantity is the independent variable.","section":"Figure 2 caption"},{"comment":"The simulation step says the core sample pool follows 'the CMF (Equation 8)', but Equation 8 is the SFR-m_max_star relation; the CMF power law is Equation 10.","section":"§4.2, step 1"},{"comment":"Footnote 3 contains a typo, 'measured in the say described in Section 3.2'; it should read 'measured in the way described in Section 3.2'.","section":"Footnote 3"},{"comment":"Footnote 1 refers to 'Footenote 3'; the spelling should be 'Footnote 3'.","section":"Footnote 1"},{"comment":"The opening sentence refers to 'Hershel'; the observatory name should be 'Herschel'.","section":"§5, Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The M_max_core-M_bound_gas correlation is a solid empirical contribution, but the manuscript's main explanatory claim is presently supported by assumptions deferred to an unpublished companion paper and by a consistency check that is partly circular. I would encourage the editor to seek a revised version that either includes the M_bound_gas≈M_dense_gas validation or clearly demotes the Gao-Solomon comparison to a predicted consequence of a model. The paper is not at the level where the derivation can be accepted as an explanation of the Gao-Solomon relation on the evidence presented."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe thing to know: the M_max_core-M_bound_gas correlation is a real, new observational result and is likely to be cited. The paper's attempt to turn it into an explanation of the Gao-Solomon relation is a chain of plausible but partly circular assumptions, not a closed argument.\n\nWhat's genuinely new: they measure the most massive core and the bound gas mass in 23 clouds from Herschel and ALMA-IMF archives, all at common 0.03 pc resolution, and find a tight correlation over three orders of magnitude. The correlation survives three different temperature estimates (slopes 0.5-0.6), so it is not a temperature artifact. The Monte Carlo test against random CMF sampling is also a useful quantitative argument: random sampling gives too much scatter and too high M_max_core for Taurus/Perseus. That is a real contribution.\n\nWhere it gets soft: the slope of Eq. 7 is quoted without uncertainty. Excluding four missing-flux fields changes it from 0.506 to 0.55, and that should be reported as a systematic range, not a single number. The bigger issue is the substitution M_bound ~ M_dense, deferred to a companion paper. Figure 3 overlays the Gao-Solomon line on SFR_theory versus M_bound; if that equivalence varies with cloud mass or environment, the comparison doesn't test the actual Gao-Solomon relation. The 30% SFE that sets the normalization is assumed, not measured. And Eq. 8 is a theoretical curve from Yan et al. 2017 (a co-author), not a direct empirical calibration; the claim that the presented data also constrain it goes through Eq. 7, the SFE, and the same M_bound ~ M_dense substitution, so that part is circular. These are not fatal to the empirical result, but they mean the paper does not close the case for explaining Gao-Solomon. The authors are explicit about the conditional nature, which helps.\n\nWho it's for: anyone working on star formation laws, core/IMF sampling, or sub-grid recipes. It deserves a serious referee: a good referee should ask for a slope uncertainty, a fuller treatment of the M_bound/M_dense equivalence (or a clear caveat), and a sharper separation between the measured correlation and the derived SFR relation. Send to peer review.","headline":"The M_max_core-M_bound_gas correlation is a real, new observational result; the Gao-Solomon explanation is a plausible but not closed chain of assumptions.","tokens_in":25667,"tokens_out":4044,"would_cite":true,"duration_ms":36845,"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":"The paper claims that the Gao-Solomon star formation law is not a statistical sampling effect but the product of three deterministic correlations, anchored by a new scaling between the most massive core and the gravitationally bound gas…","keywords":["star formation","Gao-Solomon relation","dense molecular gas","gravitationally bound gas","most massive core","core mass function","initial mass function","molecular clouds"],"falsifier":"Measure, in the same star-forming regions, both the gravitationally bound gas mass from column-density PDFs and the dense gas mass as Gao and Solomon defined it (e.g., HCN J=1-0 luminosity under their conversion assumptions): if the two are not proportional across the sample, the comparison in Figure 3 is not actually against the Gao-Solomon relation and the claimed chain collapses. As an independent check, resolve the most massive cores and weigh the star or protostar each one actually forms; an efficiency far from 30% would destroy the normalization of the predicted SFR.","tokens_in":24516,"feed_emoji":"⭐","tokens_out":23554,"duration_ms":176559,"temperature":0.7,"pith_summary":"The paper sets out to explain why the star formation rate scales linearly with the mass of dense molecular gas, the Gao-Solomon relation, without appealing to random statistical sampling of the stellar initial mass function. Using Herschel and ALMA observations of 11 nearby and 12 distant star-forming regions, it reports a power-law relation between the mass of the most massive core and the gravitationally bound gas mass of the parent cloud, with slope 0.506 over three orders of magnitude in gas mass. Chaining that relation to a constant 30% core-to-star efficiency and a theoretical link between star formation rate and the most massive star yields ${\\rm SFR} \\propto (M^{\\rm bound}_{\\rm gas})^{1.03}$, which coincides with the Gao-Solomon relation. In this picture, the linear law is a deterministic consequence of three correlated steps, not a coincidence of averaging many clouds.","feed_headline":"Three linked relations rebuild the Gao-Solomon star-formation law","feed_subtitle":"The chain gives a star formation rate proportional to bound gas mass to the 1.03 power, matching observations.","key_machinery":"The load-bearing object is the $M^{\\rm max}_{\\rm core}$-$M^{\\rm bound}_{\\rm gas}$ scaling, measured at matched 0.03 pc resolution: for nearby clouds from Herschel column density maps and for distant clouds from ALMA-IMF 1.3 mm continuum, with the same source-extraction pipeline applied to both samples. The bound gas mass is defined as the mass above the turning point where the column-density probability distribution function (N-PDF) goes from log-normal to power-law, i.e., the self-gravitating tail, identified by maximum-likelihood and MCMC fitting for each cloud. The argument's engine is the exponent chain of Eq. 9: since ${\\rm SFR} \\propto (m^{\\rm max}_{\\rm star})^{2.04}$ and $m^{\\rm max}_{\\rm star} \\propto M^{\\rm max}_{\\rm core}$ while $M^{\\rm max}_{\\rm core} \\propto (M^{\\rm bound}_{\\rm gas})^{0.506}$, multiplying the exponents gives ${\\rm SFR} \\propto (M^{\\rm bound}_{\\rm gas})^{0.506\\times2.04} \\approx (M^{\\rm bound}_{\\rm gas})^{1.03}$; the near-unity slope is what turns two unrelated power laws into the observed linear star-formation law.","core_discovery":"The central discovery is the correlation $\\log(M^{\\rm max}_{\\rm core}/M_{\\odot}) = 0.506\\,\\log(M^{\\rm bound}_{\\rm gas}/M_{\\odot}) - 0.32$ (Eq. 7), with a Spearman coefficient of 0.83, spanning $M^{\\rm bound}_{\\rm gas}$ from roughly $10^2$ to $10^5\\,M_{\\odot}$ in 23 molecular clouds; excluding four fields with known missing short-spacing flux steepens the slope to 0.55. The paper then shows that inserting this relation into two previously proposed links, a constant 30% efficiency converting the most massive core into the most massive star and the theoretical relation $\\log({\\rm SFR}/(M_{\\odot}\\,{\\rm yr}^{-1})) = 2.04\\,\\log(m^{\\rm max}_{\\rm star}/M_{\\odot}) - 5.80$ (Eq. 8), produces ${\\rm SFR} \\propto (M^{\\rm bound}_{\\rm gas})^{1.03}$, matching the Gao-Solomon relation once it is rescaled upward by the factor 2.7 used in the literature. The paper argues that the Gao-Solomon relation is therefore the combination of three non-trivial correlations, one of them new: (i) SFR versus $m^{\\rm max}_{\\rm star}$, (ii) $m^{\\rm max}_{\\rm star}$ versus $M^{\\rm max}_{\\rm core}$, and (iii) $M^{\\rm max}_{\\rm core}$ versus $M^{\\rm bound}_{\\rm gas}$. It also argues that random sampling of a canonical core mass function cannot reproduce the tightness or the normalization of the observed core-gas correlation: for Taurus and Perseus, the chance of drawing 172 cores above $3\\,M_{\\odot}$ with none exceeding $41.1\\,M_{\\odot}$ is at most about $8\\times10^{-3}$, so the stochastic picture is rejected at better than 99.1% confidence. The paper notes honestly that the galactic-scale SFR-$m^{\\rm max}_{\\rm star}$ relation turning out to hold on individual cloud scales is unexplained and could be a coincidence.","pith_inferences":["If the chain holds, the exponent product 0.506 × 2.04 ≈ 1.03 is a near-coincidence: environments that shift either exponent (temperature, metallicity, turbulence, or a different column-density break) should shift the Gao-Solomon slope away from unity, a testable prediction for galaxies with unusual star-formation efficiencies.","The core-gas correlation rests on only 23 regions; an extension the paper leaves implicit is testing whether the same 0.506 slope holds across substructures within a single giant molecular cloud, which would separate a universal deterministic law from an artifact of averaging a few regions.","If the star-formation-rate versus most-massive-star relation truly holds at cloud scale, the upper end of the stellar initial mass function encodes the current star formation rate; resolving the stellar content of the most massive cores in this sample would test that link directly."],"forward_implications":["The Gao-Solomon relation becomes a corollary of the new core-gas scaling: chaining Eq. 7 through the 30% efficiency and Eq. 8 yields a slope of 1.03, so the linear law needs no averaging over many clouds to emerge.","The same deterministic recipe works for single low-mass clouds like Taurus and massive regions like W43 and W51, implying that star formation from cloud scale to galaxy scale follows one chain.","Random sampling of the core mass function is ruled out as the origin of the most massive core: simulated clouds with canonical CMF slopes produce systematically heavier and more scattered most-massive cores than observed, so ten Taurus-mass clouds would not, in this picture, collectively form OB stars.","A single cloud measurement becomes predictive: given the gravitationally bound gas mass, the recipe yields both the expected most massive star and the expected star formation rate of the region."],"supporting_citations":[{"why":"Defines the empirical linear SFR-dense gas mass relation that the paper aims to derive, and fixes the comparison normalization in Figure 3.","marker":"Gao & Solomon 2004"},{"why":"Supplies the theoretical SFR versus most massive star relation (Eq. 8), derived from optimal sampling of the IMF, which is one of the three links in the chain.","marker":"Yan et al. 2017"},{"why":"Shows that the most massive star in a cluster is too tightly tied to cluster mass for stochastic IMF sampling, motivating the deterministic picture.","marker":"Weidner et al. 2013"},{"why":"Provides the factor 2.7 upward revision of the Gao-Solomon relation against which the derived SFRs are tested.","marker":"Lada et al. 2012"},{"why":"The ALMA-IMF large program whose 1.3 mm continuum images provide the most massive core measurements for the distant clouds.","marker":"Motte et al. 2022"},{"why":"Documents the ALMA-IMF data products and identifies the four fields with missing short-spacing flux that the paper treats as outliers.","marker":"Ginsburg et al. 2022"},{"why":"The Herschel Gould Belt survey that supplies the column density maps used for nearby cloud cores and bound gas masses.","marker":"André et al. 2010"},{"why":"Formulates the optimal sampling scheme of the IMF under which the deterministic SFR-most massive star and most massive star-cluster mass relations are derived.","marker":"Kroupa et al. 2013"},{"why":"Provides the dust opacity law that converts the 1.3 mm core fluxes into the core masses entering Eq. 7.","marker":"Ossenkopf & Henning 1994"},{"why":"Supplies the maximum-likelihood power-law fitting method used to locate the N-PDF break and thereby define the gravitationally bound gas mass.","marker":"Clauset et al. 2009"}],"fun_headline_variants":["Chain of three relations rebuilds the Gao-Solomon law","New core–gas correlation explains star formation scaling","Three linked laws reproduce SFR–dense gas proportionality","Star formation law emerges from core-mass chain","Why SFR scales with dense gas: three-step explanation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Two premises carry the whole chain: that the gravitationally bound gas mass equals the dense gas mass of the Gao-Solomon relation, an identification the paper asserts in a footnote and leaves to a companion paper, and that every most massive core converts to its most massive star at an exactly constant 30% efficiency.","fun_headline_variants_meta":{"raw":{"variants":["Chain of three relations rebuilds the Gao-Solomon law","New core–gas correlation explains star formation scaling","Three linked laws reproduce SFR–dense gas proportionality","Star formation law emerges from core-mass chain","Why SFR scales with dense gas: three-step explanation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000322,"raw_usage":{"total_tokens":2041,"prompt_tokens":1409,"completion_tokens":632,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":1025,"completion_tokens_details":{"reasoning_tokens":555}},"tokens_in":1025,"tokens_out":632,"duration_ms":5964,"temperature":1.0,"reasoning_tokens":555,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:08:48.645073+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, in the same star-forming regions, both the gravitationally bound gas mass from column-density PDFs and the dense gas mass as Gao and Solomon defined it (e.g., HCN J=1-0 luminosity under their conversion assumptions): if the two are not proportional across the sample, the comparison in Figure 3 is not actually against the Gao-Solomon relation and the claimed chain collapses. As an independent check, resolve the most massive cores and weigh the star or protostar each one actually forms; an efficiency far from 30% would destroy the normalization of the predicted SFR.","supporting_citations":[{"cited_title":"& Solomon, P","cited_arxiv_id":null,"evidence_quote":"Defines the empirical linear SFR-dense gas mass relation that the paper aims to derive, and fixes the comparison normalization in Figure 3."},{"cited_title":"2017, A&A, 607, A126","cited_arxiv_id":null,"evidence_quote":"Supplies the theoretical SFR versus most massive star relation (Eq. 8), derived from optimal sampling of the IMF, which is one of the three links in the chain."},{"cited_title":"J., Forbrich, J., Lombardi, M., & Alves, J","cited_arxiv_id":null,"evidence_quote":"Provides the factor 2.7 upward revision of the Gao-Solomon relation against which the derived SFRs are tested."},{"cited_title":"2022, A&A, 662, A9","cited_arxiv_id":null,"evidence_quote":"Documents the ALMA-IMF data products and identifies the four fields with missing short-spacing flux that the paper treats as outliers."},{"cited_title":"2013, in Planets, Stars and Stellar Systems","cited_arxiv_id":null,"evidence_quote":"Formulates the optimal sampling scheme of the IMF under which the deterministic SFR-most massive star and most massive star-cluster mass relations are derived."},{"cited_title":"& Henning, T","cited_arxiv_id":null,"evidence_quote":"Provides the dust opacity law that converts the 1.3 mm core fluxes into the core masses entering Eq. 7."},{"cited_title":"R., & Newman, M","cited_arxiv_id":null,"evidence_quote":"Supplies the maximum-likelihood power-law fitting method used to locate the N-PDF break and thereby define the gravitationally bound gas mass."}],"review_version":1}