Retrievals of Protoplanetary Disk Parameters using Thermochemical Models: I. Disk Gas Mass from Hydrogen Deuteride Spectroscopy

Abstract: We discuss statistical relationships between the mass of protoplanetary disks and the hydrogen deuteride (HD) line emission and the dust spectral energy distribution (SED) determined using 3000 ProDiMo disk models. The models have 15 free parameters describing disk physical properties, the central star, and the local radiation field. The sampling of physical parameters is done using a Monte Carlo approach to evaluate the probability density functions of observables as a function of physical parameters. We find that the HD fractional abundance is almost constant even though the UV flux varies by several orders of magnitude. Probing the statistical relation between the physical quantities and the HD flux, we find that low-mass (optically thin) disks display a tight correlation between the average disk gas temperature and HD line flux, while massive disks show no such correlation. We demonstrate that the central star luminosity, disk size, dust size distribution, and HD flux may be used to determine the disk gas mass to within a factor of three. We also find that the far-IR and sub-mm/mm SEDs and the HD flux may serve as strong constraints for determining the disk gas mass to within a factor of two. If the HD lines are fully spectrally resolved ($R\gtrsim 1.5\times10^6, \Delta v=0.2~\rm km\,s^{-1}$), the 56 $\mu$m and 112 $\mu$m HD line profiles alone may constrain the disk gas mass to within a factor of two.