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 coarse geometry of merger trees in Λ cold dark matter

We introduce the contour process to describe the geometrical properties of merger trees. The contour process translates a tree structure into a one-dimensional object: the contour walk. We characterize the contour walk measuring its length and action. The length is proportional to the number of progenitors in the tree, and the action is a proxy for the mean length of a branch in the tree.
We construct the contour walk for merger trees extracted from the public data base of the Millennium Run and also for merger trees constructed with a public Monte Carlo code which implements a Markovian algorithm. The trees correspond to haloes of final masses between 10¹¹ and  10¹⁴  h ⁻¹ M_⊙  . We study how the length and action scale with the mass of the final halo.
In all cases, except for the action measured from Markovian trees, we find a transitional scale around  1–6 × 10¹²  h ⁻¹ M_⊙  , which can be interpreted as the point where the halo growth starts to be dominated by mergers and not accretion. The transitional scale depends weakly on the mass resolution and time discretization of the merging trees.     

Constraining cold dark matter halo merger rates using the coagulation equations

We place additional constraints on the three parameters of the dark matter halo merger rate function recently proposed by Parkinson, Cole & Helly by utilizing Smoluchowski's coagulation equation, which must be obeyed by any binary merging process which conserves mass. We find that the constraints from Smoluchowski's equation are degenerate, limiting to a thin plane in the three-dimensional parameter space. This constraint is consistent with those obtained from fitting to N -body measures of progenitor mass functions, and provides a better match to the evolution of the overall dark matter halo mass function, particularly for the most massive haloes. We demonstrate that the proposed merger rate function does not permit an exact solution of Smoluchowski's equation and, therefore, the choice of parameters must reflect a compromise between fitting various parts of the mass function. The techniques described herein are applicable to more general merger rate functions, which may permit a more accurate solution of Smoluchowski's equation. The current merger rate solutions are most probably sufficiently accurate for the vast majority of applications.     

Dark matter halo creation in moving barrier models

In hierarchical models of structure formation, the time derivative of the halo mass function may be thought of as the difference of two terms – a creation term, which describes the increase in the number of haloes of mass m from mergers of less massive objects, and a destruction term, which describes the decrease in the number of m -haloes as these merge with other haloes, creating more massive haloes as a result. The first part of this paper focuses on estimating the distribution of times when these creation events take place. In models where haloes form from a spherical collapse, this distribution can be estimated from the same formalism which is used to estimate halo abundances: the constant-barrier excursion-set approach. In the excursion-set approach, moving rather than constant barriers are necessary for estimating halo abundances when the collapse is triaxial. First, we generalize the excursion-set estimate of the creation time distribution by incorporating ellipsoidal collapse. Then, we show that these moving barrier based predictions are in better agreement with measurements in numerical simulations than are the corresponding predictions of the spherical collapse model. In the second part of the paper, we link the creation time distribution to the creation term mentioned above. For this quantity, the improvement provided by the ellipsoidal collapse model is more evident. These results should be useful for studies of merger-driven star formation rates and active galactic nucleus activity. We also present a similar study of the creation of haloes conditioned on belonging to an object of a certain mass today, and reach similar conclusions – the moving barrier based estimates are in substantially better agreement with the simulations. This part of the study may be useful for understanding the tendency for the oldest stars to exist in the most massive objects, and for star formation to only occur in lower mass objects at late times.     

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.     

The assembly bias of dark matter haloes to higher orders

We use an extremely large volume  (2.4  h ⁻³ Gpc³)  , high-resolution N -body simulation to measure the higher order clustering of dark matter haloes as a function of mass and internal structure. As a result of the large simulation volume and the use of a novel 'cross-moment' counts-in-cells technique which suppresses discreteness noise, we are able to measure the clustering of haloes corresponding to rarer peaks than was possible in previous studies; the rarest haloes for which we measure the variance are 100 times more clustered than the dark matter. We are able to extract, for the first time, halo bias parameters from linear up to fourth order. For all orders measured, we find that the bias parameters are a strong function of mass for haloes more massive than the characteristic mass   M _*  . Currently, no theoretical model is able to reproduce this mass dependence closely. We find that the bias parameters also depend on the internal structure of the halo up to fourth order. For haloes more massive than   M _*  , we find that the more concentrated haloes are more weakly clustered than the less concentrated ones. We see no dependence of clustering on concentration for haloes with masses   M < M _*  ; this is contrary to the trend reported in the literature when segregating haloes by their formation time. Our results are insensitive to whether haloes are labelled by the total mass returned by the friends-of-friends group finder or by the mass of the most massive substructure. This implies that our conclusions are not an artefact of the particular choice of group finding algorithm. Our results will provide important input to theoretical models of galaxy clustering.     

A comparison of simulated and analytic major merger counts

We use large volume, high resolution, N -body simulations of three different ΛCDM models, with different clustering strengths, to generate dark-matter halo merging histories. Over the reliable range of halo masses, roughly galaxy groups to rich clusters of galaxies, we quantify the number density of major mergers for two different time intervals and compare them with analytic predictions based on the extended Press–Schechter theory.     

Warm dark matter versus bumpy power spectra

