{"id":"5197dc4b-c30d-4777-a994-f96b169e345a","arxiv_id":"2505.07562","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"JDISCS shows that mid-infrared emission from C2H2, HCN, and CO2 in 31 protoplanetary disks is generally optically thin, so line luminosities trace molecular mass, and that HCN-to-cold-water ratios are smaller in smooth disks.","lead":"This survey combines JWST mid-infrared spectra of 31 protoplanetary disks with ALMA millimeter images to connect inner disk molecular chemistry to outer disk structure. It finds that organic gas emission is mostly optically thin and that compact smooth disks show a relative surplus of cold water vapor, pointing to inward pebble drift.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Organic T/M/opt-depth claims rest on a single-temperature water subtraction that is acknowledged but never quantitatively tested; residual multi-temperature or non-LTE water could bias the C2H2/HCN/CO2 retrievals.","rationale":"The reader's weakest-assumption identifies exactly the same load-bearing point: the organic retrievals depend on subtracting a single-temperature LTE water model, and residual water could bias the organic temperatures and optical-depth conclusions. My reading of the manuscript confirms this is the most consequential vulnerability. The optical-thinness evidence based on the luminosity–mass correlation is partly circular because both quantities come from the same slab fits, but the paper has independent support from 13C-isotopologue non-detections and the pedestal analysis, so the circularity alone is not fatal. The water-subtraction issue is broader: it affects essentially every organic retrieval (T, N, M, τ) and is acknowledged by the authors but never quantitatively tested. The CO2 detection-threshold ambiguity (visual inspection of four disks) and the non-significant HCN/cold-H2O pebble-drift trend (p=0.11) are real but secondary; they do not undermine the central optical-thinness claim as directly. Given the explicit caveats in Appendix A and Section 5.2, the paper is not misleading, but the central claims would be on firmer footing with a targeted test of water-subtraction systematics. Hence the CONDITIONAL verdict stands unchanged.","tokens_in":49140,"tokens_out":5669,"duration_ms":57228,"concrete_test":"Re-fit the organic molecules for a representative subset (e.g., AS 205N, CI Tau, and the seven disks with strong residual fringing/non-LTE water noted in Section 3.3.2) using a two-temperature water model, or by simultaneously fitting water (with hot-band suppression factors as in Banzatti et al. 2025) together with C2H2, HCN, and CO2. If the retrieved median C2H2/HCN/CO2 temperatures shift by more than the quoted 1σ uncertainties, or if the L–M slopes in Figure 13 change significantly, then the single-temperature water subtraction is biasing the central optical-thinness and temperature claims. Alternatively, mask all channels within ±0.02 µm of known residual water lines after single-temperature subtraction and re-run the organic fits on AS 205N; a >1σ change in the retrieved C2H2/HCN column densities or temperatures would confirm water contamination.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that C2H2 and HCN emission is optically thin and that these molecules are hotter than CO2 is derived entirely from slab-model fits to the 12–16 µm region after subtracting a single-temperature LTE water model (Section 3.3.1, Appendix A). The paper explicitly states that the residual water emission includes non-LTE signatures and multiple temperature components that the single-temperature model cannot reproduce (Appendix A), and that CO2 'may be particularly dependent on the number of temperature components used to model the H2O' (Section 3.3.5). The C2H2 (13.7 µm) and HCN (14.0 µm) Q branches sit in a wavelength range dense with rotational water lines, so any unsubtracted water could be absorbed into the organic fits, changing the retrieved temperatures, column densities, and inferred optical depths. Because the optical-thinness conclusion and the temperature ordering (920/820 K vs 600 K) are both taken from these same fits, a systematic water-subtraction bias would propagate directly into the headline results. The paper offers a qualitative justification for using a single temperature but does not demonstrate stability, e.g., by comparing against two-temperature water fits or by masking residual water lines. Independent support from 13C-isotopologue non-detections partially mitigates the concern for C2H2 and CO2, but HCN lacks a comparable isotopologue constraint in most disks, and no such test addresses the temperature retrievals.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first Cycle 1 sample analysis of the JDISCS survey: MIRI-MRS spectra of 31 protoplanetary disks with ALMA millimeter continuum imaging. The authors report near-ubiquitous H2O, OH, and CO emission, higher molecular detection rates than Spitzer-IRS, and LTE slab-model retrievals of temperature, column density, and emitting area for C2H2, HCN, and CO2 in the 12–16 µm region after subtracting a single-temperature water model. They argue that C2H2 and HCN (and possibly CO2) emission is optically thin, so the observations constrain total molecular mass rather than column density or radius separately, and that C2H2 and HCN typically arise in hotter gas than CO2. They also examine correlations with outer disk substructure, finding a tentative trend between HCN/cold-H2O luminosity ratios and dust disk size, and they show that compact-disk 'pedestals' under the organic Q branches are not explained by optically thick C2H2 models. The paper closes with atomic lines ([Ne II], [Ne III], [Ar II]) and their relation to accretion and radiation fields.","tokens_in":49430,"tokens_out":4658,"duration_ms":48680,"significance":"If the conclusions hold, this is a valuable demographic study: it provides a homogeneous, publicly released MIRI-MRS dataset for a sample with high-quality ALMA imaging, demonstrates sensitivity gains over Spitzer, and makes a physically important claim that the inner-disk organic emission is mainly optically thin, thereby lifting the column-density/area degeneracy for total-mass measurements. The paper is commendably explicit about its limitations: it openly discusses the single-temperature water subtraction, the exclusion of flux errors from the loss function, the partial degeneracy of N and r_slab, and the moderate statistical significance of some trends. The isotopologue non-detections and the direct comparison with optically thick models in Section 5.2 provide independent support for parts of the central argument. However, as detailed below, some load-bearing points need additional quantitative support before the headline claims can be considered established.","major_comments":[{"comment":"The single-temperature LTE water subtraction is acknowledged to leave residuals that include non-LTE signatures and multiple temperature components (Appendix A), and Section 3.3.5 states that CO2 may be particularly dependent on the number of temperature components used to model the H2O. Yet no quantitative test of the impact of this subtraction on the organic retrievals is provided. Because the C2H2 (13.7 µm) and HCN (14.0 µm) Q branches sit in a wavelength region dense with rotational water lines, unsubtracted or mis-subtracted water could be absorbed into the organic slab fits and shift the retrieved temperatures, emitting masses, and inferred optical depths; this would directly propagate into the headline claims of optical thinness and of C2H2/HCN being hotter than CO2. Please add a sensitivity test for at least a subset of sources, such as fitting a two-temperature water model, masking water-dominated spectral regions before fitting the organics, or adding a non-LTE water correction factor, and report how T, M, L, and the derived optical depths change under those alternatives.","section":"Section 3.3.1 and Appendix A"},{"comment":"The near-unity slopes of log-L versus log-M for C2H2 and HCN are presented as evidence that the emission is optically thin, but both quantities are outputs of the same slab-model fits and are proportional by construction in the optically thin regime (L is integrated from the model spectrum, and M is derived from the same column density and area, M ∝ N A and L ∝ N A). The set of degenerate solutions within 5% of the minimum L2 norm preserves these products, so the observed slope near 1 is at least partly a consequence of the model construction rather than an independent physical measurement. The isotopologue non-detections and the τ versus N behavior in Figure 15 are more probative for optical thinness. Please either remove the L-M correlation as an argument for thinness, or test it against a null model (e.g., random draws of N and r_slab over the degeneracy set) and state explicitly what the correlation does and does not establish.","section":"Section 4.1, Equation (2), and Figure 13"},{"comment":"The abstract states that the HCN to cold-H2O luminosity ratios are generally smaller in smooth disks, but the reported Spearman correlation is ρ = 0.36 with p = 0.11, which is not significant at the p < 0.05 threshold used elsewhere in the paper (Section 4.1). The paper itself describes the trend as only marginally statistically significant, but this caveat is not carried into the abstract or the conclusions. Please report this as a tentative trend in the abstract and conclusions, or provide additional statistical support, such as a permutation test or a comparison restricted to the smooth versus structured subsamples, and state the resulting p-value.","section":"Section 5.1 and Figure 17"},{"comment":"The detection rates are a headline result, but the threshold for a CO2 detection is not fully reproducible: Section 3.3.5 states that four disks are added on the basis of visual inspection of the Q branch, because any luminosity or peak-to-continuum criterion that includes them also captures non-detections. Since the reported detection rates for CO2 (58% overall) are compared to Spitzer rates and used in subsequent demographic statistics, please provide an explicit decision rule for the visual-inspection step (e.g., independent inspection by two authors, or a blind test) and quantify how the detection rate changes under alternative objective thresholds. A similar concern applies to the 'rough detection threshold' for C2H2 and HCN based on log(L) ~ −5.1 L_sun, which is not tied to a noise or model-significance estimate.","section":"Sections 3.3.3 and 3.3.5, and Tables 5–7"}],"minor_comments":[{"comment":"The L2 loss function omits flux uncertainties, and the 5% threshold on the L2 norm is described as 'consistent to within roughly 2σ' without a derivation; please clarify that this is not a formal goodness-of-fit or detection significance and state the number of slab models used per source for the medians reported in Tables 5–7.","section":"Equation (1) and Section 3.3.1"},{"comment":"The text referencing 'Equation 1' for the linear regressions in Figure 13 should refer to Equation (2), and the earlier cross-reference to 'Section 3.2.2' for the slab-model fits should be Section 3.3.2.","section":"Section 5.1 and Figure 13"},{"comment":"The appendix would benefit from a single summary figure or table showing, for a few representative disks, the residuals before and after the single-temperature water subtraction in the 13.4–14.1 µm region; this would make the magnitude of the residual water contamination more transparent than the current qualitative discussion.","section":"Appendix A"},{"comment":"The footnote for Tables 5–7 says 'Since logN and r_slab are degenerate, we do not include them in the analysis,' but the tables do list representative logN and r_slab values; please state explicitly which values are shown (e.g., one member of the degenerate family used for the figures) and how the plotted model was chosen.","section":"Tables 5–7"},{"comment":"The ordering of sources by sub-mm dust disk size in Figure 9 is useful, but the two-panel layout with temperature scales repeated on the left and right makes the figure dense; a single aligned panel with shared source labels would improve readability.","section":"Figure 9"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for the journal and is likely to be an important dataset paper, but the central physical claims currently rest on a modeling chain whose weakest link—the single-temperature water subtraction—is acknowledged but not quantitatively stressed. The requested sensitivity tests and a clearer separation of model-internal correlations from independent evidence for optical thinness are essential before publication. I do not see concerns about data provenance or duplication; the authors have been commendably transparent about limitations, but transparency alone does not establish that the systematic biases are small."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — the short version: this is the first large uniform MRS sample that ties inner disk molecular emission to ALMA outer disk structure, and as a reference dataset it will get cited. The detection rates (higher than Spitzer), the median slab temperatures for C2H2 (920 K) and HCN (820 K) versus CO2 (600 K), and the sample-wide optically thin conclusion for the organics are all new. Credit where it's due: uniform reduction, public code and data, honest handling of the N–r degeneracy by working with emitting mass and luminosity, and independent checks (13C12CH2 non-detections, failure of optically thick models to reproduce the pedestals). The paper openly acknowledges its main limitations.\n\nThe soft spots are real but not fatal. First, the L–M correlation that supports optical thinness is partly model-internal: both axes come from the same slab fits, so a slope near unity is partly built in. The isotopologue non-detections and the weaker H2O correlation provide independent support, so I think the conclusion is probably right, but the paper should own the circularity more directly. Second, the single-temperature LTE water subtraction is the bigger issue. The C2H2 and HCN Q branches sit in a water-dense region, and Appendix A admits the residuals contain non-LTE and multi-temperature signatures and that CO2 may be especially sensitive. The paper gives a qualitative justification but no quantitative test, such as two-temperature water fits or masking the worst water lines. This could systematically shift the retrieved temperatures and, in principle, the optical depth conclusions. The isotopologue constraints help for C2H2 and CO2 but not for HCN temperatures. Third, the headline HCN/cold-H2O versus disk-size trend has p=0.11; the paper calls it marginal, which is fair, but the abstract's \"consistent with\" framing is a bit generous.\n\nWho is this for? Anyone working on disk chemistry demographics or JWST MIRI surveys. It deserves a serious referee. My recommendation: send it to review, and require a robustness test on the water subtraction (even a simple two-component comparison for a few disks) or a softened claim on the temperature ordering and optical thinness. The dataset stands on its own.","headline":"Valuable legacy dataset for JWST disk chemistry, but the headline pebble-drift trend is weak (p=0.11) and the optical-thinness argument carries a built-in correlation plus a water-subtraction caveat; still warrants serious peer review.","tokens_in":50108,"tokens_out":2466,"would_cite":true,"duration_ms":23190,"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":"This paper claims that emission from the organic molecules C2H2, HCN, and possibly CO2 in the inner regions of planet-forming disks is optically thin, so the observed mid-infrared line luminosity directly measures the total molecular mass…","keywords":["protoplanetary disks","JWST MIRI spectroscopy","mid-infrared molecular emission","disk chemistry","optically thin gas","pebble drift","mass accretion rate","ALMA disk substructures"],"falsifier":"Detect the 13C12CH2 isotopologue Q-branch near 13.7 microns in a C2H2-bright disk such as AS 205N: at the column densities the optically thin models require, an interstellar 12C/13C ratio near 70 predicts a readily detectable isotopologue line, so a clear non-detection would falsify the large-column-density solutions and the derived masses, while a detection at the predicted strength would confirm the mass interpretation.","tokens_in":48898,"feed_emoji":"🔭","tokens_out":5899,"duration_ms":53152,"temperature":0.7,"pith_summary":"This paper analyzes JWST mid-infrared spectra of 31 protoplanetary disks to ask what drives the chemical diversity of the gas within a few astronomical units of the star, where planets assemble. It finds that emission from acetylene (C2H2), hydrogen cyanide (HCN), and likely carbon dioxide (CO2) is optically thin, which breaks the usual degeneracy between gas amount and emitting area: the observed line luminosity is a direct measure of the total molecular mass. On that basis, it reports that C2H2 and HCN sit in hotter gas (about 920 K and 820 K) than CO2 (about 600 K), that molecular line strengths track mass accretion rate and inner-disk clearing, and that small smooth disks have less HCN relative to cold water, a signature of icy pebble drift delivering water. A sympathetic reader would care because it converts a varied collection of spectra into a demographic map linking inner-disk chemistry to outer-disk structure, the raw material for predicting planetary compositions.","feed_headline":"JWST finds inner-disk organic gas is mostly transparent","feed_subtitle":"In 31 planet-forming disks, C2H2, HCN and CO2 brightness now measures total molecular mass.","key_machinery":"The load-bearing tool is the LTE slab model, a homogeneous gas layer described by temperature, column density, and projected emitting area, fitted to the 12-16 micron spectral region with Markov chain Monte Carlo after first subtracting a single-temperature water model. Its power here is that the fits return families of degenerate (column density, radius) solutions whose integrated luminosity and total mass are nearly invariant; since the luminosity-mass correlation is close to linear, the emission is optically thin, and the degeneracy converts from a flaw into a mass measurement. The ALMA sub-millimeter continuum images supply the outer-disk context, namely dust disk sizes and ring or cavity substructures, against which the inner-disk chemistry is compared.","core_discovery":"The paper's central claim is that the organic molecules C2H2, HCN, and probably CO2 emit as optically thin gas in the inner disks of T Tauri stars, so the strong degeneracy between column density and emitting radius in slab-model fits is not a nuisance but a route to a physical quantity: total molecular mass. Fits to the 12-16 micron spectra show luminosity scaling nearly linearly with emitting mass (slope near 1 for C2H2 and HCN), as expected for optically thin emission; CO2 has a shallower slope and a cooler median temperature, suggesting it is marginally thicker. The survey also establishes higher detection rates than Spitzer for all molecules, including first detections in transition disks, and reports demographic trends: molecular luminosities correlate with accretion rate, anti-correlate with the 13-26 micron infrared spectral index, and HCN-to-cold-water ratios are low in compact, smooth disks.","pith_inferences":["Because optical thinness makes line luminosity a mass tracer, future surveys could use C2H2 and HCN luminosities to map the carbon and nitrogen budgets of inner disks as functions of stellar mass and age, not just accretion rate.","The paper's interpretation of the pedestal features as dust rather than optically thick gas predicts that high-resolution 10 micron silicate spectroscopy of compact disks should show matching dust composition features, which is a testable extension.","If CO2 is only marginally optically thin, its 13CO2/12CO2 isotopologue ratio, currently detected in only one JDISCS source, offers a direct route to break the residual CO2 column-density degeneracy.","The authors' single-temperature water subtraction leaves open the possibility that a second water component biases organic retrievals; extending the fitting to two water components simultaneously is a natural next step."],"forward_implications":["Mid-infrared organic line luminosities can be used directly as relative measures of inner-disk molecular mass, removing the need to know the emitting radius.","Demographic comparisons with ALMA substructures become meaningful, because disks with dust rings, cavities, or spiral arms can be placed on an equal mass footing with smooth disks.","The correlation of C2H2, HCN, and CO2 luminosity with mass accretion rate implies that accretion-generated ultraviolet heating is a primary controller of organic line brightness.","The low HCN-to-cold-water ratios in compact smooth disks support pebble drift as an efficient water delivery mechanism in the absence of dust traps.","All measured [Ne III]/[Ne II] ratios below unity favor X-ray, not extreme-ultraviolet, irradiation as the driver of disk winds in this sample."],"supporting_citations":[{"why":"Supplies the Spitzer-IRS detection-rate baseline against which the JDISCS detection rates are compared.","marker":"Pontoppidan et al. 2010"},{"why":"Establishes the slab-model retrieval approach and the column-density/radius degeneracy that this paper reinterprets as a mass measurement.","marker":"Carr & Najita 2011"},{"why":"Provides the IRS slab-model parameter space and the prior treatment of the same degeneracy via flux scaling.","marker":"Salyk et al. 2011b"},{"why":"Defines the line opacity-overlap treatment and the optically thick C2H2 pseudo-continuum benchmark against which JDISCS disks are tested.","marker":"Tabone et al. 2023"},{"why":"Supplies the two-temperature water picture and pebble-drift interpretation that the HCN-to-cold-water trend relies on.","marker":"Banzatti et al. 2023a"},{"why":"Provides the water-line models, Keplerian broadening treatment, and cold-water luminosities used in the ratio analysis.","marker":"Banzatti et al. 2025"},{"why":"Provides the DSHARP ALMA images from which the outer-disk dust substructures and disk sizes are taken.","marker":"Andrews et al. 2018"},{"why":"Gives the Spitzer-era correlations between molecular luminosity, accretion rate, and infrared spectral index that this paper recovers with JWST.","marker":"Banzatti et al. 2020"}],"fun_headline_variants":["Transparent organics unveil true gas mass in planet-forming disks","Optically thin C2H2 and HCN reveal disk gas mass","See-through organic gas in disk cores reveals its true mass","Inner-disk organics are see-through: mass measured"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The single-temperature LTE water model is subtracted before the organic fits, and the residuals are assumed not to bias the retrieved organic temperatures and optical depths; if hot-band or non-LTE water emission leaks through, the optically thin conclusion and median temperatures could shift.","fun_headline_variants_meta":{"raw":{"variants":["Transparent organics unveil true gas mass in planet-forming disks","Optically thin C2H2 and HCN reveal disk gas mass","See-through organic gas in disk cores reveals its true mass","Inner-disk organics are see-through: mass measured"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001302,"raw_usage":{"total_tokens":5390,"prompt_tokens":1106,"completion_tokens":4284,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":722,"completion_tokens_details":{"reasoning_tokens":4221}},"tokens_in":722,"tokens_out":4284,"duration_ms":25877,"temperature":1.0,"reasoning_tokens":4221,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:13:31.195362+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Detect the 13C12CH2 isotopologue Q-branch near 13.7 microns in a C2H2-bright disk such as AS 205N: at the column densities the optically thin models require, an interstellar 12C/13C ratio near 70 predicts a readily detectable isotopologue line, so a clear non-detection would falsify the large-column-density solutions and the derived masses, while a detection at the predicted strength would confirm the mass interpretation.","supporting_citations":[],"review_version":1}