{"id":"d01470c7-44a4-448d-b675-b8b646575ac7","arxiv_id":"2502.05061","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":16,"one_line_summary":"The first detections of C18O16O and H13CN emission in an inner protoplanetary disk are reported for the high-inclination disk around MY Lup, along with tentative C17O16O and HC15N.","lead":"JWST's MIRI instrument observed the edge-on disk around the young star MY Lup and found strong emission from rare isotopic versions of CO2 and HCN, including first-time detections of C18O16O and H13CN in an inner disk. This opens a new way to measure isotope ratios in the planet-forming region, though the derived ratios are still uncertain.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Detections are likely real, but the derived isotope ratios and column densities rest on a single-slab LTE assumption that the paper itself shows is inconsistent at the CO2 red edge; the quantitative fractionation interpretation needs re-fitting.","rationale":"The paper's primary contribution is the detection of multiple isotopologues in an inner disk, and the spectral evidence for the two secure detections (C18O16O and H13CN) is strong: the features are clear, at expected wavelengths, and not claimed to be more than tentative for the other two species. The reader's conditional verdict is appropriate because the quantitative interpretation—high column densities and isotope ratios—depends on the single-slab LTE same-reservoir assumption. I agree with the reader's weakest_assumption, and I sharpen it by pointing to the paper's own Appendix B, which provides direct evidence that a single LTE slab cannot fit all CO2 features, and to the absence of any test where isotopologue temperatures are allowed to vary independently. The paper is commendably transparent about these limitations (e.g., Section 5.1 and the self-absorption caveat in Section 5.2), but transparency does not remove the concern: the reported 16O/18O and 12C/13C values are used to discuss fractionation, and they are insufficiently robust to support a definitive claim. The concrete check of relaxing the temperature independence would settle whether the derived ratios change significantly; if they do not, the interpretation is strengthened. The central detection claim, however, remains secure, so the verdict should remain CONDITIONAL rather than shift to REJECT or ACCEPT.","tokens_in":22410,"tokens_out":7828,"duration_ms":75850,"concrete_test":"Re-run the CO2 grid fitting with the 13CO2 and C18O16O temperatures and emitting areas as free parameters independent of 12CO2, then map the derived 16O/18O ratio as a function of T_13CO2 (e.g., vary T_13CO2 by ±100 K around the best-fit 325 K). If the resulting 16O/18O ratio shifts by more than the quoted 1σ range (281–513), the claimed 18O enhancement is not robust and the same-reservoir assumption is load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of first detections is spectroscopically credible: the C18O16O and H13CN Q-branch features are distinct and appear at expected wavelengths. The load-bearing weakness is in the quantitative analysis that follows. The paper derives isotope ratios and column densities by assuming all isotopologues of a given molecule emit from a single LTE slab with the same temperature and projected area (Sections 4.2 and 4.3). This assumption is explicitly acknowledged as 'certainly not strictly true' in Section 5.1, but its impact is not quantified. The paper's own evidence shows the assumption is already strained: Appendix B demonstrates that a single LTE slab cannot simultaneously fit the 13.88 µm and 16.2 µm CO2 features (which prefer N ~ 1e18 cm-2) and the red edge of the 12CO2 Q branch (which prefers low N), suggesting non-LTE excitation or a multi-component emitting region. Furthermore, the two analysis paths (fixed 12C/13C and fixed 16O/18O) each require the other ratio to be fixed, and the resulting isotopic ratios (e.g., 16O/18O = 381) are consistent with ISM values at only ~2σ, with error bars that do not include the systematic uncertainty from the same-reservoir assumption. If the isotopologues trace different disk layers with different temperatures, the Q-branch peak ratios change, biasing the derived column densities and isotope ratios. The paper's claim that 'no clear evidence' requires different temperatures is not a test of that hypothesis. Because the interpretation of MY Lup's unique spectrum as evidence for isotopic fractionation depends on these model-dependent ratios, the quantitative conclusions should be treated as conditional until the single-slab assumption is relaxed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents JWST MIRI-MRS spectroscopy of the high-inclination protoplanetary disk around MY Lup. The spectrum shows weak H2O emission but strong CO2 and HCN emission, and the authors report detections of the rare isotopologues C18O16O and H13CN, with tentative detections of C17O16O and HC15N. LTE slab modeling is used to derive temperatures, column densities, and emitting areas for H2O, CO2, and HCN, and isotopologue ratios are estimated by fixing either the carbon or oxygen isotope ratio to ISM values. The paper attributes the unusual spectral appearance to a combination of inner disk clearing and the near-edge-on viewing geometry, and discusses the implications for isotopic fractionation studies in inner disks.","tokens_in":22843,"tokens_out":4373,"duration_ms":46765,"significance":"If the line identifications are correct, the C18O16O and H13CN detections are the first in an inner protoplanetary disk and demonstrate that JWST MIRI-MRS can access trace isotopologues in the terrestrial-planet-forming region. This is a valuable observational result with clear follow-up potential. The paper is transparent about its modeling degeneracies, makes the reduced data publicly available, and places the results in the context of other MIRI-MRS disk spectra. The main quantitative conclusions on isotope ratios and column densities are, however, conditional on a single-slab LTE assumption and on fixing one isotope ratio to ISM values; the quoted uncertainties do not include the resulting systematic errors, so the fractionation interpretation is not yet supported at the claimed confidence.","major_comments":[{"comment":"The central quantitative analysis assumes that all CO2 isotopologues arise from a single LTE slab with one temperature and column density. The paper's own residual analysis shows this assumption is already strained: the 13.88 micron and 16.2 micron features prefer high column densities (~1e18 cm^-2) while the red edge of the main CO2 Q branch near 14.98 micron prefers low column densities. Since the isotope-ratio estimates are derived from Q-branch peak ratios under the same single-slab assumption, the inconsistency introduces an unquantified systematic bias into the derived N and isotope ratios. The authors should either fit a multi-component or non-LTE model, or explicitly present the isotope ratios as conditional on the single-slab assumption with an estimate of the resulting systematic uncertainty.","section":"Section 4.2 and Appendix B"},{"comment":"The derivation of the isotope ratios is circular in a practical sense: the 16O/18O value of 381(+132,-100) is obtained after fixing 12C/13C to the ISM value of 68, while the 12C/13C value of 77(+84,-25) is obtained after fixing 16O/18O to the ISM value of 557. The paper explicitly acknowledges the N-versus-ratio degeneracy, but the quoted error bars nevertheless include only the statistical scatter under each fixed-ratio assumption. The systematic uncertainty from the fixed ratio and from the same-reservoir assumption is not propagated. Because the derived ratios are consistent with ISM values at only about 2 sigma, the statements in Section 5.1 about possible 18O enhancement and 17O depletion should be reframed as conditional constraints rather than measurements.","section":"Section 4.2, Table 1"},{"comment":"The HCN analysis has a similar load-bearing degeneracy. The HCN-only fit produces a family of models within the 1-sigma contour with 12C/13C values spanning roughly 21 to 373, and the combined HCN plus H13CN fit only tightly constrains the parameters after fixing 12C/13C to the ISM value. Consequently, the quoted NHCN = 1.4e19 cm^-2 and NH13CN = 2.0e17 cm^-2 are conditional on that fixed ratio and on the same-temperature, same-area assumption. The claim that the high HCN column density is supported by the isotopologue detection is therefore only as strong as the fixed-ratio assumption; this should be stated more prominently and the systematic dependence of NHCN on the assumed ratio should be quantified.","section":"Section 4.3 and Figure 12"},{"comment":"The geometric enhancement factor invoked to explain the high line-of-sight column densities assumes a plane-parallel slab with no self-absorption. The paper itself notes that the high inclination could produce self-absorption in the molecular lines, which would affect the measured fluxes and hence the derived column densities and ratios. Since the inclination argument is one of the two main explanations for MY Lup's unique spectrum, the authors should estimate the magnitude of possible self-absorption effects or clearly state that the column densities and ratios are upper/lower limits under the adopted geometry.","section":"Section 5.2 and Figure 15"}],"minor_comments":[{"comment":"The phrase 'observations at higher spectral resolving power is needed' should be 'observations at higher spectral resolving power are needed'.","section":"Abstract"},{"comment":"The caption refers to '12CO and 13CO' and to '12CO/C18O16O'; these should be '12CO2 and 13CO2' and '12CO2/C18O16O' to match the plotted species.","section":"Figure 7 caption"},{"comment":"The notation for isotope ratios is confusing: column (a) defines the ratio as main isotope divided by heavy isotope, while note (b) defines Ri as the heavy-to-main ratio. Please reconcile the definitions so the reader does not have to invert ratios when reading the table.","section":"Table 1"},{"comment":"The sentence describing C17O16O says 'a12CO2 P-branch line' with a missing space after the article; this should read 'a 12CO2 P-branch line'.","section":"Section 4.2"},{"comment":"The statement that 'It is not clear that there is a mechanism to explain enhancement in one heavy isotope with simultaneous depletion in another' would benefit from a reference to relevant photodissociation or chemical fractionation models, if one exists.","section":"Section 5.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is in scope for an astrophysics journal and the central observational result, the detection of rare CO2 and HCN isotopologues, is likely to be of significant interest. The main concern is that the quantitative isotope ratios and column densities are presented as measurements while resting on a single-slab LTE assumption that the paper itself shows to be imperfect, and on fixing one isotope ratio to ISM values. I believe the detections are credible enough that a revision, rather than rejection, is appropriate; the authors should either provide systematic-error estimates or soften the quantitative claims accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid observational paper with genuinely new detections—C18O16O and H13CN in an inner disk—and they look real. The modeling is careful and the paper is unusually candid about its own degeneracies. The main weakness is that the isotope ratios, which are the scientific payoff, rest on a single-slab LTE assumption that the paper's own residuals show isn't quite right. That doesn't kill the paper; it means the quantitative ratios are conditional, which the authors basically say.\n\nWhat's new: first MIRI-MRS detections of C18O16O and H13CN in an inner disk, plus tentative C17O16O and HC15N. The Q-branch features are distinct and at expected wavelengths. The comparison with other CO2-rich disks (Sz 98, GW Lup, DF Tau) is useful, and the raw data and reduction code are publicly available. The paper also reports weak water, cold temperatures, high CO2 and HCN columns, and connects this to inclination and inner clearing. That interpretation is reasonable but not unique.\n\nWhere it's soft: the single-slab LTE model. The fits to 12CO2 alone give N ~1e17, but including 13CO2 or C18O16O pushes N to ~1e18–1e19. The paper shows residuals where one feature wants low N and another wants high N (Figure 8 and Appendix B). The claim that there's 'no clear evidence' for different temperatures between isotopologues is not a strong test. The isotope ratios are derived by fixing one ratio to the ISM value and fitting the other, then swapping—so the two numbers are not independent. The error bars don't include the systematic from the same-reservoir assumption. The authors acknowledge all of this, which is to their credit, but it means the fractionation numbers (18O enhancement ~1.5, 17O depletion) should be read as illustrative, not measured. The detections stand; the ratios don't yet.\n\nCitation pattern is fine. The paper cites the relevant JDISCS/MINDS literature and the GW Lup 13CO2 work. No obvious gaps.\n\nWho's it for: disk chemistry and JWST disk observers. It's a useful data point in the growing MIRI-MRS disk sample, and it identifies a target for high-resolution follow-up.\n\nRecommendation: send to review. The referee should push on the single-slab assumption and ask for a multi-component or non-LTE exploration before the ratios can be taken at face value. But the detections themselves deserve publication, and the paper is honest enough that a light revision may be enough.","headline":"Genuinely new inner-disk isotopologue detections that hold up; the derived isotope ratios are honestly labeled as model-dependent and need a closer look, but the paper deserves refereeing.","tokens_in":23442,"tokens_out":3531,"would_cite":true,"duration_ms":31447,"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":"A JWST MIRI-MRS spectrum of the edge-on disk MY Lup yields the first inner-disk detections of the rare isotopologues C18O16O and H13CN, opening isotopic fractionation in planet-forming regions to observation.","keywords":["protoplanetary disks","isotopologues","JWST MIRI-MRS","CO2 emission","HCN emission","isotopic fractionation","inner disk chemistry","MY Lup"],"falsifier":"Take a deeper, higher-resolution MIRI spectrum of MY Lup that resolves the C$^{17}$O$^{16}$O and HC$^{15}$N Q branches: if the C$^{17}$O$^{16}$O feature disappears at higher signal-to-noise, the tentative $^{17}$O depletion is not real; if it persists with a peak ratio implying $^{16}$O/$^{17}$O below roughly 2000, the depletion is confirmed and simple mass-dependent fractionation alone cannot explain the isotope pattern.","tokens_in":22215,"feed_emoji":"🔭","tokens_out":10292,"duration_ms":74386,"temperature":0.7,"pith_summary":"This paper reports a JWST MIRI-MRS spectrum of the nearly edge-on protoplanetary disk around the young star MY Lup and claims the first detections in an inner disk of the rare isotopologues C$^{18}$O$^{16}$O and H$^{13}$CN, along with tentative detections of C$^{17}$O$^{16}$O and HC$^{15}$N. A sympathetic reader would care because these molecules carry the isotope ratios used to test how volatiles are fractionated in the terrestrial-planet-forming region. The authors use LTE slab models to derive molecular temperatures, column densities, and emitting areas, finding cold, high-column CO$_2$ and HCN and unusually weak water. They interpret the spectrum as a combination of inner disk clearing and high inclination, and argue that MIRI-MRS can open inner-disk isotopologue astronomy even though the current ratios are too uncertain to distinguish fractionation scenarios.","feed_headline":"First inner-disk detections of C18O16O and H13CN","feed_subtitle":"JWST MIRI spectra of edge-on MY Lup open isotope ratios in the planet-forming zone to direct measurement.","key_machinery":"The load-bearing tool is a single-temperature, single-column LTE slab emission model (a 'slab model') applied to the blended Q-branch spectra: for each molecular species, a grid of column densities and temperatures is generated, the projected emitting area is scaled to match peak fluxes, and chi-squared contours locate the best fit. The key trick for isotopologues is using the Q-branch peak ratios of the main and rare species, with one isotope ratio fixed to an ISM value, to break the degeneracy between column density and abundance ratio. A second geometric ingredient is the high inclination: in a plane-parallel atmosphere, the path through the gas grows roughly as the secant of the inclination angle (a factor of 3–4 at MY Lup's 77° inclination), amplifying the line-of-sight column without requiring unusual abundances.","core_discovery":"The central observational claim is that MY Lup's inner disk emits from multiple isotopologues of CO$_2$ and HCN, including the first clear detections of C$^{18}$O$^{16}$O and H$^{13}$CN in any inner disk and tentative detections of C$^{17}$O$^{16}$O and HC$^{15}$N. When fit with single-temperature, single-column LTE slabs, the emission requires CO$_2$ column densities of $3.5\\times10^{18}$ to $5.6\\times10^{18}$ cm$^{-2}$ at $T \\approx 300$–$325$ K and a high HCN column of $1.4\\times10^{19}$ cm$^{-2}$ at $T = 250$ K, with small emitting radii near 0.4–0.65 AU. The isotopologue Q-branch peak ratios, combined with assumed ISM carbon or oxygen ratios, yield $^{16}$O/$^{18}$O $\\approx 381^{+132}_{-100}$, $^{12}$C/$^{13}$C $\\approx 77^{+84}_{-25}$, and a marginal $^{16}$O/$^{17}$O that suggests mild depletion of $^{17}$O relative to ISM in both modeling approaches; all are consistent with ISM values within about 2σ. The paper proposes that MY Lup's unique spectrum is caused by a cleared inner disk that suppresses warm water, combined with the edge-on viewing geometry that lengthens the observed gas column.","pith_inferences":["A natural extension is to relax the single-slab assumption and fit the main and rare isotopologues with separate temperature layers; if the apparent $^{18}$O enhancement and $^{17}$O depletion shift together or vanish, the current isotope ratios are artifacts of vertical stratification.","A survey of high-inclination and transition disks with MIRI-MRS could map $^{12}$C/$^{13}$C and $^{16}$O/$^{18}$O/$^{17}$O in inner disks and compare directly with ALMA measurements of outer-disk CO isotopologues, exposing radial fractionation gradients.","The apparent link between high column density and isotopologue detectability implies a selection effect: any statistical sample of inner-disk isotope ratios will be biased toward the most column-rich, often edge-on, disks, which should be accounted for in population interpretations."],"forward_implications":["Trace isotopologues of CO$_2$ and HCN are detectable in inner disks with JWST MIRI-MRS, making inner-disk isotope ratios an observable quantity rather than an extrapolation from the solar system.","Nearly edge-on disks and disks with inner clearings are the most promising targets for isotopologue searches, since geometry and reduced dust opacity raise the observed gas column.","The derived ratios, while uncertain, are consistent with ISM values within about 2σ, so this dataset does not yet demand exotic fractionation; it sets the stage for higher-resolution tests.","If follow-up spectroscopy confirms the high CO$_2$ and HCN columns at cold temperatures, models invoking inner clearing and pebble drift of volatile-rich material become testable for MY Lup specifically."],"supporting_citations":[{"why":"reported the first inner-disk detection of 13CO2 in GW Lup, establishing the sensitivity baseline this work extends to rarer isotopologues.","marker":"Grant et al. (2023)"},{"why":"describes the MIRI-MRS data reduction and asteroid-based fringe removal used to produce the calibrated spectrum.","marker":"Pontoppidan et al. (2024)"},{"why":"supplies the ISM 12C/13C = 68 ratio used to fix the carbon isotope scale in the CO2 and HCN fits.","marker":"Milam et al. (2005)"},{"why":"supplies the ISM 16O/18O = 557 and 16O/17O = 2005 ratios used as the oxygen reference.","marker":"Wilson (1999)"},{"why":"HITRAN line lists used to identify isotopologue transitions and upper-level energies for CO2, HCN and water.","marker":"Gordon et al. (2022)"},{"why":"near-infrared disk inclination of 77 degrees used to convert emitting areas to radii and to estimate the secant-of-inclination column enhancement.","marker":"Avenhaus et al. (2018)"},{"why":"models inner clearing beyond the water snowline to explain strong CO2 with weak water, the scenario applied to MY Lup.","marker":"Vlasblom et al. (2024)"},{"why":"provides the spectools-ir slab fitting code used for the LTE emission models and line flux extraction.","marker":"Salyk (2022)"}],"fun_headline_variants":["First inner-disk C18O16O and H13CN","JWST reveals rare CO2 and HCN isotopologues in MY Lup","Edge-on disk exposes first C18O16O and H13CN","C18O16O and H13CN: first inner-disk sightings","MY Lup's edge-on disk: new CO2 and HCN isotopes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything derived about isotope ratios assumes the main and rare isotopologues of each molecule share a single well-mixed slab with one temperature and column density, an assumption the paper acknowledges is 'certainly not strictly true' for a disk with vertical and radial structure.","fun_headline_variants_meta":{"raw":{"variants":["First inner-disk C18O16O and H13CN","JWST reveals rare CO2 and HCN isotopologues in MY Lup","Edge-on disk exposes first C18O16O and H13CN","C18O16O and H13CN: first inner-disk sightings","MY Lup's edge-on disk: new CO2 and HCN isotopes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000501,"raw_usage":{"total_tokens":2547,"prompt_tokens":1138,"completion_tokens":1409,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":754,"completion_tokens_details":{"reasoning_tokens":1312}},"tokens_in":754,"tokens_out":1409,"duration_ms":12179,"temperature":1.0,"reasoning_tokens":1312,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T20:22:04.393542+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a deeper, higher-resolution MIRI spectrum of MY Lup that resolves the C$^{17}$O$^{16}$O and HC$^{15}$N Q branches: if the C$^{17}$O$^{16}$O feature disappears at higher signal-to-noise, the tentative $^{17}$O depletion is not real; if it persists with a peak ratio implying $^{16}$O/$^{17}$O below roughly 2000, the depletion is confirmed and simple mass-dependent fractionation alone cannot explain the isotope pattern.","supporting_citations":[{"cited_title":"E., Rothman, L","cited_arxiv_id":null,"evidence_quote":"HITRAN line lists used to identify isotopologue transitions and upper-level energies for CO2, HCN and water."}],"review_version":1}