The two-dimensional pressure structure of the HD 163296 protoplanetary disk as probed by multi-molecule kinematics

Gas kinematics is a new and unique way to study planet-forming environments by an accurate characterization of disk velocity fields. High angular resolution ALMA observations allow deep kinematical analysis of disks, by observing molecular line emission at high spectral resolution. In particular, rotation curves are key tools for studying the disk pressure structure and efficiently estimating fundamental disk parameters, such as mass and radius. In this work we explore the potential of a multi-molecule approach to gas kinematics to provide a 2D characterization of the HD 163296 disk. From the high quality data of the MAPS Large Program we extracted the rotation curves of rotational lines from seven distinct molecular species ((CO)-C-12, (CO)-C-13, (CO)-O-18, HCN, H2CO, HCO+, C2H), spanning a wide range in the disk radial and vertical extents. To obtain reliable rotation curves for the HCN and C2H hyperfine lines, we extended standard methodologies to fit multi-component line profiles. We then sampled the likelihood of a thermally stratified model that reproduces all the rotation curves simultaneously, taking into account the molecular emitting layers z(R) and disk thermal structure T(R, z). From this exploration, we obtained dynamical estimates of three fundamental parameters: the stellar mass M-star = 1.89 M-circle dot, the disk mass M-d = 0.12 M-circle dot, and the scale radius R-c = 143 au. We also explore how rotation curves, and consequently the parameter estimates, depend on the adopted emitting layers: the disk mass proves to be the most affected by these systematics, yet the main trends we find do not depend on the adopted parameterization. Finally, we investigated the impact of thermal structure on gas kinematics, and show that the thermal stratification can efficiently explain the measured rotation velocity discrepancies between tracers at different heights. Our results show that such a multi-molecule approach, tracing a large range of emission layers, can provide unique constraints on the (R, z) pressure structure of protoplanetary disks.