Global Simulations of Protoplanetary Disk Outflows with Coupled Non-ideal Magnetohydrodynamics and Consistent Thermochemistry

• Published in: The Astrophysical journal Vol. 874; no. 1; pp. 90 - 110
• Main Authors: Wang, Lile, Bai, Xue-Ning, Goodman, Jeremy
• Format: Journal Article
• Language: English
• Published: Philadelphia The American Astronomical Society 20.03.2019; IOP Publishing
• Subjects: Accretion disks; accretion, accretion disks; Astrochemistry; Astrophysics; circumstellar matter; Computational fluid dynamics; Computer simulation; Deposition; Fluid flow; Flux density; Ionization; Irradiation; Magnetic fields; Magnetohydrodynamic simulation; Magnetohydrodynamic turbulence; Magnetohydrodynamics; magnetohydrodynamics (MHD); methods: numerical; Outflow; Planet formation; planets and satellites: formation; Protoplanetary disks; Protoplanets; Radiation; Radiative transfer; Stellar winds; Thermochemistry; Wind
• ISSN: 0004-637X; 1538-4357
• DOI: 10.3847/1538-4357/ab06fd

Magnetized winds may be important in dispersing protoplanetary disks and influencing planet formation. We carry out global magnetohydrodynamic simulations in axisymmetry, coupled with ray-tracing radiative transfer, consistent thermochemistry, and non-ideal MHD diffusivities. Magnetized models lacking EUV photons ( ) feature warm molecular outflows that have typical poloidal speeds . When the magnetization is sufficient to drive accretion rates , the wind mass-loss rate is comparable. Such outflows are driven not centrifugally but by the pressure of toroidal magnetic fields produced by bending the poloidal field. Both the accretion and outflow rates increase with the poloidal field energy density, the former almost linearly. The mass-loss rate is also strongly affected by ionization due to UV and X-ray radiation near the wind base. Adding EUV irradiation to the system heats, ionizes, and accelerates the part of the outflow nearest the symmetry axis, but reduces the overall mass-loss rate by exerting pressure on the wind base. Most of our models are non-turbulent, but some with reduced dust abundance and therefore higher ionization fractions exhibit magnetorotational instabilities near the base of the wind.