{"id":"6dfef02c-7b42-40e0-91c9-faafa18bb418","arxiv_id":"2411.09011","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"Cold water vapor mass should jump at a transitional CO mass, providing a new diagnostic for whether debris disk gas is primordial or collisionally produced.","lead":"This paper proposes measuring cold water vapor in debris disks to tell whether the gas there is leftover from the protoplanetary disk or produced by collisions of icy bodies. It derives order-of-magnitude predictions and argues that a future far-infrared space telescope is needed to test them.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 10's local balance is marginal at the secondary branch's transition point and fails at the high-CO end of the primordial branch; the Appendix A 'or MCO > 1e-3' criterion conflates the two branches' mass ratios, so the predicted H2O jump may be suppressed.","rationale":"The reader's weakest assumption identifies exactly the condition on which the predicted jump rests: the local balance Eq. (10) must hold in the parameter space where the jump is drawn. I agree with that assessment and add a sharper technical point: the Appendix A validity criterion is stated as 'MZ,ref > 0.1 M⊕ (or MCO > 1e-3 M⊕)', but the equivalence between the two limits depends on the MCO/Msolid ratio, which differs by an order of magnitude between the two hypotheses. For the primary branch the relevant threshold is MCO ~ 1e-2 M⊕, so the red curve is over-extended in Figure 3; for the secondary branch the threshold is MCO ~ 1e-3 M⊕, placing the endpoint of the orange curve exactly at the boundary where freezeout can start to dominate. This does not make me reject the paper: the qualitative idea that cold water vapor may differentiate the two gas origins remains plausible, and the authors are transparent about the order-of-magnitude nature of the estimates. It does, however, strengthen the case for a conditional verdict, because the quantitative height of the jump, and hence the required FIR sensitivity, depends on the validity of Eq. (10) near the transition. The concrete test using existing observations directly probes that boundary. If the test shows the targets near the transition satisfy the timescale inequalities, the concern would be alleviated; if not, the predicted jump would be diminished and the mission sensitivity estimates in Appendix C would need revision.","tokens_in":14837,"tokens_out":11886,"duration_ms":120025,"concrete_test":"Using the observed targets in Figure 2 (compiled in Rebollido et al. 2022 and references) that have both CO and dust mass measurements and MCO within a factor of 3 of 1e-3 M⊕, compute MZ,ref from the compiled solid masses. For each target, evaluate tau_coll and tau_free using Eqs. A.2 and A.5 with the target's MZ,ref and compare to tau_deso ~ tau_diss ~ 3.2e1 yr. If any target with MCO below the transition has tau_free < tau_diss or tau_coll < tau_deso, then the orange-branch prediction at the jump location is not valid for that target, and the Monte Carlo in Appendix C should be rerun with those targets' H2O masses suppressed to test whether the non-monotonic signature survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central diagnostic is the factor-100 step in MH2O at MCO ~ 1e-3 M⊕, produced by switching between Eq. (12) and Eq. (14). Both branches rely on Eq. (10), the steady-state balance between photodesorption and photodissociation, which requires freezeout and collision timescales to be slower than the destruction and desorption timescales. Appendix A shows this holds for MZ,ref < 0.1 M⊕, but the paper phrases the limit as 'MZ,ref > 0.1 M⊕ (or MCO > 1e-3 M⊕)'. The 'or' is only valid for the secondary branch, where MCO/Msolid ~ 1e-2 gives MZ,ref ~ 100 MCO, so MZ,ref crosses 0.1 M⊕ at MCO ~ 1e-3 M⊕. For the primordial branch, MCO/Msolid ~ 0.1 (Eq. 5) gives MZ,ref ~ 10 MCO, so the threshold is MCO ~ 1e-2 M⊕, not 1e-3 M⊕. Thus the red branch in Figure 3 is drawn with the local-balance assumption up to MCO ~ 0.1 M⊕, but for MCO > 1e-2 M⊕ the assumption fails and cold water vapor could be strongly depleted. More importantly, the orange branch is drawn right up to the transition point MCO = 1e-3 M⊕, where MZ,ref = 0.1 M⊕; a modest increase in solid mass would make freezeout faster than photodissociation (Eq. A.5 vs A.3), suppressing the H2O mass and reducing the height of the very jump the method is designed to detect. The paper itself flags this limitation in the Summary: 'if debris disks are massive enough (i.e., MZ,ref > 0.1 M⊕), ... the resulting water vapor abundance could be very low.' Because the jump is the only non-monotonic feature carrying the statistical signal, this boundary sensitivity is load-bearing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the mass of cold water vapor, plotted against CO gas mass, can diagnose whether gas in debris disks originated as residual primordial protoplanetary gas (Hypothesis 1) or as collisionally produced secondary gas (Hypothesis 2). The authors perform order-of-magnitude calculations showing that for both hypotheses the cold water abundance is set by a photodesorption-photodissociation balance (Eq. 10), and from this they derive H2O/CO conversion factors of about 1e-5 for Hypothesis 1 (Eq. 12) and 1e-3 for Hypothesis 2 (Eq. 14). Switching between these conversions at a transition CO mass near 1e-3 M_Earth yields a non-monotonic jump in predicted H2O mass (Fig. 3). The paper compares this prediction tentatively with TW Hya and beta Pic, argues that current facilities cannot detect cold water vapor in the far-IR, and uses a Monte Carlo exercise (Appendix C) to estimate required sample sizes and sensitivities for a future FIR observatory.","tokens_in":15353,"tokens_out":3760,"duration_ms":33878,"significance":"If the predicted non-monotonic jump is robust, the paper offers a genuinely useful observational strategy: it identifies a specific molecule, a specific wavelength range, and a quantitative statistical target for distinguishing primordial from secondary gas in debris disks. The analytical transparency is a strength: the assumptions are stated, the order-of-magnitude algebra is straightforward, and the main caveat (validity of the local balance at high solid mass) is explicitly acknowledged in the Summary. The proposal also connects to a concrete mission context, which makes it actionable for future instrument design. However, the significance is conditional on the survival of the jump under the stated validity limits, and on the independence of the calibration from the test, both of which need strengthening before the central claim is fully supported.","major_comments":[{"comment":"The predicted non-monotonic jump relies on Eq. (10) being valid over the entire CO-mass range shown for both branches. Appendix A states the validity condition as MZ,ref > 0.1 M_Earth (or MCO > 1e-3 M_Earth), but this 'or' is only correct for Hypothesis 2. For Hypothesis 1, Eq. (5) gives MCO/Msolid ~ 0.1, so MZ,ref ~ 10 MCO and the local-balance assumption breaks down at MCO ~ 1e-2 M_Earth, not at 1e-3 M_Earth. Consequently, the red branch in Fig. 3 is drawn using Eq. (12) in a regime where freezeout and collisional timescales may be shorter than photodesorption and photodissociation, and the true water vapor abundance there could be much lower. At the same time, the orange branch is evaluated right at the transition MCO = 1e-3 M_Earth, where MZ,ref = 0.1 M_Earth, so the jump height is marginal under the paper's own criterion. This does not invalidate the qualitative idea, but the quantitative factor-of-100 jump presented as the diagnostic signal is not robust within the stated validity limits. The authors should recompute the predicted H2O mass in the parameter space where Eq. (10) holds, or clearly mark the breakdown region and discuss how the jump would be modified.","section":"Appendix A and Fig. 3"},{"comment":"The normalization fdust = 1e-2 is chosen to better reproduce the result of Du & Bergin (2014) and the cold water vapor detection toward TW Hya. Presenting the TW Hya point in Fig. 3 as a tentative test of the prediction is therefore partially circular: the same observation is used to set the free parameter that fixes the vertical normalization of the Hypothesis 1 branch. To make the comparison a genuine validation, the authors should either treat fdust as a fitted parameter and exclude TW Hya from the test, or provide an independent constraint on fdust (e.g., from dust growth and settling models) so that the TW Hya comparison is not simultaneously calibrating and checking the model.","section":"Section 2.4, Eq. (11) and Fig. 3 left"},{"comment":"The switch between the two hypotheses is placed at MCO = 1e-3 M_Earth, which the text introduces as 'as an example'. The true transition CO mass is an unknown function of disk age, stellar type, and collisional history. The Monte Carlo in Appendix C adopts this same value (with smooth transitions up to 5e-4 M_Earth), but the principal observable signature—the location and even the existence of the jump—depends on this choice. The paper should justify the transition mass more carefully and show how the predicted jump would appear if the transition occurs over a range or at a different mass; otherwise the non-monotonicity is partly assumed rather than derived.","section":"Section 2.4, Eq. (1), and Appendix C"}],"minor_comments":[{"comment":"The sentence 'it should be noted that in protoplanetary disks, solids may be in the form of dust as well as much larger bodies such as planetesimals and (proto)planets. Accordingly, our assumption (i.e., Msolid ≃ Msolid) in equation (3)...' contains an apparent typo: 'Msolid ≃ Msolid' should presumably be 'Msolid ≃ Mdust' or 'Msolid ≃ MZ,ref'.","section":"Section 2.2, Eq. (5)"},{"comment":"The phrase 'tentative constrains' should be 'tentative constraints'.","section":"Fig. 3 caption"},{"comment":"The manuscript alternates between 'dust' and 'solids' in defining Msolid; for clarity, specify explicitly whether Msolid in Eqs. (3), (5), and (9) includes only dust or also larger bodies, since the derived upper limit in Eq. (5) depends on this distinction.","section":"General notation"},{"comment":"The phrase 'where MZ,ref is estimated at MCO,tran (Figure 2)' is ambiguous: it should state that MZ,ref is evaluated using the mass ratio implied by each hypothesis at the quoted CO mass, since the conversion differs by an order of magnitude between Hypotheses 1 and 2.","section":"Appendix A, after Eq. (A.2)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a clearly written order-of-magnitude proposal with a sound qualitative core, but the central quantitative claim—the factor-of-100 jump in H2O mass at MCO ~ 1e-3 M_Earth—is weakened by the boundary of the local-balance assumption occurring exactly at the transition for Hypothesis 2 and at MCO ~ 1e-2 M_Earth for Hypothesis 1. The authors already acknowledge the catastrophic effect of massive disks in the Summary, which is commendable, but the acknowledgment is not incorporated into the figures or the Monte Carlo forecasts. I would ask the authors to re-derive the predicted jump in the valid parameter space and to make the calibration of fdust explicit and non-circular. These are fixable within the scope of the manuscript, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the paper gives the field something it didn't have—a concrete, physically motivated discriminant between the two leading explanations for gas in debris disks—and it does it with transparent order-of-magnitude arithmetic. The H2O-to-CO ratios differ by about two orders of magnitude between the primordial and collisional branches, and the authors correctly note that CO alone cannot separate the hypotheses. That is a genuine contribution, and it gives a future FIR observatory a focused science case.\n\nWhat the paper does well: the balance in Eq. (10) is standard photodesorption versus photodissociation, and the authors are unusually honest about their parameter choices. Table 1 shows how wide the possible ranges are, and the Monte-Carlo appendix at least gestures at the sample size needed. For an order-of-magnitude proposal, the reasoning is clear and reproducible.\n\nThe soft spots are real but not fatal. The normalization for Hypothesis 1 is set by fdust = 1e-2, chosen to match the TW Hya water detection; later the same detection is listed as a point of agreement. That is circular in the strict sense, and the paper should say it is calibrating rather than testing. The stress-test note on Appendix A also holds: the statement that the local balance fails when MZ,ref > 0.1 M⊕ (or MCO > 1e-3 M⊕) applies the \"or\" only to the secondary branch, where MCO/Msolid ~ 1e-2. On the primordial branch the same solid-mass threshold corresponds to MCO ~ 1e-2, not 1e-3, so the red curve in Fig. 3 extends beyond its validity domain. At the transition point itself the balance is marginal, so the height of the jump is sensitive to exactly the parameter regime the diagnostic targets. The authors flag this in the Summary but do not carry it back into the figures.\n\nThe qualitative conclusion—that cold water vapor should be far more abundant per unit CO in secondary gas—likely survives. But the sharp non-monotonic jump is, to some degree, built in by the assumed switching mass and the calibrated normalization. I would like the revised version to mark the validity boundaries in Fig. 3, separate calibration from verification, and frame the diagnostic as statistical rather than per-disk.\n\nRecommendation: send to a serious referee. This paper deserves careful review, not a desk reject. With those edits it could be a useful reference for both modelers and mission planners.","headline":"A fresh, plausible order-of-magnitude proposal for distinguishing primordial from collisional gas in debris disks via cold water vapor; the sharp predicted jump is partly built in, but the core idea deserves a serious review.","tokens_in":15914,"tokens_out":3603,"would_cite":true,"duration_ms":33605,"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":"Cold water vapor can tell whether a debris disk's gas is primordial or collisionally produced.","keywords":["debris disks","cold water vapor","CO gas","photodesorption","primordial gas","secondary gas","far-infrared observatory","protoplanetary disks"],"falsifier":"Measure $\\mathrm{H_2O}$ and CO masses in roughly two dozen debris disks spanning CO masses from below $10^{-5}$ to above $10^{-1}\\,\\mathrm{M}_\\oplus$ with a future far-infrared observatory; the diagnostic fails if the $\\mathrm{H_2O}$-mass-versus-CO-mass relation is smoothly monotonic with no jump near $10^{-3}\\,\\mathrm{M}_\\oplus$, or if disks with $\\sim 10^{-1}\\,\\mathrm{M}_\\oplus$ of CO show water masses below roughly $10^{-7}\\,\\mathrm{M}_\\oplus$, well under the primordial-branch prediction.","tokens_in":14630,"feed_emoji":"💧","tokens_out":12157,"duration_ms":95345,"temperature":0.7,"pith_summary":"Debris disks around young stars sometimes contain gas, but CO observations alone cannot tell whether that gas is leftover from the protoplanetary disk or was generated recently by collisions between icy bodies. This paper argues that cold water vapor breaks the degeneracy. Under both hypotheses, cold water vapor is produced by the same mechanism—photodesorption of water ice from dust grains by attenuated ultraviolet light—so the underlying physics is unified, yet the predicted $\\mathrm{H_2O}$-to-CO mass ratio differs by two orders of magnitude: about $10^{-5}$ for primordial gas and $10^{-3}$ for collisionally produced gas. Switching between the two branches at a CO mass near $10^{-3}\\,\\mathrm{M}_\\oplus$ turns the water-mass-versus-CO-mass relation into a non-monotonic jump that statistical surveys can identify. The obstacle is observational: cold water emits in the far-infrared, so the test requires a future space observatory with roughly a hundred times Herschel's sensitivity.","feed_headline":"Water vapor can reveal the origin of debris-disk gas","feed_subtitle":"The two gas origins give water-to-CO ratios a hundredfold apart, so a far-infrared survey could see the jump.","key_machinery":"The load-bearing object is the photodesorption–photodissociation balance, Eq. (10): $\\sigma_{\\rm photo} N_{\\rm H_2O,vap} = \\sigma_{\\rm dust} Y N_{\\rm dust}$, where photodesorption (the ejection of water molecules from dust by ultraviolet photons) is balanced against photodissociation (their destruction by the same radiation). Because both rates scale with the UV flux, the column density of cold water vapor becomes independent of the UV intensity and is set instead by the dust surface area available—which is exactly where the two hypotheses differ. The argument then reduces to the ratio identities (12) and (14) with a fixed switch point at $\\mathrm{M_{CO}} \\sim 10^{-3}\\,\\mathrm{M}_\\oplus$, producing the non-monotonic jump that serves as the observable signature.","core_discovery":"The paper's central claim is that the mass of cold water vapor in a debris disk, plotted against the disk's CO mass, encodes the origin of the gas. For residual primordial gas, Eq. (12) gives $\\mathrm{M_{H_2O}/M_{CO}} \\simeq 10^{-5}$; for collisionally produced gas, Eq. (14) gives $\\simeq 10^{-3}$. Both estimates come from Eq. (10), the local balance between photodesorption and photodissociation, evaluated with the dust population appropriate to each hypothesis: 0.1 $\\mu$m grains at 1% of the solid mass in the primordial case, and blowout-limited 1 $\\mu$m grains at full cometary abundance in the collisional case. Because the transition CO mass is fixed near $10^{-3}\\,\\mathrm{M}_\\oplus$ by photoevaporation (Eq. 1), the predicted $\\mathrm{H_2O}$ mass and $\\mathrm{H_2O}$-to-solid ratio both show a jump at that point. The paper argues this non-monotonic feature can be detected statistically with roughly ten to thirty targets, and shows that the one available detection (TW Hydrae) and the one strong upper limit ($\\beta$ Pictoris) are consistent with the predicted range.","pith_inferences":["Inference: if the predicted jump is observed, the same data could calibrate the transition CO mass empirically rather than taking it from the photoevaporation estimate, linking gas origin to disk dispersal history.","Inference: in systems where stellar UV radiation dominates and water's lifetime shortens to days, OH emission—water's photodissociation product—could serve as an indirect probe, an extension the paper mentions in its outlook.","Inference: because the primordial branch's normalization rests on a 1% dust fraction calibrated to a single disk, a multi-target survey would test whether that fraction is universal or evolves with disk age and dust processing.","Inference: a null result—no water and no jump—would shift attention to alternative secondary-gas mechanisms such as cometary evaporation or radiogenic degassing, which the paper lists but does not model."],"forward_implications":["With a future far-infrared observatory, surveying roughly 10–30 gas-bearing debris disks at $\\mathrm{H_2O}$ mass sensitivity $10^{-9}{-}10^{-8}\\,\\mathrm{M}_\\oplus$ and CO sensitivity $10^{-4}\\,\\mathrm{M}_\\oplus$ should be enough to see the predicted jump.","Disks with CO mass above about $10^{-3}\\,\\mathrm{M}_\\oplus$ should populate the primordial branch, while lower-CO disks should populate the collisional branch, so the water-vapor census doubles as an evolutionary tracer for planet-forming disks.","Detected cold water vapor in CO-rich debris disks is expected at masses around $10^{-7}{-}10^{-6}\\,\\mathrm{M}_\\oplus$, comparable to the TW Hydrae detection, so non-detection at that level in such disks would contradict the primordial branch.","For debris disks with solid masses above roughly $0.1\\,\\mathrm{M}_\\oplus$, the photodesorption balance assumption breaks down and the water abundance could be much lower, marking a boundary on where the diagnostic applies."],"supporting_citations":[{"why":"Provides the only detection of cold water vapor in a protoplanetary disk (TW Hya) and the mass estimate used to calibrate the primordial-branch dust fraction.","marker":"Hogerheijde et al. 2011"},{"why":"Supplies the outer-disk water-vapor model whose results the normalization of Eq. (11) is chosen to match.","marker":"Du & Bergin 2014"},{"why":"Establishes that cold water vapor abundance is set by photodesorption and photodissociation and is insensitive to UV flux, the basis of Eq. (10).","marker":"Dominik et al. 2005"},{"why":"Justifies focusing on interstellar UV radiation as the photodissociation driver in gas-bearing debris disks.","marker":"Kral et al. 2017"},{"why":"Supplies the compiled CO and dust masses for 21 targets used in the Fig. 2 and Fig. 3 comparisons.","marker":"Rebollido et al. 2022"},{"why":"Supplies observed CO-gas and dust masses for debris disks that anchor the data comparison.","marker":"Moór et al. 2017"},{"why":"Gives the characteristic cometary CO-ice-to-water-ice ratio used in Eq. (8) for the collisional branch.","marker":"Mumma & Charnley 2011"},{"why":"Supplies the ISM CO-to-H2 abundance ratio used to set the primordial CO fraction in Eq. (2).","marker":"Bolatto et al. 2013"},{"why":"Gives the beta Pictoris cold-water non-detection upper limits used as a tentative test of the predicted range.","marker":"Cavallius et al. 2019"}],"fun_headline_variants":["Water vapor mass jump reveals debris disk gas origin","Cold water vapor fingerprints debris disk gas source","Water-to-CO ratio splits debris disk gas origins","Far-IR water lines could trace debris disk history","Debris disk gas origin encoded in water vapor"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation rests on the assumption that cold water vapor abundance is set solely by the photodesorption–photodissociation balance with a single characteristic grain radius and a fixed dust fraction; Appendix A shows this holds only for disks with solid mass below about $0.1\\,\\mathrm{M}_\\oplus$ (CO below about $10^{-3}\\,\\mathrm{M}_\\oplus$), and for more massive, CO-rich disks the water abundance could be much lower, erasing the very jump the diagnostic relies on.","fun_headline_variants_meta":{"raw":{"variants":["Water vapor mass jump reveals debris disk gas origin","Cold water vapor fingerprints debris disk gas source","Water-to-CO ratio splits debris disk gas origins","Far-IR water lines could trace debris disk history","Debris disk gas origin encoded in water vapor"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00075,"raw_usage":{"total_tokens":3369,"prompt_tokens":1002,"completion_tokens":2367,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":618,"completion_tokens_details":{"reasoning_tokens":2297}},"tokens_in":618,"tokens_out":2367,"duration_ms":15208,"temperature":1.0,"reasoning_tokens":2297,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:10:19.404997+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $\\mathrm{H_2O}$ and CO masses in roughly two dozen debris disks spanning CO masses from below $10^{-5}$ to above $10^{-1}\\,\\mathrm{M}_\\oplus$ with a future far-infrared observatory; the diagnostic fails if the $\\mathrm{H_2O}$-mass-versus-CO-mass relation is smoothly monotonic with no jump near $10^{-3}\\,\\mathrm{M}_\\oplus$, or if disks with $\\sim 10^{-1}\\,\\mathrm{M}_\\oplus$ of CO show water masses below roughly $10^{-7}\\,\\mathrm{M}_\\oplus$, well under the primordial-branch prediction.","supporting_citations":[{"cited_title":"R., Bergin, E","cited_arxiv_id":null,"evidence_quote":"Provides the only detection of cold water vapor in a protoplanetary disk (TW Hya) and the mass estimate used to calibrate the primordial-branch dust fraction."},{"cited_title":"2005, ApJ, 635, L85","cited_arxiv_id":null,"evidence_quote":"Establishes that cold water vapor abundance is set by photodesorption and photodissociation and is insensitive to UV flux, the basis of Eq. (10)."},{"cited_title":"2022, MNRAS, 509, 693","cited_arxiv_id":null,"evidence_quote":"Supplies the compiled CO and dust masses for 21 targets used in the Fig. 2 and Fig. 3 comparisons."},{"cited_title":"2019, A&A, 628, A127","cited_arxiv_id":null,"evidence_quote":"Gives the beta Pictoris cold-water non-detection upper limits used as a tentative test of the predicted range."}],"review_version":1}