A Method for Measuring (Slopes of) the Mass Profiles of Dwarf Spheroidal Galaxies

We introduce a method for measuring the slopes of mass profiles within dwarf spheroidal (dSph) galaxies directly from stellar spectroscopic data and without adopting a dark matter halo model. Our method combines two recent results: 1) spherically symmetric, equilibrium Jeans models imply that the product of halflight radius and (squared) stellar velocity dispersion provides an estimate of the mass enclosed within the halflight radius of a dSph stellar component, and 2) some dSphs have chemo-dynamically distinct stellar subcomponents that independently trace the same gravitational potential. We devise a statistical method that uses measurements of stellar positions, velocities and spectral indices to distinguish two dSph stellar subcomponents and to estimate their individual halflight radii and velocity dispersions. For a dSph with two detected stellar subcomponents, we obtain estimates of masses enclosed at two discrete points in the same mass profile, immediately defining a slope. Applied to published spectroscopic data, our method distinguishes stellar subcomponents in the Fornax and Sculptor dSphs, for which we measure slopes $Γ\equiv Δ\log M / Δ\log r=2.61_{-0.37}^{+0.43}$ and $Γ=2.95_{-0.39}^{+0.51}$, respectively. These values are consistent with 'cores' of constant density within the central few-hundred parsecs of each galaxy and rule out `cuspy' Navarro-Frenk-White (NFW) profiles ($d\log M/d\log r \leq 2$ at all radii) with significance \ga 96% and \ga 99%, respectively. Tests with synthetic data indicate that our method tends systematically to overestimate the mass of the inner stellar subcomponent to a greater degree than that of the outer stellar subcomponent, and therefore to underestimate the slope $Γ$ (implying that the stated NFW exclusion levels are conservative).