The cold dark matter model has become the leading theoretical picture for the formation of structure in the Universe. This model, together with the theory of cosmic inflation, makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a simulation of the growth of dark matter structure using 2,160(3) particles, following them from redshift z = 127 to the present in a cube-shaped region 2.230 billion lightyears on a side. In postprocessing, we also follow the formation and evolution of the galaxies and quasars. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with future generations of observational surveys of galaxies.

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 the Millennium Simulation (MS) to study the statistics of Λ cold dark matter (ΛCDM) halo concentrations at z= 0. Our results confirm that the average halo concentration declines monotonically with mass; the concentration–mass relation is well fitted by a power law over three decades in mass, up to the most massive objects that form in a ΛCDM universe (∼ 1015 h−1 M⊙). This is in clear disagreement with the predictions of the model proposed by Bullock et al. for these rare objects, and agrees better with the original predictions of Navarro, Frenk & White. The large volume surveyed, together with the unprecedented numerical resolution of the MS, allows us to estimate with confidence the distribution of concentrations and, consequently, the abundance of systems with unusual properties. About one in a hundred cluster haloes (M200≳ 3 × 1014 h−1 M⊙) have concentrations exceeding c200= 7.5, a result that may be useful in interpreting the likelihood of unusually strong massive gravitational lenses, such as Abell 1689, in the ΛCDM cosmogony. A similar fraction of about 1 per cent of galaxy‐sized haloes (M200∼ 1012 h−1 M⊙) have c200 < 4.5 and this could be relevant to models that attempt to reconcile the ΛCDM cosmology with rotation curves of low surface brightness galaxies by appealing to haloes of unexpectedly low concentration. We find that halo concentrations are independent of spin once haloes manifestly out of equilibrium have been removed from the sample. Compared to their relaxed brethren, the concentrations of out‐of‐equilibrium haloes tend to be lower and have more scatter, while their spins tend to be higher. A number of previously noted trends within the halo population are induced primarily by these properties of unrelaxed systems. Finally, we compare the result of predicting halo concentrations using the mass assembly history of the main progenitor with predictions based on simple arguments regarding the assembly time of all progenitors. The latter are typically as good or better than the former, suggesting that halo concentration depends not only on the evolutionary path of a halo's main progenitor, but on how and when all of its constituents collapsed to form non‐linear objects.

We investigate the subhalo populations of dark matter haloes in the concordance Λ cold dark matter (ΛCDM) cosmology. We use a large cosmological simulation and a variety of high‐resolution resimulations of individual cluster and galaxy haloes to study the systematics of subhalo populations over ranges of 1000 in halo mass and 1000 in the ratio of subhalo to parent halo mass. The subhalo populations of different haloes are not scaled copies of each other, but vary systematically with halo properties. On average, the amount of substructure increases with halo mass. At fixed mass, it decreases with halo concentration and with halo formation redshift. These trends are comparable in size to the scatter in subhalo abundance between similar haloes. Averaged over all haloes of given mass, the abundance of low‐mass subhaloes per unit parent halo mass is independent of parent mass. It is very similar to the abundance per unit mass of low‐mass haloes in the Universe as a whole, once differing boundary definitions for subhaloes and haloes are accounted for. The radial distribution of subhaloes within their parent haloes is substantially less centrally concentrated than that of the dark matter. It varies at most weakly with the mass (or concentration) of the parent halo and not at all with subhalo mass. It does depend on the criteria used to define the subhalo population considered. About 90 per cent of present‐day subhaloes were accreted after z= 1 and about 70 per cent after z= 0.5. Only about 8 per cent of the total mass of all haloes accreted at z= 1 survives as bound subhaloes at z= 0. For haloes accreted at z= 2, the survival mass fraction is just 2 per cent. Subhaloes seen near the centre of their parent typically were accreted earlier and retain less of their original mass than those seen near the edge. These strong systematics mean that comparison with galaxies in real clusters is only possible if the formation of the luminous component is modelled appropriately.

We use two very large cosmological simulations to study how the density profiles of relaxed cold dark matter dark haloes depend on redshift and on halo mass. We confirm that these profiles deviate slightly but systematically from the NFW form and are better approximated by the empirical formula, d log ρ/d log r ∝ r α , first used by Einasto to fit star counts in the Milky Way. The best-fitting value of the additional shape parameter, α, increases gradually with mass, from α ∼ 0.16 for present-day galaxy haloes to α ∼ 0.3 for the rarest and most massive clusters. Halo concentrations depend only weakly on mass at z = 0, and this dependence weakens further at earlier times. At z ∼ 3 the average concentration of relaxed haloes does not vary appreciably over the mass range accessible to our simulations (M 3 × 10 11 h −1 M ). Furthermore, in our biggest simulation, the average concentration of the most massive, relaxed haloes is constant at c 200 ∼ 3.5-4 for 0 z 3. These results agree well with those of Zhao et al. and support the idea that halo densities reflect the density of the universe at the time they formed, as proposed by Navarro, Frenk & White. With their original parameters, the NFW prescription overpredicts halo concentrations at high redshift. This shortcoming can be reduced by modifying the definition of halo formation time, although the evolution of the concentrations of Milky Way mass haloes is still not reproduced well. In contrast, the muchused revisions of the NFW prescription by Bullock et al. and Eke, Navarro & Steinmetz predict a steeper drop in concentration at the highest masses and stronger evolution with redshift than are compatible with our numerical data. Modifying the parameters of these models can reduce the discrepancy at high masses, but the overly rapid redshift evolution remains. These results have important implications for currently planned surveys of distant clusters.

We use a very large simulation of the concordance LCDM cosmogony to study the clustering of dark matter haloes. For haloes less massive than about 1e13Msun/h, the amplitude of the two-point correlation function on large scales depends strongly on halo formation time. Haloes that assembled at high redshift are substantially more clustered than those that assembled more recently. The effect is a smooth function of halo formation time and its amplitude increases with decreasing halo mass. At 1e11 Msun/h the ``oldest'' 10% of haloes are more than 5 times more strongly correlated than the ``youngest'' 10%. This unexpected result is incompatible with the standard excursion set theory for structure growth, and it contradicts a fundamental assumption of the halo occupation distribution models often used to study galaxy clustering, namely that the galaxy content of a halo of given mass is statistically independent of its larger scale environment.Comment: 5 pages, 5 figures, MNRAS in press. Full resolution pdf file is avaliable at http://www.mpa-garching.mpg.de/~gaoliang/GSW.pd

We perform cosmological hydrodynamic simulations to determine to what extent galaxies lose their gas due to photoheating from an ionizing background. We find that the characteristic mass at which haloes on average have lost half of their baryons is Mc∼ 6.5 × 109 h−1 M⊙ at z= 0, which corresponds to a circular velocity of 25 km s−1. This is significantly lower than the filtering mass obtained by the linear theory, which is often used in semi‐analytical models of galaxy formation. We demonstrate it is the gas temperature at the virial radius which determines whether a halo can accrete gas. A simple model that follows the merger history of the dark matter progenitors, and where gas accretion is not allowed when this temperature is higher than the virial temperature of the halo, reproduces the results from the simulation remarkably well. This model can be applied to any reionization history, and is easy to incorporate in semi‐analytical models.

No abstract