Angular momentum and clustering properties of early dark matter haloes

In this paper, we study the angular momentum properties of simulated dark matter haloes at high redshifts that likely host the first stars in the Universe. Calculating the spin distributions of these  10⁶– 10⁷ M_⊙  haloes in redshift slices from   z = 15  to 6, we find that they are well fit by a lognormal distribution as is found for lower redshift and more massive haloes in earlier work. We find that both the mean value of the spin and dispersion are largely unchanged with redshift for all haloes. Our key result is that subsamples of low- and high-spin, 10⁶ and  10⁷ M_⊙  , haloes show difference in clustering strength. In both mass bins, higher spin haloes are more strongly clustered in concordance with a tidal torquing picture for the growth of angular momentum in dark matter haloes in the cold dark matter paradigm.     

Evolution of the phase-space density of dark matter haloes and mixing effects in merger events

Cosmological N -body simulations were performed to study the evolution of the phase-space density   Q =ρ/σ³  of dark matter haloes. No significant differences in the scale relations   Q ∝σ^−2.1  or   Q ∝ M ^−0.82  are seen for the 'cold' or 'warm' dark matter models. The follow-up of individual haloes from   z = 10  up to the present time indicate the existence of two main evolutionary phases: an early and fast one  (10 > z > 6.5)  , in which Q decreases on the average by a factor of 40 as a consequence of the randomization of bulk motions, and a late and long one  (6.5 > z ≥ 0)  , in which Q decreases by a factor of 20 because of mixing induced by merger events. The study of these haloes has also evidenced that rapid and positive variations of the velocity dispersion, induced by merger episodes, are related to a fast decrease of the phase-space density Q .     

Resolving cosmic structure formation with the Millennium-II Simulation

We present the Millennium-II Simulation (MS-II), a very large N -body simulation of dark matter evolution in the concordance Λ cold dark matter (ΛCDM) cosmology. The MS-II assumes the same cosmological parameters and uses the same particle number and output data structure as the original Millennium Simulation (MS), but was carried out in a periodic cube one-fifth the size  (100  h ⁻¹ Mpc)  with five times better spatial resolution (a Plummer equivalent softening of  1.0  h ⁻¹ kpc  ) and with 125 times better mass resolution (a particle mass of  6.9 × 10⁶  h ⁻¹ M_⊙  ). By comparing results at MS and MS-II resolution, we demonstrate excellent convergence in dark matter statistics such as the halo mass function, the subhalo abundance distribution, the mass dependence of halo formation times, the linear and non-linear autocorrelations and power spectra, and halo assembly bias. Together, the two simulations provide precise results for such statistics over an unprecedented range of scales, from haloes similar to those hosting Local Group dwarf spheroidal galaxies to haloes corresponding to the richest galaxy clusters. The 'Milky Way' haloes of the Aquarius Project were selected from a lower resolution version of the MS-II and were then resimulated at much higher resolution. As a result, they are present in the MS-II along with thousands of other similar mass haloes. A comparison of their assembly histories in the MS-II and in resimulations of 1000 times better resolution shows detailed agreement over a factor of 100 in mass growth. We publicly release halo catalogues and assembly trees for the MS-II in the same format within the same archive as those already released for the MS.     

Figure rotation of dark haloes in cold dark matter simulations

We investigate the figure rotation of dark matter haloes identified in Λ cold dark matter (CDM) simulations. We find that when strict criteria are used to select suitable haloes for study, five of the 222 haloes identified in our   z = 0  simulation output undergo coherent figure rotation over a  5 h ⁻¹ Gyr  period. We discuss the effects of varying the selection criteria and find that pattern speeds for a much larger fraction of the haloes can be measured when the criteria are relaxed. Pattern speeds measured over a  1 h ⁻¹ Gyr  period follow a lognormal distribution, centred at  Ω_p= 0.2 h rad Gyr⁻¹  with a maximum value of 0.94 h rad Gyr⁻¹. Over a  5 h ⁻¹ Gyr  period, the average pattern speed of a halo is about  0.1 h rad Gyr⁻¹  and the largest pattern speed found is  0.24 h rad Gyr⁻¹  . Less than half of the selected haloes showed alignment between their figure rotation axis and minor axis, the exact fraction being somewhat dependent on how one defines a halo. While the pattern speeds observed are lower than those generally thought capable of causing spiral structure, we note that coherent figure rotation is found over very long periods and argue that further simulations would be required before strong conclusions about spiral structure in all galaxies could be drawn. We find no correlation between halo properties such as total mass and the pattern speed.     

