{"id":"1789c598-e3e5-43f8-a0ac-ed3e3a556356","arxiv_id":"1908.05045","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"First detection of 13C17O in a protoplanetary disk indicates the HD 163296 disk's CO gas mass is 3.5 times larger than previously benchmarked models, supporting 13C17O as an optically thin mass tracer.","lead":"Using ALMA, astronomers detected the rarest stable carbon monoxide molecule, 13C17O, in a planet-forming disk for the first time. Because its signal is not hidden by thick gas, it gives a cleaner measurement of disk gas mass, suggesting the HD 163296 disk holds several times more gas than earlier estimates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The factor-3.5 disk mass increase is degenerate with the assumed 17O/16O ratio, whose 1-sigma error from Qi et al. (2011) alone can reduce the required mass boost to roughly 1.3–1.9.","rationale":"The reader's CONDITIONAL verdict is appropriate: the first detection of 13C17O and the argument that C18O is optically thick within the snowline rest on the same isotope-ratio assumptions, but the LTE column-density comparison in Section 3 (observed 12C18O/13C17O ≈ 2.5 vs. predicted ≈ 250) is a robust, parameter-light argument for C18O optical thickness that does not depend on the detailed disk structure. The weakest link is the conversion from the 13C17O line to disk gas mass, where the assumed abundance of 5.39e-10 relative to H2 is the pivotal input. The paper treats this abundance as fixed by interstellar ratios, but the Qi et al. (2011) derived ratios, especially 12C18O/12C17O = 3.8 ± 1.7, have a 45% 1-sigma uncertainty, and the resulting 17O/16O uncertainty is large enough to absorb most or all of the claimed 3.5x mass increase at only moderate significance. The paper's by-eye fitting and lack of error propagation compound this, but the central physical concern is the unquantified degeneracy between gas mass and 13C17O abundance. My analysis sharpens the reader's weakest_assumption by identifying the specific isotope-ratio error that could resolve the discrepancy without a massive disk. Because the detection and the C18O-optical-thickness conclusion remain plausible, the verdict stays CONDITIONAL rather than moving to REJECT; however, the paper should be required to propagate the isotope-ratio uncertainties and perform a quantitative fit to justify the factor 3.5.","tokens_in":10108,"tokens_out":17365,"duration_ms":160104,"concrete_test":"Re-run the Model 1 and Model 2 radiative transfer calculations with the 13C17O/H2 abundance set to the 1-sigma upper end allowed by Qi et al. (2011), i.e., 16O/17O at the lower bound of 16O/18O = 356 and 18O/17O = 2.1 (and 13C/12C = 1/59), giving 13C17O/H2 ≈ 1.3e-9. Then determine the disk mass required to match the observed 13C17O integrated intensity using a quantitative least-squares fit over the radial profile (not by eye). If the required mass increase at this abundance is less than a factor of ~1.5, the claimed factor of 3.5 is not secure; if it remains above ~2, the qualitative conclusion of a higher disk mass survives but with a revised factor.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—that matching the 13C17O J=3-2 line requires a global CO gas mass increase by a factor of 3.5—is directly proportional to the assumed 13C17O fractional abundance of 5.39e-10 relative to H2 (Section 3, Model 1). That abundance is derived from interstellar isotope ratios combined with the Qi et al. (2011) disk model. However, the Qi et al. (2011) constraints on the relevant ratios carry substantial uncertainty, in particular n(12C18O)/n(12C17O) = 3.8 ± 1.7 (Section 3). Propagating 1-sigma errors through 16O/17O = (16O/18O) × (18O/17O), with 16O/18O = 444 ± 88, gives a range of 17O/16O that is roughly a factor of 1.5–1.8 higher at the upper end. This alone reduces the required mass increase from 3.5 to about 1.9; combining with the allowed 12C/13C error yields a required boost as low as ~1.3. At ~2 sigma in 18O/17O, the inferred mass increase can vanish entirely. The paper does not propagate these isotope-ratio errors into the mass estimate, nor does it explore the degeneracy between gas mass and 13C17O abundance. Its 'lower limit' argument in Section 4.2 only considers depletion or isotope-selective photodissociation, which would lower the abundance and raise the mass, but it does not account for the possibility that the adopted 17O/16O is too low, which would lower the mass. Because the 13C17O line is the only tracer of this abundance and the detection is the sole driver of the 3.5 multiplying factor, the headline mass increase is not robust until the isotope-ratio uncertainty is folded in.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the first detection of the 13C17O J=3-2 line in a protoplanetary disk, HD 163296, using ALMA Cycle 3 data. The detection is claimed at matched-filter S/N ≈ 3.5 and integrated-intensity S/N ≈ 7. The authors compare the new line with existing 12CO, 13CO, C18O, and C17O observations, using a radiative transfer model of the disk (Qi et al. 2011) with interstellar isotopic abundances. This model underpredicts the 13C17O integrated intensity by a factor ~2.5. Scaling the disk gas mass upward by a factor 3.5 (chosen by eye) reproduces the observation, yielding a disk mass of 0.31 Msun (or 0.21 Msun at the revised Gaia distance) and a gas-to-dust ratio of ~260. The paper argues that C18O emission is optically thick inside the CO snowline while 13C17O is optically thin, and uses the new mass to evaluate the disk's gravitational stability.","tokens_in":10512,"tokens_out":9280,"duration_ms":82894,"significance":"If the result holds, it would be the first detection of the rarest stable CO isotopologue in a protoplanetary disk, and it would support the emerging picture that CO-based gas masses may be underestimated because the more abundant isotopologues are optically thick. The application of a single physical disk model to five CO isotopologues, including the new optically thin probe, is a valuable step. The paper also computes optical depth maps and considers an alternative snowline location. The modeling makes use of publicly available codes (LIME, VISIBLE), which aids reproducibility. However, the key quantitative claim—the factor-3.5 mass increase—is not accompanied by an uncertainty estimate and is strongly degenerate with the assumed isotopic abundances, so the significance of the central result is currently limited.","major_comments":[{"comment":"The factor of 3.5 mass increase is determined by an unquantified by-eye fit: the text states that the gas mass was increased from 1.5 in steps of 0.5 until the 'best by-eye fit' of 3.5 was found, with no uncertainty quoted on this factor. Since this factor is the central quantitative result of the paper, the authors should replace the by-eye procedure with a formal chi-square or likelihood analysis over the mass scaling factor (and, ideally, the model parameters) and report a confidence interval or equivalent. As it stands, the headline disk mass of 0.31 Msun has no statistical error bar.","section":"3 (Model 2)"},{"comment":"The assumed 13C17O fractional abundance of 5.39e-10 relative to H2 is derived from the Qi et al. (2011) isotope ratios, which carry substantial uncertainties: n(12C16O)/n(12C18O) = 444±88 and n(12C18O)/n(12C17O) = 3.8±1.7. Propagating these 1-sigma errors into the 17O/16O ratio changes the required gas-mass scaling by roughly a factor of 1.5–1.8; at 2 sigma in 18O/17O the inferred mass increase can be as low as ~1.3 or vanish entirely. The paper does not propagate these errors, nor does it explore the degeneracy between gas mass and 13C17O abundance. This is load-bearing because the factor of 3.5 is the basis for all subsequent conclusions (disk mass, gas-to-dust ratio, Toomre Q analysis).","section":"3 (Model 1) and 4.2"},{"comment":"The argument that the derived gas mass is a lower limit is one-sided. The authors consider only chemical effects (freeze-out, isotope-selective photodissociation, carbon depletion) that would lower the 13C17O abundance and hence raise the inferred mass; they do not consider that the random uncertainty in the adopted interstellar isotope ratios could make the true 13C17O abundance higher than assumed, which would lower the inferred mass. Therefore the statement in Section 5 that the work 'provides robust evidence that disks are more massive than previously assumed' is not supported until the isotope-ratio uncertainty is propagated into the mass estimate.","section":"4.2 and Conclusions"}],"minor_comments":[{"comment":"The phrase 'The resulting S/N ratio is ≈ 3.5' is redundant; use 'signal-to-noise ratio' or 'S/N' instead. Also, the beam size '0.87 x 0.51 (100◦)' should specify 'degrees' for the position angle to avoid ambiguity.","section":"2"},{"comment":"The term 'nH column density' is ambiguous because nH denotes the number density elsewhere in the paper; use 'column density of hydrogen nuclei, NH' instead.","section":"3"},{"comment":"The dust mass from Isella et al. (2007) that yields g/d ≈ 260 is not stated; including it (and its distance assumption) would allow the reader to assess the sensitivity of g/d to the revised Gaia distance.","section":"4.1"},{"comment":"The statement 'We note this does not account for any change in the stellar mass over this time period, which would decrease Q further' is ambiguous; specify whether the stellar mass is assumed constant or specify how it would evolve.","section":"4.4"},{"comment":"The reference entry 'Williams J. P. and Best W. M. J. 2014' should be formatted as 'Williams, J. P., & Best, W. M. J. 2014' to match the journal style.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The detection appears genuine and the paper addresses an important question, but the central mass-increase claim is under-supported because of the absence of error propagation and the by-eye fit. The authors should be encouraged to add a quantitative parameter search and to reframe the conclusions as conditional on the assumed isotopologue abundances. I would be willing to re-review a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First things first: this is a real detection and the right kind of result to push. 13C17O J=3-2 is detected in HD 163296 via matched filter at S/N ~3.5 and integrated at ~7, and the paper is careful about beam, channels, and comparison to the other CO isotopologues. The modelling uses the Qi et al. (2011) structure, which is reasonable, and it demonstrates that C18O is optically thick in the inner disk while 13C17O stays thin. That is a useful step for disk gas mass work.\n\nWhere I part ways with the abstract is the 'factor of 3.5' and the 'robust' language. The factor is a by-eye scaling of the whole disk gas mass to match one radially averaged profile, no uncertainty quoted. More worrying, the 13C17O fractional abundance is fixed at 5.39e-10 from interstellar ratios, and the isotope ratios carry sizable errors. Propagating the quoted 1-sigma errors from Qi et al. (2011) alone shrinks the required boost to about 1.3-1.9. At 2 sigma on 18O/17O, the boost basically disappears. So the headline mass increase is not robust until the abundance degeneracy is folded in. The paper acknowledges depletion and isotope-selective photodissociation as reasons the mass could be higher, but it never considers that the adopted 17O/16O could be too low, which would make the mass lower. That one-sided treatment matters.\n\nOther soft spots are minor: the fit is not quantified, the disk structure is smooth rather than the observed ringed structure, and the 'requires' language in the abstract overstates what is actually a scaled match. The Toomre Q section is clearly speculative and labelled as such; it is fine as discussion.\n\nStill, the detection stands. The observation is new, the molecular data handling is careful, and the conclusion that C18O is optically thick in this disk is supported by the comparison. This deserves a serious referee, but the referee should ask for a real uncertainty analysis on the mass, including isotope-ratio errors and some exploration of abundance-mass degeneracy, before the factor of 3.5 is taken at face value.","headline":"A genuine first detection and a useful demonstration that 13C17O is optically thin, but the factor-3.5 disk mass increase is not robust until isotope-ratio uncertainties are propagated.","tokens_in":11127,"tokens_out":3284,"would_cite":true,"duration_ms":30217,"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":"The first detection of the rarest form of carbon monoxide in a planet-forming disk shows the disk holds 3.5 times more gas than previous estimates.","keywords":["protoplanetary disks","CO isotopologues","disk gas mass","13C17O","HD 163296","radiative transfer","optically thin emission","astrochemistry"],"falsifier":"Measure a second $^{13}$C$^{17}$O rotational transition (for example $J=2$--$1$ or $J=4$--$3$) in HD 163296: optically thin LTE emission predicts a fixed intensity ratio set by the excitation temperature, so a measured ratio that departs from that prediction, or a line profile showing self-absorption, would mean the emission is not fully optically thin and the inferred factor-of-3.5 mass increase would need revision.","tokens_in":9934,"feed_emoji":"🪐","tokens_out":14859,"duration_ms":122161,"temperature":0.7,"pith_summary":"This paper reports the first detection of $^{13}$C$^{17}$O, the rarest stable form of carbon monoxide, in a planet-forming disk, and uses it to weigh the gas in the disk around the star HD 163296. Because $^{13}$C$^{17}$O is so scarce, its emission is optically thin, so its brightness scales with the total amount of CO gas. By contrast, the more abundant C$^{18}$O is shown to be optically thick inside the CO snow line, meaning it hides part of the disk's gas. To match the observed $^{13}$C$^{17}$O brightness with a radiative-transfer model, the CO gas mass must be raised by a factor of 3.5, giving a disk mass near $0.31\\,M_\\odot$ ($0.21\\,M_\\odot$ at the revised Gaia distance) and a gas-to-dust ratio near 260. If correct, this implies many disk gas masses derived from C$^{18}$O have been underestimated, with direct consequences for how much material is available for planet formation.","feed_headline":"Detecting the rarest CO reveals 3.5× more disk gas","feed_subtitle":"A newly detected rare molecule sees through the opaque gas that hid the disk's true mass from older tracers.","key_machinery":"The central object is the $^{13}$C$^{17}$O $J=3$--$2$ line, the rarest stable CO isotopologue, whose low abundance guarantees optically thin emission that scales directly with CO column density. The argument is carried by a forward radiative-transfer model: the disk density and temperature structure from Qi et al. (2011), a constant $^{13}$C$^{17}$O fractional abundance of $5.39\\times 10^{-10}$ relative to H$_2$, and the radiative transfer code described by Brinch & Hogerheijde (2010) to synthesize the hyperfine-blended line, which is then smoothed to the observed beam and compared with deprojected, azimuthally averaged radial profiles of five isotopologues. The key diagnostic is the optical depth: computed maps show $\\tau(\\mathrm{C}^{18}\\mathrm{O}\\,J=2$--$1)>1$ inside the 155 au snowline while $\\tau({}^{13}\\mathrm{C}^{17}\\mathrm{O}\\,J=3$--$2)<1$ everywhere, so the $^{13}$C$^{17}$O brightness requires the global CO gas mass to be increased by a factor of 3.5.","core_discovery":"Using ALMA Band 7 observations, the authors detect the $^{13}$C$^{17}$O $J=3$--$2$ line from the HD 163296 disk and compare its radial intensity profile with archival maps of $^{12}$CO, $^{13}$CO, C$^{18}$O, and C$^{17}$O. Adopting the disk density and temperature structure of Qi et al. (2011) and interstellar isotope ratios, which set the $^{13}$C$^{17}$O abundance to $5.39\\times 10^{-10}$ relative to H$_2$, their radiative-transfer model underpredicts the observed $^{13}$C$^{17}$O integrated intensity by a factor of 2.5. Raising the total gas mass by a factor of 3.5 reproduces all the isotopologue profiles and yields a disk gas mass of $0.31\\,M_\\odot$ ($0.21\\,M_\\odot$ at the 101.5 pc Gaia distance). The paper argues that C$^{18}$O is optically thick within the CO snow line at 155 au, so it cannot trace the bulk disk gas, whereas $^{13}$C$^{17}$O stays optically thin across the disk and therefore does trace that gas.","pith_inferences":["If isotope-selective photodissociation lowers the $^{13}$C$^{17}$O abundance in the line-forming region below the interstellar value, the required gas mass would be even higher than the factor of 3.5 quoted here, making disk masses from CO an even stronger lower limit.","The same measurement in the TW Hya disk, where $^{13}$C$^{18}$O has already been detected, would show whether the optical-depth bias in C$^{18}$O is universal or specific to massive, warm disks around Herbig Ae/Be stars such as HD 163296.","A clean test of the paper's optically thin assumption is to observe two $^{13}$C$^{17}$O transitions: their line ratio should follow the LTE optically thin prediction, and any excess ratio would reveal either residual opacity or a non-thermal excitation component not in the adopted model."],"forward_implications":["C$^{18}$O-based disk gas masses are lower limits wherever C$^{18}$O is optically thick; for HD 163296 the underestimate relative to $^{13}$C$^{17}$O is a factor of 2–6.","The HD 163296 disk gas mass is $0.31\\,M_\\odot$ ($0.21\\,M_\\odot$ at 101.5 pc), with a gas-to-dust ratio near 260, placing it at the high end of previous estimates.","With higher angular resolution and sensitivity, $^{13}$C$^{17}$O emission can directly map the midplane CO snow line, a measurement C$^{18}$O cannot provide because of its opacity.","The disk is gravitationally stable today (minimum Toomre Q near 6 at about 110 au), but extrapolated back to an age of 0.1 Myr the same mass gives Q between 1.3 and 0.7, implying early gravitational instability is plausible.","Detecting $^{13}$C$^{17}$O in additional disks could help close the gap between measured disk masses and the masses inferred for the exoplanet population."],"supporting_citations":[{"why":"Provides the benchmarked disk density and temperature structure and the isotope ratios used to set the model and the $^{13}$C$^{17}$O abundance.","marker":"Qi et al. 2011"},{"why":"Supplies the archival $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2-1 maps and disk constraints against which the models are compared.","marker":"Isella et al. 2016"},{"why":"Supplies the $^{13}$C$^{17}$O molecular line data (hyperfine components) used for the column density estimate and the line modelling.","marker":"Klapper et al. 2003"},{"why":"Provides the molecular data catalogue from which the $^{13}$C$^{17}$O transition parameters are taken.","marker":"Müller et al. 2005"},{"why":"Provides the radiative transfer code used to synthesize the hyperfine-blended $^{13}$C$^{17}$O line profiles.","marker":"Brinch & Hogerheijde 2010"},{"why":"Documents the calibration and self-calibration of the ALMA Band 7 observations from which the $^{13}$C$^{17}$O detection is made.","marker":"Notsu et al. 2019"},{"why":"Gives the interstellar carbon and oxygen isotope ratios used to set the fractional abundance of $^{13}$C$^{17}$O.","marker":"Wilson 1999"},{"why":"Provides the 122 pc source distance used for the fiducial mass estimates and literature comparison.","marker":"van den Ancker et al. 1998"},{"why":"Provides the revised 101.5 pc distance, which scales the disk mass down to 0.21 solar masses.","marker":"Gaia Collaboration et al. 2018"},{"why":"Provides the dust mass used to derive the gas-to-dust ratio of about 260.","marker":"Isella et al. 2007"}],"fun_headline_variants":["First 13C17O detection: disk gas 3.5× heavier","Rarest CO tracer reveals disk mass 3.5× larger","13C17O sees through opaque gas to weigh disk","Rare CO isotope uncovers 3.5× more disk gas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 3.5-fold mass increase rests on assuming that the amount of $^{13}$C$^{17}$O per hydrogen molecule in the disk equals the interstellar value of $5.39\\times 10^{-10}$ and that the adopted density and temperature structure of the disk is correct; if either assumption is wrong, the required scaling changes.","fun_headline_variants_meta":{"raw":{"variants":["First 13C17O detection: disk gas 3.5× heavier","Rarest CO tracer reveals disk mass 3.5× larger","13C17O sees through opaque gas to weigh disk","Rare CO isotope uncovers 3.5× more disk gas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000357,"raw_usage":{"total_tokens":1983,"prompt_tokens":1038,"completion_tokens":945,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":654,"completion_tokens_details":{"reasoning_tokens":868}},"tokens_in":654,"tokens_out":945,"duration_ms":8323,"temperature":1.0,"reasoning_tokens":868,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:24:34.195640+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure a second $^{13}$C$^{17}$O rotational transition (for example $J=2$--$1$ or $J=4$--$3$) in HD 163296: optically thin LTE emission predicts a fixed intensity ratio set by the excitation temperature, so a measured ratio that departs from that prediction, or a line profile showing self-absorption, would mean the emission is not fully optically thin and the inferred factor-of-3.5 mass increase would need revision.","supporting_citations":[{"cited_title":"2016, Physical Review Letters, 117, 251101","cited_arxiv_id":null,"evidence_quote":"Supplies the archival $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2-1 maps and disk constraints against which the models are compared."},{"cited_title":"2003, The Astrophysical Journal, 582, 262","cited_arxiv_id":null,"evidence_quote":"Supplies the $^{13}$C$^{17}$O molecular line data (hyperfine components) used for the column density estimate and the line modelling."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the radiative transfer code used to synthesize the hyperfine-blended $^{13}$C$^{17}$O line profiles."},{"cited_title":"2019, ApJ, 875, 96","cited_arxiv_id":null,"evidence_quote":"Documents the calibration and self-calibration of the ALMA Band 7 observations from which the $^{13}$C$^{17}$O detection is made."}],"review_version":1}