We present a new method for studying tracers in gravitational systems where higher moments of the line-of-sight velocity distribution are introduced via Virial equations rather than the Jeans equations. Unlike the fourth order Jeans equations, the fourth order Virial equations can simply be added to the standard second order Jeans equation without introducing a new anisotropy parameter $\beta^{\prime}$. We introduce two new global shape parameters $\zeta_A$ and $\zeta_B$ which replace the kurtosis as a more statistically robust measure of the shape of the line of sight velocity distribution. We show that in the case of stars in dwarf spheroidal galaxies these new parameters can significantly reduce the range of density profiles that are otherwise consistent with the observed stellar kinematics (a problem sometimes known as the $\beta$ degeneracy). Specifically, we find that $\zeta_A$ focuses tightly on a subset of solutions where cusped density profiles are degenerate with more concentrated cored dark matter halos. If the number density of stars $\nu(r)$ is fixed, then introducing $\zeta_B$ can further reduce the space of solutions by constraining the outer slope of the dark matter density profile. Assuming a Plummer profile for $\nu(r)$ we recover the surprising result that the dark matter in Sculptor may be cuspy after all, in contrast to the conclusions of other approaches.

The inner profile of Dark Matter (DM) halos remains one of the central problems in small-scale cosmology. At present, the problem can not be resolved in dwarf spheroidal galaxies due to a degeneracy between the DM profile and the velocity anisotropy beta of the stellar population. We discuss a method which can break the degeneracy by exploiting 3D positions and 1D line-of-sight (LOS) velocities. With the full 3D spatial information, we can determine precisely what fraction of each stars LOS motion is in the radial and tangential direction. This enables us to infer the anisotropy parameter beta directly from the data. The method is particularly effective if the galaxy is highly anisotropic. Finally, we argue that such a test could be applied to Sagittarius and potentially other dwarfs with RR Lyrae providing the necessary depth information.

The Jeans analysis is often used to infer the total density of a system by relating the velocity moments of an observable tracer population to the underlying gravitational potential. This technique has recently been applied in the search for Dark Matter in objects such as dwarf spheroidal galaxies where the presence of Dark Matter is inferred via stellar velocities. A precise account of the density is needed to constrain the expected gamma ray flux from DM self-annihilation and to distinguish between cold and warm dark matter models. Unfortunately the traditional method of fitting the second order Jeans equation to the tracer dispersion suffers from an unbreakable degeneracy of solutions due to the unknown velocity anisotropy of the projected system. To tackle this degeneracy one can appeal to higher moments of the Jeans equation. By introducing an analog to the Binney anisotropy parameter at fourth order, beta' we create a framework that encompasses all solutions to the fourth order Jeans equations rather than those in the literature that impose unnecessary correlations between anisotropy of second and fourth order moments. The condition beta' = f(beta) ensures that the degeneracy is lifted and we interpret the separable augmented density system as the order-independent case beta'= beta. For a generic choice of beta' we present the line of sight projection of the fourth moment and how it could be incorporated into a joint likelihood analysis of the dispersion and kurtosis. Having presented the mathematical framework, we then use it to develop a statistical method for the purpose of placing constraints on dark matter density parameters from discrete velocity data. The method is tested on simulated dwarf spheroidal data sets leading to results which motivate study of real dwarf spheroidal data sets.

We attempt to measure the density of dark matter in the two Dwarf Spheroidal Galaxies Fornax and Sculptor using a new method which employs Jeans equations based on both the second and fourth moment of the Collisionless Boltzmann Equation (i.e. variance and kurtosis of line of sight stellar velocities). Unlike previous related efforts, we allow the anisotropy of the radial and tangential second and fourth order moments to vary independently of each other. We apply the method to simulated data and establish that to some degree it appears to be able to break the degeneracies associated with second order only Jeans analyses. When we apply the technique to real data, we see no huge improvement in our understanding of the inner density of Fornax, which can still be fit by either a quite cuspy or cored density profile. For Sculptor however we find that the technique suggests that the data is incompatible with a steep profile and a cored profile seems more consistent. As well as presenting these new results and comparing them to other estimates in the literature, we try to understand why the technique is more effective in one system than the other.