2019ApJ...875..159E

We present new mass estimates and cumulative mass profiles (CMPs) with Bayesian credible regions for the Milky Way (MW) Galaxy, given the kinematic data of globular clusters (GCs) as provided by (1) the Gaia DR2 collaboration and the HSTPROMO team, and (2) the new catalog in Vasiliev (2019). We use GCs beyond 15 kpc to estimate the CMP of the MW, assuming a total gravitational potential model Φ(r)=Φ_○r^-γ^, which approximates an NFW-type potential at large distances when γ = 0.5. We compare the resulting CMPs given data sets (1) and (2), and find the results to be nearly identical. The median estimate for the total mass is M₂₀₀ = 0.70 x 10¹² M_☉ and the 50% Bayesian credible interval is (0.62,0.81)×10¹²M_☉. However, because the Vasiliev catalog contains more complete data at large r, the MW total mass is slightly more constrained by these data. In this work, we also supply instructions for how to create a CMP for the MW with Bayesian credible regions, given a model for M(<r) and samples drawn from a posterior distribution. With the CMP, we can report median estimates and 50% Bayesian credible regions for the MW mass within any distance (e.g., M(r=25kpc)=0.26(0.20,0.36)×10¹²M_☉, M(r=50kpc)=0.37(0.29,0.51)×10¹²M_☉, M(r=100kpc)=0.53(0.41,0.74)×10¹²M_☉, etc.), making it easy to compare our results directly to other studies.