Hall-magnetohydrodynamic simulations of X-ray photoevaporative protoplanetary disc winds

Abstract

ABSTRACT Understanding the evolution and dispersal via energetic stellar radiation of protoplanetary discs (PPDs) is a prominent challenge in astrophysics. It has been established that X-ray luminosity from the central protostar can significantly heat the surface of the disc, causing powerful photoevaporative winds that eject a considerable fraction of the disc’s mass. Recent work in the field has moreover shown the importance of global PPD simulations that simultaneously take into account non-ideal magnetohydrodynamic (MHD) effects and detailed thermochemistry. In this paper, we combine these two aspects and figure out how they interact. Focus is put on the Hall Effect (HE) and the impact it has on the overall field topology and mass-loss/accretion rates. Utilizing a novel X-ray temperature parametrization, we perform 2D-axisymmetric MHD simulations with the nirvana-iii fluid code, covering all non-ideal effects. We find that, in the aligned orientation, the HE causes prominent inward displacement of the poloidal field lines that increase the accretion rate through a laminar Maxwell stress. We find that outflows are mainly driven by photoevaporation – unless the magnetic field strength is considerable (i.e. βp ≤ 103) or the X-ray luminosity low enough (i.e. log LX ≤ 29.3). Inferred mass-loss rate are in the range of the expected values 10−8–$10^{-7}\, M_\odot \, {\rm yr}^{-1}$. Moreover, we performed pure hydrodynamic (HD) runs and compared them with the equivalent MHD runs. We concluded that the magnetic field does indeed contribute to the mass-loss rate, albeit only discernibly so for low enough LX (i.e. log LX ≤ 30.8).