Evolution of the mass function of dark matter haloes

We use a high-resolution ΛCDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of 15, in a  50 h ⁻¹ Mpc  volume containing 80 million particles. Our low-redshift results allow us to probe low-σ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth & Tormen mass function provides an excellent match to all of our data except for redshifts of 10 and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 per cent for the  10¹⁰–10¹¹ M_⊙  haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth & Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high-redshift  ( z ≃ 6) QSOs  , which are thought to be hosted by haloes at 5σ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50 per cent. We test the sensitivity of our results to force accuracy, starting redshift and halo-finding algorithm.     

Universal profile of dark matter haloes and the spherical infall model

I propose a modification of the spherical infall model for the evolution of density fluctuations with initially Gaussian probability distribution and scale-free power spectra in the Einsteinde Sitter universe as developed by Hoffman & Shaham. I introduce a generalized form of the initial density distribution around an overdense region and cut it off at half the interpeak separation, accounting in this way for the presence of the neighbouring fluctuations. Contrary to the original predictions of Hoffman & Shaham, the resulting density profiles within virial radii no longer have a power-law shape, but their steepness increases with distance. The profiles of haloes of galactic mass are well fitted by the universal profile formula of changing slope obtained as a result of N -body simulations by Navarro, Frenk & White. The trend of steeper profiles for smaller masses and higher spectral indices is also reproduced. The agreement between the model and simulations is better for smaller masses and lower spectral indices, which suggests that galaxies form mainly by accretion, while formation of clusters involves merging.     

The distribution function of dark matter in massive haloes

We study the distribution function (DF) of dark matter particles in haloes of mass range  10¹⁴–10¹⁵ M_⊙  . In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Λ cold dark matter (ΛCDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L , a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form     . This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (β₀ and  β_∞  ) and the scale of transition between them ( L ₀). The energy part   f _E ( E )  is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density   Q =ρ( r )/σ³( r )  .     

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.     

Excursion set approach to the clustering of dark matter haloes in Lagrangian space

We present a stochastic approach to the spatial clustering of dark matter haloes in Lagrangian space. Our formalism is based on a local formulation of the 'excursion set' approach by Bond et al., which automatically accounts for the 'cloud-in-cloud' problem in the identification of bound systems. Our method allows us to calculate correlation functions of haloes in Lagrangian space using either a multidimensional Fokker–Planck equation with suitable boundary conditions, or an array of Langevin equations with spatially correlated random forces. We compare the results of our method with theoretical predictions for the halo autocorrelation function considered in the literature, and find good agreement with the results recently obtained within a treatment of halo clustering in terms of 'counting fields' by Catelan et al. Finally, the possible effect of spatial correlations on numerical simulations of halo merger trees is discussed.     

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.     

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.     

An improved model for contraction of dark matter haloes in response to condensation of baryons

The cooling of gas in the centres of dark matter haloes is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. Following Gnedin et al., we improve this model by modifying the assumed invariant from M ( r ) r to     , where r and     are the current and orbit-averaged particle positions. We explore the effect of the bulge in the inner regions of the halo for different values of the bulge-to-disc mass ratio. We find that the bulge makes the velocity curve rise faster in the inner regions of the halo. We present an analytical fitting curve that describes the velocity curve of the halo after dissipation. The results should be useful for dark matter detection studies.     

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.     

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.     

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.