{"id":"7bab1be0-27d2-4d12-a9be-937300706f60","arxiv_id":"2608.12473","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Combining 12CO, 13CO, and C18O line data with optical-depth corrections recovers the self-gravitating gas in molecular clouds, matching dust-based estimates within a factor of roughly two.","lead":"Astronomers combined three forms of carbon monoxide (CO) to measure how much gas in a cloud is dense enough to form stars. The new method agrees with dust-based measurements and can separate overlapping clouds using velocity, which dust cannot do.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"CO-based N-PDFs show a high-column-density decrement relative to dust in Orion A and Aquila; its impact on M_bound is asserted but never quantified.","rationale":"The reader's weakest assumption was C18O optical thickness. That concern is plausible but not exactly the load-bearing failure: Equation 7 formally allows C18O to be optically thick, since it solves for tau_13 from the ratio using tau_C18O = tau_13/chi. The observable failure is the high-column-density decrement in the CO-based N-PDFs, which is acknowledged in Section 4 and visible in Figures 5 and 14. The origin may indeed be opacity or depletion, but the load-bearing question is quantitative: how much mass is lost from the power-law tail, and does it change M_bound? The paper gives no such estimate. The reader did mention this unquantified decrement as a secondary point, so there is partial agreement, but I do not identify C18O thinness itself as the central vulnerability. The concern does not overturn the conditional verdict; it strengthens the conditions under which the method should be accepted. If the proposed clipping/reassignment test shows a large change in M_bound, the headline claim of robustly tracing bound gas would need to be weakened. The code for N-PDF fitting is public, which is positive, but it does not address this data-level systematic. I recommend keeping the reader's CONDITIONAL verdict rather than changing it to ACCEPT or REJECT, because the method is promising and independently testable but the high-density behavior of the CO-based tail remains unquantified.","tokens_in":30214,"tokens_out":6290,"duration_ms":65432,"concrete_test":"For Orion A and Aquila, recompute M_CO_bound after clipping or reassigning the high-density pixels where the CO-based N-PDF first deviates downward from the dust-based N-PDF by more than 2σ: replace those pixel column densities with the dust-based values, refit the N-PDF and re-evaluate Equation 14, and compare the new M_bound with the published value. Also repeat the fit with the N-PDF restricted to the range below the divergence to see whether the power-law slope changes. If the fractional change in M_bound is below 10%, the claim that the decrement is negligible is supported; if it is above 30%, the robustness claim fails and the quoted uncertainties must be enlarged to include this systematic effect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the CO isotopologue combination recovers the power-law tail of the N-PDF and hence the gravitationally bound mass. Section 4 and Figures 5/14 show that in Orion A and Aquila the CO-based N-PDFs fall below the dust-based N-PDFs at the highest column densities. The paper states that this discrepancy 'does not significantly affect the identification of gravitationally bound structures,' but no quantitative support is given. Because M_bound is defined as the mass above the fitted threshold (Equation 14), a decrement in the high-density tail directly reduces the mass assigned to the densest pixels. Whether the cause is C18O opacity, CO depletion, excitation effects, or abundance variations, the observable consequence is an incomplete recovery of the very gas that the power-law tail is meant to trace. The authors' assertion that the effect is negligible needs a quantitative bound; otherwise the good slope of 0.97 ± 0.11 in Figure 7a could reflect compensating errors between threshold placement and tail shape. Table 2 also shows that the CO-based transition column density for Orion A is 42.5 × 10^21 cm^-2 versus 10.6 × 10^21 cm^-2 from dust, a factor of four difference that is hard to reconcile with the claim of broadly consistent transition densities. This is the most load-bearing soft spot: an acknowledged, visible discrepancy whose effect on the headline quantity, M_bound, is left unquantified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a method for constructing H2 column density maps from 12CO, 13CO, and C18O J=1-0 data by using optically thin 13CO where C18O is undetected and an optical-depth-corrected 13CO column density derived from the 13CO/C18O intensity ratio where both lines are detected. The column densities are converted to H2 using a Galactic 12C/13C gradient and a metallicity-dependent CO abundance. The authors fit log-normal plus power-law N-PDFs and define the bound gas mass as the mass above the transition column density, following Jiao et al. (2025). They compare CO-based and dust-based N-PDFs, bound structures, and bound masses for 16 Milky Way clouds, finding a fitted slope of 0.97 ± 0.11 in the bound-mass comparison, mean IoU around 0.55 for bound structures, and consistency with the M_bound-SFR relation. The paper concludes that the multi-line CO method is a reliable, scalable, velocity-resolved alternative to dust emission for tracing self-gravitating gas.","tokens_in":30363,"tokens_out":11244,"duration_ms":104509,"significance":"If validated, the method is significant: it would allow N-PDF and bound-mass analyses to be carried out with existing large-area CO surveys, including in the Galactic plane where velocity information can separate LOS confusion that dust cannot. The comparison against independent dust-based measurements, the public release of the fitting code, and the extension to the 12CO+13CO pair are concrete strengths. The central claim is conditional on quantifying the acknowledged high-column-density decrement in the CO-based N-PDFs and on resolving the internal inconsistency about LOS confusion in the distant-cloud sample. These issues are addressable and do not undermine the overall approach, but they need to be fixed before the headline claim is fully supported.","major_comments":[{"comment":"The high-column-density decrement of the CO-based N-PDFs relative to dust in Orion A and Aquila is acknowledged but not quantified. Because M_bound is defined as the mass above the fitted threshold (Eq. 14), a deficit in the power-law tail directly reduces the mass assigned to the densest pixels, and the statement that this discrepancy 'does not significantly affect the identification of gravitationally bound structures' needs a quantitative test. Table 2 shows that the absolute CO-based transition column density for Orion A is 42.5 × 10^21 cm^-2 versus 10.6 × 10^21 cm^-2 from dust, and for Orion B is 36.9 versus 9.0 × 10^21 cm^-2; a factor of about four in the integration threshold is hard to reconcile with the claim of broadly consistent transition densities. Please recompute M_bound excluding or correcting the affected high-density pixels (for example, replacing them with the dust-based tail or applying a conservative C18O opacity/depletion correction) and show that the fitted slope and the individual mass ratios in Figure 7a survive. Without this, the 0.97 ± 0.11 slope could reflect compensating errors between threshold placement and tail shape.","section":"§4, Figs. 5/14, Table 2"},{"comment":"The optical-depth correction assumes τ_C18O = τ_13CO/χ and identical excitation temperatures for both isotopologues at every voxel. The uncertainties quoted in Table 1 propagate the scatter in the stacked abundance ratio and the range of optical-depth thresholds, but they do not capture systematic failures of this assumption in the densest gas. The observed high-density decrement in Figures 5 and 14 shows that the assumption does not fully recover the dense gas. I ask for a sensitivity analysis that varies χ over the plausible range (including the scatter in Figure 10) and allows for C18O opacity or depletion in the power-law tail, reporting how N_thres and M_bound respond. This is load-bearing because the absolute scale of N_H2 enters M_bound linearly.","section":"§3.1.3, Eqs. (7)–(8)"},{"comment":"The claim of a controlled test with minimal LOS confusion is internally inconsistent. Section 4 states that the distant clouds have 'more substantial LOS complexity' and that the main 13CO velocity component contributes only about 40%–90% of the total 13CO integrated flux, while Section 6 says the sample 'has been shown to suffer minimal LOS confusion.' Because the dust-based reference integrates all LOS components and the CO method isolates one velocity component, the good M_bound agreement for distant clouds could partly reflect comparing different physical gas, weakening the scalability claim. Please either quantify the bias introduced by the 40%–90% flux fractions or soften the Section 6 wording and show the mass comparison with and without the most confused sources.","section":"§4 and §6"}],"minor_comments":[{"comment":"The text reports a CO-based M_bound–SFR slope of 1.08^{+0.11}_{-0.10}, while the Figure 9 caption and the adjacent text give 0.97^{+0.09}_{-0.09} and 0.98 ± 0.08; please reconcile these values.","section":"§5.2, Fig. 9"},{"comment":"The caption says the mean IoU is 0.52, whereas the in-panel label reads mean = 0.55; the text also uses different IoU thresholds in different places, so please make the numbers consistent.","section":"Fig. 6"},{"comment":"The caption does not state whether the two 12CO+13CO test clouds (red squares) are included in the fitted slope; if they are, the slope mixes two different tracer combinations and should be refit or the test points should be shown as open symbols.","section":"Fig. 7"},{"comment":"The section title contains a typo: 'Galatic' should be 'Galactic'.","section":"§5.3"},{"comment":"The facilities line spells 'Hersechel' instead of 'Herschel'.","section":"Facilities"},{"comment":"The paper does not describe the regression method used for the fitted slope in Figure 7a or state whether uncertainties on both axes are accounted for; please clarify the fitting procedure and report the scatter in log space as well as the slope.","section":"Fig. 7a"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the reference list is appropriate. The main issue is the unquantified dense-gas decrement and the LOS-confusion inconsistency; both are fixable with additional analysis and clearer presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely useful: it combines 12CO, 13CO, and C18O J=1-0 to build column density maps that recover the power-law tail of the N-PDF, which single-isotopologue CO studies have repeatedly failed to do. That alone is a step forward. The demonstration on 16 clouds, with bound masses matching independent dust-based masses at a slope of 0.97 ± 0.11, is a credible validation. The velocity-resolved nature of CO is a real advantage for separating overlapping clouds along the line of sight, and the public N-PDF fitting code is a plus. The paper is also honest about the distant-cloud limitations from dust line-of-sight contamination.\n\nNow the soft spots, in proportion. The biggest one is the high-density decrement in the CO-based N-PDFs for Orion A and Aquila. The paper states this does not significantly affect the identification of bound structures, but no quantitative support is given. Since M_bound is the mass above the fitted threshold, a decrement at the highest column densities directly trims the densest pixels, and the excellent slope in Figure 7a could in principle hide compensating errors between threshold placement and tail shape. The factor-of-four difference in the absolute transition column density for Orion A (42.5 vs 10.6 × 10^21 cm^-2) makes this harder to wave away, even if the normalized transition densities look more consistent. The related C18O optical-depth assumption is acknowledged but not tested; that is the load-bearing premise and it deserves a quantitative check.\n\nA second, more minor issue: the quoted error bars on M_bound do not propagate the abundance gradient, CO/H2 ratio, or metallicity uncertainties. The paper shows one dust-opacity alternative but does not fold systematics into the main result. The IoU threshold of 0.4 rather than 0.5 is a small blemish; the justification is reasonable, though it slightly weakens the claim of unanimous spatial agreement.\n\nNone of this sinks the paper. The central argument holds up — the method works about as well as dust for tracing gravitationally bound gas in relatively clean sightlines. What is missing is a quantitative bound on the tail decrement's effect on M_bound, plus a more honest error budget. Both are fixable in revision. I would send this to peer review, and I would ask the authors to recompute M_bound with a truncated or corrected high-density treatment as a robustness test.","headline":"A credible optical-depth-corrected CO isotopologue method that recovers N-PDF power-law tails and gives bound masses in good agreement with dust, though one acknowledged high-density discrepancy is left unquantified.","tokens_in":31129,"tokens_out":2601,"would_cite":true,"duration_ms":24743,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Combining three CO isotopologue lines recovers the same gravitationally bound gas mass as dust emission, with a fitted slope of 0.97 ± 0.11 across 16 Milky Way clouds.","keywords":["stars: formation","molecular clouds","column-density probability distribution function","CO isotopologues","gravitationally bound gas","optical depth correction","N-PDF power-law tail","Galactic plane surveys"],"falsifier":"On a benchmark cloud such as Orion A, measure the C$^{18}$O optical depth in the highest-column-density pixels that dominate the power-law tail, using an optically thin comparison line such as C$^{17}$O or a higher-$J$ C$^{18}$O transition. If the inferred $\\tau_{18}$ exceeds roughly 0.2–0.3 in those pixels, the correction in Equation 8 underpredicts the $^{13}$CO column density, so the CO-based bound mass is biased low and the agreement with dust would degrade.","tokens_in":29886,"feed_emoji":"☁️","tokens_out":11288,"duration_ms":89241,"temperature":0.7,"pith_summary":"Astronomers find that the rate of star formation tracks the mass of gravitationally bound gas, which is read off from the power-law tail of the column-density probability distribution function (N-PDF), normally measured with dust emission. This paper argues that the same bound gas can be traced with CO isotopologue lines alone: combining $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1$–0 with an optical-depth correction extends the column-density range that any single line covers. On 16 Milky Way clouds, the CO-based N-PDFs reproduce the characteristic log-normal-plus-power-law shape, and the derived bound masses agree with dust-based values to within a factor of about two, with a best-fit slope of $0.97\\pm0.11$. If this holds, it provides a scalable and velocity-resolved way to identify star-forming gas in the Galactic plane, where overlapping clouds along the line of sight contaminate dust-based maps.","feed_headline":"Three CO lines recover the same star-forming gas mass as dust","feed_subtitle":"The velocity-resolved method can map bound gas in crowded Galactic-plane sightlines where dust emission blends clouds.","key_machinery":"The central mechanism is an optical-depth-aware combination of $^{13}$CO and C$^{18}$O $J=1$–0 lines. Where C$^{18}$O is detected, the $^{13}$CO/C$^{18}$O intensity ratio is inverted through the radiative-transfer relation $T_{13}/T_{18} = (1-e^{-\\tau_{13}})/(1-e^{-\\tau_{13}/\\chi})$ to obtain $\\tau_{13}$, and the $^{13}$CO column density is boosted by the factor $\\tau_{13}/(1-e^{-\\tau_{13}})$. Where C$^{18}$O is absent, $^{13}$CO is assumed optically thin. The $^{13}$CO/C$^{18}$O abundance ratio $\\chi$ comes from stacking spectra in optically thin pixels, and conversion to H$_2$ uses a Galactocentric $^{12}$C/$^{13}$C gradient and a metallicity-dependent CO abundance. The resulting maps are fit with a piecewise log-normal-plus-power-law N-PDF (Equation 12) using a Bayesian MCMC, and the mass above the log-normal/power-law transition is the reported bound mass.","core_discovery":"The paper's central claim is that the combination of three CO isotopologue lines, processed through an optical-depth-aware column-density reconstruction, identifies the same gravitationally bound gas that dust emission identifies via the N-PDF power-law tail. The demonstration spans 16 molecular clouds from $5\\times10^3$ to $10^6$ solar masses and distances of 0.4 to 11 kpc. The bound masses from the two tracers are consistent at a fitted slope of $0.97\\pm0.11$ with most sources within a factor of two; the spatial overlap of the identified bound structures averages an IoU of about 0.52, with every source above 0.4; and the CO-based bound masses reproduce the linear bound-mass–star-formation-rate relation with a slope near unity. The paper concludes that CO isotopologues can serve as a reliable, scalable, velocity-resolved alternative to dust emission for tracing the self-gravitating component of molecular clouds.","pith_inferences":["The authors leave implicit that a large Galactic-plane survey application is now possible: applying this method to survey data would map the fraction of gravitationally bound gas across environments, and one testable prediction is that the Central Molecular Zone's low star formation efficiency appears as a low bound fraction despite abundant dense gas.","The slight deficit of very high column density gas in the CO-based N-PDFs of Orion A and Aquila suggests the method may underestimate the most extreme cores; a quantitative comparison with an optically thin dense-gas tracer such as N$_2$H$^+$ would show whether that bias matters.","At galaxy scales, the same N-PDF machinery could be applied to CO isotopologue observations of external galaxies to measure bound gas fractions, provided metallicity and isotope gradients are known, making the $M_{\\rm bound}$–SFR relation testable beyond the Milky Way.","The velocity resolution suggests a cleaner test of the $M_{\\rm bound}$–SFR relation using velocity-resolved star formation tracers (e.g., H$\\alpha$ or radio recombination lines) instead of infrared luminosity, avoiding the line-of-sight mismatch the paper notes."],"forward_implications":["The CO-based method recovers the dust-based bound gas mass to within a factor of about two across two orders of magnitude in cloud mass, with no systematic over- or underestimate at either end.","Because the CO lines are velocity-resolved, the method can isolate one cloud among overlapping line-of-sight components, making it usable in crowded Galactic-plane regions where dust-based N-PDFs blend multiple clouds.","The CO-based $M_{\\rm bound}$ reproduces the roughly linear $M_{\\rm bound}$–SFR relation, so star-formation-rate studies can proceed without dust-derived column density maps.","The same optical-depth-correction framework works with the $^{12}$CO+$^{13}$CO pair when C$^{18}$O is undetected, as demonstrated for Ophiuchus and S287, extending the method to clouds where the rare isotopologue is missing.","The derived $^{13}$CO/C$^{18}$O abundance ratios trace a Galactic gradient, giving a calibration that can be used by other isotopologue studies."],"supporting_citations":[{"why":"Supplies the parent sample of clouds with dust-based N-PDFs, the bound-mass definition via the power-law tail, and the M_bound–SFR correlation the paper reproduces.","marker":"S. Jiao et al. 2025"},{"why":"Supplies the radiative-transfer relation linking the 13CO/C18O line ratio to 13CO optical depth, which is the core of the optical-depth correction.","marker":"R. Galván-Madrid et al. 2013"},{"why":"Supplies the Galactocentric 12C/13C gradient model used to convert 13CO column density to H2.","marker":"A. M. Jacob et al. 2020"},{"why":"Supplies the oxygen-based metallicity gradient used to scale CO abundance with Galactocentric radius.","marker":"J. E. Méndez-Delgado et al. 2022"},{"why":"Supplies the local CO/H2 abundance value that anchors the abundance scaling.","marker":"J. H. Lacy et al. 2017"},{"why":"Supplies FUGIN 12CO/13CO/C18O J=1-0 data for many of the distant clouds.","marker":"T. Umemoto et al. 2017"},{"why":"Supplies MWISP data for Aquila and two of the disk clouds.","marker":"Y. Su et al. 2019"},{"why":"Supplies the constant column-density screen subtraction used to correct dust maps for line-of-sight contamination.","marker":"N. Schneider et al. 2015b"},{"why":"Supplies the lower dust opacity alternative used to test the systematic offset in the mass comparison.","marker":"K. M. Pontoppidan et al. 2024"}],"fun_headline_variants":["Three CO lines match dust for star-forming gas mass","CO isotopologues trace bound gas as well as dust","Velocity-resolved CO method recovers dust-like gas mass","Triple CO technique measures self-gravitating gas without dust","16 clouds show CO lines rival dust for bound gas mass"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes the rarer CO variant C$^{18}$O stays optically thin in exactly the dense gas that forms the power-law tail, so the correction for $^{13}$CO opacity is valid; if C$^{18}$O becomes opaque there, the computed bound masses come out too low.","fun_headline_variants_meta":{"raw":{"variants":["Three CO lines match dust for star-forming gas mass","CO isotopologues trace bound gas as well as dust","Velocity-resolved CO method recovers dust-like gas mass","Triple CO technique measures self-gravitating gas without dust","16 clouds show CO lines rival dust for bound gas mass"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000219,"raw_usage":{"total_tokens":1469,"prompt_tokens":995,"completion_tokens":474,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":611,"completion_tokens_details":{"reasoning_tokens":393}},"tokens_in":611,"tokens_out":474,"duration_ms":4482,"temperature":1.0,"reasoning_tokens":393,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:07:43.327674+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On a benchmark cloud such as Orion A, measure the C$^{18}$O optical depth in the highest-column-density pixels that dominate the power-law tail, using an optically thin comparison line such as C$^{17}$O or a higher-$J$ C$^{18}$O transition. If the inferred $\\tau_{18}$ exceeds roughly 0.2–0.3 in those pixels, the correction in Equation 8 underpredicts the $^{13}$CO column density, so the CO-based bound mass is biased low and the agreement with dust would degrade.","supporting_citations":[],"review_version":1}