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.     

Survival of substructure within dark matter haloes

Using high-resolution cosmological N -body simulations, we investigate the survival of dark matter satellites falling into larger haloes. Satellites preserve their identity for some time after merging. We compute their loss of mass, energy and angular momentum as they are dissolved by dynamical friction, tidal forces and collisions with other satellites. We also analyse the evolution of their internal structure. Satellites with less than a few per cent of the mass of the main halo may survive for several billion years, whereas larger satellites rapidly sink into the centre of the main halo potential well and lose their identity. Penetrating encounters between satellites are frequent and may lead to significant mass loss and disruption. Only a minor fraction of cluster mass (10–15 per cent on average) is bound to substructure at most redshifts of interest. We discuss the application of these results to the survival and extent of dark matter haloes associated with galaxies in clusters, and to their interactions. We find that a minor fraction of galaxy-size dark matter haloes are disrupted by redshift z  = 0. The fraction of satellites undergoing close encounters is similar to the observed fraction of interacting or merging galaxies in clusters at moderate redshift.     

On the dynamics of the satellite galaxies in NGC 5044

The NGC 5044 galaxy group is dominated by a luminous elliptical galaxy that is surrounded by ∼160 dwarf satellites. The projected number density profile of this dwarf population deviates within ∼1/3 of the virial radius from a projected Navarro, Frenk and White (NFW) profile, which is assumed to approximate the underlying total matter distribution. By means of a semi-analytic model, we demonstrate that the interplay between gravitation, dynamical friction and tidal mass loss and destruction can explain the observed number density profile. We use only two parameters in our models: the total to stellar mass fraction of the satellite haloes and the disruption efficiency. The disruption efficiency is expressed by a minimum radius. If the tidal radius of a galaxy (halo) falls below this radius, it is assumed to become unobservable. The preferred parameters are an initial total to stellar mass fraction of ∼20 and a disruption radius of  4 kpc  . In that model, about 20 per cent of all the satellites are totally disrupted on their orbits within the group environment. Dynamical friction is less important in shaping the inner slope of the number density profile because the reduction in mass by tidal forces lowers the impact of the friction term. The main destruction mechanism is tide. In the preferred model, the total B -band luminosity of all disrupted galaxies is about twice the observed luminosity of the central elliptical galaxy, indicating that a significant fraction of stars are scattered into the intragroup medium. Dwarf galaxy satellites closer to the centre of the NGC 5044 group may exhibit optical evidence of partial tidal disruption. If dynamical friction forces the satellite to merge with the central elliptical, the angular momentum of the satellite tends to be removed at the apocentre passage. Afterwards, the satellite drops radially towards the centre.     

Infall of substructures on to a Milky Way-like dark halo

We analyse the dynamical properties of substructures in a high-resolution dark matter simulation of the formation of a Milky Way-like halo in a Λcold dark matter cosmology. Our goal is to shed light on the dynamical peculiarities of the Milky Way satellites. Our simulations show that about one-third of the subhaloes have been accreted in groups. We quantify this clustering by measuring the alignment of the angular momentum of subhaloes in a group. We find that this signal is visible even for objects accreted up to z ∼ 1, i.e. 8 Gyr ago, and long after the spatial coherence of the groups has been lost due the host tidal field. This group infall may well explain the ghostly streams proposed by Lynden-Bell & Lynden-Bell to orbit the Milky Way. Our analyses also show that if most satellites originate in a few groups, the disc-like distribution of the Milky Way satellites would be almost inevitable. This non-random assignment of satellites to subhaloes implies an environmental dependence on whether these low-mass objects are able to form stars, possibly related to the nature of reionization in the early Universe. With this picture, both the 'ghostly streams' and the 'disc-like configuration' are manifestations of the same phenomenon: the hierarchical growth of structure down to the smallest scales.     

Merging Timescales and Merger Rates of Star Clusters in Dense Star Cluster Complexes

