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.     

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.     

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.     

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.     

Dark matter haloes within clusters

We examine the properties of dark matter haloes within a rich galaxy cluster using a high-resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 haloes with circular velocities larger than 80 km s⁻¹ within the cluster virial radius of 2 Mpc (with Hubble constant H ₀ = 50 km s⁻¹ Mpc⁻¹). This enables an unprecedented study of the statistical properties of a large sample of dark matter haloes evolving in a dense environment. The cumulative fraction of mass attached to these haloes varies from close to zero per cent at 200 kpc to 13 per cent at the virial radius. Even at this resolution the overmerging problem persists; haloes that pass within 100–200 kpc of the cluster centre are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor haloes. The median ratio of apocentric to pericentric radii is 6:1, so that the orbital distribution is close to isotropic, circular orbits are rare and radial orbits are common. The orbits of haloes are unbiased with respect to both position within the cluster and the orbits of the smooth dark matter background, and no velocity bias is detected. The tidal radii of surviving haloes are generally well-fitted using the simple analytic prediction applied to their orbital pericentres. Haloes within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other haloes that cause significant mass loss —'galaxy harassment'. Mergers between haloes do not occur inside the cluster virial radius.     

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.