We present and apply rigorous dynamical modeling with which we infer unprecedented constraints on the stellar and dark matter mass distribution within our Milky Way (MW), based on large sets of phase-space data on individual stars. Specifically, we model the dynamics of 16,269 G-type dwarfs from SEGUE, which sample 5 < R_GC/kpc < 12 and 0.3 < |Z|/kpc < 3. We independently fit a parameterized MW potential and a three-integral, action-based distribution function (DF) to the phase-space data of 43 separate abundance-selected sub-populations (MAPs), accounting for the complex selection effects affecting the data. We robustly measure the total surface density within 1.1 kpc of the mid-plane to 5% over 4.5 < R_GC/kpc < 9. Using metal-poor MAPs with small radial scale lengths as dynamical tracers probes 4.5 < R_GC/kpc < 7, while MAPs with longer radial scale lengths sample 7 < R_GC/kpc < 9. We measure the mass-weighted Galactic disk scale length to be R_d = 2.15+/-0.14 kpc, in agreement with the photometrically inferred spatial distribution of stellar mass. We thereby measure dynamically the mass of the Galactic stellar disk to unprecedented accuracy: M_* = 4.6+/-0.3+3.0x(R_0/kpc-8)x10^{10}Msun and a total local surface density of \Sigma_{R_0}(Z=1.1 kpc) = 68+/-4 Msun/pc^2 of which 38+/-4 Msun/pc^2 is contributed by stars and stellar remnants. By combining our surface density measurements with the terminal velocity curve, we find that the MW's disk is maximal in that V_{c,disk} / V_{c,total} = 0.83+/-0.04 at R=2.2 R_d. We also constrain for the first time the radial profile of the dark halo at such small Galactocentric radii, finding that \rho_{DM} (r;near R_0) \propto 1 / r^\alpha with \alpha < 1.53 at 95% confidence. Our results show that action-based distribution-function modeling of complex stellar data sets is now a feasible approach that will be fruitful for interpreting Gaia data.
Read full abstract