Interacting galaxies like the famous Antennae (NGC 4038/4039) or Stephan's Quintet (HCG 92) show considerable star forming activity in their tidal arms. High resolution images (e.g. from HST-observations) indicate that these regions consist of up to hundreds of massive stellar clusters or tidal dwarf galaxies (TDG). In this paper we want to investigate the future fate of these clusters of massive star clusters (in this work called super-clusters). We simulate compact super-clusters in the tidal field of a host-galaxy and investigate the influence of orbital and internal parameters on the rate and timescale of the merging process. We show that it is possible that such configurations merge and build a dwarf galaxy, which could be an important mechanism of how long-lived dwarf satellite galaxies form. A detailed study of the merger object will appear in a follow-up paper.This revised version was published online in October 2005 with corrections to the Cover Date.     

The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

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.     

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.     

Sinking, tidally stripped, galactic satellites

Tidally interacting galaxies offer an interesting field for the investigation of chaotic phenomena in stellar systems. When the galaxies are gravitationally bound, and one of them is much larger than the other, the latter can be regarded as a satellite of the former. The study of their dynamics is somewhat simplified in this case, which presents well observed examples in nature (e.g., globular clusters). Galactic satellites suffer orbital decay due to dynamical friction, a process that may be greatly enhanced in the presence of chaotic motions. Besides, the satellite is stripped by the field of tidal forces and, in the long run, it will disintegrate completely. Modern observations are able to show the signature of these processes taking place at present.     

Investigating the Andromeda stream – III. A young shell system in M31

Published maps of red giant stars in the halo region of M31 exhibit a giant stellar stream to the south of this galaxy, as well as a giant 'shelf' to the northeast of M31's centre. Using these maps, we find that there is a fainter shelf of comparable size on the western side as well. By choosing appropriate structural and orbital parameters for an accreting dwarf satellite within the accurate M31 potential model of Geehan et al., we produce a very similar structure in an N -body simulation. In this scenario, the tidal stream produced at pericentre of the satellite's orbit matches the observed southern stream, while the forward continuation of this tidal stream makes up two orbital loops, broadened into fan-like structures by successive pericentric passages; these loops correspond to the north-eastern and western shelves. The tidal debris from the satellite also reproduces a previously observed 'stream' of counterrotating planetary nebulae and a related stream seen in red giant stars. The debris pattern in our simulation resembles the shell systems detected around many elliptical galaxies, though this is the first identification of a shell system in a spiral galaxy and the first in any galaxy close enough to allow measurements of stellar velocities and relative distances. We discuss the physics of these partial shells, highlighting the role played by spatial and velocity caustics in the observations. We show that kinematic surveys of the tidal debris will provide a sensitive measurement of M31's halo potential, while quantifying the surface density of debris in the shelves will let us reconstruct the original mass and time of disruption of the progenitor satellite.     

Orbital parameters of infalling dark matter substructures

We present distributions of the orbital parameters of dark matter substructures at the time of merging into their host halo. Accurate knowledge of the orbits of dark matter substructures is a crucial input to studies which aim to assess the effects of the cluster environment on galaxies, the heating of galaxy discs and many other topics. Orbits are measured for satellites in a large number of N -body simulations. We focus on the distribution of radial and tangential velocities, but consider also distributions of orbital eccentricity and semimajor axis. We show that the distribution of radial and tangential velocities has a simple form and provide a fitting formula for this distribution. We also search for possible correlations between the infall directions of pairs of satellites, finding evidence for positive correlation at small angular separations as expected if some infall occurs along filaments. We also find (weak) evidence for correlations between the direction of the infall and infall velocity and the spin of the host halo.     

Active galactic nuclei and the minor merger hypothesis

We have investigated the dynamics of the merging process in the minor merger hypothesis for active galactic nuclei. We find that for a satellite galaxy to be able to merge directly with the nucleus of the host galaxy (for example, to give rise to the compact dust discs which are seen in early-type active galaxies), the initial orbit of the satellite is required to be well aimed. For the case of the host galaxy being a disc galaxy, if the initial orbits of the satellites are randomly oriented with respect to the host galaxy, then the orbits of those which reach the host nuclear regions in a reasonable time are also fairly randomly oriented once they reach the nucleus. We note that this result might be able to provide an explanation of why the jet directions in the nuclei of Seyfert galaxies are apparently unrelated to the plane of the galaxy discs.     

