Simulating subhaloes at high redshift: merger rates, counts and types

Galaxies are believed to be in one-to-one correspondence with simulated dark matter subhaloes. We use high-resolution N -body simulations of cosmological volumes to calculate the statistical properties of subhalo (galaxy) major mergers at high redshift ( z = 0.6–5). We measure the evolution of the galaxy merger rate, finding that it is much shallower than the merger rate of dark matter host haloes at   z > 2.5  , but roughly parallels that of haloes at   z < 1.6  . We also track the detailed merger histories of individual galaxies and measure the likelihood of multiple mergers per halo or subhalo. We examine satellite merger statistics in detail: 15–35 per cent of all recently merged galaxies are satellites, and satellites are twice as likely as centrals to have had a recent major merger. Finally, we show how the differing evolution of the merger rates of haloes and galaxies leads to the evolution of the average satellite occupation per halo, noting that for a fixed halo mass, the satellite halo occupation peaks at   z ∼ 2.5  .     

Global environmental effects versus galaxy interactions

We explore properties of close galaxy pairs and merging systems selected from the Sloan Digital Sky Survey Data Release 4 in different environments with the aim to assess the relative importance of the role of interactions over global environmental processes. For this purpose, we perform a comparative study of galaxies with and without close companions as a function of local density and host halo mass, carefully removing sources of possible biases. We find that at low- and high-local-density environments, colours and concentration indices of close galaxy pairs are very similar to those of isolated galaxies. At intermediate densities, we detect significant differences, indicating that close pairs could have experienced a more rapid transition on to the red sequence than isolated galaxies. The presence of a correlation between concentration index and colours indicates that the physical mechanism responsible for the colour transformation also operates in the transformation of the luminous matter distribution. At fixed local densities, we find a dependence of the red galaxy fraction on dark matter halo mass for galaxies with or without a close companion. This suggests the action of host halo mass related effects. Regardless of dark matter halo mass, we show that the percentage of red galaxies in close pairs and in the control sample are comparable at low- and high-local-density environments. However, at intermediate local densities, the gap in the red fraction between close pairs and the control galaxies increases from ∼10 per cent in low-mass haloes up to ∼50 per cent in the most massive ones. Interestingly, we also detect that 50 per cent of merging systems populate the intermediate local environments, with a large fraction of them being extremely red and bulge dominated. Our findings suggest that in intermediate-density environments galaxies are efficiently pre-processed by close encounters and mergers before entering higher local density regions.     

The growth of central and satellite galaxies in cosmological smoothed particle hydrodynamics simulations

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to  7 × 10⁹ M_⊙  . Most friends-of-friends haloes in the simulation have a distinct central galaxy, typically 2–5 times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive haloes form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5–1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Typically, satellites in our simulation are 0.1–0.2 mag bluer than in models that assume no gas accretion on to satellites after a halo merger. Since   z = 1  , 27 per cent of central galaxies (above  3 × 10¹⁰ M_⊙  ) and 22 per cent of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain 'central' objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.     

Dynamical friction and galaxy merging time-scales

The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N -body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M _sat/ M _host≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M _sat/ M _host≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halo's centre for a wide range of M _sat/ M _host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≲10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.     

Merging history trees for dark matter haloes: tests of the Merging Cell Model in a CDM cosmology

The merging history of dark matter haloes is computed with the Merging Cell Model proposed by Rodrigues & Thomas. While originally discussed in the case of scale-free power spectra, it is developed and tested here in the framework of the cold dark matter cosmology. The halo mass function, the mass distribution of progenitors and child haloes, as well as the probability distribution of formation times, have been computed and compared with the available analytic predictions. The halo autocorrelation function has also been obtained (a first for a semi-analytic merging tree), and tested against analytic formulae. An overall good agreement is found between results of the model, and the predictions derived from the Press & Schechter theory and its extensions. More severe discrepancies appear when formulae that better describe N -body simulations are used for comparison. In many instances, the model can be a useful tool for following the hierarchical growth of structures. In particular, it is suitable for addressing the issue of the formation and evolution of galaxy clusters, as well as the population of Lyman-break galaxies at high redshift, and their clustering properties.     

