If "seed" central black holes were common in the subgalactic building blocks that merged to form present-day massive galaxies, then relic intermediate-mass black holes (IMBHs) should be present in the Galactic bulge and halo. We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with black holes, and assess the size, properties, and detectability of the leftover population. The method assigns a black hole to the most tightly bound central particle of each subhalo at infall according to an extrapolation of the M_BH-sigma_* relation, and self-consistently follows the accretion and disruption of Milky Way progenitor dwarfs and their holes in a cosmological "live" host from high redshift to today. We show that, depending on the minimum stellar velocity dispersion, sigma_m, below which central black holes are assumed to be increasingly rare, as many as ~2000 (sigma_m=3 km/s) or as few as ~70 (sigma_m=12 km/s) IMBHs may be left wandering in the halo of the Milky Way today. The fraction of IMBHs kicked out of their host by gravitational recoil is < 20%. We identify two main Galactic subpopulations, "naked" IMBHs, whose host subhalos were totally destroyed after infall, and "clothed" IMBHs residing in dark matter satellites that survived tidal stripping. Naked IMBHs typically constitute 40-50% of the total and are more centrally concentrated. We show that, in the sigma_m=12 km/s scenario, the clusters of tightly bound stars that should accompany naked IMBHs would be fainter than m_V=16 mag, spatially resolvable, and have proper motions of 0.1-10 milliarcsec per year. Their detection may provide an observational tool to constrain the formation history of massive black holes in the early Universe.

We discuss a novel approach to "weighing" the Milky Way dark matter halo, one that combines the latest samples of halo stars selected from the Sloan Digital Sky Survey (SDSS) with state-of-the-art numerical simulations of Milky Way analogs. The fully cosmological runs employed in the present study include "Eris", one of the highest-resolution hydrodynamical simulations of the formation of a M_vir=8e11 M_sun late-type spiral, and the dark-matter only M_vir=1.7e12 M_sun "Via Lactea II" simulation. Eris provides an excellent laboratory for creating mock SDSS samples of tracer halo stars, and we successfully compare their density, velocity anisotropy, and radial velocity dispersion profiles with the observational data. Most mock SDSS realizations show the same "cold veil" recently observed in the distant stellar halo of the Milky Way, with tracers as cold as sigma_los ~ 50 km/s between 100-150 kpc. Controlled experiments based on the integration of the spherical Jeans equation as well as a particle tagging technique applied to Via Lactea II show that a "heavy" M_vir 2e12 M_sun realistic host produces a poor fit to the kinematic SDSS data. We argue that these results offer added evidence for a "light," centrally-concentrated Milky Way halo.

A galaxy halo may contain a large number of intermediate mass black holes (IMBHs) with masses in the range of 10^{2-6} solar mass. We propose to directly detect these IMBHs by observing multiply imaged QSO-galaxy or galaxy-galaxy strong lens systems in the submillimeter bands with high angular resolution. The silhouette of an IMBH in the lensing galaxy halo would appear as either a monopole-like or a dipole-like variation at the scale of the Einstein radius against the Einstein ring of the dust-emitting region surrounding the QSO. We use a particle tagging technique to dynamically populate a Milky Way-sized dark matter halo with black holes, and show that the surface mass density and number density of IMBHs have power-law dependences on the distance from the center of the host halo if smoothed on a scale of ~ 1 kpc. Most of the black holes orbiting close to the center are freely roaming as they have lost their dark matter hosts during infall due to tidal stripping. Next generation submillimeter telescopes with high angular resolution (< 0.3 mas) will be capable of directly mapping such off-nuclear freely roaming IMBHs with a mass of ~ 10^6 solar mass in a lensing galaxy that harbours a O(10^9) solar mass supermassive black hole in its nucleus.

