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

On the nature of superoutbursts in dwarf novae

We determine a crucial feature of the dark halo density distribution from the fact that the luminous matter dominates the gravitational potential at about one disc scalelength R _d, but at the optical edge     the dark matter has already become the main component of the galaxy density. From the kinematics of 137 spirals we find that the dark matter halo density profiles are self-similar at least out to R _opt and show core radii much larger than the corresponding disc scalelengths. The luminous regions of spirals consist of stellar discs embedded in dark haloes with roughly constant density. This invariant dark matter profile is very difficult to reconcile with the fundamental properties of the density distribution of cold dark matter haloes. With respect to previous work, the present evidence is obtained by means of a robust method and for a large and complete sample of normal spirals.     

Dark matter in elliptical galaxies – I. Is the total mass density profile of the NFW form or even steeper?

Elliptical galaxies are modelled as Sérsic luminosity distributions with density profiles (DPs) for the total mass adopted from the DPs of haloes within dissipationless ΛCDM (cold dark matter) N -body simulations. Ellipticals turn out to be inconsistent with cuspy low-concentration NFW models representing the total mass distribution, neither are they consistent with a steeper −1.5 inner slope, nor with the shallower models proposed by Navarro et al., nor with NFW models 10 times more concentrated than predicted, as deduced from several X-ray observations – the mass models, extrapolated inwards, lead to local mass-to-light ratios that are smaller than the stellar value inside an effective radius ( R _e), and to central aperture velocity dispersions that are much smaller than observed. This conclusion remains true as long as there is no sharp steepening (slope < −2) of the dark matter DPs just inside 0.01 virial radii.
The very low total mass and velocity dispersion produced within R _e by an NFW-like total mass profile suggests that the stellar component should dominate the dark matter component out to at least R _e. It should then be difficult to kinematically constrain the inner slope of the DP of ellipticals. The high-concentration parameters deduced from X-ray observations appear to be a consequence of fitting an NFW model to the total mass DP made up of a stellar component that dominates inside and a dark matter component that dominates outwards.
An appendix gives the virial mass dependence of the concentration parameter, central density and total mass of the Navarro et al. model. In a second appendix are given single integral expressions for the velocity dispersions averaged along the line of sight, in circular apertures and in thin slits, for general luminosity density and mass distributions, with isotropic orbits.     

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.     

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.     

An excursion set model for the distribution of dark matter and dark matter haloes

A model of the gravitationally evolved dark matter distribution, in the Eulerian space, is developed. It is a simple extension of the excursion set model that is commonly used to estimate the mass function of collapsed dark matter haloes. In addition to describing the evolution of the Eulerian space distribution of the haloes, the model allows one to describe the evolution of the dark matter itself. It can also be used to describe density profiles, on scales larger than the virial radius of these haloes, and to quantify the way in which matter flows in and out of Eulerian cells. When the initial Lagrangian space distribution is white noise Gaussian, the model suggests that the Inverse Gaussian distribution should provide a reasonably good approximation to the evolved Eulerian density field, in agreement with numerical simulations. Application of this model to clustering from more general Gaussian initial conditions is discussed at the end.     

Velocity dispersions of dwarf spheroidal galaxies: dark matter versus MOND

We present predictions for the line-of-sight velocity dispersion profiles of dwarf spheroidal galaxies and compare them to observations in the case of the Fornax dwarf. The predictions are made in the framework of standard dynamical theory of spherical systems with different velocity distributions. The stars are assumed to be distributed according to Sérsic laws with parameters fitted to observations. We compare predictions obtained assuming the presence of dark matter haloes (with density profiles adopted from N -body simulations) with those resulting from Modified Newtonian Dynamics (MOND). If the anisotropy of velocity distribution is treated as a free parameter, observational data for Fornax are reproduced equally well by models with dark matter and with MOND. If stellar mass-to-light ratio of 1 M_⊙/L_⊙ is assumed, the required mass of the dark halo is     , two orders of magnitude larger than the mass in stars. The derived MOND acceleration scale is     . In both cases a certain amount of tangential anisotropy in the velocity distribution is needed to reproduce the shape of the velocity dispersion profile in Fornax.     

The tightening of wide binaries in dSph galaxies through dynamical friction as a test of the dark matter hypothesis

We estimate the time-scales for orbital decay of wide binaries embedded within dark matter haloes, due to dynamical friction against the dark matter particles. We derive analytical scalings for this decay and calibrate and test them through the extensive use of N -body simulations, which accurately confirm the predicted temporal evolution. For density and velocity dispersion parameters as inferred for the dark matter haloes of local dSph galaxies, we show that the decay time-scales become shorter than the ages of the dSph stellar populations for binary stars composed of  1 M_⊙  stars, for initial separations larger than 0.1 pc. Such wide binaries are conspicuous and have been well measured in the solar neighbourhood. The prediction of the dark matter hypothesis is that they should now be absent from stellar populations embedded within low velocity dispersion, high-density dark mater haloes, as currently inferred for the local dSph galaxies, having since evolved into tighter binaries. Relevant empirical determinations of this will become technically feasible in the near future, and could provide evidence to discriminate between dark matter particle haloes or modified gravitational theories, to account for the high dispersion velocities measured for stars in local dSph galaxies.     

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.     

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.     

