pith. machine review for the scientific record. sign in

arxiv: 2605.07016 · v1 · submitted 2026-05-07 · 🌌 astro-ph.SR

Recognition: no theorem link

Characterizing the Extended Molecular Hydrogen Winds in Protoplanetary Disks from the JWST Disk Infrared Spectroscopic Chemistry Survey

Andrea Banzatti, Colette Salyk, Emma Dahl, Geoffrey A. Blake, Giovanni Rosotti, Ilaria Pascucci, Jane Huang, Joan Najita, Joel Green, Karin Oberg, Ke Zhang, Klaus M. Pontoppidan, L. Ilsedore Cleeves, Mayank Narang, Nicole Arulanantham, Sebastiaan Krijt, Till Kaeufer

Pith reviewed 2026-05-11 01:09 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords protoplanetary disksmolecular hydrogendisk windsJWST spectroscopymass loss ratesdisk dispersalH2 emissionMHD winds
0
0 comments X

The pith

Molecular winds traced by H2 emission can disperse typical protoplanetary disks in 2-3 million years.

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

The paper examines extended pure rotational H2 lines across 34 protoplanetary disks using JWST data and finds clear wind signatures in 16 of them, including monopolar, bipolar, and ring-like structures. These winds are modeled as slow, wide-angle MHD outflows with a median half-opening angle of 45 degrees and a power-law index near 1.6, yielding gas temperatures around 624 K. The derived median total mass-loss rate of 10 to the power of -9 solar masses per year implies that a typical disk of 2-3 Jupiter masses would lose its mass on a 2-3 Myr timescale. This matches observed disk lifetimes and positions the H2 emission as a common tracer for molecular disk winds that operate independently from stellar accretion.

Core claim

Extended H2 emission is common in protoplanetary disks and traces wide-angle molecular winds driven by MHD processes at velocities of about 4 km/s. For ten disks, morphological modeling gives a median half-opening angle of 45 degrees and density power-law index of 1.6, with excitation conditions showing median temperatures of 624 K and total column densities of 10^18.6 cm^{-2}. The resulting wind mass-loss rates cluster tightly around a median log10 value of -9 solar masses per year, indicating that these winds alone could account for the dispersal of typical 2-3 Jupiter-mass disks within 2-3 Myr, consistent with empirical disk lifetimes, while showing little direct correlation with stellar

What carries the argument

Spatially extended pure rotational H2 emission modeled as wide-angle disk winds using half-opening angles near 45 degrees and power-law density profiles to convert line fluxes into total mass-loss rates.

If this is right

  • Molecular winds traced by H2 represent a widespread disk dispersal channel operating across many systems.
  • Typical disks reach dispersal in 2-3 Myr when wind mass loss dominates.
  • Wind mass-loss rates occupy a narrow range of roughly 2 dex and remain decoupled from stellar accretion rates.
  • H2 lines serve as reliable tracers for both inclined and face-on wind geometries.

Where Pith is reading between the lines

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

  • Disk evolution models may need to incorporate MHD wind mass loss as a primary channel alongside photoevaporation.
  • Higher-resolution spectroscopy could map wind launching radii and test whether the narrow rate distribution holds for larger samples.
  • The decoupling from accretion rates suggests that wind-driven dispersal and stellar accretion probe distinct evolutionary phases or episodic behaviors.

Load-bearing premise

The spatially extended H2 emission originates from disk winds rather than other kinematic components and the adopted wind morphology models accurately convert the observed fluxes into total mass-loss rates.

What would settle it

A direct measurement or alternative tracer showing that the H2 emission arises from bound disk gas or non-wind flows, or a median wind mass-loss rate differing by more than an order of magnitude from 10^{-9} solar masses per year.

Figures

Figures reproduced from arXiv: 2605.07016 by Andrea Banzatti, Colette Salyk, Emma Dahl, Geoffrey A. Blake, Giovanni Rosotti, Ilaria Pascucci, Jane Huang, Joan Najita, Joel Green, Karin Oberg, Ke Zhang, Klaus M. Pontoppidan, L. Ilsedore Cleeves, Mayank Narang, Nicole Arulanantham, Sebastiaan Krijt, Till Kaeufer.

