Encounters between spherical galaxies – I. Systems without a dark halo

We report here on a survey of N -body simulations of encounters between spherical galaxies. Initial systems are isotropic Jaffe models. Different sets of mass ratios, impact parameters and orbital energies are studied. Both merger remnants and systems perturbed after a non-merging encounter are analysed and compared to real-life elliptical galaxies. The properties of merger remnants show a large variety. Merger remnants resulting from head-on encounters are mainly non-rotating prolate spheroids. Merger remnants from models with   J _orb≠ 0  are tri-axial or mildly oblate spheroids, supported in part by rotation. The velocity distributions are biased towards the radial direction in the prolate case and the tangential direction in the oblate case. Non-mergers are affected in various ways, depending on the orbital characteristics. We conclude that many of the global properties of real-life ellipticals can, in principle, be attributed to a merger of spherical progenitors.     

The elliptical galaxy colour–magnitude relation as a discriminant between the monolithic and merger paradigms

The colour–magnitude relation (CMR) of cluster elliptical galaxies has been widely used to constrain their star formation histories (SFHs) and to discriminate between the monolithic collapse and merger paradigms of elliptical galaxy formation. We use a Λ cold dark matter hierarchical merger model of galaxy formation to investigate the existence and redshift evolution of the elliptical galaxy CMR in the merger paradigm. We show that the SFH of cluster ellipticals predicted by the model is quasi-monolithic , with only ∼10 per cent of the total stellar mass forming after   z ∼ 1  . The quasi-monolithic SFH results in a predicted CMR that agrees well with its observed counterpart in the redshift range  0 < z < 1.27  . We use our analysis to argue that the elliptical-only CMR can be used to constrain the SFHs of present-day cluster ellipticals only if we believe a priori in the monolithic collapse model. It is not a meaningful tool for constraining the SFH in the merger paradigm, since a progressively larger fraction of the progenitor set of present-day cluster ellipticals is contained in late-type star-forming systems at higher redshift, which cannot be ignored when deriving the SFHs. Hence, the elliptical-only CMR is not a useful discriminant between the two competing theories of elliptical galaxy evolution.     

Evidence for merger-driven activity in the clustering of high-redshift quasars

Recently, a very large clustering length has been measured for quasars at a redshift of   z ∼ 4  . In combination with the observed quasar luminosity function, we assess the implications of this clustering for the relationship between quasar luminosity and dark matter halo mass. Our analysis allows for non-linearity and finite scatter in the relation between quasar luminosity and halo mass, as well as a luminosity dependent quasar lifetime. The additional novel ingredient in our modelling is the allowance for an excess in the observed bias over the underlying halo bias owing to the merger driven nature of quasar activity. We find that the observations of clustering and luminosity function can be explained only if both of the following conditions hold: (i) the luminosity to halo mass ratio increases with halo mass; (ii) the observed clustering amplitude is in excess of that expected solely from halo bias. The latter result is statistically significant at the 99 per cent level. Taken together, the observations provide compelling evidence for merger driven quasar activity, with a black-hole growth that is limited by feedback. In difference from previous analyses, we show that there could be scatter in the luminosity–halo mass relation of up to 1 dex, and that quasar clustering cannot be used to estimate the quasar lifetime.     

HST/WFPC2 imaging of the QDOT ultraluminous infrared galaxy sample

We present HST WFPC2 V -band imaging for 23 ultraluminous infrared galaxies (ULIRGs) taken from the QDOT redshift survey. The fraction of sources observed to be interacting is 87 per cent. Most of the merging systems show a number of compact 'knots', whose colour and brightness differ substantially from their immediate surroundings. Colour maps for nine of the objects show a non-uniform colour structure. Features include blue regions located towards the centres of merging systems which are likely to be areas of enhanced star formation, and compact red regions which are likely to be dust shrouded starbursts or active galactic nuclei. The host galaxies of the quasi-stellar objects (QSOs) in the sample were found to be either interacting systems or ellipticals. Our data show no evidence that ULIRGs are a simple transition stage between galaxy mergers and QSOs. We propose an alternative model for ULIRGs based on the morphologies in our sample and previous N -body simulations. Under this model ULIRGs as a class are much more diverse than a simple transition between galaxy merger and QSO. The evolution of IR power source and merger morphology in ULIRGs is driven solely by the local environment and the morphologies of the merger progenitors.     

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  .     

