Using a newly-developed modeling technique, we present orbit-based dynamical models of the Carina, Draco, Fornax, Sculptor, and Sextans dwarf spheroidal (dSph) galaxies. These models calculate the dark matter profiles non-parametrically without requiring any assumptions to be made about their profile shapes. By lifting this restriction we discover a host of dark matter profiles in the dSphs that are different from the typical profiles suggested by both theorists and observers. However, when we scale these profiles appropriately and plot them on a common axis they appear to follow an approximate r^-1 power law with considerable scatter.

We introduce a novel implementation of orbit-based (or Schwarzschild) modeling that allows dark matter density profiles to be calculated non-parametrically in nearby galaxies. Our models require no assumptions to be made about velocity anisotropy or the dark matter profile. The technique can be applied to any dispersion-supported stellar system, and we demonstrate its use by studying the Local Group dwarf spheroidal (dSph) galaxy Draco. We use existing kinematic data at larger radii and also present 12 new radial velocities within the central 13 pc obtained with the VIRUS-W integral field spectrograph on the 2.7m telescope at McDonald Observatory. Our non-parametric Schwarzschild models find strong evidence that the dark matter profile in Draco is cuspy for 20 < r < 700 pc. The profile for r > 20 pc is well-fit by a power law with slope \alpha=-1.0 +/- 0.2, consistent with predictions from Cold Dark Matter (CDM) simulations. Our models confirm that, despite its low baryon content relative to other dSphs, Draco lives in a massive halo.

We construct axisymmetric Schwarzschild models to measure the mass profile of the local group dwarf galaxy Fornax. These models require no assumptions to be made about the orbital anisotropy of the stars, as is the case for commonly used Jeans models. We test a variety of parameterizations of dark matter density profiles and find cored models with uniform density rho_c = (1.6 +/- 0.1) x 10^-2 M_sun pc^-3 fit significantly better than the cuspy halos predicted by cold dark matter simulations. We also construct models with an intermediate-mass black hole, but are unable to make a detection. We place a 1-sigma upper limit on the mass of a potential intermediate-mass black hole at M_BH < 3.2 x 10^4 M_sun.

We present axisymmetric, orbit-based models to study the central black hole, stellar mass-to-light ratio, and dark matter halo of NGC 4594 (M104, the Sombrero Galaxy). For stellar kinematics, we use published high-resolution kinematics of the central region taken with the Hubble Space Telescope, newly obtained Gemini long-slit spectra of the major axis, and integral field kinematics from the SAURON instrument. At large radii, we use globular cluster kinematics to trace the mass profile and apply extra leverage to recovering the dark matter halo parameters. We find a black hole of mass M_{\bullet}=(6.6 +/- 0.4) x 10^8 M_{\odot}, and determine the stellar M/L_I=3.4 +/- 0.05 (uncertainties are the 68% confidence band marginalized over the other parameters). Our best fit dark matter halo is a cored logarithmic model with asymptotic circular speed V_c=376 +/- 12 km/s and core radius r_c= 4.7 +/- 0.6 kpc. The fraction of dark to total mass contained within the half-light radius is 0.52. Taking the bulge and disk components into account in our calculation of \sigma_e puts NGC 4594 squarely on the M-\sigma relation. We also determine that NGC 4594 lies directly on the M-L relation.

We test the proposal by El-Zant et al that the dark matter density of halos could be reduced through dynamical friction acting on heavy baryonic clumps in the early stages of galaxy formation. Using N-body simulations, we confirm that the inner halo density cusp is flattened to 0.2 of the halo break radius by the settling of a single clump of mass \ga 0.5% of the halo mass. We also find that an ensemble of 50 clumps each having masses \ga 0.2% can flatten the cusp to almost the halo break radius on a time scale of \sim9 Gyr, for an NFW halo of concentration 15. We summarize some of the difficulties that need to be overcome if this mechanism is to resolve the apparent conflict between the observed inner densities of galaxy halos and the predictions of LCDM.