We simulate the growth of galaxies and their central supermassive black holes by implementing a suite of semi‐analytic models on the output of the Millennium Run, a very large simulation of the concordance Λ cold dark matter cosmogony. Our procedures follow the detailed assembly history of each object and are able to track the evolution of all galaxies more massive than the Small Magellanic Cloud throughout a volume comparable to that of large modern redshift surveys. In this first paper we supplement previous treatments of the growth and activity of central black holes with a new model for ‘radio’ feedback from those active galactic nuclei that lie at the centre of a quasi‐static X‐ray‐emitting atmosphere in a galaxy group or cluster. We show that for energetically and observationally plausible parameters such a model can simultaneously explain: (i) the low observed mass drop‐out rate in cooling flows; (ii) the exponential cut‐off at the bright end of the galaxy luminosity function; and (iii) the fact that the most massive galaxies tend to be bulge‐dominated systems in clusters and to contain systematically older stars than lower mass galaxies. This success occurs because static hot atmospheres form only in the most massive structures, and radio feedback (in contrast, for example, to supernova or starburst feedback) can suppress further cooling and star formation without itself requiring star formation. We discuss possible physical models that might explain the accretion rate scalings required for our phenomenological ‘radio mode’ model to be successful.

We use semi-analytic techniques to study the formation and evolution of brightest cluster galaxies (BCGs). We show the extreme hierarchical nature of these objects and discuss the limits of simple ways to capture their evolution. In a model where cooling flows are suppressed at late times by AGN activity, the stars of BCGs are formed very early (50 per cent at z~5, 80 per cent at z~3) and in many small galaxies. The high star formation rates in these high-z progenitors are fuelled by rapid cooling, not by merger-triggered starbursts. We find that model BCGs assemble surprisingly late: half their final mass is typically locked-up in a single galaxy after z~0.5. Because most of the galaxies accreted onto BCGs have little gas content and red colours, late mergers do not change the apparent age of BCGs. It is this accumulation of a large number of old stellar populations -- driven mainly by the merging history of the dark matter halo itself -- that yields the observed homogeneity of BCG properties. In the second part of the paper, we discuss the evolution of BCGs to high redshifts, from both observational and theoretical viewpoints. We show that our model BCGs are in qualitative agreement with high-z observations. We discuss the hierarchical link between high-z BCGs and their local counter-parts. We show that high-z BCGs belong to the same population as the massive end of local BCG progenitors, although they are not in general the same galaxies. Similarly, high-z BCGs end-up as massive galaxies in the local Universe, although only a fraction of them are actually BCGs of massive clusters.Comment: 13 pages, 17 figures, MNRAS accepted versio

We have updated and extended our semi-analytic galaxy formation modelling capabilities and applied them simultaneously to the stored halo/subhalo merger trees of the Millennium and Millennium-II Simulations (MS and MS-II, respectively). These differ by a factor of 125 in mass resolution, allowing explicit testing of resolution effects on predicted galaxy properties. We have revised the treatment of the transition between the rapid infall and cooling flow regimes of gas accretion, of the sizes of bulges, and of gaseous and stellar discs, of supernova feedback, of the transition between central and satellite status as galaxies fall into larger systems, and of gas and star stripping, once they become satellites. Plausible values of efficiency and scaling parameters yield an excellent fit not only to the observed abundance of low-redshift galaxies over five orders of magnitude in stellar mass and 9 mag in luminosity, but also to the observed abundance of Milky Way satellites. This suggests that reionization effects may not be needed to solve the 'missing-satellite' problem, except, perhaps, for the faintest objects. The same model matches the observed large-scale clustering of galaxies as a function of stellar mass and colour. The fit remains excellent down to ∼30 kpc for massive galaxies. For M * < 6 × 10 10 M , however, the model overpredicts clustering at scales below ∼1 Mpc, suggesting that the assumed fluctuation amplitude, σ 8 = 0.9, is too high. The observed difference in clustering between active and passive galaxies is matched quite well for all masses. Galaxy distributions within rich clusters agree between the simulations and match those observed, but only if galaxies without dark matter subhaloes (so-called orphans) are included. Even at MS-II resolution, schemes which assign galaxies only to resolved dark matter subhaloes cannot match observed clusters. Our model predicts a larger passive fraction among low-mass galaxies than is observed, as well as an overabundance of ∼10 10 M galaxies beyond z ∼ 0.6. (The abundance of ∼10 11 M galaxies is matched out to z ∼ 3.) These discrepancies appear to reflect deficiencies in the way star formation rates are modelled.

We take advantage of the largest high-resolution simulation of cosmic structure growth ever carried out -- the Millennium Simulation of the concordance LambdaCDM cosmogony -- to study how the star formation histories, ages and metallicities of elliptical galaxies depend on environment and on stellar mass. We concentrate on a galaxy formation model which is tuned to fit the joint luminosity/colour/morphology distribution of low redshift galaxies. Massive ellipticals in this model have higher metal abundances, older luminosity-weighted ages, shorter star formation timescales, but lower assembly redshifts than less massive systems. Within clusters the typical masses, ages and metal abundances of ellipticals are predicted to decrease, on average, with increasing distance from the cluster centre. We also quantify the effective number of progenitors of ellipticals as a function of present stellar mass, finding typical numbers below 2 for M* < 10^{11} Msun, rising to about 5 for the most massive systems. These findings are consistent with recent observational results that suggest ``down-sizing'' or ``anti-hierarchical'' behaviour for the star formation history of the elliptical galaxy population, despite the fact that our model includes all the standard elements of hierarchical galaxy formation and is implemented on the standard, LambdaCDM cosmogony.Comment: 12 pages, 11 figures, minor revisions, MNRAS accepte

We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1