{"id":"2d298c80-17a8-4d52-a341-11ba28074441","arxiv_id":"2411.12916","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"A single-dish survey of DR21(OH) resolves multiple warm and cold gas components and finds that grain chemistry plus non-thermal desorption, not warm carbon chain chemistry, explains the observed methanol abundance.","lead":"This paper maps three molecules, CH3CCH, CH3OH, and H2CO, across the high-mass star-forming clump DR21(OH) using two single-dish telescopes, and separates the emission into warm, cold, and outflow components. It then runs chemical models to argue that methanol needs a shock-like, non-thermal desorption mechanism, while the other two molecules can form thermally.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Methanol conclusion rests on hand-fixed LTE excitation temperatures; RADEX column densities differ enough that the claimed need for non-thermal desorption may be an artifact of assumed abundances.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the fixed CH3OH LTE excitation temperatures set the absolute column densities, hence the abundances that the chemical model is tuned to reproduce. This is the most vulnerable link in the chain from observations to the abstract's central claim that non-thermal mechanisms are needed to produce CH3OH. The observational component decomposition and maps are a solid contribution, and the paper is transparent about its own limitations, including the timescale mismatch and the borrowed shock profile. The concern is not an outside-consensus disagreement but a propagation gap: no uncertainties are attached to the fixed Tex values, and the RADEX results already in the paper indicate that the LTE column densities are not uniquely determined. A single refit with free Tex, plus a rerun of the warm-up model with revised abundances, would settle whether the non-thermal conclusion is robust. Because the qualitative result is plausible and the manuscript explicitly acknowledges several of the uncertainties, I would keep the conditional verdict rather than escalate to reject or de-escalate to accept.","tokens_in":20738,"tokens_out":7237,"duration_ms":80113,"concrete_test":"Re-fit the CH3OH spectra with Tex as a free parameter, or with priors from the RADEX Tkin values in Table 3, and propagate the resulting column-density uncertainties. Then rerun the NAUTILUS warm-up model for components N44 3, 4, and 5 using the revised abundances with both the default and the RADEX post-shock densities. If the warm-up model reaches the revised observed abundance within 1.6e4 yr, the non-thermal-desorption conclusion fails; if it still falls short by more than an order of magnitude, the central claim survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All CH3OH excitation temperatures in Table 2 are fixed by hand (e.g., 82, 75, 22, 15, 25, 30, 6, 6, 30, 8, 8, 20, 8 K) with no quoted uncertainties. The NAUTILUS 'observed abundances' are ratios of these column densities to the H2 column, so the chemical conclusion inherits whatever Tex was assumed. Table 3 shows the same spectra can be fit by RADEX with substantially different column densities: across the N44 components the reported A-only Ntot values are factors of about 1.7-3 below the LTE totals, and the paper does not specify how the A/E splitting should be applied to restore the total. If the true CH3OH column is at the lower end, the warm-up model that already matches H2CO and CH3CCH may also match CH3OH, and the need for non-thermal desorption would disappear. The shock model does not rescue the inference because it is parameterized from IRAS 20126 (Palau et al. 2017), not from DR21(OH), and the fiducial shock timescale exceeds the 1.6e4 yr outflow age the paper itself derives. The abstract states the non-thermal requirement as a result, but it is conditional on unpropagated Tex assumptions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a multi-line, single-dish study of the high-mass star-forming clump DR21(OH)/N44 using IRAM 30-m and GBT observations of CH3CCH, CH3OH, and H2CO. A multi-component LTE model is used to decompose the spectra toward N44 and the neighboring clumps N40/N36/N41, N38, and N48, identifying warm (20-80 K), cold (6-15 K), and broad outflow components. RADEX is used to test non-LTE effects in CH3OH at two positions. The NAUTILUS gas-grain chemical code is then run in two-phase (cold-cloud plus warm-up) and three-phase (cold-cloud plus shock plus post-shock) configurations to reproduce the observed abundances. The central conclusion is that thermal (warm-up) chemistry can produce the observed H2CO and CH3CCH abundances, while CH3OH requires non-thermal desorption, modelled as a C-type shock. The paper also maps the spatial distribution of the components and argues for grain-surface formation of CH3CCH and CH3OH.","tokens_in":21036,"tokens_out":11600,"duration_ms":105926,"significance":"The paper's qualitative component decomposition is a useful contribution: it connects single-dish clump-scale observations with interferometric core-scale results and identifies multiple temperature/velocity components in a generally careful way. The RADEX follow-up for CH3OH is a reasonable check on LTE assumptions, and the extensive appendix material (model spectra for every component, production/destruction route figures) is transparent and helpful. The chemical modelling addresses a current question - whether WCCC or grain-surface chemistry dominates in high-mass clumps - and the conclusion that H2CO and CH3CCH can be made by thermal mechanisms is consistent with previous work. However, the quantitative claim that CH3OH needs non-thermal desorption is fragile: the CH3OH excitation temperatures are all fixed, the RADEX column densities are not reconciled with the LTE totals, and the fiducial shock model cannot match the outflow age derived in the same paper. These issues need to be resolved before the abstract's strong statement can be accepted.","major_comments":[{"comment":"All CH3OH excitation temperatures in Table 2 are fixed to integer values (82, 75, 22, 15, 25, 30, 6, 6, 30, 8, 8, 20, 8 K) and the table caption states that values without errors were fixed in the model. The LTE column densities, and therefore the observed fractional abundances used in §4, depend on these assumed temperatures; for an optically thin transition the column density scales roughly as exp(E_up/T_ex) once the partition function is accounted for, so a change from 15 K to 25 K can alter the column density by a factor of order two. The paper does not present a sensitivity test over T_ex for any CH3OH component. Because the central conclusion that CH3OH requires non-thermal desorption is based on these abundances, I request a grid over T_ex (or a fit where the lines permit) and a discussion of how the chemical-model conclusion would change if the RADEX-derived column densities of Table 3 were used instead.","section":"§3.1, Table 2"},{"comment":"The RADEX column densities in Table 3 are explicitly labelled as A-CH3OH only, while the LTE model does not differentiate A and E forms (Sect. 3). The paper never states the A/E ratio used to convert the RADEX values to total CH3OH column densities. In the centre pixel, the LTE totals are a factor of 1.4-3.1 larger than the RADEX A-only values (e.g., N44 4: 2.08e15 vs 7.84e14). Without a specified conversion, the quantitative comparison between LTE and RADEX is ambiguous, and the claim in §3.1.1 that CH3OH is 'near LTE' is not backed by a well-defined column-density comparison. Please provide total RADEX column densities or an explicit A/E ratio, and clarify how the non-LTE column densities would modify the input abundances to the chemical model.","section":"§3.1.1, Table 3"},{"comment":"The paper derives an outflow dynamical age of 1.6e4 yr from the outflow component's 15 km/s width and 10 arcsec extent, and then states that this age is less than the time needed for the fiducial shock model to match the observed CH3OH abundances in the post-shock stage. The only model variant that matches on a timescale of a few 1e3 yr uses the RADEX densities, but this variant is rejected because it requires an unexplained density jump between the shock and post-shock phases. Thus the fiducial model that supports the conclusion 'CH3OH needs some form of non-thermal desorption' is not dynamically consistent with the source's own outflow age. Either adopt the RADEX-density model and justify the density evolution, or explicitly present the conclusion as conditional on an unresolved timescale discrepancy.","section":"§4, Fig. 10"},{"comment":"The CH3CCH model invokes a cloud-edge environment with A_V = 1, and the text acknowledges 'we do not know if the extent of CH3CCH is truly on the envelope edge' (Sect. 4). The paper also states that A_V = 50 requires a shock to reproduce the CH3CCH abundances. Therefore the conclusion that CH3CCH can be produced by thermal mechanisms depends on an unverified geometrical assumption. The abstract's statement that thermal mechanisms are adequate for CH3CCH should be qualified unless independent evidence for the low-A_V location is provided, or the model should be presented as one of two admissible scenarios.","section":"§4, CH3CCH"}],"minor_comments":[{"comment":"In the description of the post-shock phase, 'set the density at 10^5 cm^-2' should read '10^5 cm^-3'; earlier in the same paragraph, 'stimulates the abundances over time' should be 'simulates'. These typos should be corrected.","section":"§4 (chemical modelling)"},{"comment":"The 135 GHz IRAM data are stated to be excluded from the analysis, yet their channel width, beam size, and pixel size are given in detail. Consider moving this information to a footnote or removing it, to avoid confusing the reader about which data are used.","section":"§2.1"},{"comment":"The statement that DR21(OH) was the first HMSFR found with evidence for WCCC (citing Mookerjea et al. 2012) should be accompanied by a brief description of what constitutes WCCC and a reference to the original WCCC papers (e.g., Sakai et al. 2008; Aikawa et al. 2008), which are currently cited later but not in this context.","section":"§1"},{"comment":"The caption says components are shown by 'dotted, dashed, or solid lines' but does not identify which component has which line style. A legend or explicit mapping from line style to component (C1-C5) would make the figure interpretable.","section":"Fig. 10 caption"},{"comment":"The sentence 'the chemical and physical environments of star-forming regions is revealed' has a subject-verb agreement error ('environments ... is'); it should be 'are revealed.'","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid, observationally rich study, and the component decomposition is likely correct. The main weakness is the overstatement of the chemical-model conclusion relative to the unpropagated T_ex assumptions and the acknowledged timescale problem. I believe the paper is suitable for A&A after major revision, provided the authors add a sensitivity analysis and soften the abstract accordingly. The citation of 'Freeman et al. (submitted)' for line identification is a reviewability concern; the authors should ensure the companion paper is available or summarize the line IDs. I do not see grounds for rejection, as the central qualitative result (non-thermal desorption needed for CH3OH) may well be correct but is not yet demonstrated with the required robustness."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this if you work on high-mass star-forming clumps. The mapping and multi-component LTE analysis are the real contribution: new IRAM 30m and GBT maps of CH3CCH, CH3OH, and H2CO across DR21(OH)/N44, and a transparent decomposition into warm and cold components toward MM1 and MM2, a broad blue-shifted outflow, plus components toward N48, N38, N40, N36, N41. The fits look careful, and the comparisons to prior interferometric work give the component separation credibility. This is honest observational progress on a benchmark source.\n\nThe chemistry section is the soft spot. The NAUTILUS models take the LTE column densities as observed abundances and then adjust environmental parameters until the model matches them. For CH3CCH they pick a low AV; for CH3OH they graft on a shock parameterized from IRAS 20126, not from DR21(OH). The paper flags the timescale mismatch itself—the outflow age is 1.6e4 yr and the CH3OH model needs longer—so it is not misleading, but the abstract's 'non-thermal mechanisms are needed' reads stronger than the evidence supports.\n\nThe stress-test note lands. All CH3OH excitation temperatures in Table 2 are fixed by hand with no uncertainties, and the RADEX results in Table 3 give A-only column densities about a factor of 2-3 below the LTE totals. Since the chemical inference is anchored to column density ratios, you cannot rule out that a lower true methanol column, combined with the warm-up chemistry that already works for H2CO and CH3CCH, would remove the need for a shock. The paper does not propagate the Tex uncertainty through the chemistry.\n\nNo code or data are released, which limits reuse of otherwise useful maps. The LTE components with the 84.5 GHz maser excluded are a valuable reference.\n\nI would send this to review. The observational component analysis justifies publication, and the chemistry claims can be tested by varying the fixed Tex values and re-running NAUTILUS. A good referee report should ask for that. The paper is honest on its own terms; it just overreaches in the abstract.","headline":"Worth reviewing for the component decomposition and maps; the methanol non-thermal claim is real but rests on fixed excitation temperatures and a tuned shock.","tokens_in":21609,"tokens_out":2755,"would_cite":false,"duration_ms":28027,"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":"The DR21(OH) clump's methanol requires non-thermal desorption, while H2CO and CH3CCH can form thermally.","keywords":["astrochemistry","high-mass star formation","DR21(OH)","methanol","complex organic molecules","carbon-chain molecules","gas-grain chemical modelling","LTE modelling"],"falsifier":"Measure the methanol excitation temperature independently in each component, for example with interferometric observations of a ladder of CH3OH lines that resolve MM1 and MM2; if the fitted temperatures differ from the fixed values enough to bring the thermal warm-up model into agreement with the observed abundances, the need for non-thermal desorption would be overturned.","tokens_in":20473,"feed_emoji":"🧪","tokens_out":5649,"duration_ms":54980,"temperature":0.7,"pith_summary":"Using wide-band single-dish spectra of CH3CCH, CH3OH, and H2CO toward the high-mass star-forming clump DR21(OH), the paper disentangles multiple gas components at clump scale: warm components of roughly 20 to 80 K toward the dusty cores MM1 and MM2, cooler extended envelopes, and a broad blue-shifted outflow. It then feeds the observed column densities into a gas-grain chemical model to ask which formation routes reproduce them. The central finding is that grain-surface chemistry plus thermal desorption can reproduce H2CO and CH3CCH, but CH3OH requires an extra non-thermal desorption mechanism, simulated as a shock, to reach its observed abundance. A sympathetic reader would take this as evidence that clump-scale observations, modelled component by component, can connect core-scale chemistry to parent-cloud structure.","feed_headline":"Methanol needs shocks: DR21(OH) chemistry disentangled","feed_subtitle":"Warm gas can make H2CO and CH3CCH; methanol demands non-thermal desorption.","key_machinery":"The argument is carried by a multi-component local thermodynamic equilibrium (LTE) model that fits every observed line of each species as a sum of Gaussian components, producing a column density, excitation temperature, line width, and velocity for each component; the methanol fits are checked against a non-LTE radiative transfer calculation. The component column densities are converted to fractional abundances using an H2 column density map and then fed into a three-phase gas-grain chemical model that evolves the gas, grain surface, and grain mantle under cold-cloud, warm-up, and simulated-shock stages. The load-bearing comparison is between the model-predicted and observed abundances at matched temperature and density conditions.","core_discovery":"The paper claims that in DR21(OH) the observed methanol cannot be made by warm gas alone: a gas-grain model with a cold cloud, a warm-up phase, and a shock-like burst of non-thermal desorption reproduces all five methanol components, while a plain warm-up model fails. For H2CO the same warm-up model succeeds in about 100 to over 1000 years depending on the component, and for CH3CCH it succeeds only in a low-extinction, cloud-edge environment with visual extinction around 1 to 5 magnitudes, where grain-surface hydrogenation dominates rather than warm carbon chain chemistry. The LTE decomposition that underlies these abundances finds roughly 80 K inner-core components and roughly 20 to 40 K envelope components near MM1 and MM2, plus cold 6 to 15 K gas toward the neighbouring clumps N40, N36, N41, N38, and N48.","pith_inferences":["A sharper test would be to map CH3OH with an interferometer and fit excitation temperatures per core; if methanol's true excitation temperatures are lower than the fixed values, the required abundance would drop and thermal routes might suffice.","The same LTE-plus-chemical-model approach could be applied to other clumps in the region, with each LTE component compared against interferometric maps to see whether it corresponds to a spatially distinct structure.","If non-thermal desorption is required, methanol column densities should be enhanced where outflow shocks interact with dense gas, so a targeted map of CH3OH against shock tracers could test this prediction.","The model's success for CH3CCH only at low extinction implies the molecule traces the cloud edge; higher-resolution observations of CH3CCH emission morphology could confirm that it is limb-brightened."],"forward_implications":["In this source, methanol emission from clump scales should be interpreted as shock- or outflow-liberated grain-surface ice, not as a thermal hot-core product.","CH3CCH is not a reliable tracer of warm carbon chain chemistry here; its abundance is set by grain-surface hydrogenation in a low-extinction envelope.","The component decomposition shows that single-dish clump-scale spectra can separate core, envelope, and outflow gas when a multi-component LTE fit is used.","H2CO can be used as a thermal probe of warm gas near the embedded protostars, since its production is reproduced by thermal mechanisms.","The inferred outflow age of about 1.6 times 10^4 years is shorter than the chemical post-shock timescale, so the model requires the methanol to have been liberated recently."],"supporting_citations":[{"why":"Supplies the gas-grain chemical model used to compute abundances under cold-cloud, warm-up, and shock stages.","marker":"[Ruaud et al. 2016]"},{"why":"Provides the non-LTE radiative transfer code used to check whether methanol behaves non-thermally.","marker":"[van der Tak et al. 2007]"},{"why":"Provides the H2 column-density map used to convert LTE column densities into fractional abundances.","marker":"[Cao et al. 2022]"},{"why":"Defines the dense clumps N44, N48, N38, N40, N41, and N36 whose components are modelled.","marker":"[Motte et al. 2007]"},{"why":"Supplies the temperature and density evolution of a C-type shock used to simulate non-thermal desorption for methanol.","marker":"[Palau et al. 2017]"},{"why":"Provides the molecular line catalogues used to identify transitions and run the LTE fits.","marker":"[Müller et al. 2005]"},{"why":"Details the observations and the LTE model formalism that this analysis builds on.","marker":"[Freeman et al. 2023]"}],"fun_headline_variants":["Methanol in DR21(OH) needs shocks, not just warm gas","DR21(OH): Warm gas makes H2CO and CH3CCH, but methanol needs shocks","Non-thermal desorption drives methanol in DR21(OH) clump","Shock-driven methanol: DR21(OH) chemistry disentangled","Methanol's shock requirement: New model of DR21(OH) chemistry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The methanol temperatures used to calculate how much methanol is present were set by hand rather than measured, and the methanol abundance the chemical model must match shifts directly if those temperatures are wrong.","fun_headline_variants_meta":{"raw":{"variants":["Methanol in DR21(OH) needs shocks, not just warm gas","DR21(OH): Warm gas makes H2CO and CH3CCH, but methanol needs shocks","Non-thermal desorption drives methanol in DR21(OH) clump","Shock-driven methanol: DR21(OH) chemistry disentangled","Methanol's shock requirement: New model of DR21(OH) chemistry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000247,"raw_usage":{"total_tokens":1614,"prompt_tokens":1090,"completion_tokens":524,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":706,"completion_tokens_details":{"reasoning_tokens":424}},"tokens_in":706,"tokens_out":524,"duration_ms":5179,"temperature":1.0,"reasoning_tokens":424,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:03:58.113780+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the methanol excitation temperature independently in each component, for example with interferometric observations of a ladder of CH3OH lines that resolve MM1 and MM2; if the fitted temperatures differ from the fixed values enough to bring the thermal warm-up model into agreement with the observed abundances, the need for non-thermal desorption would be overturned.","supporting_citations":[{"cited_title":"2022, ApJ, 927, 106","cited_arxiv_id":null,"evidence_quote":"Provides the H2 column-density map used to convert LTE column densities into fractional abundances."},{"cited_title":"2007, A&A, 476, 1243 Müller, H","cited_arxiv_id":null,"evidence_quote":"Defines the dense clumps N44, N48, N38, N40, N41, and N36 whose components are modelled."},{"cited_title":"2017, MNRAS, 467, 2723","cited_arxiv_id":null,"evidence_quote":"Supplies the temperature and density evolution of a C-type shock used to simulate non-thermal desorption for methanol."},{"cited_title":"Parsec scales of carbon chain and complex organic molecules in AFGL 2591 and IRAS 20126","cited_arxiv_id":"2308.09584","evidence_quote":"Details the observations and the LTE model formalism that this analysis builds on."}],"review_version":1}