The population of Milky Way satellites in the Λ cold dark matter cosmology

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic haloes of the Aquarius project. The model, calculated using the galform code, incorporates improved treatments of the relevant physics in the ? cold dark matter cosmogony, particularly a self-consistent calculation of reionization by ultraviolet (UV) photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhaloes is suppressed primarily by supernova feedback, while for smaller subhaloes, it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present-day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order to completely reionize the universe by redshift z? 8. In the most successful model, the local sources photoionize the pre-galactic region completely by z? 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminositymetallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of similar to 2 among the six simulations.