Testing tidal-torque theory – I. Spin amplitude and direction

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N -body simulation with thousands of well-resolved haloes. The protohaloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the protohaloes match, on average, the spin amplitudes of the virialized haloes of today, if linear growth is assumed until ∼ t ₀/3, or  55–70  per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ∼50° in the TTT prediction at t ₀, such that the linear spatial correlations of spins on scales ≥1  h ⁻¹ Mpc are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second-order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.     

Constraints on dark matter physics from dwarf galaxies through galaxy cluster haloes

One of the predictions of the standard cold dark matter model is that dark haloes have centrally divergent density profiles. An extensive body of rotation curve observations of dwarf and low surface brightness galaxies shows the dark haloes of those systems to be characterized by soft constant-density central cores. Several physical processes have been proposed to produce soft cores in dark haloes, each one with different scaling properties. With the aim of discriminating among them we have examined the rotation curves of dark-matter-dominated dwarf and low surface brightness galaxies and the inner mass profiles of two clusters of galaxies lacking a central cD galaxy and with evidence of soft cores in the centre. The core radii and central densities of these haloes scale in a well-defined manner with the depth of their potential wells, as measured through the maximum circular velocity. As a result of our analysis we identify self-interacting cold dark matter as a viable solution to the core problem, where a non-singular isothermal core is formed in the halo centre surrounded by a Navarro, Frenk & White profile in the outer parts. We show that this particular physical situation predicts core radii in agreement with observations. Furthermore, using the observed scalings, we derive an expression for the minimum cross-section ( σ ) which has an explicit dependence with the halo dispersion velocity ( v ). If m _x is the mass of the dark matter particle: σ m _x ≈4×10⁻²⁵ (100 km s⁻¹  v ⁻¹) cm² GeV⁻¹.     

Dark matter subhaloes in numerical simulations

We use cosmological Λ cold dark matter (CDM) numerical simulations to model the evolution of the substructure population in 16 dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10¹¹ to  10¹⁴ M_⊙  at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation  d n /d v = (1/8)( v _cmax/ v _cmax,host)⁻⁴  . Halo to halo variance in the VDF is a factor of roughly 2 to 4. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for   r ≳ 0.3 r _vir  , where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease towards the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias towards more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power-law cosmologies of   P ∝ kⁿ   , and demonstrate that a steeper spectral index, n , results in significantly less substructure.     

The distribution of ejected subhaloes and its implication for halo assembly bias

Using a high-resolution cosmological N -body simulation, we identify the ejected population of subhaloes, which are haloes at redshift   z = 0  but were once contained in more massive 'host' haloes at high redshifts. The fraction of the ejected subhaloes in the total halo population of the same mass ranges from 9 to 4 per cent for halo masses from  ∼10¹¹  to  ∼10¹²  h ⁻¹ M_⊙  . Most of the ejected subhaloes are distributed within four times the virial radius of their hosts. These ejected subhaloes have distinct velocity distribution around their hosts in comparison to normal haloes. The number of subhaloes ejected from a host of given mass increases with the assembly redshift of the host. Ejected subhaloes in general reside in high-density regions, and have a much higher bias parameter than normal haloes of the same mass. They also have earlier assembly times, so that they contribute to the assembly bias of dark matter haloes seen in cosmological simulations. However, the assembly bias is not dominated by the ejected population, indicating that large-scale environmental effects on normal haloes are the main source for the assembly bias.     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

Scale-dependent bias induced by local non-Gaussianity: a comparison to N-body simulations

