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.     

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.     

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 fate of substructures in cold dark matter haloes

We use the Millennium Simulation, a large, high-resolution N -body simulation of the evolution of structure in a Λ cold dark matter cosmology, to study the properties and fate of substructures within a large sample of dark matter haloes. We find that the subhalo mass function departs significantly from a power law at the high-mass end. We also find that the radial and angular distributions of substructures depend on subhalo mass. In particular, high-mass subhaloes tend to be less radially concentrated and to have angular distributions closer to the direction perpendicular to the spin of the host halo than their less massive counterparts. We find that mergers between subhaloes occur. These tend to be between substructures that were already dynamically associated before accretion into the main halo. For subhaloes larger than 0.001 times the mass of the host halo, it is more likely that the subhalo will merge with the central or main subhalo than with another subhalo larger than itself. For lower masses, subhalo–subhalo mergers become equally likely to mergers with the main subhalo. Our results have implications for the variation of galaxy properties with environment and for the treatment of mergers in galaxy formation 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.     

An excursion set model for the distribution of dark matter and dark matter haloes

A model of the gravitationally evolved dark matter distribution, in the Eulerian space, is developed. It is a simple extension of the excursion set model that is commonly used to estimate the mass function of collapsed dark matter haloes. In addition to describing the evolution of the Eulerian space distribution of the haloes, the model allows one to describe the evolution of the dark matter itself. It can also be used to describe density profiles, on scales larger than the virial radius of these haloes, and to quantify the way in which matter flows in and out of Eulerian cells. When the initial Lagrangian space distribution is white noise Gaussian, the model suggests that the Inverse Gaussian distribution should provide a reasonably good approximation to the evolved Eulerian density field, in agreement with numerical simulations. Application of this model to clustering from more general Gaussian initial conditions is discussed at the end.     

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.     

Phase-space structure in the local dark matter distribution and its signature in direct detection experiments

We study predictions for dark matter (DM) phase-space structure near the Sun based on high-resolution simulations of six galaxy haloes taken from the Aquarius project. The local DM density distribution is predicted to be remarkably smooth; the density at the Sun differs from the mean over a best-fitting ellipsoidal equidensity contour by less than 15 per cent at the 99.9 per cent confidence level. The local velocity distribution is also very smooth, but it differs systematically from a (multivariate) Gaussian distribution. This is not due to the presence of individual clumps or streams, but to broad features in the velocity modulus and energy distributions that are stable in both space and time and reflect the detailed assembly history of each halo. These features have a significant impact on the signals predicted for weakly interacting massive particle and axion searches. For example, weakly interacting massive particles recoil rates can deviate by ∼10 per cent from those expected from the best-fitting multivariate Gaussian models. The axion spectra in our simulations typically peak at lower frequencies than in the case of multivariate Gaussian velocity distributions. Also in this case, the spectra show significant imprints of the formation of the halo. This implies that once direct DM detection has become routine, features in the detector signal will allow us to study the DM assembly history of the Milky Way. A new field, 'DM astronomy', will then emerge.     

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.     

Scalar fields as dark matter in spiral galaxies: comparison with experiments

We present an exact solution for a static and axially symmetric spacetime, which is obtained from a scalar-tensor theory that comes from unification theories. As an attempt to model the dark matter (DM) in spiral galaxies we find that an exponential scalar potential is enough to explain the rotation curves in such galaxies. We also present the fitting to the rotation curve of six spiral galaxies and we find an excellent agreement between observational data and the results of our model.     

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.     

Generating dark matter halo merger trees

We present a new Monte Carlo algorithm to generate merger trees describing the formation history of dark matter haloes. The algorithm is a modification of the algorithm of Cole et al. used in the galform semi-analytic galaxy formation model. As such, it is based on the Extended Press–Schechter theory and so should be applicable to hierarchical models with a wide range of power spectra and cosmological models. It is tuned to be in accurate agreement with the conditional mass functions found in the analysis of merger trees extracted from the Λ cold dark matter Millennium N -body simulation. We present a comparison of its predictions not only with these conditional mass functions, but also with additional statistics of the Millennium Simulation halo merger histories. In all cases, we find it to be in good agreement with the Millennium Simulation and thus it should prove to be a very useful tool for semi-analytic models of galaxy formation and for modelling hierarchical structure formation in general. We have made our merger tree generation code and code to navigate the trees available at http://star-www.dur.ac.uk/~cole/merger_trees .