Galactic halo cusp–core: tidal compression in mergers

We explain in simple terms how the build-up of dark haloes by merging compact satellites, as in the cold dark matter (CDM) cosmology, inevitably leads to an inner cusp of density profile  ρ∝ r ^−α  with  α≳ 1  , as seen in cosmological N -body simulations. A flatter halo core with  α < 1  exerts on the satellites tidal compression in all directions, which prevents the deposit of stripped satellite material in the core region. This makes the satellite orbits decay from the radius where  α∼ 1  to the halo centre with no local tidal mass transfer, and thus causes a rapid steepening of the inner profile to  α > 1  . These tidal effects, the resultant steepening of the profile to a cusp, and the stability of this cusp to tandem mergers with compact satellites are demonstrated using N -body simulations. The transition at  α∼ 1  is then addressed using toy models in the limiting cases of impulse and adiabatic approximations and using tidal radii for satellites on radial and circular orbits. In an associated paper, we address the subsequent slow convergence from either side to an asymptotic stable cusp with  α≳ 1  . Our analysis thus implies that an inner cusp is enforced when small haloes are typically more compact than larger haloes, as in the CDM scenario, such that enough satellite material makes it intact into the inner halo and is deposited there. We conclude that a necessary condition for maintaining a flat core, as indicated by observations, is that the inner regions of the CDM satellite haloes be puffed up by about 50 per cent such that when they merge into a larger halo they would be disrupted outside the halo core. This puffing up could be due to baryonic feedback processes in small haloes, which may be stimulated by the tidal compression in the halo cores.     

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.     

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.     

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.     

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.     

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.     

Measuring the distribution of galaxies between haloes

We develop a method to measure the probability, P ( N;   M ), of finding N galaxies in a dark matter halo of mass M from the theoretically determined clustering properties of dark matter haloes and the observationally measured clustering properties of galaxies. Knowledge of this function and the distribution of the dark matter completely specifies all clustering properties of galaxies on scales larger than the size of dark matter haloes. Furthermore, P ( N;   M ) provides strong constraints on models of galaxy formation, since it depends upon the merger history of dark matter haloes and the galaxy–galaxy merger rate within haloes. We show that measurements from a combination of the Two Micron All Sky Survey and Sloan Digital Sky Survey or Two-degree Field Galaxy Redshift Survey data sets will allow P ( N;   M ) averaged over haloes occupied by bright galaxies to be accurately measured for N =0–2 .     

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.     

The formation and survival of discs in a ΛCDM universe

We study the formation of galaxies in a Λ cold dark matter (ΛCDM) universe using high-resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight isolated haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at   z = 0  , none of our galaxies contains a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The   z = 0  spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at   z = 0  nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at   z ≳ 2  , regardless of their   z = 0  morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic ΛCDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.     

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.     

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.     

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

The importance of interactions for mass loss from satellite galaxies in cold dark matter haloes

We investigate the importance of interactions between dark matter substructures for the mass loss they suffer whilst orbiting within a sample of high-resolution galaxy cluster mass cold dark matter (CDM) haloes formed in cosmological N -body simulations. We have defined a quantitative measure that gauges the degree to which interactions are responsible for mass loss from substructures. This measure indicates that interactions are more prominent in younger systems when compared to older more relaxed systems. We show that this is due to the increased number of encounters a satellite experiences and a higher mass fraction in satellites. This is in spite of the uniformity in the distributions of relative distances and velocities of encounters between substructures within the different host systems in our sample.
Using a simple model to relate the net force felt by a single satellite to the mass loss it suffers, we show that interactions with other satellites account for ∼30 per cent of the total mass loss experienced over its lifetime. The relation between the age of the host and the importance of interactions increases the scatter about this mean value from ∼25 per cent for the oldest to ∼45 per cent for the youngest system we have studied. We conclude that satellite interactions play a vital role in the evolution of substructure in dark matter haloes and that a significant fraction of the tidally stripped material can be attributed to these interactions.     

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.