We investigate the effect of primordial non-Gaussianity of the local f _NL type on the auto- and cross-power spectra of dark matter haloes using simulations of the Λ cold dark matter cosmology. We perform a series of large N -body simulations of both positive and negative f _NL, spanning the range between 10 and 100. Theoretical models predict a scale-dependent bias correction  Δ b ( k , f _NL)  that depends on the linear halo bias   b ( M )  . We measure the power spectra for a range of halo mass and redshifts covering the relevant range of existing galaxy and quasar populations. We show that auto- and cross-correlation analyses of bias are consistent with each other. We find that for low wavenumbers with   k < 0.03  h  Mpc⁻¹  the theory and the simulations agree well with each other for biased haloes with   b ( M ) > 1.5  . We show that a scale-independent bias correction improves the comparison between theory and simulations on smaller scales, where the scale-dependent effect rapidly becomes negligible. The current limits on f _NL from Slosar et al. come mostly from very large scales   k < 0.01  h  Mpc⁻¹  and, therefore, remain valid. For the halo samples with   b ( M ) < 1.5 − 2  , we find that the scale-dependent bias from non-Gaussianity actually exceeds the theoretical predictions. Our results are consistent with the bias correction scaling linearly with f _NL.     

Galaxy–dark matter correlations applied to galaxy–galaxy lensing: predictions from the semi-analytic galaxy formation models

We use semi-analytic models of galaxy formation combined with high-resolution N -body simulations to make predictions for galaxy–dark matter correlations and apply them to galaxy–galaxy lensing. We analyse cross-power spectra between the dark matter and different galaxy samples selected by luminosity, colour or star formation rate. We compare the predictions with the recent detection by the Sloan Digital Sky Survey (SDSS). We show that the correlation amplitude and the mean tangential shear depend strongly on the luminosity of the sample on scales below 1  h ⁻¹ Mpc, reflecting the correlation between the galaxy luminosity and the halo mass. The cross-correlation cannot, however, be used to infer the halo profile directly because different halo masses dominate on different scales and because not all galaxies are at the centres of the corresponding haloes. We compute the redshift evolution of the cross-correlation amplitude and compare it with those of galaxies and dark matter. We also compute the galaxy–dark matter correlation coefficient and show that it is close to unity on scales above 1  h ⁻¹ Mpc for all considered galaxy types. This would allow one to extract the bias and the dark matter power spectrum on large scales from the galaxy and galaxy–dark matter correlations.     

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.     

Numerical simulations of hot halo gas in galaxy mergers

Galaxy merger simulations have explored the behaviour of gas within the galactic disc, yet the dynamics of hot gas within the galaxy halo have been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with metal-free hot halo gas. To isolate the effect of the halo gas, we simulate only the dark matter halo and the hot halo gas over a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that (i) a strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to   T ∼ 10^6.3 K  . (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is  ≳10³⁹ erg s⁻¹  for halo gas fractions larger than 10 per cent. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as  10⁴² erg s⁻¹  . This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) ∼10–20 per cent of the initial gas mass is unbound from the galaxies for equal-mass mergers, while 3–5 per cent of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the intergalactic medium with metals.     

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

Evaluating approximations for halo merging histories

We study the merging history of dark matter haloes in N -body simulations and semi-analytical 'merger trees' based on the extended Press–Schechter (EPS) formalism. The main focus of our study is the joint distribution of progenitor number and mass as a function of redshift and parent halo mass. We begin by investigating the mean quantities predicted directly by the Press–Schechter (PS) and EPS formalism, such as the halo mass and conditional mass functions, and compare these predictions with the results of the simulations. The higher moments of this distribution are not predicted by the EPS formalism alone and must be obtained from the merger trees. We find that the Press–Schechter model deviates from the simulations at the level of 30–50 per cent on certain mass scales, and that the sense of the discrepancy changes as a function of redshift. We show that this discrepancy is reflected in the higher moments of the distribution of progenitor mass and number. We investigate some related statistics such as the accretion rate and the mass ratio of the largest two progenitors. For galaxy sized haloes ( M ∼10¹² M_⊙), we find that the merging history of haloes, as represented by these statistics, is well reproduced in the merger trees compared with the simulations. The agreement deteriorates for larger mass haloes. We conclude that merger trees based on the extended Press–Schechter formalism provide a reasonably reliable framework for semi-analytical models of galaxy formation.     