Figure 1
Figure 1. Figure 1: The H2 and [Ne II] emission detected towards DoAr 25 on top of the ALMA 240 GHz continuum emission (shown in the background as a reddish-yellow image). The white central circle shows the inner working angle from JWST. In this case, the disk itself is seen in absorption against line emission from the Ophiuchus PDR for H2 S(1)-S(3) (17.03-9.66 µm) [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The observed and median subtracted H2 S(3) lime maps towards WSB 52 on top of the ALMA 240 GHz continuum emission (shown in the background as a reddish-yellow image). The white central circle delineates the inner working angle of the observation. The median subtraction removes the ambient fore/background H2 emission. The orange square represents the 250 au × 250 au aperture from which the line intensities … view at source ↗
Figure 3
Figure 3. Figure 3: (Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram of WSB 52, presented as a representative example. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range d… view at source ↗
Figure 4
Figure 4. Figure 4: The velocity distribution for H2 S(1) line in WSB 52 from the spaxels within the 250 au square aperture used to measure the line fluxes in Section 3.2. Also shown is the best fit gaussian to the velocity distribution. The systemic velocity of the star as measured using H2O lines is shown as the dashed vertical line with the error bar showing the uncertainty in determining the systemic velocity. In [PITH_F… view at source ↗
Figure 5
Figure 5. Figure 5: (a) Schematic illustrating the fitting procedure used to derive the opening angle of the H2 wind. (b) The wide-angle H2 wind emission as observed for DoAr 25. (c) The image rotated to align the wind axis vertically. The outer contour tracing the edge of the flow is shown in lime green, and the magenta points indicate the data selected for model fitting. (d) The resulting best-fit power-law model is shown i… view at source ↗
Figure 6
Figure 6. Figure 6: The ten sources with detected wide angle flows. The color scale is the H2 line image. The dashed line is the best fit curve and the grey lines show the range of fits from the MCMC. The red circle delineates the inner working angle [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The measured total wind mass-loss rates (M˙ tot wind) compared against stellar mass accretion rates obtained from (a) previous literature and (b) HI lines simultaneously ob￾served by MIRI MRS. Sources are color-coded according to their morphology: the ten sources that show clear wide-angle wind structures are shown as red triangles, the six sources that exhibit ring- or bubble-like morphologies are shown a… view at source ↗
Figure 8
Figure 8. Figure 8: 252.3149 252.3141 252.3133 252.3126 RA (J2000) [degrees] -14.3703 -14.3696 -14.3688 -14.3681 Dec (J2000) [degrees] AS 209 / H2 S(1) 0 2.5e-19 5e-19 7.5e-19 1e-18 1.2e-18 1.5e-18 Lin e flu x [W/m2 /s q. a rc s e c] 252.3149 252.3141 252.3133 252.3126 RA (J2000) [degrees] -14.3703 -14.3696 -14.3688 -14.3681 Dec (J2000) [degrees] AS 209 / H2 S(2) 0 9.9e-19 2e-18 3e-18 4e-18 4.9e-18 5.9e-18 Lin e flu x [W/m2 /… view at source ↗
Figure 9
Figure 9. Figure 9: The H2 and [Ne II] emission detected towards AS 209 on top of the ALMA 240 GHz continuum emission (shown in the background as a reddish-yellow image). The white central circle shows the inner working angle from JWST. The same plotting scheme is followed for all the following images in the gallery, and for AS 205 above [PITH_FULL_IMAGE:figures/full_fig_p026_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: 246.5999 246.5991 246.5983 246.5974 RA (J2000) [degrees] -24.7216 -24.7209 -24.7201 -24.7194 Dec (J2000) [degrees] DoAr 25 / H2 S(1) 0 3.3e-19 6.6e-19 9.8e-19 1.3e-18 1.6e-18 2e-18 2.3e-18 Lin e flu x [W/m2 /s q. a rc s e c] 246.5999 246.5991 246.5983 246.5974 RA (J2000) [degrees] -24.7216 -24.7209 -24.7201 -24.7194 Dec (J2000) [degrees] DoAr 25 / H2 S(2) 0 2.9e-19 5.9e-19 8.8e-19 1.2e-18 1.5e-18 1.8e-18 … view at source ↗
Figure 11
Figure 11. Figure 11 [PITH_FULL_IMAGE:figures/full_fig_p027_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: 71.7771 71.7763 71.7755 71.7748 RA (J2000) [degrees] 16.9774 16.9781 16.9788 16.9796 Dec (J2000) [degrees] DR Tau / H2 S(1) 0 2.3e-19 4.5e-19 6.8e-19 9.1e-19 1.1e-18 1.4e-18 1.6e-18 Lin e flu x [W/m2 /s q. a rc s e c] 71.7771 71.7763 71.7755 71.7748 RA (J2000) [degrees] 16.9774 16.9781 16.9788 16.9796 Dec (J2000) [degrees] DR Tau / H2 S(2) 0 6.2e-19 1.2e-18 1.8e-18 2.5e-18 3.1e-18 3.7e-18 4.3e-18 Lin e fl… view at source ↗
Figure 13
Figure 13. Figure 13 [PITH_FULL_IMAGE:figures/full_fig_p028_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: 246.6015 246.6007 246.5999 246.5991 RA (J2000) [degrees] -24.2715 -24.2708 -24.2700 -24.2693 Dec (J2000) [degrees] Elias 24 / H2 S(1) 0 3.2e-19 6.4e-19 9.6e-19 1.3e-18 1.6e-18 Lin e flu x [W/m2 /s q. a rc s e c] 246.6015 246.6007 246.5999 246.5991 RA (J2000) [degrees] -24.2715 -24.2708 -24.2700 -24.2693 Dec (J2000) [degrees] Elias 24 / H2 S(2) 0 5.4e-19 1.1e-18 1.6e-18 2.2e-18 2.7e-18 Lin e flu x [W/m2 /s… view at source ↗
Figure 15
Figure 15. Figure 15 [PITH_FULL_IMAGE:figures/full_fig_p029_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: 68.1336 68.1328 68.1320 68.1312 RA (J2000) [degrees] 24.3329 24.3337 24.3344 24.3352 Dec (J2000) [degrees] FZ Tau / H2 S(1) 0 8.1e-20 1.6e-19 2.4e-19 3.2e-19 4e-19 Lin e flu x [W/m2 /s q. a rc s e c] 68.1336 68.1328 68.1320 68.1312 RA (J2000) [degrees] 24.3329 24.3337 24.3344 24.3352 Dec (J2000) [degrees] FZ Tau / H2 S(2) 0 2.4e-19 4.8e-19 7.2e-19 9.6e-19 1.2e-18 Lin e flu x [W/m2 /s q. a rc s e c] 68.133… view at source ↗
Figure 17
Figure 17. Figure 17 [PITH_FULL_IMAGE:figures/full_fig_p030_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: 73.7971 73.7962 73.7953 73.7945 RA (J2000) [degrees] 30.3654 30.3662 30.3669 30.3677 Dec (J2000) [degrees] GM Aur / H2 S(1) 0 1.6e-20 3.3e-20 4.9e-20 6.6e-20 8.2e-20 9.9e-20 1.1e-19 1.3e-19 1.5e-19 Lin e flu x [W/m2 /s q. a rc s e c] 73.7971 73.7962 73.7953 73.7945 RA (J2000) [degrees] 30.3654 30.3662 30.3669 30.3677 Dec (J2000) [degrees] GM Aur / H2 S(2) 0 2.7e-20 5.4e-20 8e-20 1.1e-19 1.3e-19 1.6e-19 1.… view at source ↗
Figure 19
Figure 19. Figure 19 [PITH_FULL_IMAGE:figures/full_fig_p031_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: 237.3017 237.3008 237.2999 237.2990 RA (J2000) [degrees] -35.6527 -35.6519 -35.6512 -35.6504 Dec (J2000) [degrees] GQ Lup / H2 S(1) 0 8.9e-20 1.8e-19 2.7e-19 3.5e-19 4.4e-19 Lin e flu x [W/m2 /s q. a rc s e c] 237.3017 237.3008 237.2999 237.2990 RA (J2000) [degrees] -35.6527 -35.6519 -35.6512 -35.6504 Dec (J2000) [degrees] GQ Lup / H2 S(2) 0 1.4e-19 2.8e-19 4.2e-19 5.6e-19 7.1e-19 Lin e flu x [W/m2 /s q. … view at source ↗
Figure 21
Figure 21. Figure 21 [PITH_FULL_IMAGE:figures/full_fig_p032_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: 239.1679 239.1671 239.1663 239.1655 RA (J2000) [degrees] -22.0290 -22.0283 -22.0275 -22.0268 Dec (J2000) [degrees] HD 142666 / H2 S(1) 0 4.9e-19 9.7e-19 1.5e-18 1.9e-18 2.4e-18 2.9e-18 3.4e-18 Lin e flu x [W/m2 /s q. a rc s e c] 239.1679 239.1671 239.1663 239.1655 RA (J2000) [degrees] -22.0290 -22.0283 -22.0275 -22.0268 Dec (J2000) [degrees] HD 142666 / H2 S(2) 0 1.3e-18 2.7e-18 4e-18 5.4e-18 6.7e-18 8e-1… view at source ↗
Figure 23
Figure 23. Figure 23 [PITH_FULL_IMAGE:figures/full_fig_p033_23.png] view at source ↗
Figure 24
Figure 24. Figure 24: 269.0899 269.0891 269.0883 269.0875 RA (J2000) [degrees] -21.9574 -21.9566 -21.9559 -21.9552 Dec (J2000) [degrees] HD 163296 / H2 S(1) 0 8.2e-19 1.6e-18 2.5e-18 3.3e-18 4.1e-18 4.9e-18 5.8e-18 Lin e flu x [W/m2 /s q. a rc s e c] 269.0899 269.0891 269.0883 269.0875 RA (J2000) [degrees] -21.9574 -21.9566 -21.9559 -21.9552 Dec (J2000) [degrees] HD 163296 / H2 S(2) 0 6.7e-18 1.3e-17 2e-17 2.7e-17 3.4e-17 4e-1… view at source ↗
Figure 25
Figure 25. Figure 25 [PITH_FULL_IMAGE:figures/full_fig_p034_25.png] view at source ↗
Figure 26
Figure 26. Figure 26: 236.3049 236.3040 236.3031 236.3022 RA (J2000) [degrees] -34.2931 -34.2923 -34.2916 -34.2908 Dec (J2000) [degrees] HT Lup / H2 S(1) 0 1.6e-19 3.1e-19 4.7e-19 6.2e-19 7.8e-19 Lin e flu x [W/m2 /s q. a rc s e c] 236.3049 236.3040 236.3031 236.3022 RA (J2000) [degrees] -34.2931 -34.2923 -34.2916 -34.2908 Dec (J2000) [degrees] HT Lup / H2 S(2) 0 3.6e-19 7.3e-19 1.1e-18 1.5e-18 1.8e-18 Lin e flu x [W/m2 /s q. … view at source ↗
Figure 27
Figure 27. Figure 27 [PITH_FULL_IMAGE:figures/full_fig_p035_27.png] view at source ↗
Figure 28
Figure 28. Figure 28: 70.4130 70.4122 70.4114 70.4106 RA (J2000) [degrees] 25.9395 25.9403 25.9410 25.9418 Dec (J2000) [degrees] IRAS 04385 / H2 S(1) 0 6.9e-20 1.4e-19 2.1e-19 2.8e-19 3.5e-19 Lin e flu x [W/m2 /s q. a rc s e c] 70.4130 70.4122 70.4114 70.4106 RA (J2000) [degrees] 25.9395 25.9403 25.9410 25.9418 Dec (J2000) [degrees] IRAS 04385 / H2 S(2) 0 1.5e-19 3e-19 4.4e-19 5.9e-19 7.4e-19 Lin e flu x [W/m2 /s q. a rc s e c… view at source ↗
Figure 29
Figure 29. Figure 29 [PITH_FULL_IMAGE:figures/full_fig_p036_29.png] view at source ↗
Figure 30
Figure 30. Figure 30: 74.6940 74.6932 74.6923 74.6915 RA (J2000) [degrees] 29.8425 29.8432 29.8440 29.8447 Dec (J2000) [degrees] MWC 480 / H2 S(1) 0 3.6e-19 7.3e-19 1.1e-18 1.5e-18 1.8e-18 Lin e flu x [W/m2 /s q. a rc s e c] 74.6940 74.6932 74.6923 74.6915 RA (J2000) [degrees] 29.8425 29.8432 29.8440 29.8447 Dec (J2000) [degrees] MWC 480 / H2 S(2) 0 1.2e-18 2.5e-18 3.7e-18 4.9e-18 6.1e-18 Lin e flu x [W/m2 /s q. a rc s e c] 74… view at source ↗
Figure 31
Figure 31. Figure 31 [PITH_FULL_IMAGE:figures/full_fig_p037_31.png] view at source ↗
Figure 32
Figure 32. Figure 32: 239.8697 239.8687 239.8677 239.8668 RA (J2000) [degrees] -40.3655 -40.3647 -40.3640 -40.3632 Dec (J2000) [degrees] RY Lup / H2 S(1) 0 1.2e-19 2.3e-19 3.5e-19 4.6e-19 5.8e-19 Lin e flu x [W/m2 /s q. a rc s e c] 239.8697 239.8687 239.8677 239.8668 RA (J2000) [degrees] -40.3655 -40.3647 -40.3640 -40.3632 Dec (J2000) [degrees] RY Lup / H2 S(2) 0 2.8e-19 5.7e-19 8.5e-19 1.1e-18 1.4e-18 Lin e flu x [W/m2 /s q. … view at source ↗
Figure 33
Figure 33. Figure 33 [PITH_FULL_IMAGE:figures/full_fig_p038_33.png] view at source ↗
Figure 34
Figure 34. Figure 34: 164.1312 164.1278 164.1245 164.1211 RA (J2000) [degrees] -77.1953 -77.1946 -77.1939 -77.1931 Dec (J2000) [degrees] SY Cha / H2 S(1) 0 4e-20 8.1e-20 1.2e-19 1.6e-19 2e-19 Lin e flu x [W/m2 /s q. a rc s e c] 164.1312 164.1278 164.1245 164.1211 RA (J2000) [degrees] -77.1953 -77.1946 -77.1939 -77.1931 Dec (J2000) [degrees] SY Cha / H2 S(2) 0 6.3e-20 1.3e-19 1.9e-19 2.5e-19 3.2e-19 Lin e flu x [W/m2 /s q. a rc… view at source ↗
Figure 35
Figure 35. Figure 35 [PITH_FULL_IMAGE:figures/full_fig_p039_35.png] view at source ↗
Figure 36
Figure 36. Figure 36: 239.8201 239.8191 239.8181 239.8171 RA (J2000) [degrees] -41.9541 -41.9533 -41.9526 -41.9519 Dec (J2000) [degrees] Sz 129 / H2 S(1) 0 4e-20 8e-20 1.2e-19 1.6e-19 2e-19 Lin e flu x [W/m2 /s q. a rc s e c] 239.8201 239.8191 239.8181 239.8171 RA (J2000) [degrees] -41.9541 -41.9533 -41.9526 -41.9519 Dec (J2000) [degrees] Sz 129 / H2 S(2) 0 4.2e-20 8.3e-20 1.2e-19 1.7e-19 2.1e-19 Lin e flu x [W/m2 /s q. a rc s… view at source ↗
Figure 37
Figure 37. Figure 37 [PITH_FULL_IMAGE:figures/full_fig_p040_37.png] view at source ↗
Figure 38
Figure 38. Figure 38: 165.4673 165.4664 165.4655 165.4646 RA (J2000) [degrees] -34.7059 -34.7052 -34.7044 -34.7037 Dec (J2000) [degrees] TW Hya / H2 S(1) 0 2.5e-19 5e-19 7.5e-19 1e-18 1.3e-18 Lin e flu x [W/m2 /s q. a rc s e c] 165.4673 165.4664 165.4655 165.4646 RA (J2000) [degrees] -34.7059 -34.7052 -34.7044 -34.7037 Dec (J2000) [degrees] TW Hya / H2 S(2) 0 1.5e-19 3.1e-19 4.6e-19 6.1e-19 7.6e-19 Lin e flu x [W/m2 /s q. a rc… view at source ↗
Figure 39
Figure 39. Figure 39 [PITH_FULL_IMAGE:figures/full_fig_p041_39.png] view at source ↗
Figure 40
Figure 40. Figure 40: 246.9155 246.9147 246.9139 246.9131 RA (J2000) [degrees] -24.6554 -24.6547 -24.6539 -24.6532 Dec (J2000) [degrees] WSB 52 / H2 S(1) 0 2.6e-19 5.2e-19 7.8e-19 1e-18 1.3e-18 Lin e flu x [W/m2 /s q. a rc s e c] 246.9155 246.9147 246.9139 246.9131 RA (J2000) [degrees] -24.6554 -24.6547 -24.6539 -24.6532 Dec (J2000) [degrees] WSB 52 / H2 S(2) 0 5.7e-19 1.1e-18 1.7e-18 2.3e-18 2.9e-18 Lin e flu x [W/m2 /s q. a … view at source ↗
Figure 41
Figure 41. Figure 41: B. ROTATION DIAGRAMS This section of the appendix contains the rotation diagrams for each of our targets [PITH_FULL_IMAGE:figures/full_fig_p042_41.png] view at source ↗
Figure 42
Figure 42. Figure 42: (Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:f… view at source ↗
Figure 43
Figure 43. Figure 43: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
Figure 44
Figure 44. Figure 44: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
Figure 45
Figure 45. Figure 45: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
Figure 46
Figure 46. Figure 46: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
Figure 47
Figure 47. Figure 47: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
Figure 48
Figure 48. Figure 48: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
Figure 49
Figure 49. Figure 49: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
Figure 50
Figure 50. Figure 50: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
Figure 51
Figure 51. Figure 51: Left) Corner plot showing the posterior distributions and best-fit parameters from the MCMC fit to the rotation diagram. (Right) Rotation diagram with observed data points in green and extinction-corrected points in blue for ortho-H2 and red for para-H2. The black line denotes the best-fit model, and the shaded grey region represents the uncertainty range derived from the MCMC fitting [PITH_FULL_IMAGE:fi… view at source ↗
read the original abstract

We present a comprehensive analysis of extended H$_2$ emission from 34 protoplanetary disks observed with the JWST Disk Infrared Spectroscopic Chemistry Survey (JDISCS), supplemented by archival data. We investigated the morphology, kinematics, excitation conditions, and mass dynamics of H$_2$. Extended emission from pure rotational H$_2$ lines is found to be common, with 16 sources exhibiting clear signatures of disk winds. These include monopolar and bipolar structures in inclined disks and ring-like or bubble-like morphologies in face-on systems features indicative of wide-angle disk winds. Our analysis shows that the H$_2$ is consistent with slow {(4.2$^{+6.7}_{-3.0}$ km s$^{-1}$)} MHD driven winds. For ten disks, we model the wind morphology and find a median half-opening angle of $45\arcdeg^{+5}_{-4}$ and a characteristic power-law index of $\alpha \sim$ 1.6. Excitation analysis yields a median gas temperature of 624 $\pm$ 130 K and a column density of $\log(N_{\mathrm{tot}}\,[\mathrm{cm}^{-2}]) = 18.6 \pm 0.6$. The median wind mass-loss rate, ${\rm log_{10}}(\dot{\rm M}_{\rm wind}^{\rm tot}) = -9_{-0.4}^{+0.8}\,{\rm M_\odot\,yr^{-1}}$, implies that, if molecular winds are the dominant mechanism responsible for disk dispersal, a typical disk with a mass of $2-3\,M_{\rm Jup}$ would dissipate on a $\sim$2-3 Myr timescale, consistent with observed disk lifetimes. The $\dot{\rm M}_{\mathrm{\rm wind}}^{\rm tot}$ span a relatively narrow range ($\sim$2 dex) and do not correlate strongly with accretion rates onto the star, suggesting that the mass loss rate and the accretion rates are probing different timescales. Our findings demonstrate that spatially extended warm H$_2$ emission is a widespread and reliable tracer of molecular disk winds in protoplanetary systems.

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 manuscript analyzes extended pure-rotational H2 emission in 34 protoplanetary disks from the JWST JDISCS survey (plus archival data). It reports clear wind signatures in 16 sources, including monopolar/bipolar structures and ring/bubble morphologies. For ten disks, wind morphology is modeled with a median half-opening angle of 45°+5−4 and power-law index α∼1.6; excitation analysis yields median T=624±130 K and log N_tot=18.6±0.6. The resulting median log10(Ṁ_wind^tot)=−9−0.4+0.8 M⊙ yr−1 is used to argue that molecular winds can disperse a typical 2–3 MJup disk in ∼2–3 Myr, consistent with observed lifetimes. The rates span ∼2 dex and show no strong correlation with stellar accretion rates, supporting slow MHD winds as a widespread dispersal mechanism.

Significance. If the mass-loss rates prove robust, the work supplies one of the largest observational samples of molecular disk winds to date and directly links JWST-detected H2 structures to disk-evolution timescales. The narrow range of Ṁ values and the decoupling from accretion rates are potentially important constraints on wind-launching physics. The survey-scale detection statistics strengthen the case that extended warm H2 is a reliable wind tracer.

major comments (2)
  1. [Wind morphology modeling for the ten disks (values and assumptions stated in the abstract and corresponding analysis)] The headline implication (2–3 Myr dissipation for a 2–3 MJup disk) rests on the median log10(Ṁ_wind^tot)=−9 reported for the ten modeled sources. This value is obtained only after fixing the half-opening angle at 45°+5−4 and the density power-law index at α∼1.6 to convert observed line fluxes into total mass-loss rates. No sensitivity analysis quantifies how Ṁ (and therefore the timescale) changes when these parameters are varied within plausible ranges or when the wind is allowed to be narrower or steeper.
  2. [Kinematic and morphological classification of the 16 wind sources] The claim that the spatially extended H2 traces disk winds (rather than bound Keplerian emission, shocks, or other kinematic components) is central to interpreting all 16 detections and the derived Ṁ values. While morphologies are described as “indicative of wide-angle winds,” the manuscript does not present a quantitative test or upper limit on non-wind contributions to the line fluxes.
minor comments (2)
  1. [Abstract] The abstract quotes a median velocity of 4.2+6.7−3.0 km s−1 but does not indicate the section or figure from which this value and its asymmetric uncertainty are derived.
  2. [Results section on modeled sources] A summary table listing the individual half-opening angles, α values, temperatures, column densities, and Ṁ for each of the ten modeled sources would allow readers to evaluate the robustness of the reported medians.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed review. The comments have prompted us to strengthen the robustness of our mass-loss rate derivations and the interpretation of the H2 emission. We address each major comment below and have revised the manuscript accordingly.

read point-by-point responses
  1. Referee: [Wind morphology modeling for the ten disks (values and assumptions stated in the abstract and corresponding analysis)] The headline implication (2–3 Myr dissipation for a 2–3 MJup disk) rests on the median log10(Ṁ_wind^tot)=−9 reported for the ten modeled sources. This value is obtained only after fixing the half-opening angle at 45°+5−4 and the density power-law index at α∼1.6 to convert observed line fluxes into total mass-loss rates. No sensitivity analysis quantifies how Ṁ (and therefore the timescale) changes when these parameters are varied within plausible ranges or when the wind is allowed to be narrower or steeper.

    Authors: We agree that an explicit sensitivity analysis strengthens the headline result. The reported median half-opening angle and α were derived directly from morphological fits to the ten disks rather than arbitrarily fixed, but the original manuscript did not propagate plausible variations in these parameters into the Ṁ uncertainties. In the revised manuscript we have added a dedicated sensitivity subsection (Section 4.3) that re-derives Ṁ while allowing the half-opening angle to range from 30° to 60° and α from 1.0 to 2.2. The resulting median log10(Ṁ_wind^tot) shifts by at most 0.4 dex, leaving the 2–3 Myr dispersal timescale unchanged within the quoted uncertainties. We have also updated the abstract and discussion to quote the enlarged error budget that incorporates this analysis. revision: yes

  2. Referee: [Kinematic and morphological classification of the 16 wind sources] The claim that the spatially extended H2 traces disk winds (rather than bound Keplerian emission, shocks, or other kinematic components) is central to interpreting all 16 detections and the derived Ṁ values. While morphologies are described as “indicative of wide-angle winds,” the manuscript does not present a quantitative test or upper limit on non-wind contributions to the line fluxes.

    Authors: We acknowledge that the original text relied on qualitative morphological and kinematic descriptors without a formal quantitative decomposition. In the revision we have added a new paragraph in Section 3.2 that compares the observed spatial profiles and velocity fields against simple Keplerian disk models and shock templates. The extended emission is shown to be inconsistent with bound Keplerian rotation at >3σ in the majority of sources, and we now quote an estimated upper limit of ~25% on any non-wind contribution to the total line flux based on the radial surface-brightness distribution. A full multi-component radiative-transfer decomposition is not feasible with the current spectral resolution and is noted as a limitation for future work; the revised discussion is accordingly more cautious while retaining the wind interpretation as the dominant component. revision: partial

Circularity Check

0 steps flagged

No circularity: mass-loss rates derived from data-driven modeling, timescale follows by direct division

full rationale

The paper observes extended H2 lines, fits wind morphology (half-opening angle ~45°, α~1.6) and excitation parameters to the fluxes of ten sources to obtain per-source Mdot values, then reports the median and computes a dissipation timescale as typical disk mass divided by that median Mdot. This is a standard forward modeling chain from data to derived quantities; the output Mdot is not an input, the timescale is not renamed or fitted, and no self-citation, uniqueness theorem, or ansatz is invoked to force the result. The derivation remains self-contained against the observed line fluxes and standard wind assumptions.

Axiom & Free-Parameter Ledger

5 free parameters · 2 axioms · 0 invented entities

The central claim rests on the interpretation of observed line fluxes as wind emission and on several fitted morphological and excitation parameters; no new physical entities are postulated.

free parameters (5)
  • median half-opening angle = 45 deg
    Fitted from modeling 10 disks
  • power-law index alpha = 1.6
    Characteristic value from wind morphology models
  • median gas temperature = 624 K
    Derived from excitation analysis of H2 lines
  • median total column density = log N_tot = 18.6 cm^-2
    From excitation diagrams
  • median wind mass-loss rate = log Mdot = -9 M_sun/yr
    Converted from observed fluxes using wind models
axioms (2)
  • domain assumption Extended H2 emission traces disk winds
    Invoked when classifying monopolar, bipolar, ring-like, and bubble-like morphologies as wind signatures
  • domain assumption Winds are MHD-driven
    Inferred from the low observed velocities (4.2 km/s median)

pith-pipeline@v0.9.0 · 5781 in / 1553 out tokens · 38178 ms · 2026-05-11T01:09:48.996685+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

31 extracted references · 31 canonical work pages

  1. [1]

    The winds appears to originate primarily from the northern component, AS 205N

    AS 205 shows prominent and spatially extended H 2 emission, particularly in the strong ortho H2 lines, consistent with a wide-angle wind. The winds appears to originate primarily from the northern component, AS 205N. Given the close separation of the two components, the wind may also interact with the secondary star, potentially influencing the observed s...

  2. [2]

    The observed H 2 morphology is consistent with this viewing geometry, with emission symmetric around the central source

    AS 209 is a moderately inclined disk. The observed H 2 morphology is consistent with this viewing geometry, with emission symmetric around the central source. Residual features associated with PSF subtraction are present in the immediate vicinity of the central source, but the extended emission beyond this region remains clearly detectable mainly in S(1),...

  3. [3]

    Only the blue-shifted side of the wind is detected, likely because the red-shifted emission is obscured by the disk

    CI Tau exhibits a clear wide-angle wind. Only the blue-shifted side of the wind is detected, likely because the red-shifted emission is obscured by the disk. The wind is detected in the S(1), S(2), S(3), and S(5) lines. The S(4) line is faint and is dominated by PSF artifacts close to the center. The monopolar wind shows has a large opening angle. Complem...

  4. [4]

    The uniform background line emission originates from the photon-dominated region (PDR) near the Ophiuchus core

    DoAr 25 stands out in the sample because its dust disk is observed in absorption against bright background H 2 emission, particularly in the S(1), S(2), and S(3) lines. The uniform background line emission originates from the photon-dominated region (PDR) near the Ophiuchus core. Superimposed on the background emission, the H2 data show clear signatures o...

  5. [5]

    We detect extended H 2 emission associated with the source in the S(1), S(2), and S(3) lines

    DoAr 33 shows strong, uniform background H 2 emission in the S(1) and S(2) lines. We detect extended H 2 emission associated with the source in the S(1), S(2), and S(3) lines. The S(4) line shows little to no evidence of wind-like emission, while the S(5) line is affected by residal instrumental artifacts and PSF subtraction effects, making detailed inter...

  6. [6]

    This is consistent with the face-on orientation of the system

    DR Tau is a nearly face-on disk, which exhibits faint extended H2 emission, and a symmetric morphology around the central source. This is consistent with the face-on orientation of the system. However, the faint nature of the extended component and limited S/N prevents further characterization of the spatial structure of the emission

  7. [7]

    The background H 2 emission exhibits a gradient that increases toward the northeast, reflecting the structured nature of the PDR environment

    Elias 20 is superimposed on bright background H 2 line emission associated with the the Ophiuchus core PDR. The background H 2 emission exhibits a gradient that increases toward the northeast, reflecting the structured nature of the PDR environment. The source displays a monopolar outflow, particularly prominent in the higher-J transitions (S(3)–S(5)), wh...

  8. [8]

    The ring is asymmetric, extending more along the northeast–southwest direction and less along the east–west direction

    Elias 24 shows a ring-like morphology with a bright knot to the northeast in all H 2 lines. The ring is asymmetric, extending more along the northeast–southwest direction and less along the east–west direction. The bright northeastern knot stands out as a localized enhancement in H2 emission, suggesting a region of higher excitation, density, or possible ...

  9. [9]

    In the S(3) and S(5) lines, a monopolar outflow extending from the disk along the northeast side is observed

    Elias 27 is observed in silhouette against strong background H 2 line emission from the Ophiuchus PDR in the S(1) and S(2) lines due to the dust in the disk. In the S(3) and S(5) lines, a monopolar outflow extending from the disk along the northeast side is observed. Extended [Ne II] emission is observed, indicating the presence of a collimated, ionized j...

  10. [10]

    FZ Tau is a moderately inclined system. It represents one of the clearest examples of ring-like H 2 emission in our sample, with an average radius of∼2.4 ′′ (310 au) and a width of about 1.15 ′′ (148 au) (Pontoppidan et al. 2024b). A localized region of enhanced emission is visible toward the southwest portion of the ring, suggesting possible asymmetry in...

  11. [11]

    The emission is diffuse, without a clear preferred orientation

    GK Tau shows extended H 2 emission in the S(1), S(2), S(3), and S(5) lines. The emission is diffuse, without a clear preferred orientation. There are no obvious morphological signatures indicative of a collimated wind. The lack of structure suggests that any molecular wind, if present, is intrinsically weak with a wide opening angle

  12. [12]

    GM Aur exhibits a monopolar outflow oriented perpendicular to the dust disk plane detected in all observed H 2 lines. The emission traces a blue-shifted wide-angle structure extending away from the central source, consistent with a molecular wind emerging from the disk surface, while the red-shifted emission is obscured by the dust disk. GM Aur one of the...

  13. [13]

    In contrast, the S(1) line reveals more discernible structure, suggestive of a wide-angle outflow and a disk silhouette with major axis oriented north-east/south-west

    GO Tau shows weak extended emission in the S(3), S(4), and S(5) lines. In contrast, the S(1) line reveals more discernible structure, suggestive of a wide-angle outflow and a disk silhouette with major axis oriented north-east/south-west

  14. [14]

    Some PSF subtraction artifacts are present, particularly in the S(5) line, where a bright artifact is visible to the north of the source

    GQ Lup shows tentative extended emission with a wide-angle morphology in the S(1), S(2), and S(3) lines. Some PSF subtraction artifacts are present, particularly in the S(5) line, where a bright artifact is visible to the north of the source

  15. [15]

    The extended emission is concentrated around the star

    GW Lup shows bright H 2 emission surrounding the central source in the S(1) and S(2) lines. The extended emission is concentrated around the star. Some of the bright features seen in the S(2) line near the central source may be caused by PSF subtraction artifacts

  16. [16]

    The apparent bright emission surrounding the star is consistent with PSF subtraction artifacts

    HD 142666, HD 143006, HD 163296, and MWC 480 are Herbig Ae stars and show no extended H 2 emission in any of the observed lines. The apparent bright emission surrounding the star is consistent with PSF subtraction artifacts

  17. [17]

    The tail is prominent and coherent, suggesting material flowing or being entrained along this direction

    HP Tau exhibits a striking ring-like H 2 morphology with a tail extending toward the west. The tail is prominent and coherent, suggesting material flowing or being entrained along this direction. In addition, the [Ne II] emission is also spatially extended along the western tail, indicating that both warm molecular and ionized gas trace the same outflow o...

  18. [18]

    The primary component is associated with extended H2 emission, although not with an classical wind-like morphology

    HT Lup is a triple system with bright continuum emission. The primary component is associated with extended H2 emission, although not with an classical wind-like morphology. A tentative wide-angle wind component becomes more discernible in the S(3), S(4), and S(5) lines. This suggests that the primary hosts a molecular wind that is partially disrupted by ...

  19. [19]

    The northwestern lobe of the outflow has lower S/N in lines higher than S(1) and S(2)

    IQ Tau is an inclined source that exhibits a clear bipolar H 2 outflow oriented perpendicular to the dust disk plane. The northwestern lobe of the outflow has lower S/N in lines higher than S(1) and S(2). The [Ne II] emission appears slightly extended along the outflow direction, suggesting that both warm molecular and ionized gas trace the outflow structure. 25

  20. [20]

    The bipolar structure is well-defined across the observed H 2 lines

    IRAS 04385+2550 is moderately inclined (i= 60 ◦) and presents an unambiguous example of wide-angle bipolar H2 emission, reminiscent of emission seen in protostars, though with a significantly wider opening angle. The bipolar structure is well-defined across the observed H 2 lines. This demonstrates that bipolar winds can be visible even in sources that ar...

  21. [21]

    The higher-Jtransitions have lower S/N

    MY Lup exhibits a clear bipolar morphology in the H 2 S(1) line, with tentative evidence for the S(2) line. The higher-Jtransitions have lower S/N. Weak bipolar emission is detected in the S(3) line, while the S(4) line shows little to no detectable emission. The [Ne II] emission is marginally extended

  22. [22]

    The symmetric structure is consistent with the nearly face-on orientation of this system

    RU Lup shows bright, but compact, emission close to the central source in the S(1), S(2), and S(3) transitions, suggesting that any extended molecular wind is confined to small spatial scales. The symmetric structure is consistent with the nearly face-on orientation of this system

  23. [23]

    RY Lup is unique within the sample, exhibiting H 2 emission both parallel and perpendicular to the disk plane. The perpendicular component likely traces a molecular outflow or wind emerging from the disk surface, whereas the emission aligned with the disk may originate from cooler, quiescent disk gas. This dual morphology provides a rare view of both the ...

  24. [24]

    No definitive sig- natures of outflows or winds are detected in the H 2 emission, except in the S(5) transition

    SR 4 is dominated by bright background H 2 emission originating from a nearby PDR in the Ophiuchus core, showing a clear spatial gradient increasing from the southwestern to the northeastern side. No definitive sig- natures of outflows or winds are detected in the H 2 emission, except in the S(5) transition. Extended [Ne II] emission is centered on the so...

  25. [25]

    SY Cha exhibits a conical monopolar H 2 outflow detected in all observed transitions, oriented toward the east (Schwarz et al. 2025a). Emission from the western side is observed, albeit significantly weaker. Additionally, SY Cha shows an extended [Ne II] jet on both sides of the disk, aligned along the wide-angle H 2 wind, indicating the presence of ioniz...

  26. [26]

    The extended emission is concentrated around the central source, consistent with a face-on disk wind, and is detected in all observed H 2 lines

    Sz 114 is a nearly face-on source (i= 20 ◦). The extended emission is concentrated around the central source, consistent with a face-on disk wind, and is detected in all observed H 2 lines

  27. [27]

    Extended emission from S(1), S(2), and S(3) H 2 lines is concentrated around the central source, consistent with a face-on disk wind

    Sz 129 is a low-inclination source (i= 34 ◦). Extended emission from S(1), S(2), and S(3) H 2 lines is concentrated around the central source, consistent with a face-on disk wind. The [Ne II] emission is also marginally extended

  28. [28]

    The [Ne II] emission also appears extended, albeit more compact

    TW Cha displays a centrally concentrated emission from the S(1), S(2), S(3), and S(5) lines, distributed sym- metrically around the central source, consistent with a face-on disk wind. The [Ne II] emission also appears extended, albeit more compact

  29. [29]

    The emission appears symmetric around the central source and likely traces the wide-angle disk wind

    TW Hya shows a clear, symmetric, extended H2 morphology in all observed lines, consistent with a nearly face-on disk wind. The emission appears symmetric around the central source and likely traces the wide-angle disk wind. The [Ne II] emission also appears extended, albeit more compact

  30. [30]

    The bubble has a brighter rim and in the inner region is less bright

    VZ Cha exhibits an unusual H 2 outflow morphology characterized by a bubble-like structure on the northwestern side of the source. The bubble has a brighter rim and in the inner region is less bright. The bubble component appears offset from the central source, possibly indicating a past episodic ejection event. The [Ne II] emission clearly reveals a jet ...

  31. [31]

    The S(1) transition is dominated by bright, extended emission, while the higher-excitation H 2 lines reveal a weaker but discernible bipolar morphology

    WSB 52 exhibits a wide-angle H 2 outflow. The S(1) transition is dominated by bright, extended emission, while the higher-excitation H 2 lines reveal a weaker but discernible bipolar morphology. The southwestern side of the outflow appears brighter than the northeastern side, indicating asymmetry in the emission or excitation conditions. The source also s...