Evaluating approximations for halo merging histories

We study the merging history of dark matter haloes in N -body simulations and semi-analytical 'merger trees' based on the extended Press–Schechter (EPS) formalism. The main focus of our study is the joint distribution of progenitor number and mass as a function of redshift and parent halo mass. We begin by investigating the mean quantities predicted directly by the Press–Schechter (PS) and EPS formalism, such as the halo mass and conditional mass functions, and compare these predictions with the results of the simulations. The higher moments of this distribution are not predicted by the EPS formalism alone and must be obtained from the merger trees. We find that the Press–Schechter model deviates from the simulations at the level of 30–50 per cent on certain mass scales, and that the sense of the discrepancy changes as a function of redshift. We show that this discrepancy is reflected in the higher moments of the distribution of progenitor mass and number. We investigate some related statistics such as the accretion rate and the mass ratio of the largest two progenitors. For galaxy sized haloes ( M ∼10¹² M_⊙), we find that the merging history of haloes, as represented by these statistics, is well reproduced in the merger trees compared with the simulations. The agreement deteriorates for larger mass haloes. We conclude that merger trees based on the extended Press–Schechter formalism provide a reasonably reliable framework for semi-analytical models of galaxy formation.     

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.     

A marked correlation function analysis of halo formation times in the Millennium Simulation

We study the environmental dependence of the formation epoch of dark matter haloes in the Millennium Simulation: a ten billion particle N -body simulation of standard Lambda cold dark matter cosmology. A sensitive test of this dependence – the marked correlation function – reveals highly significant evidence that haloes of a given mass form earlier in denser regions. We define a marked cross-correlation function, which helps quantify how this effect depends upon the choice of the halo population used to define the environment. The mean halo formation redshift as a function of the local overdensity in dark matter is also well determined, and we see an especially clear dependence for galaxy-sized haloes. This contradicts one of the basic predictions of the excursion set model of structure formation, even though we see that this theory predicts other features of the distribution of halo formation epochs rather well. It also invalidates an assumption usually employed in the popular halo, or halo occupation distribution, models of galaxy clustering, namely that the distribution of halo properties is a function of halo mass but not of halo environment.     

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.     

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.     

On the origin of cold dark matter halo density profiles

N -body simulations predict that cold dark matter (CDM) halo-assembly occurs in two phases: (i) a fast-accretion phase with a rapidly deepening potential well; and (ii) a slow-accretion phase characterized by a gentle addition of mass to the outer halo with little change in the inner potential well. We demonstrate, using one-dimensional simulations, that this two-phase accretion leads to CDM haloes of the Navarro, Frenk & White (NFW) form and provides physical insight into the properties of the mass-accretion history that influence the final profile. Assuming that the velocities of CDM particles are effectively isotropized by fluctuations in the gravitational potential during the fast-accretion phase, we show that gravitational collapse in this phase leads to an inner profile  ρ( r ) ∝ r ⁻¹  . Slow accretion on to an established potential well leads to an outer profile with  ρ( r ) ∝ r ⁻³  . The concentration of a halo is determined by the fraction of mass that is accreted during the fast-accretion phase. Using an ensemble of realistic mass-accretion histories, we show that the model predictions of the dependence of halo concentration on halo formation time and, hence, the dependence of halo concentration on halo mass, and the distribution of halo concentrations all match those found in cosmological N -body simulations. Using a simple analytic model that captures much of the important physics, we show that the inner   r ⁻¹  profile of CDM haloes is a natural result of hierarchical mass assembly with an initial phase of rapid accretion.     

Analytical approach to subhalo population in dark matter haloes

In the standard model of cosmic structure formation, dark matter haloes form by gravitational instability. The process is hierarchical: smaller systems collapse earlier, and later merge to form larger haloes. The galaxy clusters, hosted by the largest dark matter haloes, are at the top of this hierarchy and representing the largest as well as the last structures formed in the Universe, while the smaller and first haloes are those Earth-sized dark subhaloes that have been both predicted by theoretical considerations and found in numerical simulations, though there do not exist any observational hints of their existence. The probability that a halo of mass m at redshift z will be part of a larger halo of mass M at the present time can be described in the frame of the extended Press & Schecter theory making use of the progenitor (conditional) mass function. Using the progenitor mass function, we calculate analytically, at redshift zero, the distribution of subhaloes in mass, formation epoch and rarity of the peak of the density field at the formation epoch. That is done for a Milky Way size system, assuming both a spherical and an ellipsoidal collapse model. Our calculation assumes that small progenitors do not lose mass due to dynamical processes after entering the parent halo, and that they do not interact with other subhaloes. For a Λ cold dark matter power spectrum, we obtain a subhalo mass function  d n /d m   proportional to   m ^−α  with a model-independent  α∼ 2  . Assuming that the dark matter is a weakly interacting massive particle, the inferred distributions are used to test the feasibility of an indirect detection in the γ-ray energy band of such a population of subhaloes with a Gamma-ray Large Area Space Telescope like satellite.     