Elliptical galaxies from mergers of discs

We analyse N -body galaxy merger experiments involving disc galaxies. Mergers of disc–bulge–halo models are compared to those of bulgeless, disc–halo models to quantify the effects of the central bulge on merger dynamics and the structure of the remnant. Our models explore galaxy mass ratios 1:1 through 3:1, and use higher bulge mass fractions than previous studies. A full comparison of the structural and dynamical properties with our observations is carried out. The presence of central bulges results in longer tidal tails, oblate final intrinsic shapes, surface brightness profiles with a higher Sérsic index, steeper rotation curves and oblate-rotator internal dynamics. Mergers of bulgeless galaxies do not generate long-lasting tidal tails, and their strong triaxiality seems inconsistent with observations; these remnants show shells, which we do not find in models including central bulges. Giant ellipticals with boxy isophotes and anisotropic dynamics cannot be produced by the mergers modelled here; they could be the result of mergers between lower luminosity ellipticals, themselves plausibly formed in disc-disc mergers.     

On the physical connections between galaxies of different types

Galaxies can be classified in two broad sequences which are likely to reflect their formation mechanism. The 'main sequence', consisting of spirals, irregulars and all dwarf galaxies, is probably produced by gas settling within dark matter haloes. We show that the sizes and surface densities along this sequence are primarily determined by the distributions of the angular momentum and formation time of dark haloes. They are well reproduced by current cosmogonies provided that galaxies form late, at z  ≲ 2. In this scenario, dwarf ellipticals were small 'discs' at z  ∼ 1 and become 'ellipticals' after they fall into cluster environments. The strong clustering of dwarf ellipticals is then a natural by-product of the merging and transformation process. The number of dwarf galaxies predicted in a cluster such as Virgo is in good agreement with the observed number. On the other hand, the 'giant branch', consisting of giant ellipticals and bulges, is probably produced by the merging of disc galaxies. Based on the observed phase-space densities of galaxies, we show that the main bodies of all giant ellipticals can be produced by dissipationless mergers of high-redshift discs. However, high-redshift discs, although denser than present-day ones, are still not compact enough to produce the high central phase-space density of some low-luminosity ellipticals. Dissipation must have occurred in the central parts of these galaxies during the merger which formed them.     

An improved model for contraction of dark matter haloes in response to condensation of baryons

The cooling of gas in the centres of dark matter haloes is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. Following Gnedin et al., we improve this model by modifying the assumed invariant from M ( r ) r to     , where r and     are the current and orbit-averaged particle positions. We explore the effect of the bulge in the inner regions of the halo for different values of the bulge-to-disc mass ratio. We find that the bulge makes the velocity curve rise faster in the inner regions of the halo. We present an analytical fitting curve that describes the velocity curve of the halo after dissipation. The results should be useful for dark matter detection studies.     

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.     

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

Constraints on the merger models of elliptical galaxies from their globular cluster systems

The discovery of protoglobular cluster candidates in many present-day mergers allows us to understand better the possible effects of a merger event on the globular cluster system of a galaxy, and to foresee the properties of the end-product. By comparing these expectations with the properties of globular cluster systems of elliptical galaxies at the present time we can constrain merger models. The observational data indicate that (i) every gaseous merger induces the formation of new star clusters, and (ii) the number of new clusters formed in such a merger increases with the gas content of the progenitor galaxies. Low-luminosity (about M _V  > −21), discy ellipticals are generally thought to be the result of a gaseous merger. As such, new globular clusters are expected to form but have not been detected to date. We investigate various reasons for the non-detection of subpopulations in low-luminosity ellipticals, i.e. absence of an old population, absence of a new population, destruction of one of the populations and, finally, an age–metallicity conspiracy that allows old and new globular clusters to appear indistinguishable at the present epoch. All of these possibilities lead us to a similar conclusion, namely that low-luminosity ellipticals did not form recently ( z  < 1) in a gas-rich merger, and might not have formed in a major merger of stellar systems at all. High-luminosity ellipticals do reveal globular cluster subpopulations. However, it is difficult to account for the two populations in terms of mergers alone and, in particular, we can rule out scenarios in which the second subpopulation is the product of a recent, gas-poor merger.     