In this paper we explore the differences between a warm dark matter (WDM) model and a cold dark matter (CDM) model where the power on a certain scale is reduced by introducing a narrow negative feature ('dip'). This dip is placed in a way so as to mimic the loss of power in the WDM model: both models have the same integrated power out to the scale where the power of the dip model rises to the level of the unperturbed CDM spectrum again.
Using N -body simulations we show that some of the large-scale clustering patterns of this new model follow more closely the usual CDM scenario while simultaneously suppressing small-scale structures (within galactic haloes) even more efficiently than WDM. The analysis in the paper shows that the new Dip model appears to be a viable alternative to WDM, but it is based on different physics. Where WDM requires the introduction of a new particle species, the Dip model is based on a non-standard inflationary period. If we are looking for an alternative to the currently challenged standard ΛCDM structure formation scenario, neither the ΛWDM nor the new Dip model can be ruled out based on the analysis presented in this paper. They both make very similar predictions and the degeneracy between them can only be broken with observations yet to come.     

Dark matter halo profiles in scale-free cosmologies

We explore the dependence of the central logarithmic slope of dark matter halo density profiles α on the spectral index n of the linear matter power spectrum P ( k ) using cosmological N -body simulations of scale-free models [i.e. P ( k ) ∝ k ⁿ ]. These simulations are based on a set of clear, reproducible and physically motivated criteria that fix the appropriate starting and stopping times for runs, and allow one to compare haloes across models with different spectral indices and mass resolutions. For each of our simulations we identify samples of well-resolved haloes in dynamical equilibrium and we analyse their mass profiles. By parametrizing the mass profile using a 'generalized' Navarro, Frenk & White profile in which the central logarithmic slope α is allowed to vary while preserving the r ⁻³ asymptotic form at large radii, we obtain preferred central slopes for haloes in each of our models. There is a strong correlation between α and n , such that α becomes shallower as n becomes steeper. However, if we normalize our mass profiles by r ₋₂, the radius at which the logarithmic slope of the density profile is −2, we find that these differences are no longer present. This is apparent if we plot the maximum slope     as a function of r / r ₋₂– we find that the profiles are similar for haloes forming in different n models. This reflects the importance of concentration, and reveals that the concentrations of haloes forming in steep- n cosmologies tend to be smaller than those of haloes forming in shallow- n cosmologies. We conclude that there is no evidence for convergence to a unique central asymptotic slope, at least on the scales that we can resolve.     

Tidal tails in cold dark matter cosmologies

We study the formation of tidal tails in pairs of merging disc galaxies with structural properties motivated by current theories of cold dark matter (CDM) cosmologies. In a recent study, Dubinski, Mihos & Hernquist showed that the formation of prominent tidal tails can be strongly suppressed by massive and extended dark haloes. For the large halo-to-disc mass ratio expected in CDM cosmologies their sequence of models failed to produce strong tails like those observed in many well-known pairs of interacting galaxies. In order to test whether this effect can constrain the viability of CDM cosmologies, we construct N ‐body models of disc galaxies with structural properties derived in analogy to the recent analytical work of Mo, Mao & White. With a series of self-consistent collisionless simulations of galaxy–galaxy mergers we demonstrate that even the discs of very massive dark haloes have no problems developing long tidal tails, provided the halo spin parameter is large enough. For our class of models, the halo-to-disc mass ratio is not a good indicator of the ability to produce tails. Instead, the relative size of disc and halo or, alternatively, the ratio of circular velocity to local escape speed at the half mass radius of the disc is a more useful criterion. This result holds in all CDM models. While tidal tails can provide useful information on the structure of galaxies, it thus appears unlikely that they are able to constrain the values of the cosmological parameters within these models.     

Clustering of galaxy clusters in cold dark matter universes

We use very large cosmological N -body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ( τ CDM) and a flat low density model with a cosmological constant (ΛCDM). Our simulations each use 10⁹ particles to follow the mass distribution within cubes of side 2  h ⁻¹ Gpc ( τ CDM) and 3  h ⁻¹ Gpc (ΛCDM) with a force resolution better than 10⁻⁴ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from     for samples with mean separation     to     for samples with     The lower value here corresponds to τ CDM and the upper to ΛCDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10 per cent using the simple model of Sheth, Mo & Tormen. This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the automated plate measurement (APM) cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for ΛCDM. The observed power spectrum (for a cluster sample of mean separation     ) lies significantly above the predictions of both models.     

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.     

Orbital parameters of infalling dark matter substructures

We present distributions of the orbital parameters of dark matter substructures at the time of merging into their host halo. Accurate knowledge of the orbits of dark matter substructures is a crucial input to studies which aim to assess the effects of the cluster environment on galaxies, the heating of galaxy discs and many other topics. Orbits are measured for satellites in a large number of N -body simulations. We focus on the distribution of radial and tangential velocities, but consider also distributions of orbital eccentricity and semimajor axis. We show that the distribution of radial and tangential velocities has a simple form and provide a fitting formula for this distribution. We also search for possible correlations between the infall directions of pairs of satellites, finding evidence for positive correlation at small angular separations as expected if some infall occurs along filaments. We also find (weak) evidence for correlations between the direction of the infall and infall velocity and the spin of the host halo.