Disk Masses in the Embedded and T Tauri Phases of Stellar Evolution

Motivated by a growing concern that masses of circumstellar disks may have been systematically underestimated by conventional observational methods, we present a numerical hydrodynamics study of time-averaged disk masses (langM _drang) around low-mass Class 0, Class I, and Class II objects. Mean disk masses (\overline{M}_d) are then calculated by weighting the time-averaged disk masses according to the corresponding stellar masses using a power-law weight function with a slope typical for the Kroupa initial mass function of stars. Two distinct types of disks are considered: self-gravitating disks, in which mass and angular momentum are redistributed exclusively by gravitational torques, and viscous disks, in which both the gravitational and viscous torques are at work. We find that self-gravitating disks have mean masses that are slowly increasing along the sequence of stellar evolution phases. More specifically, Class 0/I/II self-gravitating disks have mean masses \overline{M}_d=0.09, 0.10, and 0.12 M _sun, respectively. Viscous disks have similar mean masses (\overline{M}_d=0.10-0.11 M _sun) in the Class 0/I phases but almost a factor of 2 lower mean mass in the Class II phase (\overline{M}_d,CII=0.06 M_⊙). In each evolution phase, time-averaged disk masses show a large scatter around the mean value. Our obtained mean disk masses are larger than those recently derived by Andrews & Williams and Brown et al., regardless of the physical mechanisms of mass transport in the disk. The difference is especially large for Class II disks, for which we find \overline{M}_d,CII=0.06-0.12 M _sun but Andrews & Williams report median masses of the order 3 × 10^-3 M _sun. When Class 0/I/II systems are considered altogether, a least-squares best fit yields the following relation between the time-averaged disk and stellar masses, langM _drang = (0.2 ± 0.05)langM _*rang^1.3±0.15. The dependence of langM _drang on langM _*rang becomes progressively steeper along the sequence of stellar evolution phases, with exponents 0.7 ± 0.2, 1.3 ± 0.15, and 2.2 ± 0.2 for Class 0, Class I, and Class II systems, respectively.