An analytic model for the epoch of halo creation

In this paper we describe the Bayesian link between the cosmological mass function and the distribution of times at which isolated haloes of a given mass exist. By assuming that clumps of dark matter undergo monotonic growth on the time-scales of interest, this distribution of times is also the distribution of 'creation' times of the haloes. This monotonic growth is an inevitable aspect of gravitational instability. The spherical top-hat collapse model is used to estimate the rate at which clumps of dark matter collapse. This gives the prior for the creation time given no information about halo mass. Applying Bayes' theorem then allows any mass function to be converted into a distribution of times at which haloes of a given mass are created. This general result covers both Gaussian and non-Gaussian models. We also demonstrate how the mass function and the creation time distribution can be combined to give a joint density function, and discuss the relation between the time distribution of major merger events and the formula calculated. Finally, we determine the creation time of haloes within three N -body simulations, and compare the link between the mass function and creation rate with the analytic theory.     

Survival of substructure within dark matter haloes

Using high-resolution cosmological N -body simulations, we investigate the survival of dark matter satellites falling into larger haloes. Satellites preserve their identity for some time after merging. We compute their loss of mass, energy and angular momentum as they are dissolved by dynamical friction, tidal forces and collisions with other satellites. We also analyse the evolution of their internal structure. Satellites with less than a few per cent of the mass of the main halo may survive for several billion years, whereas larger satellites rapidly sink into the centre of the main halo potential well and lose their identity. Penetrating encounters between satellites are frequent and may lead to significant mass loss and disruption. Only a minor fraction of cluster mass (10–15 per cent on average) is bound to substructure at most redshifts of interest. We discuss the application of these results to the survival and extent of dark matter haloes associated with galaxies in clusters, and to their interactions. We find that a minor fraction of galaxy-size dark matter haloes are disrupted by redshift z  = 0. The fraction of satellites undergoing close encounters is similar to the observed fraction of interacting or merging galaxies in clusters at moderate redshift.     

A simple model for the evolution of supermassive black holes and the quasar population

An empirically motivated model is presented for accretion-dominated growth of supermassive black holes (SMBH) in galaxies, and the implications are studied for the evolution of the quasar population in the Universe. We investigate the core aspects of the quasar population, including space density evolution, evolution of the characteristic luminosity, plausible minimum masses of quasars, the mass function of SMBH and their formation epoch distribution. Our model suggests that the characteristic luminosity in the quasar luminosity function arises primarily as a consequence of a characteristic mass scale above which there is a systematic separation between the black hole and the halo merging rates. At lower mass scales, black hole merging closely tracks the merging of dark haloes. When combined with a declining efficiency of black hole formation with redshift, the model can reproduce the quasar luminosity function over a wide range of redshifts. The observed space density evolution of quasars is well described by formation rates of SMBH above  ∼10⁸  M_⊙  . The inferred mass density of SMBH agrees with that found independently from estimates of the SMBH mass function derived empirically from the quasar luminosity function.     

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.     

Dynamical friction in dwarf galaxies

We present a simplified analytic approach to the problem of the spiralling of a massive body orbiting within the dark halo of a dwarf galaxy. This dark halo is treated as the core region of a King distribution of dark matter particles, in consistency with the observational result of dwarf galaxies having solid-body rotation curves. Thus we derive a simple formula which provides a reliable and general first-order solution to the problem, totally analogous to the one corresponding to the dynamical friction problem in an isothermal halo. This analytic approach allows a clear handling and a transparent understanding of the physics and the scaling of the problem. A comparison with the isothermal case shows that in the core regions of a King sphere, dynamical friction proceeds at a different rate, and is sensitive to the total core radius. Thus, in principle, observable consequences may result. In order to illustrate the possible effects, we apply this formula to the spiralling of globular cluster orbits in dwarf galaxies, and show how present-day globular cluster systems could, in principle, be used to derive better limits on the structure of dark haloes around dwarf galaxies, when the observational situation improves. As a second application, we study the way a massive black hole population forming a fraction of these dark haloes would gradually concentrate towards the centre, with a consequent deformation of an originally solid-body rotation curve. This effect allows us to set limits on the fraction/mass of any massive black hole minority component of the dark haloes of dwarf galaxies. In essence, we take advantage of the way the global matter distribution fixes the local distribution function for the dark matter particles, which in turn determines the dynamical friction problem.     

