pith. machine review for the scientific record. sign in

arxiv: 2604.12242 · v1 · submitted 2026-04-14 · 🌌 astro-ph.EP · astro-ph.GA· astro-ph.SR

Recognition: unknown

JWST/MIRI Hydrocarbon and Water Absorption in the Wind of a Young Disk: Signatures of Pebble Drift and Carbon Grain Sublimation

Mar\'ia Jos\'e Colmenares , Edwin A. Bergin , Ke Zhang , Geoffrey A. Blake , Klaus M. Pontoppidan , Alexa R. Anderson , John Carr , Emma Dahl
show 22 more authors
Joan Najita Jonathan P. Williams Colette Salyk Till Kaeufer Mayank Narang Ilaria Pascucci Beno\^it Tabone Lucas Cieza Miguel Vioque Adrien Houge Sebastiaan Krijt Aditya M. Arabhavi Giovanni Rosotti John Carpenter Feng Long Paola Pinilla Jayatee Kanwar Eshan Raul Karina Mauco James Miley Abygail Waggoner the JDISCS collaboration
Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:28 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.GAastro-ph.SR
keywords protoplanetary disksmolecular absorptiondisk windspebble drifthydrocarbon chemistrywater vaporcarbon grainsJWST MIRI
0
0 comments X

The pith

Blueshifted absorption in a young disk wind reveals enhanced hydrocarbons and water from inward pebble drift and carbon grain processing at the soot line.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper uses JWST mid-infrared spectra of an inclined young disk to detect strong blueshifted absorption from HCN, C2H2, CH4, CO2, H2O and CO, with no ice features. These lines arise from warm gas in the inner few AU and show higher hydrocarbon-to-CO and water-to-CO ratios than seen in other similar systems. The authors interpret the ratios as the direct chemical signature of icy pebbles drifting inward and sublimating, combined with high-temperature processing of carbon grains. A reader would care because this indicates that distinctive carbon-rich compositions can appear in the inner disk within the first million years and can be observed in real time through outflowing material. The finding links observable wind chemistry to the processes that set the inventory of solids and volatiles available for planet formation.

Core claim

The absorption lines are significantly blueshifted relative to systemic velocity, with mid-IR lines showing larger shifts than near-IR CO, indicating a velocity- and temperature-stratified molecular disk wind. In this framework the lines directly sample disk material lifted from the inner disk surface. ISO-Oph 37 exhibits unusually high (C2H2 + CH4)/HCN, (C2H2 + CH4)/CO and H2O/CO column-density ratios while CO and HCN columns remain typical. These ratios are best explained by simultaneous enhancement of hydrocarbons and water driven by inward drift and sublimation of icy pebbles together with thermal processing of carbonaceous grains at the soot line. The source therefore shows that carbon-

What carries the argument

The velocity- and temperature-stratified molecular disk wind, whose blueshifted absorption lines preserve the chemical imprint of the wind-launching region without significant contributions from unrelated components.

If this is right

  • Carbon-rich inner-disk chemistry can be established within the first million years of disk evolution.
  • Molecular absorption in disk winds provides a direct probe of the chemical composition at the wind-launching region.
  • Inward pebble drift and sublimation plus carbon-grain processing at the soot line simultaneously boost both hydrocarbon and water abundances.
  • Similar hydrocarbon-rich signatures may appear in other young, inclined disks with outflows.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The observed wind chemistry could seed the carbon inventory of terrestrial planets forming in the inner disk.
  • Multi-epoch observations of additional inclined sources would test whether this enhancement phase is common or transient.
  • Linking these inner-disk wind compositions to atmospheric measurements of close-in exoplanets could constrain how early disk processing affects final planetary building blocks.

Load-bearing premise

The blueshifted absorption lines directly sample the unaltered chemical composition of the wind-launching region in the inner few au without significant contributions from other velocity or temperature components along the line of sight.

What would settle it

High-resolution spectra showing that the absorption profiles contain multiple velocity components unrelated to a single stratified wind, or molecular ratios in other young wind sources that do not match the pebble-drift and soot-line predictions.

