Modeling the Evolution of Protoplanetary Disks: Two Pathways from Gravitational Instability to MHD Wind-Driven Accretion
arXiv:2607.12007
The paper introduces a semi‑two‑dimensional model of protoplanetary disk evolution that includes gravitational instability, MHD winds, MRI, and radiative effects, and shows that the radial distribution of large‑scale magnetic flux leads to two distinct pathways: a compact, shadowed, massive disk or an extended, flared disk.
Abstract
The global evolution of protoplanetary disks sets the initial conditions for planet formation. However, most models focus on individual evolutionary phases, with idealized initial conditions and oversimplified prescriptions for angular momentum transport and thermodynamics. We present a more realistic semi-two-dimensional (D) model incorporating gravitational instability (GI), magnetohydrodynamic (MHD) winds, magneto-rotational instability (MRI), stellar irradiation, self-shadowing, and radiation transport. The radial distribution of large-scale magnetic flux drives two different pathways of disk evolution. When the vertical field is spatially uniform, a puffed-up, MRI-heated inner rim shadows the disk beyond it, sustaining a massive, gravitationally unstable region for Myr and, for several Myr, a compact ( AU), cold ( K), low-turbulence (), high-density (), optically thick reservoir, so that the disk mass inferred from mm-continuum emission can be greatly underestimated. When the field instead scales with midplane gas pressure, it drives stronger transport in the inner disk and eventually strips the shadow, leaving an extended, flared disk whose observable mass closely traces the true mass. Our results connect GI-dominated Class~0/I disks to MHD wind-driven Class~II disks, and point to three broader conclusions: (i) disk physics is strongly inhomogeneous in space and time, so constant- treatments miss essential physics; (ii) thermodynamics plays an active role, with self-shadowing simultaneously preserving GI and weakening MHD winds; and (iii) the distribution of large-scale magnetic flux is the key uncertainty, closely linked to whether the shadow is maintained. The two pathways align, respectively, with observations of compact, shadowed disks and extended, irradiated disks.
22 pages, 5 figures; submitted to ApJ