The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase-space distribution of dynamical tracers. We test this approach using realistic mock stellar halos constructed from the Aquarius N-body simula
tions of dark matter halos in the $Lambda$CDM cosmology. We extend the standard treatment to include a Navarro-Frenk-White (NFW) potential and use a maximum likelihood method to recover the parameters describing the simulated halos from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is highly correlated with the estimate of halo concentration. The best-fit halo masses within the virial radius, $R_{200}$, are biased, ranging from a 40% underestimate to a 5% overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass within the half-mass radius of tracers even when restricted to tracers within 60kpc. The recovered velocity anisotropies of tracers, $beta$, are biased systematically, but this does not affect other parameters if tangential velocity data are used as constraints.
We investigate the angular and kinematic distributions of satellite galaxies around a large sample of bright isolated primaries in the spectroscopic and photometric catalogues of the Sloan Digital Sky Survey (SDSS). We detect significant anisotropy i
n the spatial distribution of satellites. To test whether this anisotropy could be related to the rotating disks of satellites recently found by Ibata et al. in a sample of SDSS galaxies, we repeat and extend their analysis. Ibata et al. found an excess of satellites on opposite sides of their primaries having anticorrelated radial velocities. We find that this excess is sensitive to small changes in the sample selection criteria which can greatly reduce its significance. In addition, we find no evidence for correspondingly correlated velocities for satellites observed on the same side of their primaries, which would be expected for rotating disks of satellites. We conclude that the detection of coherent rotation in the satellite population in current observational samples is not robust. We compare our data to the $Lambda$CDM Millennium simulations populated with galaxies according to the semi-analytic model of Guo et al. We find excellent agreement with the spatial distribution of satellites in the SDSS data and the lack of a strong signal from coherent rotation.
We present a new technique for creating mock catalogues of the individual stars that make up the accreted component of stellar haloes in cosmological simulations and show how the catalogues can be used to test and interpret observational data. The ca
talogues are constructed from a combination of methods. A semi-analytic galaxy formation model is used to calculate the star formation history in haloes in an N-body simulation and dark matter particles are tagged with this stellar mass. The tags are converted into individual stars using a stellar population synthesis model to obtain the number density and evolutionary stage of the stars, together with a phase-space sampling method that distributes the stars while ensuring that the phase-space structure of the original N-body simulation is maintained. A set of catalogues based on the $Lambda$CDM Aquarius simulations of Milky Way mass haloes have been created and made publicly available on a website. Two example applications are discussed that demonstrate the power and flexibility of the mock catalogues. We show how the rich stellar substructure that survives in the stellar halo precludes a simple measurement of its density profile and demonstrate explicitly how pencil-beam surveys can return almost any value for the slope of the profile. We also show that localized variations in the abundance of particular types of stars, a signature of differences in the composition of stellar populations, allow streams to be easily identified.
We use the SDSS/DR8 galaxy sample to study the radial distribution of satellite galaxies around isolated primaries, comparing to semi-analytic models of galaxy formation based on the Millennium and Millennium-II simulations. SDSS satellites behave di
fferently around high- and low-mass primaries: those orbiting objects with $M_*>10^{11}M_odot$ are mostly red and are less concentrated towards their host than the inferred dark matter halo, an effect that is very pronounced for the few blue satellites. On the other hand, less massive primaries have steeper satellite profiles that agree quite well with the expected dark matter distribution and are dominated by blue satellites, even in the inner regions where strong environmental effects are expected. In fact, such effects appear to be strong only for primaries with $M_* > 10^{11}M_odot$. This behaviour is not reproduced by current semi-analytic simulations, where satellite profiles always parallel those of the dark matter and satellite populations are predominantly red for primaries of all masses. The disagreement with SDSS suggests that environmental effects are too efficient in the models. Modifying the treatment of environmental and star formation processes can substantially increase the fraction of blue satellites, but their radial distribution remains significantly shallower than observed. It seems that most satellites of low-mass primaries can continue to form stars even after orbiting within their joint halo for 5 Gyr or more.
