Galaxy structure from multiple tracers - II. M87 from parsec to megaparsec scales

Following a number of conflicting studies of M87's mass profile, we undertake a dynamical analysis of multiple tracer populations to constrain its mass over a large radius range. We combine stellar kinematics in the central regions with the dynamics of 612 globular clusters out to 200 kpc and satellite galaxies extending to scales comparable with the virial radius. Using a spherical Jeans analysis, we are able to disentangle the mass contributions from the dark and baryonic components and set constraints on the structure of each. Assuming isotropy, we explore four different models for the dark matter halo and find that a centrally cored dark matter distribution is preferred. We infer a stellar mass-to-light ratio Upsilon(star,nu) = 6.9 +/- 0.1 - consistent with a Salpeter-like initial mass function (IMF) - and a core radius r(c) = 67 +/- 20 kpc. We then introduce anisotropy and find that, while the halo remains clearly cored, the radial stellar anisotropy has a strong impact on both Upsilon(star,nu) and the core's radius; here we find Upsilon(star,nu) = 3.50(-0.36)(+0.32) - consistent with a Chabrier-like IMF - and r(c) = 19.00(-8.34)(+8.38) kpc. Thus, the presence of a core at the centre of the dark halo is robust against anisotropy assumptions, while the stellar mass and core size are not. We are able to reconcile previously discrepant studies by showing that modelling the globular cluster data alone leads to the very different inference of a super-NFW cusp, thus highlighting the value of multiple-population modelling, and we point to the possible role of M87's AGN and the cluster environment in forming the central dark matter core.