Profiles of dark matter haloes at high redshift

I study the evolution of halo density profiles as a function of time in the SCDM and ΛCDM cosmologies. Following Del Popolo, I calculate the concentration parameter c = r _v / a and study its time evolution. For a given halo mass, I find that c ( z ) ∝ 1/(1+ z ) in both the ΛCDM and SCDM cosmology, in agreement with the analytic model of Bullock et al. and N -body simulations. In both models, a ( z ) is roughly constant. The present model predicts a stronger evolution of c ( z ) with respect to the Navarro, Frenk & White model. Finally I show some consequences of the results on galaxy modelling.     

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.     

Decaying dark matter and the deficit of dwarf haloes

The hierarchical clustering inherent in Λcold dark matter cosmology seems to produce many of the observed characteristics of large-scale structure. But some glaring problems still remain, including the overprediction (by a factor of 10) of the number of dwarf galaxies within the virialized population of the local group. Several secondary effects have already been proposed to resolve this problem. It is still not clear, however, whether the principal solution rests with astrophysical processes, such as early feedback from supernovae, or possibly with as yet undetermined properties of the dark matter itself. In this paper, we carry out a detailed calculation of the dwarf halo evolution incorporating the effects of a hypothesized dark matter decay, D → D'+ l , where D is the unstable particle, D ' is the more massive daughter particle and l is the other, lighter (or possibly massless) daughter particle. This process preferentially heats the smaller haloes, expanding them during their evolution and reducing their present-day circular velocity. We find that this mechanism can account very well for the factor of 4 deficit in the observed number of systems with velocity  10–20 km s⁻¹  compared to those predicted by the numerical simulations, if     , where Δ m is the mass difference between the initial and final states. The corresponding lifetime τ cannot be longer than ∼30 Gyr, but may be as short as just a few Gyr.     

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.     

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.     

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.     

The orientation of galaxy dark matter haloes around cosmic voids

Using the Millennium N -body Simulation we explore how the shape and angular momentum of galaxy dark matter haloes surrounding the largest cosmological voids are oriented. We find that the major and intermediate axes of the haloes tend to lie parallel to the surface of the voids, whereas the minor axis points preferentially in the radial direction. We have quantified the strength of these alignments at different radial distances from the void centres. The effect of these orientations is still detected at distances as large as 2.2 R _void from the void centre. Taking a subsample of haloes expected to contain disc-dominated galaxies at their centres we detect, at the 99.9 per cent confidence level, a signal that the angular momentum of those haloes tends to lie parallel to the surface of the voids. Contrary to the alignments of the inertia axes, this signal is only detected in shells at the void surface  (1 < R < 1.07  R _void)  and disappears at larger distances. This signal, together with the similar alignment observed using real spiral galaxies, strongly supports the prediction of the Tidal Torque theory that both dark matter haloes and baryonic matter have acquired, conjointly, their angular momentum before the moment of turnaround.     

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 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.