Testing the modified Press–Schechter model against N-body simulations

A modified version of the extended Press–Schechter model for the growth of dark-matter haloes was introduced in two previous papers, with the aim of explaining the mass–density relation shown by haloes in high-resolution cosmological simulations. In this model, major mergers are well separated from accretion, thereby allowing a natural definition of halo formation and destruction. This makes it possible to derive analytic expressions for halo formation and destruction rates, the mass accretion rate and the probability distribution functions of halo formation times and progenitor masses. The stochastic merger histories of haloes can be readily derived and easily incorporated into semi-analytical models of galaxy formation, thus avoiding the usual problems encountered in the construction of Monte Carlo merger trees from the original extended Press–Schechter formalism. Here we show that the predictions of the modified Press–Schechter model are in good agreement with the results of N -body simulations for several scale-free cosmologies.     

Merger rates of dark matter haloes

We derive analytic merger rates for dark matter haloes within the framework of the extended Press–Schechter (EPS) formalism. These rates become self-consistent within EPS once we realize that the typical merger in the limit of a small time-step involves more than two progenitors, contrary to the assumption of binary mergers adopted in earlier studies. We present a general method for computing merger rates that span the range of solutions permitted by the EPS conditional mass function, and focus on a specific solution that attempts to match the merger rates in N -body simulations. The corrected EPS merger rates are more accurate than the earlier estimates of Lacey & Cole by ∼20 per cent for major mergers and by up to a factor of ∼3 for minor mergers of mass ratio 1:10⁴. Based on the revised merger rates, we provide a new algorithm for constructing Monte Carlo EPS merger trees, which could be useful in semi-analytic modelling. We provide analytic expressions and plot numerical results for several quantities that are very useful in studies of galaxy formation. This includes (i) the rate of mergers of a given mass ratio per given final halo, (ii) the fraction of mass added by mergers to a halo and (iii) the rate of mergers per given main progenitor. The creation and destruction rates of haloes serve for a self-consistency check. Our method for computing merger rates can be applied to conditional mass functions beyond EPS, such as those obtained by the ellipsoidal collapse model or extracted from N -body simulations.     

On the assembly history of dark matter haloes

We study the mass assembly history (MAH) of dark matter haloes. We compare MAHs obtained using (i) merger trees constructed with the extended Press–Schechter (EPS) formalism, (ii) numerical simulations and (iii) the Lagrangian perturbation code pinocchio . We show that the pinocchio MAHs are in excellent agreement with those obtained using numerical simulations, while the EPS formalism predicts MAHs that occur too late. pinocchio , which is much less CPU intensive than N -body simulation, can be run on a simple personal computer, and does not require any labour intensive post-simulation analysis, therefore provides a unique and powerful tool to investigate the growth history of dark matter haloes. Using a suite of 55 pinocchio simulations, with 256³ particles each, we study the MAHs of 12 924 cold dark matter (CDM) haloes in a ΛCDM concordance cosmology. This is by far the largest set of haloes used for any such analysis. For each MAH we derive four different formation redshifts, which characterize different epochs during the assembly history of a dark matter halo. We show that haloes less massive than the characteristic non-linear mass scale establish their potential wells much before they acquire most of their mass. The time when a halo reaches its maximum virial velocity roughly divides its mass assembly into two phases, a fast-accretion phase which is dominated by major mergers, and a slow-accretion phase dominated by minor mergers. Each halo experiences about 3 ± 2 major mergers since its main progenitor had a mass equal to 1 per cent of the final halo mass. This major merger statistic is found to be virtually independent of halo mass. However, the average redshift at which these major mergers occur is strongly mass dependent, with more massive haloes experiencing their major mergers later.     

The statistical properties of Λ cold dark matter halo formation

We present a comparison of the statistical properties of dark matter halo merger trees extracted from the Millennium Simulation with Extended Press–Schechter (EPS) formalism and the related galform Monte Carlo method for generating ensembles of merger trees. The volume, mass resolution and output frequency make the Millennium Simulation a unique resource for the study of the hierarchical growth of structure. We construct the merger trees of present-day friends-of-friends groups and calculate a variety of statistics that quantify the masses of their progenitors as a function of redshift, accretion rates, and the redshift distribution of their most recent major merger. We also look in the forward direction and quantify the present-day mass distribution of haloes into which high-redshift progenitors of a specific mass become incorporated. We find that the EPS formalism and its Monte Carlo extension capture the qualitative behaviour of all these statistics, but as redshift increases they systematically underestimate the masses of the most massive progenitors. This shortcoming is worst for the Monte Carlo algorithm. We present a fitting function to a scaled version of the progenitor mass distribution and show how it can be used to make more accurate predictions of both progenitor and final halo mass distributions.     

