{"id":"5e775dd8-0841-4566-b665-a509851bf71a","arxiv_id":"1909.00209","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Deuterium fractionation of six molecules mapped toward two high-mass clumps shows N2H+ deuteration enhanced by an order of magnitude over HCN, HNC, and HCO+, while methanol deuteration exceeds 10% only at the CO-depletion peak.","lead":"This paper maps six deuterated molecules across two young high-mass star-forming clumps in the same infrared dark cloud, comparing the chemistry of a colder, younger clump with a more evolved one. The maps show that deuteration of N2H+ responds to temperature and density while methanol deuteration is high only where carbon monoxide is frozen onto dust grains.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The CH3OH deuteration headline is not supported at the reported precision: at Soff the quoted D(CH2DOH/CH3OH) = 33.3±26.7% has a 1σ lower bound of 6.6%, below the '>10%' in the abstract; the high D-fraction rests on a single weak line and needs error propagation before publication.","rationale":"The reader's formal 'weakest_assumption' was the uniform 12C/13C ratio of 42.9, which mainly affects absolute D-fractions of HCO+, HCN, and HNC rather than the robust N2H+ gradient or the order-of-magnitude separation from N2H+. The reader's rationale, however, already flagged the CH3OH detection as the reason for CONDITIONAL, and my stress-test identifies that as the more load-bearing issue: the paper's third headline finding, single-deuterated methanol >10%, is not supported at the 1σ level by the values in Table 4. The line ratio at Soff is detected, but the column-density ratio has an uncertainty of ±26.7 percentage points, so the central value of 33.3% cannot be cleanly separated from 10% or from zero. The single-line, LTE-based derivation, combined with known maser contamination in some CH3OH transitions and uncertain beam filling for compact CH2DOH emission, adds systematic risk. I do not see a reason to reject the paper: the multi-species maps, the N2H+ and HCO+ gradients, and the careful error discussion are genuine strengths. The appropriate outcome remains CONDITIONAL, with the CH3OH claim softened or reworded as tentative.","tokens_in":37792,"tokens_out":13078,"duration_ms":119034,"concrete_test":"Recompute D(CH2DOH/CH3OH) at Soff from the published line intensities, Einstein A coefficients, partition functions, and T_rot map, propagating all 1σ uncertainties (line ratio, T_rot, column densities, flux calibration). If the 1σ lower bound remains below 10%, the abstract and Section 5.5 should be revised to report a tentative value or an upper limit. As a second check, verify the factor ~12 amplification from the measured line ratio (2.8±0.8%) to the column ratio (33.3±26.7%) by reproducing the partition-function and excitation correction; a discrepancy would indicate an error in the D(CH3OH) map.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Table 4 lists D(CH2DOH/CH3OH) = 33.3±26.7% at Soff and 10.3±8.0% at P1, with a line ratio at Soff of only 2.8±0.8% (CH2DOH 2(0,2)-1(0,1) / CH3OH 2(0,2)-1(0,1) A). The 1σ uncertainties overlap zero (lower bounds 6.6% and 2.3%, respectively), so the abstract's assertion of a '>10%' single-deuterated methanol fraction and Section 5.5's '10%-39%' are not established by the quoted errors. The D-fraction is derived from a single CH2DOH transition at 89.408 GHz under LTE, with T_rot from CH3OH lines that include known Class I maser transitions (e.g., 84.521, 95.169, and 218.440 GHz; Table 1). Section 5.4.3 further notes that CH2DOH emission is compact (~10'' wide), so beam-filling and excitation corrections are not independently constrained. The paper's own caveat in Section 5.5 that an underestimated CH3OH column could inflate D(CH3OH) makes the high value vulnerable until the conversion is checked.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents IRAM-30 m on-the-fly spectral line maps of two parsec-scale clumps, P1 and S, in the infrared dark cloud G28.34+0.06, together with archival Herschel, CSO, SCUBA-2, and VLA+Effelsberg NH3 data. The authors derive dust and gas temperature maps, H2 and NH3 column density maps, and deuterium fraction maps for six species: NH3, N2H+, HCN, HNC, HCO+, and CH3OH. The main claims are that low-J lines of the dense gas tracers are subthermally excited; that D(N2H+) is several percent and exceeds D(HCO+), D(HCN), and D(HNC) by about an order of magnitude; that these D-fractions are enriched by a factor of 2-3 in the colder, denser, younger clump S relative to P1; that D(NH3) is essentially uniform at the few x 10^-3 level; and that single-deuterated methanol reaches >10%, with a peak of ~30-40% toward the CO-depletion peak Soff. The paper argues that these chemical variations trace evolutionary differences between the two clumps and that the deuterium fraction of different species can serve as an evolutionary diagnostic for high-mass star formation.","tokens_in":38103,"tokens_out":6960,"duration_ms":69709,"significance":"If the main claims hold, this is a valuable comparative dataset: it provides simultaneous, spatially resolved deuterium fraction maps for six species in two clumps of different evolutionary stages in the same cloud, with physical properties derived from independent gas and dust temperature diagnostics. The paper's strengths are its mapping approach over single-pointing studies, the explicit error budget in Section 5.4, the use of isotopologue ratios to test optical depth assumptions, and the cross-check of LTE rotational diagrams against LVG/RADEX fits. The qualitative conclusion that N2H+ deuteration is higher than the carbon-bearing species and that the D-fractions of these species increase toward the colder, denser clump S is credible. However, the headline quantitative claim of >10% single-deuterated methanol is not supported at the stated precision, and the absolute D-fractions for HCN, HNC, and HCO+ carry an unquantified factor-1.5-2 systematic uncertainty from the adopted 12C/13C ratio.","major_comments":[{"comment":"The stated detection of single-deuterated methanol with D(CH2DOH/CH3OH) > 10% is not supported by the quoted errors. At Soff, Table 4 gives D = 33.3 ± 26.7%, so the 1σ lower bound is 6.6%, below the >10% threshold asserted in the Abstract and Conclusions; the underlying line ratio in Table 4 is only 2.8 ± 0.8% and rests on a single CH2DOH transition at 89.408 GHz. The '10%-39%' range in §5.5 also appears to be a scatter range rather than an uncertainty-weighted interval. Please propagate all relevant uncertainties (CH2DOH integrated intensity, T_rot from the CH3OH rotational diagram, beam-filling, and column density uncertainties) and either report the CH3OH deuteration as a tentative detection with a firm upper limit, or present additional checks (e.g., a second CH2DOH transition) before claiming >10%.","section":"§5.5, Table 4, Abstract"},{"comment":"The absolute D-fractions for HCN, HNC, and HCO+ are linearly scaled by a single assumed 12C/13C ratio of 42.9, while §5.4 item 4 explicitly notes that isotopic exchange in cold, dense gas can raise the ratio to 47-80 for these species. This factor of 1.5-2 is not folded into the quoted uncertainties in Table 4, and because the correction differs per species (e.g., ~80 for HCN/HNC versus ~47 for HCO+ according to Furuya et al. 2011), the relative ordering of D(HCN), D(HNC), and D(HCO+) and the comparison with literature values are not robust at the stated precision. Please either quote D-fractions with a systematic error bar reflecting the plausible 12C/13C range, or present the raw isotopologue line ratios (already listed in Table 4) as the primary quantitative result.","section":"§5.2, §5.4, Table 4"},{"comment":"The D(NH3) map combines NH3 column densities at ~5'' resolution (VLA+Effelsberg) with NH2D column densities from a single ~30'' IRAM beam. Unless the NH3 map is explicitly smoothed to the NH2D beam before forming the ratio, the claimed uniformity of D(NH3) at (5 ± 3) x 10^-3 could be a beam-mismatch artifact: compact NH2D emission would be beam-diluted while the NH3 reference would not, and the local value toward Soff could be underestimated. Please state the smoothing/resolution-matching procedure used for Figure 5 and Table 4, or add a resolution-matched analysis.","section":"§5.1, Table 1, Figure 5"}],"minor_comments":[{"comment":"The sentence 'CO starts to catastrophy called freezes out' contains a typo; it should read something like 'CO starts to freeze out catastrophically.'","section":"§7"},{"comment":"There are small language errors, including 'an line imaging survey' and 'determine the the extremely young protostellar objects'; these should be corrected in a language pass.","section":"§2.1, §7"},{"comment":"The H13CO+ and DCO+ rows appear in both subtable I and subtable II with different fit methods; consider clarifying in the table notes which rows correspond to GAUSS fits and which to HFS fits to avoid apparent duplication.","section":"Table A1"},{"comment":"The axis label 'Log1 N_r' appears garbled and should read 'Log10 N_rot' for clarity.","section":"Figure A4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid observational study with a useful multi-species dataset and mostly careful error accounting. The main obstacle is the CH3OH deuteration claim, which is overstated relative to the quoted uncertainties; the 12C/13C systematic also needs to be incorporated into the quantitative values. I expect the authors can address these with a revision that reframes the CH3OH result as tentative and adds the requested systematic error bars."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core of this paper is good: it delivers the first comparative, parsec-scale maps of six deuterated species toward the P1–S pair in G28.34, and the qualitative result—higher D-fraction for N2H+, HCO+, HCN, and HNC toward the colder, denser clump S, with N2H+ deuteration an order of magnitude above the others—looks robust. The detection of compact CH2DOH emission at the CO-depletion peak (Soff) is genuinely new and worth reporting. The authors also do several things right: they check isotopologue ratios to argue the 13C lines are optically thin, compare LTE and LVG excitation, map the gas temperature from NH3 independently of the dust temperature, and write an unusually honest error budget. Most of the systematic caveats a referee would raise are already in Section 5.4.\n\nThat said, the headline CH3OH result is not supported at the reported precision. Table 4 gives D(CH2DOH/CH3OH) = 33.3±26.7% at Soff and 10.3±8.0% at P1, with 1σ lower bounds of 6.6% and 2.3%—both below the abstract's flat \">10%\". The ratio itself rests on one weak CH2DOH line at 89.4 GHz under LTE, with rotation temperature taken from CH3OH lines that include known Class I maser transitions, and the paper concedes the CH2DOH emission is compact enough that beam-filling is unconstrained. The stress-test note is right: the authors should either propagate the errors properly, soften the abstract and Section 5.5, or both. The same applies to the \"up to ~40%\" language in the conclusions. This is a fixable problem, not a fatal one—the underlying detection and the qualitative enhancement at Soff are believable.\n\nThe uniform 12C/13C = 42.9 conversion is a genuine second-order uncertainty for the absolute D-fractions of HCN, HNC, and HCO+, and the paper's own acknowledgment that fractionation could raise it to 47–80 means the absolute numbers carry an extra factor of 1.5–2. That does not kill the P1-versus-S gradient, which is a line-intensity effect largely insensitive to the conversion, but it does mean the absolute values quoted in Table 4 should be treated as indicative, not precise. The NH2D/NH3 beam mismatch (5\" vs 30\") is also a real limitation, properly noted.\n\nOverall: this deserves serious peer review. It is an observational paper with new maps, a solid qualitative finding, and honest caveats—it just needs a revision where the CH3OH claims are scaled back to match the statistics. I would send it out, and I would cite it for the N2H+/HCO+/HCN/HNC gradients.","headline":"The mapped D-fraction gradients for N2H+, HCO+, HCN, and HNC are a solid new result, but the paper's >10% CH3OH deuteration claim outruns its own error bars.","tokens_in":38699,"tokens_out":1155,"would_cite":true,"duration_ms":65194,"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":"By mapping deuterium fractionation across two parsec-scale clumps in the same infrared dark cloud, this paper establishes that colder, younger clump S is 2-3 times more deuterated in N2H+, HCO+, HCN, and HNC than the more evolved P1, with…","keywords":["deuteration","deuterium fractionation","high-mass star formation","infrared dark clouds","molecular line mapping","N2H+","CH3OH","CO depletion"],"falsifier":"A direct measurement of the local 12C/13C ratio toward P1, S, and Soff, from optically thin double-isotopologue lines, would settle the main systematic: if the ratio is 47-80 rather than 42.9, the reported HCN, HNC, and HCO+ D-fractions drop by a factor of 1.5-2 and the species ordering could change.","tokens_in":37576,"feed_emoji":"🌌","tokens_out":6500,"duration_ms":54763,"temperature":0.7,"pith_summary":"This paper maps the ratio of deuterated to normal molecules for six species across two neighboring high-mass protostellar clumps, P1 and S, that sit in the same infrared dark cloud but at different evolutionary stages. The authors show that N2H+ is about ten times more deuterated than HCO+, HCN, and HNC, and that all four species are 2-3 times more deuterated in the colder, denser, younger clump S than in the warmer, more evolved clump P1. Single-deuterated methanol reaches 10-40% at a spot offset from S where CO is most depleted, a detection that single-pointing surveys had missed. The study argues that deuterium fractionation tracks environmental differences (temperature, density, CO freeze-out) and can serve as a chemical clock for the earliest phases of high-mass star formation.","feed_headline":"Colder, younger clump shows 2-3x more deuteration","feed_subtitle":"Mapped D-fractions of six molecules in IRDC G28.34 tie chemistry to clump age and temperature.","key_machinery":"The analysis rests on column-density ratio maps between deuterated and hydrogenated isotopologues (D/hyd), constructed from optically thin lines of 13C or 15N/18O isotopologues. For HCN, HNC, and HCO+, the hydrogenated column densities are recovered from their 13C isotopologues using a fixed 12C/13C ratio of 42.9; excitation temperatures are derived either from hyperfine structure fitting or from the kinetic temperature map from para-NH3. Gas temperature and density maps, built from Herschel/SED dust fits, NH3 hyperfine fitting, and LVG modeling of H2CO and CH3OH, provide the physical backdrop against which the D-fraction gradients are interpreted.","core_discovery":"The paper reports measured D-fractions for six molecules toward two high-mass protostellar clumps, P1 and S, in the filamentary infrared dark cloud G28.34+0.06. It finds that D(N2H+) is 2.1% toward P1 and 6.1% toward S, roughly an order of magnitude higher than D(HCO+), D(HCN), and D(HNC), which lie in the 0.4-1.5% range, and that these four D-fractions are enriched toward the colder (Tkin ~14 K), denser clump S by a factor of 2-3. By contrast, D(NH3) stays at $(5\\pm 3)\\times 10^{-3}$ across the region, and D(CH3OH) peaks at 10-40% at the location Soff, 20 arcseconds northeast of S, where CO depletion reaches about a factor of 10. The authors interpret these patterns as the chemical signature of an evolutionary sequence: deuterium fractionation favors cold, dense, CO-freeze-out gas, and different species respond at different rates.","pith_inferences":["One testable extension is to observe the Soff region with sub-arcsecond resolution (e.g., with ALMA): the model that CH3OH forms on grains as CO freezes out predicts a compact, cold, dense object at that position, which single-dish mapping cannot resolve.","If the 12C/13C exchange models are correct, pixel-by-pixel isotopic ratios would steepen the D-fraction gradients reported here, making the deuteration maps lower limits; this could be tested with higher-J 13C lines or by observing additional rare isotopologues.","The comparative two-clump design could be scaled to a larger sample of 70-micron-dark/bright clump pairs to convert deuterium fractionation from a chemical tracer into a quantitative evolutionary clock for high-mass star formation."],"forward_implications":["D(N2H+) of 2-6% places P1 and S between high-mass prestellar and protostellar objects, consistent with embedded but very young protostars.","D(CH3OH) above 10% toward Soff, comparable to Class 0 low-mass protostars, suggests a deeply embedded, dense object at the CO-depletion peak.","Maps reveal peaks and gradients that pointing observations miss, for example D(CH3OH) below 0.3% at the pointed location versus 10-40% at Soff.","The gradient from P1 to S in D(N2H+), D(HCO+), D(HCN), and D(HNC) supports temperature and density, rather than external UV, as the controlling factors.","Because P1 and S share the same natal cloud, the chemical differences are attributable to evolutionary stage rather than to environmental differences."],"supporting_citations":[{"why":"Introduces the CTEX method and the deuterium fractionation chemistry that the column density analysis builds on.","marker":"Caselli et al. 2002b"},{"why":"Provides the D(N2H+) versus CO depletion correlation used to compare P1 and S with low-mass prestellar cores.","marker":"Crapsi et al. 2005"},{"why":"The pointing survey of D(CH3OH) and D(NH3) in high-mass clumps that this mapping study extends and partly reproduces.","marker":"Fontani et al. 2015"},{"why":"Supplies the Galactic 12C/13C ratio R=7.5 DGC+7.6 = 42.9 used to convert 13C lines to 12C column densities.","marker":"Giannetti et al. 2014"},{"why":"Previous line observations establishing CO depletion and line profiles toward P1 and S.","marker":"Feng et al. 2016b"},{"why":"The iterative SED fitting method that yields dust temperature and H2 column density maps.","marker":"Lin et al. 2017"},{"why":"Model predictions of 12C/13C fractionation in cold gas (values of 47-80) used in the error budget for the D-fraction estimates.","marker":"Furuya et al. 2011"},{"why":"Provides the RADEX code used for the LVG fitting of CH3OH and H2CO lines.","marker":"van der Tak et al. 2007"}],"fun_headline_variants":["Deuteration rises 3x in colder, younger high-mass clump","N2H+ deuteration 10x higher than HCO+ and HCN","Deuteration maps unveil chemical youth in high-mass clumps","CH3OH deuteration spikes where CO freezes out"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The maps assume a single uniform 12C/13C ratio of 42.9 across the whole region; if isotopic exchange raises the true ratio to 47-80, the reported D-fractions for HCN, HNC, and HCO+ would be overestimated by a factor of 1.5-2.","fun_headline_variants_meta":{"raw":{"variants":["Deuteration rises 3x in colder, younger high-mass clump","N2H+ deuteration 10x higher than HCO+ and HCN","Deuteration maps unveil chemical youth in high-mass clumps","CH3OH deuteration spikes where CO freezes out"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000564,"raw_usage":{"total_tokens":2773,"prompt_tokens":1140,"completion_tokens":1633,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":756,"completion_tokens_details":{"reasoning_tokens":1554}},"tokens_in":756,"tokens_out":1633,"duration_ms":11971,"temperature":1.0,"reasoning_tokens":1554,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:58:39.854970+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement of the local 12C/13C ratio toward P1, S, and Soff, from optically thin double-isotopologue lines, would settle the main systematic: if the ratio is 47-80 rather than 42.9, the reported HCN, HNC, and HCO+ D-fractions drop by a factor of 1.5-2 and the species ordering could change.","supporting_citations":[{"cited_title":"2014, , 570, A65, 10.1051/0004-6361/201423692","cited_arxiv_id":null,"evidence_quote":"Supplies the Galactic 12C/13C ratio R=7.5 DGC+7.6 = 42.9 used to convert 13C lines to 12C column densities."}],"review_version":1}