Effect of MHD wind-driven disk evolution on the observed sizes of protoplanetary disks

TL;DR

This study investigates how magnetic wind-driven disk evolution influences the observed sizes of protoplanetary disks, using models to compare with observed data across different star-forming regions.

Contribution

It introduces a combined thermochemical and analytical model to analyze the impact of magnetic disk winds on disk size evolution, highlighting the importance of the characteristic radius R_c.

Findings

Disk size decreases linearly with time under constant wind strength.

Larger R_c leads to a steeper decrease in disk size.

Models match observed sizes in older star-forming regions but overpredict sizes in very young disks.

Abstract

It is still unclear whether the evolution of protoplanetary disks, a key ingredient in the theory of planet formation, is driven by viscous turbulence or magnetic disk winds. As viscously evolving disks expand outward over time, the evolution of disk sizes is a discriminant test for studying disk evolution. However, it is unclear how the observed disk size changes over time if disk evolution is driven by magnetic disk winds. Combining the thermochemical code DALI with the analytical wind-driven disk evolution model presented in Tabone et al. (2021a), we study the time evolution of the observed gas outer radius as measured from CO rotational emission ($R_{\rm CO, 90\%}$). The evolution of $R_{\rm CO, 90\%}$ is driven by the evolution of the disk mass, as the physical radius stays constant over time. For a constant $\alpha_{\rm DW}$, an extension of the $\alpha-$Shakura-Sunyaev parameter to wind-driven accretion, $R_{\rm CO, 90\%}$ decreases linearly with time. Its initial size is set by the disk mass and the characteristic radius $R_c$, but only $R_c$ affects the evolution of $R_{\rm CO, 90\%}$, with a larger $R_c$ resulting in a steeper decrease of $R_{\rm CO, 90\%}$. For a time-dependent $\alpha_{\rm DW}$ $R_{\rm CO, 90\%}$ stays approximately constant during most of the disk lifetime until $R_{\rm CO, 90\%}$ rapidly shrinks as the disk dissipates. The constant $\alpha_{\rm DW}$-models are able to reproduce the observed gas disk sizes in the $\sim1-3$ Lupus and $\sim5-11$ Myr old Upper Sco star-forming regions. However, they likely overpredict the gas disk size of younger $(\lessapprox0.7\ \mathrm{Myr})$ disks.