The universal mass accretion history of cold dark matter haloes

We use the extended Press–Schechter formalism to investigate the rate at which cold dark matter haloes accrete mass. We discuss the shortcomings of previous methods that have been used to compute the mass accretion histories of dark matter haloes, and present an improved method based on the N -branch merger tree algorithm of Somerville & Kolatt. We show that this method no longer suffers from inconsistencies in halo formation times, and compare its predictions with high-resolution N -body simulations. Although the overall agreement is reasonable, there are slight inconsistencies which are most easily interpreted as a reflection of ellipsoidal collapse (as opposed to spherical collapse assumed in the Press–Schechter formalism). We show that the average mass accretion histories follow a simple, universal profile, and we present a simple recipe for computing the two scale-parameters which is applicable to a wide range of halo masses and cosmologies. Together with the universal profiles for the density and angular momentum distributions of cold dark matter haloes, these universal mass accretion histories provide a simple but accurate framework for modelling the structure and formation of dark matter haloes. In particular, they can be used as a backbone for modelling various aspects of galaxy formation where one is not interested in the detailed effects of merging. As an example we use the universal mass accretion history to compute the rate at which dark matter haloes accrete mass, which we compare with the cosmic star formation history of the Universe.     

Modelling angular-momentum history in dark-matter haloes

We model the acquisition of spin by dark-matter haloes in semi-analytic merger trees. We explore two different algorithms: one in which halo spin is acquired from the orbital angular momentum of merging satellites, and another in which halo spin is gained via tidal torquing on shells of material while still in the linear regime. We find that both scenarios produce the characteristic spin distribution of haloes found in N -body simulations, namely, a log-normal distribution with mean ≈ 0.04 and standard deviation ≈ 0.5 in the log. A perfect match requires fine-tuning of two free parameters. Both algorithms also reproduce the general insensitivity of the spin distribution to halo mass, redshift and cosmology seen in N -body simulations. The spin distribution can be made strictly constant by physically motivated scalings of the free parameters. In addition, both schemes predict that haloes that have had recent major mergers have systematically larger spin values. These algorithms can be implemented within semi-analytic models of galaxy formation based on merger trees. They yield detailed predictions of galaxy properties that strongly depend on angular momentum (such as size and surface brightness) as a function of merger history and environment.     

On the distribution of haloes, galaxies and mass

The stochasticity in the distribution of dark haloes in the cosmic density field is reflected in the distribution function   P _V ( N _h| δ _m)  , which gives the probability of finding N _h haloes in a volume V with mass density contrast δ _m. We study the properties of this function using high-resolution N -body simulations, and find that   P _V ( N _h| δ _m)  is significantly non-Poisson. The ratio between the variance and the mean goes from ∼1 (Poisson) at  1+ δ _m≪1  to <1 (sub-Poisson) at  1+ δ _m∼1  to >1 (super-Poisson) at  1+ δ _m≫1  . The mean bias relation is found to be well described by halo bias models based on the Press–Schechter formalism. The sub-Poisson variance can be explained as a result of halo exclusion, while the super-Poisson variance at high δ _m may be explained as a result of halo clustering. A simple phenomenological model is proposed to describe the behaviour of the variance as a function of δ _m. Galaxy distribution in the cosmic density field predicted by semi-analytic models of galaxy formation shows similar stochastic behaviour. We discuss the implications of the stochasticity in halo bias to the modelling of higher order moments of dark haloes and of galaxies.     

Dynamical friction and galaxy merging time-scales

The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N -body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M _sat/ M _host≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M _sat/ M _host≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halo's centre for a wide range of M _sat/ M _host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≲10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.     

Environmental dependence of dark matter halo growth – I. Halo merger rates

In an earlier paper, we quantified the mean merger rate of dark matter haloes as a function of redshift z , descendant halo mass M ₀, and progenitor halo mass ratio ξ using the Millennium simulation of the Λ cold dark matter cosmology. Here, we broaden that study and investigate the dependence of the merger rate of haloes on their surrounding environment. A number of local mass overdensity variables, both including and excluding the halo mass itself, are tested as measures of a halo's environment. The simple functional dependence on   z , M ₀  , and ξ of the merger rate found in our earlier work, is largely preserved in different environments, but we find that the overall amplitude of the merger rate has a strong positive correlation with the environmental densities. For galaxy-mass haloes, we find mergers to occur ∼2.5 times more frequently in the densest regions than in voids at both   z = 0  and higher redshifts. Higher mass haloes show similar trends. We present a fitting form for this environmental dependence that is a function of both mass and local density and valid out to   z = 2  . The amplitude of the progenitor (or conditional) mass function shows a similar correlation with local overdensity, suggesting that the extended Press–Schechter model for halo growth needs to be modified to incorporate environmental effects.     

