Abstract
Disk winds are thought to play a critical role in the evolution and dispersal of protoplanetary disks. A primary diagnostic of this physics is emission from the wind, especially in the low-velocity component of the [O I] λ6,300 line. However, the interpretation of the line is usually based on spectroscopy alone, which leads to confusion between magnetohydrodynamic winds and photoevaporative winds. Here we report that in high-resolution spectral mapping of TW Hya by the multi-unit spectroscopic explorer at the Very Large Telescope, 80% of the [O I] emission is confined to within 1 au radially from the star. A generic model of a magnetothermal wind produces [O I] emission at the base of the wind that broadly matches the flux and the observed spatial and spectral profiles. The emission at large radii is much fainter than predicted from models of photoevaporation, perhaps because the magnetothermal wind partially shields the outer disk from energetic radiation from the central star. This result calls into question the previously assessed importance of photoevaporation in disk dispersal predicted by models.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 digital issues and online access to articles
$119.00 per year
only $9.92 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
The raw MUSE data can be taken from ESO archive under programme ID no. 0104.C-0449. The combined MUSE data used in this research are available in Figshare (https://doi.org/10.6084/m9.figshare.21906306). The reduced spectral data from other ESO facilities can be downloaded through ESO Phase III archive (https://archive.eso.org/wdb/wdb/adp/phase3_spectral/form) with programme ID nos. 074.A-9021(A), 078.A-9059(A), 079.A-9006(A), 079.A-9007(A), 079.A-9017(A), 081.A-9005(A), 081.A-9023(A), 089.A-9007(A), 0101.A-9012(A), 0102.A-9008(A), 60.A-9036(A), 075.C-0202(A), 081.C-0779(A), 082.C-0390(A), 082.C-0427(C), 65.I-0404(A), 082.C-0218(A), 089.C-0299(A) and 106.20Z8.011. The reduced Keck spectral data were downloaded from the Keck Observatory Archive (https://www2.keck.hawaii.edu/koa/public/koa.php) with programme ID nos. C179Hr, C269Hr, C199Hb, N125Hr, C186Hr, N204Hr, C189Hr, C252Hr and C247Hr. The reduced ESPaDOnS spectra were downloaded from PolarBase (http://polarbase.irap.omp.eu/).
Code availability
The MUSE data are reduced with the MUSE pipeline v.2.8.5, which is publicly available. Upon request, the first author will provide code (IDL) to generate figures.
References
Pascucci, I. et al. in Protostars and Planets VII Ch. 16 (in the press).
Frank, A. et al. in Protostars and Planets VI (eds Beuther, H. et al.) 451–474 (Univ. Arizona Press, 2014).
Ferreira, J., Dougados, C. & Cabrit, S. Which jet launching mechanism(s) in T Tauri stars? Astron. Astrophys. 453, 785–796 (2006).
Bai, X.-N. Global simulations of the inner regions of protoplanetary disks with comprehensive disk microphysics. Astrophys. J. 845, 75 (2017).
Wang, L., Bai, X.-N. & Goodman, J. Global simulations of protoplanetary disk outflows with coupled non-ideal magnetohydrodynamics and consistent thermochemistry. Astrophys. J. 874, 90 (2019).
Hartigan, P., Edwards, S. & Ghandour, L. Disk accretion and mass loss from young stars. Astrophys. J. 452, 736 (1995).
Rigliaco, E. et al. Understanding the origin of the [O I] low-velocity component from T Tauri stars. Astrophys. J. 772, 60 (2013).
Natta, A. et al. X-shooter spectroscopy of young stellar objects. V. Slow winds in T Tauri stars. Astron. Astrophys. 569, A5 (2014).
Simon, M. N. et al. Tracing slow winds from T Tauri stars via low-velocity forbidden line emission. Astrophys. J. 831, 169 (2016).
Gaia Collaboration et al. Gaia early data release 3 – summary of the contents and survey properties. Astron. Astrophys. 649, A1 (2021).
Qi, C. et al. Imaging the disk around TW Hydrae with the submillimeter array. Astrophys. J. Lett. 616, L11–L14 (2004).
Hughes, A. M. et al. Empirical constraints on turbulence in protoplanetary accretion disks. Astrophys. J. 727, 85 (2011).
Teague, R. et al. Measuring turbulence in TW Hydrae with ALMA: methods and limitations. Astron. Astrophys. 592, A49 (2016).
Huang, J. et al. CO and dust properties in the TW Hya disk from high-resolution ALMA observations. Astrophys. J. 852, 122 (2018).
Andrews, S. M. et al. Ringed substructure and a gap at 1 AU in the nearest protoplanetary disk. Astrophys. J. Lett. 820, L40 (2016).
Calvet, N. et al. Evidence for a developing gap in a 10 Myr old protoplanetary disk. Astrophys. J. 568, 1008–1016 (2002).
GRAVITY Collaboration A measure of the size of the magnetospheric accretion region in TW Hydrae. Nature 584, 547–550 (2020).
Pascucci, I. & Sterzik, M. Evidence for disk photoevaporation driven by the central star. Astrophys. J. 702, 724–732 (2009).
Pascucci, I. et al. The photoevaporative wind from the disk of TW Hya. Astrophys. J. 736, 13 (2011).
Fang, M. et al. A new look at T Tauri star forbidden lines: MHD-driven winds from the inner disk. Astrophys. J. 868, 28 (2018).
Banzatti, A. et al. Kinematic links and the coevolution of MHD winds, jets, and inner disks from a high-resolution optical [O I] survey. Astrophys. J. 870, 76 (2019).
Pascucci, I. et al. The evolution of disk winds from a combined study of optical and infrared forbidden lines. Astrophys. J. 903, 78 (2020).
Yang, H., Johns-Krull, C. M. & Valenti, A. Spectropolarimetry of the classical T Tauri star TW Hydrae. Astron. J. 133, 73–80 (2007).
Alexander, R. D., Clarke, C. J. & Pringle, J. E. Photoevaporation of protoplanetary discs – I. Hydrodynamic models. Mon. Not. R. Astron. Soc. 369, 216–228 (2006).
Alexander, R. D., Clarke, C. J. & Pringle, J. E. Photoevaporation of protoplanetary discs – II. Evolutionary models and observable properties. Mon. Not. R. Astron. Soc. 369, 229–239 (2006).
Gorti, U. & Hollenbach, D. Line emission from gas in optically thick dust disks around young stars. Astrophys. J. 683, 287–303 (2008).
Gorti, U. & Hollenbach, D. Photoevaporation of circumstellar disks by far-ultraviolet, extreme-ultraviolet and X-ray radiation from the central star. Astrophys. J. 690, 1539–1552 (2009).
Nemer, A., Goodman, J. & Wang, L. L. The role of far-ultraviolet pumping in exciting the [O I] lines in protostellar disks and winds. Astrophys. J. 904, L27 (2020).
Owen, J. E., Ercolano, B., Clarke, C. J. & Alexander, R. D. Radiation-hydrodynamic models of X-ray and EUV photoevaporating protoplanetary discs. Mon. Not. R. Astron. Soc. 401, 1415–1428 (2010).
Owen, J. E., Clarke, C. J. & Ercolano, B. On the theory of disc photoevaporation. Mon. Not. R. Astron. Soc. 422, 1880–1901 (2012).
Ercolano, B. & Owen, J. E. Blueshifted [O I] lines from protoplanetary discs: the smoking gun of X-ray photoevaporation. Mon. Not. R. Astron. Soc. 460, 3472–3478 (2016).
Wang, L. & Goodman, J. Hydrodynamic photoevaporation of protoplanetary disks with consistent thermochemistry. Astrophys. J. 847, 11 (2017).
Ercolano, B. & Owen, J. E. Theoretical spectra of photoevaporating protoplanetary discs: an atlas of atomic and low-ionization emission lines. Mon. Not. R. Astron. Soc. 406, 1553–1569 (2010).
Stelzer, B. & Schmitt, J. H. M. M. X-ray emission from a metal depleted accretion shock onto the classical T Tauri star TW Hya. Astron. Astrophys. 418, 687–697 (2004).
Kastner, J. H. et al. Evidence for accretion: high-resolution X-ray spectroscopy of the classical T Tauri star TW Hydrae. Astrophys. J. 567, 434–440 (2002).
Bai, X.-N., Ye, J., Goodman, J. & Yuan, F. Magneto-thermal disk winds from protoplanetary disks. Astrophys. J. 818, 152 (2016).
Blandford, R. D. & Payne, D. G. Hydromagnetic flows from accretion disks and the production of radio jets. Mon. Not. R. Astron. Soc. 199, 883–903 (1982).
Bai, X.-N. & Stone, J. M. Wind-driven accretion in protoplanetary disks. I. Suppression of the magnetorotational instability and launching of the magnetocentrifugal wind. Astrophys. J. 769, 76 (2013).
Wang, L. L., Bai, X. N. & Goodman, J. Global simulations of protoplanetary disk outflows with coupled non-ideal magnetohydrodynamics and consistent thermochemistry. Astrophys. J. 874, 90 (2019).
Bai, X.-N. & Stone, J.-M. Hall effect-mediated magnetic flux transport in protoplanetary disks. Astrophys. J. 836, 46 (2017).
Lavalley, C. et al. Sub-arcsecond morphology and kinematics of the DG Tauri jet in the [O I]λ6300 line. Astron. Astrophys. 327, 671–680 (1997).
Schneider, P. C. et al. HST far-ultraviolet imaging of DG Tauri. Fluorescent molecular hydrogen emission from the wide opening-angle outflow. Astron. Astrophys. 557, A110 (2013).
Schneider, P. C. et al. Discovery of a jet from the single HAe/Be star HD 100546. Astron. Astrophys. 638, L3 (2020).
Xie, C. et al. A MUSE view of the asymmetric jet from HD 163296. Astron. Astrophys. 650, L6 (2021).
Flores-Rivera, L. et al. Forbidden emission lines in protostellar outflows and jets with MUSE. Astron. Astrophys. 670, A126 (2023).
Whelan, E. T. et al. Evidence for an MHD disk wind via optical forbidden line spectroastrometry. Astrophys. J. 913, 43 (2021).
Gangi, M. et al. GIARPS high-resolution observations of T Tauri stars (GHOsT). II. Connecting atomic and molecular winds in protoplanetary disks. Astron. Astrophys. 643, A32 (2020).
Bacon, R. et al. The MUSE Hubble Ultra Deep Field Survey. I. Survey description, data reduction, and source detection. Astron. Astrophys. 608, A1 (2017).
Weilbacher, P. M. et al. The data processing pipeline for the MUSE instrument. Astron. Astrophys. 641, A28 (2020).
Fang, M., Kim, J. S., Pascucci, I. & Apai, D. An improved Hertzsprung–Russell diagram for the Orion Trapezium Cluster. Astrophys. J. 908, 49 (2021).
Herczeg, G. J. & Hillenbrand, L. A. An optical spectroscopic study of T Tauri stars. I. Photospheric properties. Astrophys. J. 786, 97 (2014).
Manara, C. F. et al. X-shooter spectroscopy of young stellar objects. II. Impact of chromospheric emission on accretion rate estimates. Astron. Astrophys. 551, A107 (2013).
Petit, P. et al. PolarBase: a database of high-resolution spectropolarimetric stellar observations. Publ. Astron. Soc. Pac. 126, 469–475 (2014).
Bertaux, J. L. et al. TAPAS, a web-based service of atmospheric transmission computation for astronomy. Astron. Astrophys. 564, A46 (2014).
Noll, S. et al. An atmospheric radiation model for Cerro Paranal. I. The optical spectral range. Astron. Astrophys. 543, A92 (2012).
Husser, T. O. et al. A new extensive library of PHOENIX stellar atmospheres and synthetic spectra. Astron. Astrophys. 553, A6 (2013).
Acknowledgements
M.F. acknowledges the science research grants from the China Manned Space Project with no. CMS-CSST-2021-B06. This research is based on observations made with ESO Telescopes under programme ID nos. 0104.C-0449, 074.A-9021(A), 078.A-9059(A), 079.A-9006(A), 079.A-9007(A), 079.A-9017(A), 081.A-9005(A), 081.A-9023(A), 089.A-9007(A), 0101.A-9012(A), 0102.A-9008(A), 60.A-9036(A), 075.C-0202(A), 081.C-0779(A), 082.C-0390(A), 082.C-0427(C), 65.I-0404(A), 082.C-0218(A), 089.C-0299(A) and 106.20Z8.011. The research made use of the Keck Observatory Archive, which is operated by the W. M. Keck Observatory and the National Aeronautics and Space Administration (NASA) Exoplanet Science Institute, under contract with NASA. We acknowledge S. E. Dahm, L. Hillenbrand, J. Carpenter, W. Borucki and G. W. Marcy for datasets that were obtained through the Keck Observatory Archive.
Author information
Authors and Affiliations
Contributions
G.J.H. started the project. J.H. reduced the data. M.F. analysed the data. L.W. contributed to the MHD modelling. A.N. and L.W. performed the line synthesis. M.F., G.J.H. and L.W. were the primary writers of the paper. All the authors contributed to the scientific interpretation of the results and reviewed the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Astronomy thanks Luca Ricci, Emma Whelan and Giulia Ballabio for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 [O I] λ6300 line profiles of TW Hya.
a, Distribution of the [O I] λ6300 line centroids of TW Hya and. These spectra are taken with ESPaDOnS (gray color filled histogram), FEROS or HARPS (maroon line filled histogram), and UVES or HIRES (olive open histogram). The numbers of the used spectra from individual instruments are shown in the figure. b, Comparisons of the [O I] λ6300 line profiles of TW Hya. For the comparisons, the centroids of the lines have been shifted to the mean centroid (-0.8 km s−1) and the spectral resolutions have been degraded to 40000. The colors are the same as in a.
Extended Data Fig. 2 [O I] λ6300 and [Ne II] 12.81 μm line profiles.
Comparisons of the observed [O I] λ6300 (a) and [Ne II] 12.81 μm (b) line profiles (cyan filled circles) with the simulated ones from the MHD wind model (golden-color filled area for disk inclinations ranging from 5∘ to 7∘). In each panel, the error bars (gray color) are 1σ uncertainties, corresponding to the 68% confidence level.
Extended Data Fig. 3 Template fits to the observed MUSE spectra near [O I] λ6300 and He Iλ6678 lines.
a, Example spectra near [O I] λ6300 (gray lines) of two exposures. In each panel, the red dashed line is the best-fit template, and the blue dash-dotted line for the best-fit template correcting for the residual. The grey color filled area marks the wavelength range within which the [O I] λ6300 line is extracted, and the grey line-filed areas mark the wavelength ranges within which the uncertainties are estimated and the fluxes are integrated to construct PSF. b, same as a, but for He I λ6678 line.
Extended Data Fig. 4 Error-weighted mean He I λ6678 intensity map.
a, Error-weighted mean He I λ6678 intensity map of TW Hya. b, Two slice cuts (white lines) for He I λ6678 line map shown in a, overplotted with the PSF (dark dashed lines) near He I λ6678 line. The dotted line in each panel is the standard deviation along the slice cut. c, He I λ6678 line residual map after a subtraction of a PSF which has been normalized to line map by the peak emission. d, the signal-to-noise ratio map for the He I λ6678 line residual map shown in c. The cyan and blue lines show contours with S/N=2 and 3, respectively.
Extended Data Fig. 5 Reduced χ2 for the toy models.
a-b, Reduced χ2 from the coarse grids of the three free parameters: α, Rin and Rout. In each panel, the reduced χ2 for each pair of two parameters is lowest one by setting the remaining parameter free. c, Reduced χ2 from the fine grid of α and Rin with Rout=61 AU. In each panel, the pentagram symbol marks the minimum reduced χ2.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Fang, M., Wang, L., Herczeg, G.J. et al. High-resolution [O I] line spectral mapping of TW Hya supportive of a magnetothermal wind. Nat Astron 7, 905–912 (2023). https://doi.org/10.1038/s41550-023-02004-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41550-023-02004-x