Constraining cold dark matter halo merger rates using the coagulation equations

We place additional constraints on the three parameters of the dark matter halo merger rate function recently proposed by Parkinson, Cole & Helly by utilizing Smoluchowski's coagulation equation, which must be obeyed by any binary merging process which conserves mass. We find that the constraints from Smoluchowski's equation are degenerate, limiting to a thin plane in the three-dimensional parameter space. This constraint is consistent with those obtained from fitting to N -body measures of progenitor mass functions, and provides a better match to the evolution of the overall dark matter halo mass function, particularly for the most massive haloes. We demonstrate that the proposed merger rate function does not permit an exact solution of Smoluchowski's equation and, therefore, the choice of parameters must reflect a compromise between fitting various parts of the mass function. The techniques described herein are applicable to more general merger rate functions, which may permit a more accurate solution of Smoluchowski's equation. The current merger rate solutions are most probably sufficiently accurate for the vast majority of applications.     

Environmental dependence of dark matter halo growth – I. Halo merger rates

In an earlier paper, we quantified the mean merger rate of dark matter haloes as a function of redshift z , descendant halo mass M ₀, and progenitor halo mass ratio ξ using the Millennium simulation of the Λ cold dark matter cosmology. Here, we broaden that study and investigate the dependence of the merger rate of haloes on their surrounding environment. A number of local mass overdensity variables, both including and excluding the halo mass itself, are tested as measures of a halo's environment. The simple functional dependence on   z , M ₀  , and ξ of the merger rate found in our earlier work, is largely preserved in different environments, but we find that the overall amplitude of the merger rate has a strong positive correlation with the environmental densities. For galaxy-mass haloes, we find mergers to occur ∼2.5 times more frequently in the densest regions than in voids at both   z = 0  and higher redshifts. Higher mass haloes show similar trends. We present a fitting form for this environmental dependence that is a function of both mass and local density and valid out to   z = 2  . The amplitude of the progenitor (or conditional) mass function shows a similar correlation with local overdensity, suggesting that the extended Press–Schechter model for halo growth needs to be modified to incorporate environmental effects.     

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.     

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 .     

Separation of the visible and dark matter in the Einstein ring LBG J213512.73−010143

We model the mass distribution in the recently discovered Einstein ring LBG J213512.73−010143 (the 'Cosmic Eye') using archival Hubble Space Telescope imaging. We reconstruct the mass density profile of the z = 0.73 lens and the surface brightness distribution of the z = 3.07 source and find that the observed ring is best fitted with a dual-component lens model consisting of a baryonic Sersic component nested within a dark matter halo. The dark matter halo has an inner slope of 1.42^+0.24_−0.22, consistent with cold dark matter simulations after allowing for baryon contraction. The baryonic component has a mass-to-light ratio of  1.71^+0.28_−0.38 M_⊙/L _{B ⊙}  which when evolved to the present day is in agreement with local ellipticals. Within the Einstein radius of 0.77 arcsec (5.6 kpc), the baryons account for 46 ± 11 per cent of the projected lens mass. External shear from a nearby foreground cluster is accurately predicted by the model. The reconstructed surface brightness distribution in the source plane clearly shows two peaks. Through a generalization of our lens inversion method, we conclude that the redshifts of both peaks are consistent with each other, suggesting that we are seeing structure within a single galaxy.