Cosmological origin of the lowest metallicity halo stars

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars that are likely to be found within our Galaxy today. By PopIII we refer to stars formed at large redshift ( z >4), with low metallicities ([ Z /Z_⊙]<−2.5) and in small systems (total mass ≲ 2×10⁸ M_⊙) that are extremely sensitive to stellar feedback, and which through a prescribed merging history end up becoming part of the Milky Way today. An analytic, extended Press–Schechter formalism is used to obtain the mass functions of haloes which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modelled as a mini-galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disc instabilities, all of which determine the fraction of the baryons that are subject to star formation. The use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, and give accurate estimates of their expected numbers today and their location in L /L_⊙ versus T /K Hertzsprung–Russell (HR) diagrams. A first comparison with observational data suggests that the initial mass function (IMF) of the first stars was increasingly high-mass weighted towards high redshifts, levelling off at z ≳9 at a characteristic stellar mass scale m _s=10–15 M_⊙.     

Merging history trees for dark matter haloes: tests of the Merging Cell Model in a CDM cosmology

The merging history of dark matter haloes is computed with the Merging Cell Model proposed by Rodrigues & Thomas. While originally discussed in the case of scale-free power spectra, it is developed and tested here in the framework of the cold dark matter cosmology. The halo mass function, the mass distribution of progenitors and child haloes, as well as the probability distribution of formation times, have been computed and compared with the available analytic predictions. The halo autocorrelation function has also been obtained (a first for a semi-analytic merging tree), and tested against analytic formulae. An overall good agreement is found between results of the model, and the predictions derived from the Press & Schechter theory and its extensions. More severe discrepancies appear when formulae that better describe N -body simulations are used for comparison. In many instances, the model can be a useful tool for following the hierarchical growth of structures. In particular, it is suitable for addressing the issue of the formation and evolution of galaxy clusters, as well as the population of Lyman-break galaxies at high redshift, and their clustering properties.     

Numerical simulations of hot halo gas in galaxy mergers

Galaxy merger simulations have explored the behaviour of gas within the galactic disc, yet the dynamics of hot gas within the galaxy halo have been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with metal-free hot halo gas. To isolate the effect of the halo gas, we simulate only the dark matter halo and the hot halo gas over a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that (i) a strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to   T ∼ 10^6.3 K  . (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is  ≳10³⁹ erg s⁻¹  for halo gas fractions larger than 10 per cent. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as  10⁴² erg s⁻¹  . This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) ∼10–20 per cent of the initial gas mass is unbound from the galaxies for equal-mass mergers, while 3–5 per cent of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the intergalactic medium with metals.     

Simulating subhaloes at high redshift: merger rates, counts and types

Galaxies are believed to be in one-to-one correspondence with simulated dark matter subhaloes. We use high-resolution N -body simulations of cosmological volumes to calculate the statistical properties of subhalo (galaxy) major mergers at high redshift ( z = 0.6–5). We measure the evolution of the galaxy merger rate, finding that it is much shallower than the merger rate of dark matter host haloes at   z > 2.5  , but roughly parallels that of haloes at   z < 1.6  . We also track the detailed merger histories of individual galaxies and measure the likelihood of multiple mergers per halo or subhalo. We examine satellite merger statistics in detail: 15–35 per cent of all recently merged galaxies are satellites, and satellites are twice as likely as centrals to have had a recent major merger. Finally, we show how the differing evolution of the merger rates of haloes and galaxies leads to the evolution of the average satellite occupation per halo, noting that for a fixed halo mass, the satellite halo occupation peaks at   z ∼ 2.5  .     

The build-up of haloes within Press–Schechter theory

Modelling the build-up of haloes is important for linking the formation of galaxies with cosmological models. A simple model of halo growth is provided by Press–Schechter (PS) theory, where the initial field of density fluctuations is smoothed using spherically symmetric filters centred on a given position to obtain information about the likelihood of later collapse on varying scales. In this paper the predicted halo mass growth is compared for three filter shapes: Gaussian, top-hat and sharp k -space. Preliminary work is also presented analysing the build-up of haloes within numerical simulations using a friends-of-friends group finder. The best-fit to the simulation mass function was obtained using PS theory with a top-hat filter. By comparing both the backwards conditional mass function, which gives the distribution of halo progenitors, and the distribution of halo mergers in time, the build-up of haloes in the simulations is shown to be better fitted by PS theory with a sharp k -space filter. This strengthens previous work, which also found the build-up of haloes in simulations to be well matched to PS theory with a sharp k -space filter by providing a direct comparison of different filters and by extending the statistical tools used to analyse halo mass growth. The usefulness of this work is illustrated by showing that the cosmological evolution in the proportion of haloes that have undergone recent merger is predicted to be independent of mass and power spectrum and to only depend upon cosmology. Recent results from observations of field galaxies are shown to match the evolution expected, but are not sufficiently accurate to distinguish usefully between cosmological parameters.     

