Modeling protoplanetary disk SEDs with artificial neural networks. Revisiting the viscous disk model and updated disk masses

We model the spectral energy distributions (SEDs) of 23 protoplanetary disks in the Taurus-Auriga star-forming region using detailed disk models and a Bayesian approach. This is made possible by combining these models with artificial neural networks to drastically speed up their performance. Such a setup allows us to confront $\alpha$-disk models with observations while accounting for several uncertainties and degeneracies. Our results yield high viscosities and accretion rates for many sources, which is not consistent with recent measurements of low turbulence levels in disks. This inconsistency could imply that viscosity is not the main mechanism for angular momentum transport in disks, and that alternatives such as disk winds play an important role in this process. We also find that our SED-derived disk masses are systematically higher than those obtained solely from (sub)mm fluxes, suggesting that part of the disk emission could still be optically thick at (sub)mm wavelengths. This effect is particularly relevant for disk population studies and alleviates previous observational tensions between the masses of protoplanetary disks and exoplanetary systems.