Dwarf galaxy rotation curves and the core problem of dark matter haloes

The standard cold dark matter (CDM) model has recently been challenged by the claim that dwarf galaxies have dark matter haloes with constant-density cores, whereas CDM predicts haloes with steeply cusped density distributions. Consequently, numerous alternative dark matter candidates have recently been proposed. In this paper we scrutinize the observational evidence for the incongruity between dwarf galaxies and the CDM model. To this end, we analyse the rotation curves of 20 late-type dwarf galaxies studied by Swaters. Taking the effects of beam smearing and adiabatic contraction into account, we fit mass models to these rotation curves with dark matter haloes with different cusp slopes, ranging from constant-density cores to r ⁻² cusps. Even though the effects of beam smearing are small for these data, the uncertainties in the stellar mass-to-light ratio and the limited spatial sampling of the halo's density distribution hamper a unique mass decomposition. Consequently, the rotation curves in our sample cannot be used to discriminate between dark haloes with constant-density cores and r ⁻¹ cusps. We show that the dwarf galaxies analysed here are consistent with CDM haloes in a ΛCDM cosmology, and that there is thus no need to abandon the idea that dark matter is cold and collisionless. However, the data are also consistent with any alternative dark matter model that produces dark matter haloes with central cusps less steep than r ^−1.5. In fact, we argue that based on existing H  i rotation curves alone, at best weak limits can be obtained on cosmological parameters and/or the nature of the dark matter. In order to make progress, rotation curves with higher spatial resolution and independent measurements of the mass-to-light ratio of the disc are required.     

Dark matter halo profiles in scale-free cosmologies

We explore the dependence of the central logarithmic slope of dark matter halo density profiles α on the spectral index n of the linear matter power spectrum P ( k ) using cosmological N -body simulations of scale-free models [i.e. P ( k ) ∝ k ⁿ ]. These simulations are based on a set of clear, reproducible and physically motivated criteria that fix the appropriate starting and stopping times for runs, and allow one to compare haloes across models with different spectral indices and mass resolutions. For each of our simulations we identify samples of well-resolved haloes in dynamical equilibrium and we analyse their mass profiles. By parametrizing the mass profile using a 'generalized' Navarro, Frenk & White profile in which the central logarithmic slope α is allowed to vary while preserving the r ⁻³ asymptotic form at large radii, we obtain preferred central slopes for haloes in each of our models. There is a strong correlation between α and n , such that α becomes shallower as n becomes steeper. However, if we normalize our mass profiles by r ₋₂, the radius at which the logarithmic slope of the density profile is −2, we find that these differences are no longer present. This is apparent if we plot the maximum slope     as a function of r / r ₋₂– we find that the profiles are similar for haloes forming in different n models. This reflects the importance of concentration, and reveals that the concentrations of haloes forming in steep- n cosmologies tend to be smaller than those of haloes forming in shallow- n cosmologies. We conclude that there is no evidence for convergence to a unique central asymptotic slope, at least on the scales that we can resolve.     

Dark matter subhaloes in numerical simulations

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

Dark matter distribution in dwarf spheroidal galaxies

We study the distribution of dark matter in dwarf spheroidal galaxies by modelling the moments of their line-of-sight velocity distributions. We discuss different dark matter density profiles, both cuspy and possessing flat density cores. The predictions are made in the framework of standard dynamical theory of two-component (stars and dark matter) spherical systems with different velocity distributions. We compare the predicted velocity dispersion profiles to observations in the case of Fornax and Draco dwarfs. For isotropic models the dark haloes with cores are found to fit the data better than those with cusps. Anisotropic models are studied by fitting two parameters, dark mass and velocity anisotropy, to the data. In this case all profiles yield good fits, but the steeper the cusp of the profile, the more tangential is the velocity distribution required to fit the data. To resolve this well-known degeneracy of density profile versus velocity anisotropy, we obtain predictions for the kurtosis of the line-of-sight velocity distribution for models found to provide best fits to the velocity dispersion profiles. It turns out that profiles with cores typically yield higher values of kurtosis which decrease more steeply with distance than the cuspy profiles, which will allow us to discriminate between the profiles once the kurtosis measurements become available. We also show that with present quality of the data the alternative explanation of velocity dispersions in terms of Modified Newtonian Dynamics cannot yet be ruled out.     

The dark matter distribution in disc galaxies

We use high-quality optical rotation curves of nine low-luminosity disc galaxies to obtain the velocity profiles of the surrounding dark matter haloes. We find that they increase linearly with radius at least out to the edge of the stellar disc, implying that, over the entire stellar region, the density of the dark halo is about constant.
The properties of the mass structure of these haloes are similar to those found for a number of dwarf and low surface brightness galaxies, but provide a more substantial evidence of the discrepancy between the halo mass distribution predicted in the cold dark matter scenario and those actually detected around galaxies. We find that the density law proposed by Burkert reproduces the halo rotation curves, with halo central densities ( ρ ₀∼1–4×10⁻²⁴ g cm⁻³) and core radii ( r ₀∼5–15 kpc) scaling as ρ ₀∝ r ₀^−2/3.     

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.     

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

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

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.     

The graininess of dark matter haloes

We use the recently completed one billion particle Via Lactea II Λ cold dark matter simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as the structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions, the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo cannot be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.