Dissipationless mergers of elliptical galaxies and the evolution of the fundamental plane

We carry out numerical simulations of dissipationless major mergers of elliptical galaxies using initial galaxy models that consist of a dark matter haloes and a stellar bulge with properties consistent with the observed fundamental plane. By varying the density profile of the dark matter haloes [standard Navarro, Frenk & White (NFW) profile versus adiabatically contracted NFW profile], the global stellar to dark matter mass ratio and the orbit of the merging galaxies, we are able to assess the impact of each of these factors on the structure of the merger remnant. Our results indicate that the properties of the remnant bulge depend primarily on the angular momentum and energy of the orbit; for a cosmologically motivated orbit, the effective radius and velocity dispersion of the remnant bulge remain approximately on the fundamental plane. This indicates that the observed properties of elliptical galaxies are consistent with significant growth via late dissipationless mergers. We also find that the dark matter fraction within the effective radius of our remnants increases after the merger, consistent with the hypothesis that the tilt of the fundamental plane from the virial theorem is due to a varying dark matter fraction as a function of galaxy mass.     

Dark matter halo response to the disc growth

We consider the sensitivity of the circular-orbit adiabatic contraction approximation to the baryon condensation rate and the orbital structure of dark matter haloes in the Λ cold dark matter (ΛCDM) paradigm. Using one-dimensional hydrodynamic simulations including the dark matter halo mass accretion history and gas cooling, we demonstrate that the adiabatic approximation is approximately valid even though haloes and discs may assemble simultaneously. We further demonstrate the validity of the simple approximation for ΛCDM haloes with isotropic velocity distributions using three-dimensional N -body simulations. This result is easily understood: an isotropic velocity distribution in a cuspy halo requires more circular orbits than radial orbits. Conversely, the approximation is poor in the extreme case of a radial orbit halo. It overestimates the response of a core dark matter halo, where radial orbit fraction is larger. Because no astronomically relevant models are dominated by low angular momentum orbits in the vicinity of the disc and the growth time-scale is never shorter than a dynamical time, we conclude that the adiabatic contraction approximation is useful in modelling the response of dark matter haloes to the growth of a disc.     

Mass distribution in nearby Abell clusters

We study the mass distribution in six nearby  ( z < 0.06)  relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters, we approximate their total density distribution by the Navarro, Frenk & White (NFW) formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120–420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from the N -body simulations of dark matter haloes. We find that the best-fitting Galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters, we conclude that the trend of decreasing concentration for higher masses found in the cosmological N -body simulations is consistent with the data. By scaling the individual cluster data by mass, we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range  1.5 < c < 14  and the anisotropy parameter in the range  −1.1 < β < 0.5  at the 95 per cent confidence level.     

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  .     

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.     

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.     

Testing the modified Press–Schechter model against N-body simulations

A modified version of the extended Press–Schechter model for the growth of dark-matter haloes was introduced in two previous papers, with the aim of explaining the mass–density relation shown by haloes in high-resolution cosmological simulations. In this model, major mergers are well separated from accretion, thereby allowing a natural definition of halo formation and destruction. This makes it possible to derive analytic expressions for halo formation and destruction rates, the mass accretion rate and the probability distribution functions of halo formation times and progenitor masses. The stochastic merger histories of haloes can be readily derived and easily incorporated into semi-analytical models of galaxy formation, thus avoiding the usual problems encountered in the construction of Monte Carlo merger trees from the original extended Press–Schechter formalism. Here we show that the predictions of the modified Press–Schechter model are in good agreement with the results of N -body simulations for several scale-free cosmologies.     

Modelling angular-momentum history in dark-matter haloes

We model the acquisition of spin by dark-matter haloes in semi-analytic merger trees. We explore two different algorithms: one in which halo spin is acquired from the orbital angular momentum of merging satellites, and another in which halo spin is gained via tidal torquing on shells of material while still in the linear regime. We find that both scenarios produce the characteristic spin distribution of haloes found in N -body simulations, namely, a log-normal distribution with mean ≈ 0.04 and standard deviation ≈ 0.5 in the log. A perfect match requires fine-tuning of two free parameters. Both algorithms also reproduce the general insensitivity of the spin distribution to halo mass, redshift and cosmology seen in N -body simulations. The spin distribution can be made strictly constant by physically motivated scalings of the free parameters. In addition, both schemes predict that haloes that have had recent major mergers have systematically larger spin values. These algorithms can be implemented within semi-analytic models of galaxy formation based on merger trees. They yield detailed predictions of galaxy properties that strongly depend on angular momentum (such as size and surface brightness) as a function of merger history and environment.     

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.     

The Planetary Nebula Spectrograph elliptical galaxy survey: the dark matter in NGC 4494

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 planetary nebulae out to seven effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the planetary nebulae agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside  1 R _e  . The velocity dispersion profile declines with radius, though not very steeply, down to  ∼70 km s⁻¹  at the last data point.
We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component Λ cold dark matter (CDM) motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fitting solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model.
Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration haloes, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.     

Collisional baryonic dark matter haloes

