Dynamical model of the Milky Way using APOGEE and Gaia data

TL;DR

This paper develops a detailed dynamical model of the Milky Way using Gaia and APOGEE data, providing updated estimates of dark matter distribution, circular velocity, and mass density across a large radial range.

Contribution

It extends previous models by incorporating larger spatial data, improving estimates of the Milky Way's dark matter halo and disk properties, and analyzing effects of non-axisymmetric features and disk flaring.

Findings

Dark matter density slope: -1.602

Circular velocity at Sun: 234.7 km/s

Total surface density: 55.5 M_sun/pc^2

Abstract

We construct a dynamical model of the Milky Way disk from a data set, which combines Gaia EDR3 and APOGEE data throughout Galactocentric radii between $5.0\leq R\leq19.5$ kpc. We make use of the spherically-aligned Jeans Anisotropic Method to model the stellar velocities and their velocity dispersions. Building upon our previous work, our model now is fitted to kinematic maps that have been extended to larger Galactocentric radii due to the expansion of our data set, probing the outer regions of the Galactic disk. Our best-fitting dynamical model suggests a logarithmic density slope of $\alpha_{\rm DM}=-1.602\pm0.079_{\rm syst}$ for the dark matter halo and a dark matter density of $\rho_{\rm DM}(R_{\odot})=(8.92\pm0.56_{\rm syst})\times 10^{-3}$ M$_{\odot}$ pc$^{-3}$ ($0.339\pm0.022_{\rm syst}$ GeV cm$^{3}$). We estimate a circular velocity at the solar radius of $v_{\rm circ}=(234.7\pm1.7_{\rm syst})$ km s$^{-1}$ with a decline towards larger radii. The total mass density is $\rho_{\rm tot}(R_{\odot})$=$(0.0672\pm0.0015_{\rm syst})$ M$_{\odot}$ pc$^{-3}$ with a slope of $\alpha_{\rm tot}$=$-2.367\pm0.047_{\rm syst}$ for $5\leq R\leq19.5$ kpc and the total surface density is $\Sigma(R{_\odot}, |z|\leq$ 1.1 kpc)=$(55.5\pm1.7_{\rm syst})$ M$_{\odot}$ pc$^{-2}$. While the statistical errors are small, the error budget of the derived quantities is dominated by the 3 to 7 times larger systematic uncertainties. These values are consistent with our previous determination, but systematic uncertainties are reduced due to the extended data set covering a larger spatial extent of the Milky Way disk. Furthermore, we test the influence of non-axisymmetric features on our resulting model and analyze how a flaring disk model would change our findings.