H i observations of interacting galaxy pair NGC 4038/9

We present the results of new radio interferometer H  i line observations for the merging galaxy pair NGC 4038/9 ('The Antennae'), obtained using the Australia Telescope Compact Array. The results improve substantially with respect to those of van der Hulst and show in detail the two merging galactic discs and the two tidal tails produced by their interaction. The small edge-on spiral dwarf galaxy ESO 572–G045 is also seen near the tip of the southern tail, but distinct from it. It shows no signs of tidal interaction. The northern tidal tail of the Antennae shows no H  i connection to the discs and has an extension towards the west. The southern tidal tail is continuous, with a prominent H  i concentration at its tip, roughly at the location of the tidal dwarf galaxy observed optically by Mirabel, Dottori & Lutz. Clear velocity structure is seen along the tidal tails and in the galactic discs. Radio continuum images at 20 and 13 cm are also presented, showing the discs in detail.     

Equipotential surface and normal gravity of the Galilean satellites

In this article, expanded equations of normal gravity on the equipotential surface are proposed for a natural satellite whose orbital plane is close to its equatorial plane. Tidal effects on the normal gravity are also discussed. The authors apply these to the Galilean satellites. Calculations suggest that the tides raised by Jupiter weakly affect the Galilean satellites. The radial displacements of the gravity due to the tides are in the range between 10⁻³ and 10⁻⁵ m s⁻², which are similar to the latitudinal and longitudinal displacements. The variations along the latitude circle are larger than those along the longitude circle. We conclude that the tidal effects on most of the Galilean satellites are larger than those on the Moon.     

Tidal interactions in binary stellar systems

The tidal interactions in binary stellar systems are studied under the assumption that the orbital motion of the binary is negligible in comparison with the stellar motion. By integrating over time the tidal forces acting on the stars, the energy changes are derived. These are used to obtain simple analytical expressions for the rates of disruption and merging. This method gives appropriate value for the Roche density ρ_r and it is found that the disruption rate of a satellite of density ρ changes drastically at ρ ρ_R A comparison is made with earlier results obtained under the simplifying assumption that stellar motion is negligible in comparison with the orbital motion of the binary and its implications are discussed.     

Formation of gas discs in merging galaxies

Observations indicate that much of the interstellar gas in merging galaxies may settle into extended gaseous discs. Here, I present simulations of disc formation in mergers of gas-rich galaxies. Up to half of the total gas settles into embedded discs; the most massive instances result from encounters in which both galaxies are inclined to the orbital plane. These discs are often warped, many have rather complex kinematics, and roughly a quarter have counter-rotating or otherwise decoupled central components. Discs typically grow from the inside out; infall from tidal tails may continue disc formation over long periods of time.     

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.     

The dynamical evolution of substructure

The evolution of substructure embedded in non-dissipative dark haloes is studied through N -body simulations of isolated systems, both in and out of initial equilibrium, complementing cosmological simulations of the growth of structure. We determine by both analytic calculations and direct analysis of the N -body simulations the relative importance of various dynamical processes acting on the clumps, such as the removal of material by global tides, clump–clump heating, clump–clump merging and dynamical friction. The ratio of the internal clump velocity dispersion to that of the dark halo is an important parameter; as this ratio approaches a value of unity, heating by close encounters between clumps becomes less important, while the other dynamical processes continue to increase in importance. Our comparison between merging and disruption processes implies that spiral galaxies cannot be formed in a protosystem that contains a few large clumps, but can be formed through the accretion of many small clumps; elliptical galaxies form in a more clumpy environment than do spiral galaxies. Our results support the idea that the central cusp in the density profiles of dark haloes is the consequence of self-limiting merging of small, dense haloes. This implies that the collapse of a system of clumps/substructure is not sufficient to form a cD galaxy, with an extended envelope; plausibly, subsequent accretion of large galaxies is required. The post-collapse system is in general triaxial, with rounder systems resulting from fewer, but more massive, clumps. Persistent streams of material from disrupted clumps can be found in the outer regions of the final system, and at an overdensity of around 0.75, can cover 10 to 30 per cent of the sky.     

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.