Figures

Figures reproduced from arXiv: 2604.12242 by Abygail Waggoner, Aditya M. Arabhavi, Adrien Houge, Alexa R. Anderson, Beno\^it Tabone, Colette Salyk, Edwin A. Bergin, Emma Dahl, Eshan Raul, Feng Long, Geoffrey A. Blake, Giovanni Rosotti, Ilaria Pascucci, James Miley, Jayatee Kanwar, Joan Najita, John Carpenter, John Carr, Jonathan P. Williams, Karina Mauco, Ke Zhang, Klaus M. Pontoppidan, Lucas Cieza, Mar\'ia Jos\'e Colmenares, Mayank Narang, Miguel Vioque, Paola Pinilla, Sebastiaan Krijt, the JDISCS collaboration, Till Kaeufer.

Figure 1
Figure 1. Figure 1: Spectral energy distribution of ISO-Oph 37. Photometric points are compiled from K. Zhang et al. (2025); G. Marton et al. (2024). The X-SHOOTER spectrum is taken from C. F. Manara et al. (2015). The dotted gray lines show the best-fit SEDs, described in Appendix D. The inset shows the JWST/MIRI-MRS spectrum, and Spitzer/IRS spectrum taken from the CASSIS database (V. Lebouteiller et al. 2011), with some id… view at source ↗
Figure 2
Figure 2. Figure 2: Observed CO lines in ISO-Oph 37. The pink line shows the Keck/NIRSPEC spectrum and the purple line shows the JWST/MIRI-MRS data. The MRS spectrum is shifted down for clarity. Vertical lines and makers show the locations of the 12CO v = 1 − 0 and v = 2 − 1, and 13CO v = 1 − 0 lines. Labels denote individual 12CO v = 1 − 0 rotational quantum numbers [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of the effect of the covering frac￾tion fc on the slab model spectra for CH4 and C2H2. As the column density increases, the flux of the P branch also increases, while the Q branch flux decreases. All models as￾sume Tex = 500 K. In the case of the organics in the 12–16 µm range, the C2H2 absorption is abruptly “cut off” by the sili￾cate absorption feature. In addition, because C2H2 is the stronge… view at source ↗
Figure 4
Figure 4. Figure 4: Best-fit LTE slab models beyond 12 µm. The top panel shows the absolute flux and the fitted continuum, whereas the rest of the panels show the continuum-subtracted flux. Molecular hydrogen and atomic lines are indicated with a vertical dashed line. The contribution from warm H2O in the top panel is shown only for reference, as it was not used in the fitting procedure of the organics [PITH_FULL_IMAGE:figur… view at source ↗
Figure 5
Figure 5. Figure 5 [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Stacked 12CO profiles for the NIRSPEC observa￾tions. Light gray lines show the individual profiles while the black line shows the final stacked profile. The dotted vertical lines show the velocity centroid of the absorption component, v0 from [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Resulting CO profiles from the LTE slab fitting. The purple dashed line shows the best-fit absorption slab model for the CO v = 1–0 transitions, each line has been velocity shifted to match the absorption centroid. Transition line centers are shown for reference at the top. The 13CO lines are analyzed in Appendix B. Some additional absorption features not associated with CO are consistent with water lines,… view at source ↗
Figure 8
Figure 8. Figure 8: Comparison of molecular column densities and temperatures for ISO-Oph 37, IRS 46, and GV Tau N. When multiple literature values are available for a given source/species, the absolute column density panel shows the maximum reported value, corresponding to a best-case scenario for hydrocarbon richness. The temperature panel shows the excitation temperatures associated with these maximum column densities. The… view at source ↗
Figure 9
Figure 9. Figure 9: Stacked profiles of individual H2O, C2H2, CH4 and prompt OH transitions. The lines are in the velocity frame of the source (T. Sullivan et al. 2019). Dashed lines located at 0km/s and -50 km/s show the blueshift of the ab￾sorption component. Cieza et al. 2021). This morphology is consistent with the broader ODISEA result that embedded and flat￾spectrum disks typically lack prominent ring-gap sub￾structure … view at source ↗
Figure 10
Figure 10. Figure 10: ALMA and JWST-MIRI images of ISO-Oph 37. Top: ALMA 1.3mm continuum image from ODISEA (L. A. Cieza et al. 2019), CO isotopologues from AGE-PRO (K. Zhang et al. 2025; D. A. Ruiz-Rodriguez et al. 2025). Dotted white lines show the disk surface, traced from the H2 (S3) line emission from MIRI. Middle: H2 line images from MIRI. The black circle shows the inner working angle of JWST, and the contours show the 1… view at source ↗
Figure 11
Figure 11. Figure 11: Line velocity maps of ISO-Oph 37. The top and bottom panels show the same MIRI lines as [PITH_FULL_IMAGE:figures/full_fig_p018_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Schematic diagram of ISO-Oph 37. The thermal continuum serves as a background for some of the species in the wind (CO, C2H2, CH4, H2O, CO2) to absorb against. The thick arrows show the direction of the wind velocity. The molecular labels in the wind are placed to indicate relative velocity, temperature, and covering fraction, and should not be interpreted as a chemical progression along the wind. All spec… view at source ↗
Figure 13
Figure 13. Figure 13: Posterior distribution from the MCMC slab model fitting for the CO v = 1– absorption component. 50 0 50 Velocity (km s 1 ) 0.7 0.8 0.9 1.0 1.1 1.2 Normalized flux 13CO v = 1 0 Stacked profile 0 20 40 60 E /k [K] 32.0 32.5 33.0 33.5 34.0 34.5 ln(N / g ) Fit: Trot = 29 ± 1 K log10 Ntot = 16.30 ± 0.01 13CO v = 1 0 lines [PITH_FULL_IMAGE:figures/full_fig_p025_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Left: 13CO stacked profile. Grey lines show individual line transitions, vertical dotted line shows rest ve￾locity. Right: rotation diagram of 13CO line absorption. blended or minimally blended transitions for each species. For H2O, we adopted unblended rovibra￾tional and rotational transitions within the MIRI wave￾length range from A. Banzatti et al. (2024). For C2H2, the absorption is dominated by relat… view at source ↗
Figure 15
Figure 15. Figure 15: Best-fit SED models for ISO-Oph 37. Each panel shows a different model set, with the icon describing the contributions. Details on the specific model names are explained in T. P. Robitaille (2017). Cieza, L. A., Gonz´alez-Ruilova, C., Hales, A. S., et al. 2021, MNRAS, 501, 2934, doi: 10.1093/mnras/staa3787 Colmenares, M. J., Bergin, E. A., Salyk, C., et al. 2024, ApJ, 977, 173, doi: 10.3847/1538-4357/ad8b… view at source ↗
read the original abstract

We present JWST/MIRI-MRS observations of ISO-Oph 37, a highly inclined flat-spectrum ($\lesssim$1 Myr old) source, to investigate the chemical composition and dynamical origin of its inner-disk gas. The spectrum reveals a rich combination of molecular emission and absorption: H$_2$O, CO, and OH are detected in emission, while strong absorption is observed from CO, H$_2$O, CO$_2$, HCN, C$_2$H$_2$, and CH$_4$, with no detectable ice absorption features. LTE slab modeling of the absorption yields excitation temperatures of $T_{\rm ex}\sim400-600$ K and column densities of $\log N/{\rm cm}^{2}\sim16-19$, characteristic of warm gas located within the inner few au. The absorption lines are significantly blueshifted relative to the systemic velocity, with mid-IR lines exhibiting larger shifts than near-IR CO absorption. This velocity structure points to a velocity- and temperature-stratified molecular disk wind. In this framework, the absorption directly samples disk material lifted from the inner disk surface, preserving the chemical imprint of the wind-launching region. Along the line of sight, ISO-Oph 37 is unusually hydrocarbon-rich compared to other known absorption systems (GV Tau N and IRS 46), exhibiting high (C$_2$H$_2$+CH$_4$)/HCN, (C$_2$H$_2$+CH$_4$)/CO and H$_2$O/CO column density ratios, while the CO and HCN columns remain broadly typical. We find that these molecular ratios are best explained by enhancement of both hydrocarbons and water, driven by inward drift and sublimation of icy pebbles and by thermal processing of carbonaceous grains at the soot line. ISO-Oph 37 thus demonstrates that carbon-rich inner-disk chemistry can be established early in disk evolution and that it can be directly probed through molecular absorption in disk winds.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper presents JWST/MIRI-MRS spectra of the young, highly inclined flat-spectrum source ISO-Oph 37, detecting blueshifted absorption from CO, H2O, CO2, HCN, C2H2, and CH4 (plus emission from H2O, CO, OH) with no ice features. LTE slab modeling yields Tex of 400-600 K and log N of 16-19, interpreted as warm gas in a velocity- and temperature-stratified molecular disk wind launched from the inner few au. The source exhibits elevated (C2H2+CH4)/HCN, (C2H2+CH4)/CO, and H2O/CO column density ratios relative to GV Tau N and IRS 46; these are attributed to inward drift and sublimation of icy pebbles plus thermal processing of carbonaceous grains at the soot line, implying carbon-rich inner-disk chemistry can be established early and probed via wind absorption.

Significance. If the central interpretation holds, the work supplies direct observational evidence that pebble drift and soot-line carbon processing can enhance hydrocarbons and water in the inner disk within the first Myr, offering a new probe of the chemical environment relevant to terrestrial planet formation. The rich molecular inventory and velocity-stratified wind detection in a single source add to the growing JWST sample of disk winds.

major comments (2)
  1. [Velocity structure and LTE slab modeling] The velocity structure section notes that mid-IR absorption lines show larger blueshifts than near-IR CO, indicating temperature-velocity stratification in the wind. However, the LTE slab modeling (single Tex per species, 400-600 K) assumes the absorption samples a uniform composition from the wind-launching region. This is load-bearing for the claim that the observed ratios directly trace unaltered enhancement from pebble drift and soot-line processing; if the line of sight integrates over stratified components with differing local chemistry (e.g., re-condensation or radial variations), the slab-derived columns would not map one-to-one to the proposed mechanism. No multi-component fits or tests of this assumption are described.
  2. [Molecular column density ratios and interpretation] The ratio comparison (Table or figure showing columns for ISO-Oph 37 vs. GV Tau N and IRS 46) is central to the enhancement claim, yet the manuscript provides no quantitative sensitivity analysis on how the derived (C2H2+CH4)/HCN or H2O/CO ratios respond to plausible variations in optical depth treatment, line selection, or temperature gradients. Without this, it is unclear whether the 'best explained by' interpretation is robust or could be reproduced by alternative wind geometries.
minor comments (2)
  1. [Abstract] The abstract omits any mention of line selection criteria, optical depth handling, or goodness-of-fit metrics for the LTE slabs; these details (presumably in the methods) should be briefly summarized to allow readers to assess the column density reliability.
  2. [Results] Notation for column density ratios could be clarified (e.g., explicit definition of summed C2H2+CH4) to avoid ambiguity when comparing across papers.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments, which highlight important aspects of our modeling assumptions and the robustness of our interpretations. We address each major comment below and outline the revisions we will incorporate.

read point-by-point responses

  1. Referee: [Velocity structure and LTE slab modeling] The velocity structure section notes that mid-IR absorption lines show larger blueshifts than near-IR CO, indicating temperature-velocity stratification in the wind. However, the LTE slab modeling (single Tex per species, 400-600 K) assumes the absorption samples a uniform composition from the wind-launching region. This is load-bearing for the claim that the observed ratios directly trace unaltered enhancement from pebble drift and soot-line processing; if the line of sight integrates over stratified components with differing local chemistry (e.g., re-condensation or radial variations), the slab-derived columns would not map one-to-one to the proposed mechanism. No multi-component fits or tests of this assumption are described.

    Authors: We agree that the single-component LTE slab approach is an approximation and that the observed differences in blueshift between mid-IR and near-IR lines indicate velocity stratification within the wind. However, the derived excitation temperatures remain consistent across species (400-600 K), suggesting that the dominant absorbing gas originates from a relatively narrow range of conditions in the inner few au where the wind is launched. This supports our view that the line of sight primarily samples material whose chemistry reflects the inner-disk processing (pebble drift and soot-line sublimation) rather than significant radial or re-condensation variations. To address the concern directly, we will add a dedicated paragraph discussing the limitations of the single-slab assumption, including a simple two-component test using the highest-S/N lines to assess whether differential chemistry could alter the reported ratios. Given the magnitude of the observed enhancements relative to GV Tau N and IRS 46, we expect these tests to confirm that the core interpretation remains robust. revision: partial

  2. Referee: [Molecular column density ratios and interpretation] The ratio comparison (Table or figure showing columns for ISO-Oph 37 vs. GV Tau N and IRS 46) is central to the enhancement claim, yet the manuscript provides no quantitative sensitivity analysis on how the derived (C2H2+CH4)/HCN or H2O/CO ratios respond to plausible variations in optical depth treatment, line selection, or temperature gradients. Without this, it is unclear whether the 'best explained by' interpretation is robust or could be reproduced by alternative wind geometries.

    Authors: We acknowledge that the original manuscript did not include a formal sensitivity analysis for the column-density ratios. The values were obtained via standard LTE slab fits with line selection that avoided obvious blends and accounted for optical-depth effects by prioritizing weaker transitions. To strengthen the result, we will add a new subsection (and associated figure or table) that quantifies the impact of varying Tex by ±100 K, adjusting for possible saturation corrections, and testing alternative line subsets. These tests demonstrate that the elevated (C2H2+CH4)/HCN, (C2H2+CH4)/CO, and H2O/CO ratios persist at levels well above those in the comparison sources. We will also expand the discussion of alternative wind geometries (e.g., more extended or radially varying flows) and show why the combination of blueshifted warm absorption, lack of ice features, and specific molecular inventory favors the pebble-drift plus soot-line scenario over other explanations. revision: yes

Circularity Check

0 steps flagged

No significant circularity; ratios measured from spectra and compared to literature

full rationale

The paper extracts molecular absorption features from new JWST/MIRI-MRS spectra of ISO-Oph 37, performs LTE slab fits to derive Tex and column densities, computes ratios such as (C2H2+CH4)/HCN and H2O/CO, and contrasts them with published values for GV Tau N and IRS 46. The physical interpretation (pebble drift and soot-line processing) is offered as the best explanation for the observed enhancement but does not mathematically reduce any reported quantity to a parameter fitted from the same dataset. Self-citations to prior disk-wind studies exist but are not load-bearing for the central observational result.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The interpretation rests on standard assumptions of LTE slab modeling and the physical picture of a stratified disk wind; no new entities are postulated, but two fitted quantities are introduced by the modeling.

free parameters (2)
  • excitation temperature = 400-600 K
    Fitted per molecule in the LTE slab model to match observed absorption line depths and widths.
  • column density = log N = 16-19 cm^-2
    Fitted per molecule to reproduce the strength of the absorption features.
axioms (2)
  • domain assumption Local thermodynamic equilibrium holds for the absorbing gas slab
    Invoked to convert observed line strengths directly into Tex and column density without solving statistical equilibrium equations.
  • domain assumption Absorption lines trace the wind-launching region without significant line-of-sight contamination or non-LTE effects
    Used to link the observed velocity shifts and molecular ratios directly to inner-disk chemistry.

pith-pipeline@v0.9.0 · 5825 in / 1593 out tokens · 53190 ms · 2026-05-10T16:28:43.272545+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 1 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

  3. [3]

    4^x,m,

    thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

    work page doi:10.1051/0004-6361/201832860 2017