We find that the relative contribution of satellite galaxies accreted at high redshift to the stellar population of the Milky Way's smooth halo increases with distance, becoming observable relative to the classical smooth halo about 15 kpc from the Galactic center. In particular, we determine line-of-sight-averaged [Fe/H] and [alpha/Fe] in the metal-poor main-sequence turnoff (MPMSTO) population along every Sloan Extension for Galactic Understanding and Exploration (SEGUE) spectroscopic line of sight. Restricting our sample to those lines of sight along which we do not detect elements of cold halo substructure (ECHOS), we compile the largest spectroscopic sample of stars in the smooth component of the halo ever observed in situ beyond 10 kpc. We find significant spatial autocorrelation in [Fe/H] in the MPMSTO population in the distant half of our sample beyond about 15 kpc from the Galactic center. Inside of 15 kpc however, we find no significant spatial autocorrelation in [Fe/H]. At the same time, we perform SEGUE-like observations of N-body simulations of Milky Way analog formation. While we find that halos formed entirely by accreted satellite galaxies provide a poor match to our observations of the halo within 15 kpc of the Galactic center, we do observe spatial autocorrelation in [Fe/H] in the simulations at larger distances. This observation is an example of statistical chemical tagging and indicates that spatial autocorrelation in metallicity is a generic feature of stellar halos formed from accreted satellite galaxies.

We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with stars. The method is calibrated using the observed luminosity function of Milky Way satellites and the concentration of their stellar populations, and self-consistently follows the accretion and disruption of progenitor dwarfs and the build-up of the stellar halo in a cosmological "live host". Simple prescriptions for assigning stellar populations to collisionless particles are able to reproduce many properties of the observed Milky Way halo and its surviving dwarf satellites, like velocity dispersions, sizes, brightness profiles, metallicities, and spatial distribution. Our model predicts the existence of approximately 1,850 subhalos harboring "extremely faint" satellites (with mass-to-light ratios >5,000) lying beyond the Sloan Digital Sky Survey detection threshold. Of these, about 20 are "first galaxies", i.e. satellites that formed a stellar mass above 10 Msun before redshift 9. The ten most luminous satellites (L> 1e6 Lsun) in the simulation are hosted by subhalos with peak circular velocities today in the range V_max=10-40 km/s that have shed between 80% and 99% of their dark mass after being accreted at redshifts 1.7< z <4.6. The satellite maximum circular velocity and stellar line-of-sight velocity dispersion today follow the relation V_max=2.2 sigma_los. We apply a standard mass estimation algorithm based on Jeans modelling of the line-of-sight velocity dispersion profiles to the simulated dwarf spheroidals, and test the accuracy of this technique. The inner (within 300 pc) mass-luminosity relation for currently detectable satellites is nearly flat in our model, in qualitative agreement with the "common mass scale" found in Milky Way dwarfs. We do, however, predict a weak, but significant positive correlation for these objects: M_300 ~L^{0.088 \pm 0.024}.

We present an updated version of MegaZ-LRG (Collister et al.,(2007)) with photometric redshifts derived with the neural network method, ANNz as well as five other publicly available photo-z codes (HyperZ, SDSS, Le PHARE, BPZ and ZEBRA) for ~1.5 million Luminous Red Galaxies (LRGs) in SDSS DR6. This allows us to identify how reliable codes are relative to each other if used as described in their public release. We compare and contrast the relative merits of each code using ~13000 spectroscopic redshifts from the 2SLAQ sample. We find that the performance of each code depends on the figure of merit used to assess it. As expected, the availability of a complete training set means that the training method performs best in the intermediate redshift bins where there are plenty of training objects. Codes such as Le PHARE, which use new observed templates perform best in the lower redshift bins. All codes produce reasonable photometric redshifts, the 1-sigma scatters ranging from 0.057 to 0.097 if averaged over the entire redshift range. We also perform tests to check whether a training set from a small region of the sky such as 2SLAQ produces biases if used to train over a larger area of the sky. We conclude that this is not likely to be a problem for future wide-field surveys. The complete photometric redshift catalogue including redshift estimates and errors on these from all six methods can be found at www.star.ucl.ac.uk/~mbanerji/MegaZLRGDR6/megaz.html