{"id":"563eb696-fa8c-4dc0-956b-66e6b0c5efbe","arxiv_id":"1908.06999","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Embedded cosmic rays from accreting protostars can regulate the CO-to-H2 conversion factor to the Milky Way value for star-forming efficiencies above 2%, while the CI-to-H2 conversion factor varies by 1.2 dex depending on cosmic ray transport.","lead":"This paper models how cosmic rays from newborn stars inside gas clouds change the standard formulas astronomers use to convert carbon monoxide or neutral carbon emission into a measurement of molecular gas. The authors find that embedded cosmic ray sources can stabilize the CO conversion factor to the Milky Way value, while the neutral carbon conversion factor varies by more than tenfold depending on the assumed cosmic ray model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Headline claim depends on the unconstrained 1/r diffusive CR transport model; under the paper's own 1/r^2 rectilinear branch, embedded sources have no effect on X_CO.","rationale":"I agree with the reader's weakest assumption. The paper is careful and honest: it flags the unknown transport law, makes the code public, and presents a full grid with the r^-2 case as a null result. Those are real strengths. However, the central conclusion is drawn from the r^-1 branch, and the paper's own Section 3.3 concedes that the alternative branch removes the effect entirely. That is not a small parameter variation; it is a switch. The quoted uncertainty (2.5 +/- 1) x 10^20 is the spread across the grid under one transport law, not an estimate that includes transport uncertainty. The optically thin flux assumption and lack of freeze-out/grain chemistry are additional caveats, but they do not change the main point. A CONDITIONAL verdict remains appropriate: the mechanism is plausible and the models are internally consistent, but the headline 'regulation' claim should be presented as contingent on 1/r propagation until the transport law is constrained observationally or by higher-fidelity CR transport simulations. Since the reader already reached CONDITIONAL, no change in verdict is needed.","tokens_in":9435,"tokens_out":4036,"duration_ms":43662,"concrete_test":"Use ionization-sensitive observations toward a protocluster with known embedded protostars (e.g., HCO+/H13CO+ or DCO+/HCO+ maps, or H3+ absorption where feasible) to measure the radial CRIR profile within the cloud. Fit the derived zeta(r) to r^-beta over the relevant 0.1-10 pc scales; the 1/r and 1/r^2 branches differ by orders of magnitude at small radii, so this distinguishes them. Then rerun the six-model grid with the measured beta (or a beta = 1.5 intermediate as a quick sensitivity test). If the X_CO regulation and epsilon_g > 2% MW consistency vanish for beta >~ 1.3, the headline should be reframed as conditional on diffusive transport; if observations favor beta ~ 1, the claim stands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The positive claim that embedded protostellar cosmic rays regulate X_CO and reproduce the nearly constant Milky Way conversion factor rests entirely on the 1/r 'diffusive' transport model. The paper itself defines two transport regimes in Section 2 and states that it is not known exactly how CRs transport through molecular clouds. Section 3.3 then reports that if embedded CRs transport as r^-2 (rectilinear), there is no impact on X_CO because the CRIR is dominated by external sources. Thus the abstract-level statement that embedded sources 'regulate X_CO and decrease its variance' is not a robust conclusion of the model family; it is a property of one unvalidated branch. The quoted X_CO = (2.5 +/- 1) x 10^20 cm^-2 (K km s^-1)^-1 for LDI and the MW consistency shading in Figure 1 would not exist under the alternative branch. Because the transport exponent determines whether the mechanism operates at all, not just its magnitude, the central astrophysical claim has a binary sensitivity to an assumption the authors explicitly identify as unknown. This is a correctness risk rather than an internal inconsistency: the models are self-consistent, but the choice of 1/r is not independently constrained by the calculations presented here.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the modified 3d-pdr astrochemistry code from Paper I, which includes in-situ cosmic-ray attenuation, to compute the CO-to-H2 and CI-to-H2 conversion factors X_CO and X_CI for one-dimensional models of protoclusters embedded in molecular clouds. The models vary the external cosmic-ray spectrum (low and high), the presence of embedded protostellar cosmic rays, and the cosmic-ray transport law (diffusive 1/r vs rectilinear 1/r^2). The paper reports that for the diffusive-transport model with a low external spectrum, clouds with star formation efficiencies above 2% have X_CO = (2.5 +/- 1) x 10^20 cm^-2 (K km s^-1)^-1, consistent with Milky Way observations; that embedded sources reduce X_CO and X_CI and lower their variance; and that X_CI is highly sensitive to the assumed cosmic-ray model.","tokens_in":9708,"tokens_out":8119,"duration_ms":89257,"significance":"If the central mechanism holds, the paper offers a physically motivated explanation for the near-constancy of X_CO in the Milky Way and Local Group, and it predicts systematically lower X_CO in extreme cosmic-ray environments such as starbursts and galactic centers. The study is a forward-modeling calculation rather than a fit to the Milky Way value: it uses explicit cosmic-ray spectra, a public astrochemistry code, realistic protostellar mass sampling, and treats the observed Milky Way X_CO as an external benchmark. These are genuine strengths. However, the headline regulator claim is realized only in the 1/r diffusive transport branch, and the paper itself states that the transport law is unknown; this is the main load-bearing weakness addressed below.","major_comments":[{"comment":"The abstract-level claim that embedded protostellar cosmic rays \"regulate X_CO and decrease its variance\" is only true for the diffusive (1/r) transport branch. Section 2 states that it is not known how cosmic rays transport through molecular clouds and therefore considers both 1/r and 1/r^2 regimes, while Section 3.3 reports that if embedded cosmic rays transport as r^-2 \"there is no impact on X_CO because the CRIR is lower and dominated by the CRs originating from external sources rather than internal.\" The quoted value X_CO = (2.5 +/- 1) x 10^20 cm^-2 (K km s^-1)^-1 for clouds with SFE > 2% and the Milky Way consistency shading in Figure 1 are properties of the LDI branch only. Because the transport exponent determines whether the mechanism operates at all, the central astrophysical claim has a binary sensitivity to an assumption the authors explicitly identify as unknown. I request that the abstract and Section 3.4 be rephrased to make the LDI/diffusive-transport condition explicit, and that the paper either provide an independent physical or observational reason to prefer 1/r transport or present both branches as equally viable and report the resulting range of X_CO and X_CI.","section":"Section 2; Section 3.3; Abstract"},{"comment":"The UMIST12 network used here excludes freeze-out, desorption, and gas-grain reactions. Freeze-out can strongly reduce gas-phase CO in the cold, dense gas that dominates the H2 column in these models, and CO depletion is a first-order effect for the absolute value of X_CO. The sensitivity test reported for grain-assisted recombination of C+ and He+ does not constrain the freeze-out/desorption uncertainty. Since the paper's headline is a specific numerical agreement with the Milky Way value to within ~0.3 dex, the authors should either include a freeze-out/desorption sensitivity test or give a quantitative argument for why its omission does not bias X_CO and X_CI at that level.","section":"Section 2 (chemistry network)"},{"comment":"The text says that the integrated flux defined by Eq. (4) assumes the interstellar medium is entirely optically thin, immediately after stating that 3d-pdr uses an escape probability method to account for line opacity. These two statements are in tension. If the emissivity epsilon in Eq. (4) already includes the escape probability, then the integrated flux is not optically thin and the sentence is misleading; if it does not, the CO(1-0) flux from optically thick lines will be overestimated and X_CO underestimated. Please clarify which is the case and, if necessary, quantify the effect on the reported X_CO values.","section":"Section 2, Eq. (4)"}],"minor_comments":[{"comment":"The y-axis labels appear to render as \"log X_CI,i X_CI,LNA\" and \"log X_CO,i X_CO,LNA\" without a division sign; please correct the typesetting so the ratio is clear.","section":"Figure 3"},{"comment":"The word \"efficiences\" should be \"efficiencies\".","section":"Section 3.5"},{"comment":"The reference to \"Offner et al. 2019, sub.\" is to an unpublished manuscript; please update it or mark it as in preparation consistently.","section":"Section 2"},{"comment":"The notation \"CI\" is conventionally written as \"C I\" for neutral atomic carbon; consider using \"C I\" to avoid confusion with the cyanogen molecule CI or with the abbreviation for confidence interval.","section":"Title and throughout"}],"recommendation":"major_revision","confidential_remarks":"The main concern is the conditional nature of the headline result: the 1/r transport assumption is unvalidated and switches the mechanism off entirely, as the paper itself acknowledges in Section 3.3. This is fixable in revision by making the condition explicit and adding a plausibility or observational argument for the transport law; I do not see evidence of circularity or inappropriate fitting to the Milky Way value. The paper is otherwise a solid forward-modeling contribution, and I would be willing to review a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Brandt,\n\nQuick read on 1908.06999. The genuinely new thing here is that they put CR sources inside the cloud—protostellar accretion shocks—with in-situ attenuation, instead of just dialing up an external CRIR as prior work did. That changes the chemistry in a spatially resolved way and, in the diffusive branch, gives a mechanism that can hold X_CO near the Milky Way value across a wide range of surface densities and star formation efficiencies. X_CI turning out to be extremely CR-sensitive is a useful caution for high-redshift work.\n\nThe paper does a lot right. The model grid is transparent, the code is public from Paper I, the Milky Way comparison is an external check rather than a fit, and the limitations are flagged in the text rather than buried. The citation pattern to the CRIR-X_CO literature looks complete to me; the self-citations to Paper I and Gaches & Offner (2018b) are legitimate because the code and CR prescription live there.\n\nThe soft spot is the transport law, and it is not minor. The whole \"embedded sources regulate X_CO\" story runs through the 1/r diffusive model. The paper itself says in Section 3.3 that if the transport is 1/r^2, internal sources have no impact on X_CO because the external CR field dominates. They also note in Section 2 that it is not known how CRs transport through molecular clouds. So the abstract-level claim \"embedded sources regulate X_CO and decrease its variance\" is true only for one unvalidated branch of the model. The stress-test note gets this exactly right. I'd want the abstract and summary to say \"in the diffusive transport model\" rather than making it sound like a robust conclusion.\n\nTwo smaller things. The flux definition assumes optically thin emission, which is hard to defend for CO 1-0 in these clouds; even with escape probability, the emergent line is not what Equation 4 gives. That could bias the absolute X_CO values, though maybe not the trends. And the quoted ±1 × 10^20 is the grid spread, not a propagated uncertainty—fine as a spread, but it reads like an error bar.\n\nWho is this for? Anyone using CO or CI to estimate molecular gas in starbursts or the high-redshift universe, and people working on CR feedback. It deserves a serious referee; the mechanism is interesting and the modeling is careful. I would send it out, but the referee should insist that the transport dependence be in the abstract.","headline":"A useful modeling study of how embedded protostellar CRs affect X_CO and X_CI, but the headline 'regulation' result holds only for one of two unconstrained transport models.","tokens_in":10284,"tokens_out":2589,"would_cite":false,"duration_ms":25101,"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":"Cosmic rays from newborn stars may hold the Milky Way's CO-to-H2 conversion factor steady.","keywords":["cosmic rays","CO-to-H2 conversion factor","X_CO","X_CI","protostellar accretion shocks","molecular cloud chemistry","star formation efficiency","cosmic-ray ionization rate"],"falsifier":"Map the radial falloff of the cosmic-ray ionization rate around an embedded protostar using ${\\rm H}_3^+$ or OH$^+$ absorption toward multiple sightlines: a $1/r$ profile supports the claim, while a $1/r^2$ profile would remove the predicted regulation of $X_{\\rm CO}$.","tokens_in":9169,"feed_emoji":"🌌","tokens_out":9967,"duration_ms":82998,"temperature":0.7,"pith_summary":"This paper asks whether cosmic rays produced by accreting newborn stars inside a forming cluster change how astronomers convert carbon-monoxide or neutral-carbon emission into a molecular gas mass. Using one-dimensional cloud models whose chemistry computes cosmic-ray attenuation throughout the gas, it argues that embedded protostellar cosmic rays can regulate the conversion factors. The main result is that clouds with more than about 2% star formation efficiency have $X_{\\rm CO} = (2.5 \\pm 1)\\times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, matching the Milky Way average, even as cloud surface density and efficiency vary widely. A reader should care because this suggests the Galaxy's nearly constant CO-to-H2 conversion factor is partly a consequence of the star formation process itself, not just a fixed property of molecular gas.","feed_headline":"Protostar cosmic rays may steady the Milky Way's gas-mass yardstick","feed_subtitle":"New models find embedded cosmic rays hold X_CO near observed values once star formation passes 2%.","key_machinery":"The load-bearing object is the model of protostellar cosmic rays built from accretion-shock acceleration: each protostar produces a power-law momentum spectrum via diffuse shock acceleration, normalized by the shock energy, with maximum proton energies typically between 1 and 10 GeV. Those spectra are attenuated through the protostellar core column and then transported into the surrounding cloud either diffusively, with intensity falling as $1/r$, or rectilinearly, falling as $1/r^2$, combined with one of two external cosmic-ray spectra. Feeding these spectra into a modified photodissociation-region astrochemistry code that computes cosmic-ray attenuation in situ, together with the chemical network and CO and CI line emissivities, is what turns the cosmic-ray prescription into specific $X_{\\rm CO}$ and $X_{\\rm CI}$ predictions. The diffusive $1/r$ case is the one that lets embedded sources dominate the ionization rate and regulate the conversion factors.","core_discovery":"The paper's central claim is that cosmic rays accelerated at protostellar accretion shocks, if they diffuse through the surrounding cloud, act as an internal source that regulates the CO-to-H2 conversion factor. Once the star formation efficiency exceeds roughly 2%, these embedded cosmic rays lower $X_{\\rm CO}$ and flatten its dependence on cloud surface density and star formation efficiency, producing $X_{\\rm CO} = (2.5 \\pm 1)\\times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, consistent with the Milky Way value and its scatter. The same mechanism lowers $X_{\\rm CI}$ by over an order of magnitude and leaves a 1.2 dex dispersion, from $2\\times 10^{20}$ to $4\\times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, so the neutral-carbon conversion factor is far more sensitive to the assumed cosmic-ray model than $X_{\\rm CO}$ is.","pith_inferences":["If this is right, $X_{\\rm CO}$ may be partly self-regulating: more star formation makes more cosmic rays, which lowers $X_{\\rm CO}$ and trims the inferred gas mass, so the conversion factor could buffer against environmental changes.","A direct test would be to measure $X_{\\rm CO}$ in clouds with matched surface density but different star formation efficiencies; the model predicts a systematic drop in $X_{\\rm CO}$ with efficiency only under diffusive cosmic-ray transport.","Because $X_{\\rm CI}$ is so sensitive to the cosmic-ray model, comparing CI and CO maps in starbursts could constrain how cosmic rays actually travel through dense gas, which the paper identifies as an open question.","At high redshift, where galaxies form stars efficiently, CO-based gas masses computed with the Milky Way $X_{\\rm CO}$ could be systematically too high if embedded cosmic rays matter as modeled here."],"forward_implications":["If the diffusive-transport model is right, $X_{\\rm CO}$ should be nearly flat across clouds with star formation efficiency above roughly 2%, matching the observed Milky Way scatter.","Raising the assumed external cosmic-ray ionization rate to the observationally preferred $\\zeta \\approx 10^{-16}$ s$^{-1}$ lowers $X_{\\rm CO}$ by 0.2 dex in low-surface-density clusters with $\\Sigma_{\\rm cl} < 3$ g cm$^{-2}$, which would trim typical molecular-gas mass estimates.","Embedded cosmic rays reduce $X_{\\rm CO}$ by up to about 0.5 dex and $X_{\\rm CI}$ by up to about 1.2 dex relative to models with no internal sources, so CI-based gas masses depend strongly on the assumed cosmic-ray environment.","The models reproduce the observed tendency for $X_{\\rm CO}$ to decrease in extreme cosmic-ray environments such as starburst galaxies and the Galactic Center, implying that applying a single Milky Way conversion factor there would misestimate gas mass."],"supporting_citations":[{"why":"supplies the modified astrochemistry code with in-situ cosmic-ray attenuation and the model grid this paper analyzes","marker":"Gaches et al. 2019 (Paper I)"},{"why":"provides the fiducial Milky Way value and observed scatter of $X_{\\rm CO}$ that the models are compared against","marker":"Bolatto et al. 2013"},{"why":"defines the low and high external cosmic-ray spectra labeled L and H","marker":"Ivlev et al. 2015"},{"why":"gives the diffuse-shock-acceleration prescription for protostellar cosmic-ray spectra and their typical maximum energies","marker":"Gaches & Offner 2018b"},{"why":"supplies the attenuation of cosmic rays through the protostellar core column","marker":"Padovani et al. 2009"},{"why":"describes the photodissociation-region astrochemistry code extended to compute CO and CI line emissivities","marker":"Bisbas et al. 2012"},{"why":"provides the UMIST12 gas-phase chemical network used in the models","marker":"McElroy et al. 2013"},{"why":"establishes the prior external-CR scaling of $X_{\\rm CO}$ with star formation rate that this work contrasts with embedded sources","marker":"Clark & Glover 2015"}],"fun_headline_variants":["Protostar cosmic rays clamp CO-to-H2 conversion near Milky Way value","Embedded cosmic rays steady X_CO once star formation tops 2%","Protostar cosmic rays set CO yardstick, leave CI tracer sensitive","Embedded cosmic rays flatten CO gas-mass scatter, disrupt CI tracer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands on the premise that cosmic rays from embedded protostars diffuse through the cloud with a $1/r$ intensity profile, so internal sources dominate the ionization rate; if transport were rectilinear instead, declining as $1/r^2$, the paper finds essentially no effect on $X_{\\rm CO}$.","fun_headline_variants_meta":{"raw":{"variants":["Protostar cosmic rays clamp CO-to-H2 conversion near Milky Way value","Embedded cosmic rays steady X_CO once star formation tops 2%","Protostar cosmic rays set CO yardstick, leave CI tracer sensitive","Embedded cosmic rays flatten CO gas-mass scatter, disrupt CI tracer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000768,"raw_usage":{"total_tokens":3479,"prompt_tokens":1093,"completion_tokens":2386,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":2307}},"tokens_in":709,"tokens_out":2386,"duration_ms":18752,"temperature":1.0,"reasoning_tokens":2307,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:29:11.514008+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Map the radial falloff of the cosmic-ray ionization rate around an embedded protostar using ${\\rm H}_3^+$ or OH$^+$ absorption toward multiple sightlines: a $1/r$ profile supports the claim, while a $1/r^2$ profile would remove the predicted regulation of $X_{\\rm CO}$.","supporting_citations":[{"cited_title":"V., Padovani, M., Galli, D., & Caselli, P","cited_arxiv_id":null,"evidence_quote":"defines the low and high external cosmic-ray spectra labeled L and H"},{"cited_title":"C., & Glover, S","cited_arxiv_id":null,"evidence_quote":"establishes the prior external-CR scaling of $X_{\\rm CO}$ with star formation rate that this work contrasts with embedded sources"}],"review_version":1}