Effects of correlation between merging steps on the global halo formation

The excursion set theory of halo formation is modified by adopting the fractional Brownian motion, to account for possible correlation between merging steps. We worked out analytically the conditional mass function, halo merging rate and formation time distribution in the spherical collapse model. We also developed an approximation for the ellipsoidal collapse model and applied it to the calculation of the conditional mass function and the halo formation time distribution. For models in which the steps are positively correlated, the halo merger rate is enhanced when the accreted mass is less than  ∼25 M *  , while for the negatively correlated case this rate is reduced. Compared with the standard model in which the steps are uncorrelated, the models with positively correlated steps produce more aged population in small mass haloes and more younger population in large mass haloes, while for the models with negatively correlated steps the opposite is true. An examination of simulation results shows that a weakly positive correlation between successive merging steps appears to fit best. We have also found a systematic effect in the measured mass function due to the finite volume of simulations. In future work, this will be included in the halo model to accurately predict the three-point correlation function estimated from simulations.     

Effects of merging histories on angular momentum distribution of dark matter haloes

The effects of merging histories of proto-objects on the angular momentum distributions of the present-time dark matter haloes are analysed. An analytical approach to the analysis of the angular momentum distributions assumes that the haloes are initially homogeneous ellipsoids and that the growth of the angular momentum of the haloes halts at their maximum expansion time. However, the maximum expansion time cannot be determined uniquely, because in the hierarchical clustering scenario each progenitor, or subunit, of the halo has its own maximum expansion time. Therefore the merging history of the halo may be important in estimating its angular momentum. Using the merger tree model by Rodrigues &38; Thomas, which takes into account the spatial correlations of the density fluctuations, we have investigated the effects of the merging histories on the angular momentum distributions of dark matter haloes. It was found that the merger effects, that is, the effects of the inhomogeneity of the maximum expansion times of the progenitors which finally merge together into a halo, do not strongly affect the final angular momentum distributions, so that the homogeneous ellipsoid approximation happens to be good for the estimation of the angular momentum distribution of dark matter haloes. This is because the effect of the different directions of the angular momenta of the progenitors cancels out the effect of the inhomogeneity of the maximum expansion times of the progenitors.   The contribution of the orbital angular momentum to the total angular momentum when two or more pre-existing haloes merge together was also investigated. It is shown that this contribution is more important than that of the angular momentum of diffuse accreting matter to the total angular momentum, especially when the mergers occur many times.     

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.     

Constraints on dark matter physics from dwarf galaxies through galaxy cluster haloes

One of the predictions of the standard cold dark matter model is that dark haloes have centrally divergent density profiles. An extensive body of rotation curve observations of dwarf and low surface brightness galaxies shows the dark haloes of those systems to be characterized by soft constant-density central cores. Several physical processes have been proposed to produce soft cores in dark haloes, each one with different scaling properties. With the aim of discriminating among them we have examined the rotation curves of dark-matter-dominated dwarf and low surface brightness galaxies and the inner mass profiles of two clusters of galaxies lacking a central cD galaxy and with evidence of soft cores in the centre. The core radii and central densities of these haloes scale in a well-defined manner with the depth of their potential wells, as measured through the maximum circular velocity. As a result of our analysis we identify self-interacting cold dark matter as a viable solution to the core problem, where a non-singular isothermal core is formed in the halo centre surrounded by a Navarro, Frenk & White profile in the outer parts. We show that this particular physical situation predicts core radii in agreement with observations. Furthermore, using the observed scalings, we derive an expression for the minimum cross-section ( σ ) which has an explicit dependence with the halo dispersion velocity ( v ). If m _x is the mass of the dark matter particle: σ m _x ≈4×10⁻²⁵ (100 km s⁻¹  v ⁻¹) cm² GeV⁻¹.