The bias field of dark matter haloes

Accepted 1998 January 26. Received 1998 January 26; in original form 1997 August 13This paper presents a stochastic approach to the clustering evolution of dark matter haloes in the Universe. Haloes, identified by a Press–Schechter-type algorithm in Lagrangian space, are described in terms of 'counting fields', acting as non-linear operators on the underlying Gaussian density fluctuations. By ensemble-averaging these counting fields, the standard Press–Schechter mass function as well as analytic expressions for the halo correlation function and corresponding bias factors of linear theory are obtained, extending the recent results by Mo & White. The non-linear evolution of our halo population is then followed by solving the continuity equation, under the sole hypothesis that haloes move by the action of gravity. This leads to an exact and general formula for the bias field of dark matter haloes, defined as the local ratio between their number density contrast and the mass density fluctuation. Besides being a function of position and 'observation' redshift, this random field depends upon the mass and formation epoch of the objects and is both non-linear and non-local. The latter features are expected to leave a detectable imprint on the spatial clustering of galaxies, as described, for instance, by statistics like the bispectrum and the skewness. Our algorithm may have several interesting applications, among which is the possibility of generating mock halo catalogues from low-resolution N -body simulations.     

Generating dark matter halo merger trees

We present a new Monte Carlo algorithm to generate merger trees describing the formation history of dark matter haloes. The algorithm is a modification of the algorithm of Cole et al. used in the galform semi-analytic galaxy formation model. As such, it is based on the Extended Press–Schechter theory and so should be applicable to hierarchical models with a wide range of power spectra and cosmological models. It is tuned to be in accurate agreement with the conditional mass functions found in the analysis of merger trees extracted from the Λ cold dark matter Millennium N -body simulation. We present a comparison of its predictions not only with these conditional mass functions, but also with additional statistics of the Millennium Simulation halo merger histories. In all cases, we find it to be in good agreement with the Millennium Simulation and thus it should prove to be a very useful tool for semi-analytic models of galaxy formation and for modelling hierarchical structure formation in general. We have made our merger tree generation code and code to navigate the trees available at http://star-www.dur.ac.uk/~cole/merger_trees .     

The non-Gaussian tail of cosmic-shear statistics

Owing to gravitational instability, an initially Gaussian density field develops non-Gaussian features as the Universe evolves. The most prominent non-Gaussian features are massive haloes, visible as clusters of galaxies. The distortion of high-redshift galaxy images because of the tidal gravitational field of the large-scale matter distribution, called cosmic shear, can be used to investigate the statistical properties of the large‐scale structure (LSS) . In particular, non-Gaussian properties of the LSS will lead to a non-Gaussian distribution of cosmic-shear statistic. The aperture mass ( M _ap) statistics, recently introduced as a measure for cosmic shear, is particularly well suited for measuring these non-Gaussian properties. In this paper we calculate the highly non-Gaussian tail of the aperture mass probability distribution, assuming Press–Schechter theory for the halo abundance and the 'universal' density profile of haloes as obtained from numerical simulations. We find that for values of M _ap much larger than its dispersion, this probability distribution is closely approximated by an exponential, rather than a Gaussian. We determine the amplitude and shape of this exponential for various cosmological models and aperture sizes, and show that wide-field imaging surveys can be used to distinguish between some of the currently most popular cosmogonies. Our study here is complementary to earlier cosmic-shear investigations, which focused more on two- and three-point statistical properties.     

Extended Press–Schechter theory and the density profiles of dark matter haloes

An inside–out model for the formation of haloes in a hierarchical clustering scenario is studied. The method combines the picture of the spherical infall model and a modification of the extended Press–Schechter theory. The mass accretion rate of a halo is defined to be the rate of its mass increase due to minor mergers. The accreted mass is deposited at the outer shells without changing the density profile of the halo inside its current virial radius. We applied the method to a flat Λ-cold dark matter universe. The resulting density profiles are compared with analytical models proposed in the literature, and a very good agreement is found. A trend is found of the inner density profile to become steeper for larger halo mass, which also results from recent N -body simulations. Additionally, present-day concentrations as well as their time evolution are derived and it is shown that they reproduce the results of large cosmological N -body simulations.