If dark haloes are composed of dense gas clouds, as has recently been inferred, then collisions between clouds lead to galaxy evolution. Collisions introduce a core in an initially singular dark matter distribution, and can thus help to reconcile scale-free initial conditions – such as are found in simulations – with observed haloes, which have cores. A pseudo-Tully–Fisher relation, between halo circular speed and visible mass (not luminosity), emerges naturally from the model: M _vis∝ V ^7/2.
Published data conform astonishingly well to this theoretical prediction. For our sample of galaxies, the mass–velocity relationship has much less scatter than the Tully–Fisher relation, and holds as well for dwarf galaxies (where diffuse gas makes a sizeable contribution to the total visible mass) as it does for giants. It seems very likely that this visible-mass/velocity relationship is the underlying physical basis for the Tully–Fisher relation, and this discovery in turn suggests that the dark matter is both baryonic and collisional.     

Mass profiles and anisotropies of early-type galaxies

We discuss the problem of using stellar kinematics of early-type galaxies to constrain the orbital anisotropies and radial mass profiles of galaxies. We demonstrate that compressing the light distribution of a galaxy along the line of sight produces approximately the same signature in the line-of-sight velocity profiles as radial anisotropy. In particular, fitting spherically symmetric dynamical models to apparently round, isotropic face-on flattened galaxies leads to a spurious bias towards radial orbits in the models, especially if the galaxy has a weak face-on stellar disc. Such face-on stellar discs could plausibly be the cause of the radial anisotropy found in spherical models of intermediate luminosity ellipticals such as NGC 2434, 3379 and 6703.
In the light of this result, we use simple dynamical models to constrain the outer mass profiles of a sample of 18 round, early-type galaxies. The galaxies follow a Tully–Fisher relation parallel to that for spiral galaxies, but fainter by at least 0.8 mag ( I -band) for a given mass. The most luminous galaxies show clear evidence for the presence of a massive dark halo, but the case for dark haloes in fainter galaxies is more ambiguous. We discuss the observations that would be required to resolve this ambiguity.     

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.     

The image separation distribution of strong lenses: halo versus subhalo populations

We present a halo model prediction of the image separation distribution of strong lenses. Our model takes into account the subhalo population, which has been ignored in previous studies, as well as the conventional halo population. Haloes and subhaloes are linked to central and satellite galaxies by adopting a universal scaling relation between masses of (sub)haloes and luminosities of galaxies. Our model predicts that 10–20 per cent of lenses should be caused by the subhalo population. The fraction of lensing by satellite galaxies (subhaloes) peaks at ∼1 arcsec and decreases rapidly with increasing image separations. We compute fractions of lenses which lie in groups and clusters and find them to be ∼14 and ∼4 per cent, respectively; nearly half of such lenses are expected to be produced by satellite galaxies, rather than central parts of haloes. We also study mass distributions of lensing haloes and find that, even at image separations of ∼3 arcsec, the deviation of lens mass distributions from isothermal profiles is large; at or beyond ∼3 arcsec, image separations are enhanced significantly by surrounding haloes. Our model prediction agrees reasonably well with observed image separation distributions from galaxy to cluster scales.     

Mass-to-light ratio gradients in early-type galaxy haloes

Owing to the fact that the near future should see a rapidly expanding set of probes of the halo masses of individual early-type galaxies, we introduce a convenient parameter for characterizing the halo masses from both observational and theoretical results:  ∇_ℓϒ  , the logarithmic radial gradient of the mass-to-light ratio. Using halo density profiles from Λ-cold dark matter (CDM) simulations, we derive predictions for this gradient for various galaxy luminosities and star formation efficiencies  ε_SF  . As a pilot study, we assemble the available  ∇_ℓϒ  data from kinematics in early-type galaxies – representing the first unbiased study of halo masses in a wide range of early-type galaxy luminosities – and find a correlation between luminosity and  ∇_ℓϒ  , such that the brightest galaxies appear the most dark-matter dominated. We find that the gradients in most of the brightest galaxies may fit in well with the ΛCDM predictions, but that there is also a population of fainter galaxies whose gradients are so low as to imply an unreasonably high star formation efficiency  ε_SF > 1  . This difficulty is eased if dark haloes are not assumed to have the standard ΛCDM profiles, but lower central concentrations.     

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.     

Extended Press–Schechter theory and the density profiles of dark matter haloes

An inside–out model for the formation of haloes in a hierarchical clustering scenario is studied. The method combines the picture of the spherical infall model and a modification of the extended Press–Schechter theory. The mass accretion rate of a halo is defined to be the rate of its mass increase due to minor mergers. The accreted mass is deposited at the outer shells without changing the density profile of the halo inside its current virial radius. We applied the method to a flat Λ-cold dark matter universe. The resulting density profiles are compared with analytical models proposed in the literature, and a very good agreement is found. A trend is found of the inner density profile to become steeper for larger halo mass, which also results from recent N -body simulations. Additionally, present-day concentrations as well as their time evolution are derived and it is shown that they reproduce the results of large cosmological N -body simulations.     

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.     

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.