Halo masses for optically selected and for radio-loud AGN from clustering and galaxy–galaxy lensing

We compute two-point correlation functions and measure the shear signal due to galaxy–galaxy lensing for 80 000 optically identified and 5700 radio-loud active galactic nuclei (AGN) from Data Release 4 of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed   M _*  . Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive haloes than optical AGN: in our samples, their mean halo masses are  1.6 × 10¹³  and  8 × 10¹¹  h ⁻¹ M_⊙  , respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter haloes of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in haloes where gas cools quasi-statically.     

Numerical action reconstruction of the dynamical history of dark matter haloes in N-body simulations

We test the ability of the numerical action method (NAM) to recover the individual orbit histories of mass tracers in an expanding universe, given the masses and redshift-space coordinates at the present epoch. The mass tracers are represented by dark matter (DM) haloes identified in a region of radius  26  h ⁻¹ Mpc  of a high-resolution N -body simulation of the standard Λ cold dark matter (CDM) cosmology. Since previous tests of NAM at this scale have traced the underlying distribution of DM particles rather than extended haloes, our study offers an assessment of the accuracy of NAM in a scenario which more closely approximates the complex dynamics of actual galaxy haloes. We show that NAM can recover the present-day halo distances with typical distance errors of less than 3 per cent and radial peculiar velocities with a dispersion of  ∼130 km s⁻¹  . The accuracy of individual orbit reconstructions was limited by the inability of NAM, in some instances, to correctly model the positions of haloes at early times solely on the basis of the redshifts, angular positions and masses of the haloes at the present epoch. Improvements in the quality of NAM reconstructions may be possible using the present-day three-dimensional halo velocities and distances to further constrain the dynamics. This velocity data is expected to become available for nearby galaxies in the coming generations of observations by Space Interferometry Mission ( SIM ) and Global Astrometric Interferometer for Astrophysics ( GAIA ).     

Gas in simulations of high-redshift galaxies and minihaloes

We study the gas content of haloes in the early universe using high-resolution hydrodynamical simulations. We extract from the simulations and also predict, based on linear theory, the halo mass for which the enclosed baryon fraction equals half of the mean cosmic fraction. We find a rough agreement between the simulations and the predictions, which suggests that during the high-redshift era before stellar heating, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was  ∼3 × 10⁴ M_⊙  . We also carry out a detailed resolution analysis and show that in order to determine a halo's gas fraction even to 20 per cent accuracy, the halo must be resolved into at least 500 dark matter particles.     

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 concentration–velocity dispersion relation in galaxy groups

Based on results from cold dark matter N -body simulations, we develop a dynamical model for the evolution of subhaloes within group-sized host haloes. Only subhaloes more massive than 5 × 10⁸ M_⊙ are considered, because they are massive enough to possibly host luminous galaxies. On their orbits within a growing host potential the subhaloes are subject to tidal stripping and dynamical friction. At the present time  ( z = 0)  , all model hosts have equal mass  ( M _vir= 3.9 × 10¹³ M_⊙)  but different concentrations associated with different formation times. We investigate the variation of subhalo (or satellite galaxy) velocity dispersion with host concentration and/or formation time. In agreement with the Jeans equation, the velocity dispersion of subhaloes increases with the host concentration. Between concentrations of ∼5 and ∼20, the subhalo velocity dispersions increase by a factor of ∼1.25. By applying a simplified tidal disruption criterion, that is, rejection of all subhaloes with a tidal truncation radius below 3  kpc at   z = 0  , the central velocity dispersion of the 'surviving' subhalo sample increases substantially for all concentrations. The enhanced central velocity dispersions in the surviving subhalo samples are caused by a lack of slow tangential motions. Additionally, we present a fitting formula for the anisotropy parameter which does not depend on concentration if the group-centric distances are scaled by r _s, the characteristic radius of the Navarro, Frenk & White profile. Since the expected loss of subhaloes and galaxies due to tidal disruption increases the velocity dispersion of surviving galaxies, the observed galaxy velocity dispersion can substantially overestimate the virial mass.