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 origin and formation of cuspy density profiles through violent relaxation of stellar systems

It is shown that the cuspy density distributions observed in the cores of elliptical galaxies can be realized by dissipationless gravitational collapse. The initial models consist of power-law density spheres such as ρ ∝ r ⁻¹ with anisotropic velocity dispersions. Collapse simulations are carried out by integrating the collisionless Boltzmann equation directly, on the assumption of spherical symmetry. From the results obtained, the extent of constant density cores, formed through violent relaxation, decreases as the velocity anisotropy increases radially, and practically disappears for extremely radially anisotropic models. As a result, the relaxed density distributions become more cuspy with increasing radial velocity anisotropy. It is thus concluded that the velocity anisotropy could be a key ingredient for the formation of density cusps in a dissipationless collapse picture. The velocity dispersions increase with radius in the cores according to the nearly power-law density distributions. The power-law index, n , of the density profiles, defined as ρ ∝ r ^{− n} , changes from n ≈2.1 at intermediate radii to a shallower power than n ≈2.1 toward the centre. This density bend can be explained from our postulated local phase-space constraint that the phase-space density accessible to the relaxed state is determined at each radius by the maximum phase-space density of the initial state.     

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.     

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.     

Spherical galaxy models with power-law logarithmic slope

We present a new family of spherically symmetric models for the luminous components of elliptical and spiral galaxies and their dark matter haloes. Our starting point is a general expression for the logarithmic slope  α( r ) = d log ρ/d log  r   from which most of the cuspy models yet available in literature may be derived. We then dedicate our attention to a particular set of models where the logarithmic slope is a power-law function of the radius r investigating in detail their dynamics assuming isotropy in the velocity space. While the basic properties (such as the density profile and the gravitational potential) may be expressed analytically, both the distribution function and the observable quantities (the surface brightness and the line-of-sight velocity dispersion) have to be evaluated numerically. We also consider the extension to anisotropic models, trying two different parametrization. As the model recently proposed by Navarro et al. as the best fit to their sample of numerically simulated haloes belongs to the family we present here, analytical approximations are given for the most useful quantities.     

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

How black holes turn cusps into cores

Collapsing collisionless particle systems form gravitational bound halos with cuspy density profiles. Also hierarchical merging of these systems produce remnants with cuspy central density profiles. These results lead to the assumption of cuspy NFW (Navarro et al., 1996) profiles for the density distribution in dark matter halos. However, observed rotation curves in disk galaxies suggest dark matter halos with isothermal cores. The same kind of problem can be found for globular clusters, which show cores in their density profiles but should be cuspy if they where formed through a cold collapse(King, 1962; Jenkins et al., 1998). We are showing how small, massive, and compact objects can efficiently transform cuspy stellar systems into density profiles with an isothermal core. This revised version was published online in August 2006 with corrections to the Cover Date.     

Tailoring triaxial N-body models via a novel made-to-measure method

The made-to-measure N -body method slowly adapts the particle weights of an N -body model, whilst integrating the trajectories in an assumed static potential, until some constraints are satisfied, such as optimal fits to observational data. I propose a novel technique for this adaption procedure, which overcomes several limitations and shortcomings of the original method. The capability of the new technique is demonstrated by generating realistic N -body equilibrium models for dark matter haloes with prescribed density profile, triaxial shape and slowly outwardly growing radial velocity anisotropy.     

Density–potential pairs for spherical stellar systems with Sérsic light profiles and (optional) power-law cores

Popular models for describing the luminosity-density profiles of dynamically hot stellar systems (e.g. Jaffe, Hernquist, Dehnen) were constructed with the desire to match the deprojected form of an   R ^1/4  light profile. Real galaxies, however, are now known to have a range of different light-profile shapes that scale with mass. Consequently, although highly useful, the above models have implicit limitations, and this is illustrated here through their application to a number of real galaxy density profiles. On the other hand, the analytical density profile given by Prugniel & Simien closely matches the deprojected form of Sérsic   R ^{1/ n}   light profiles – including deprojected exponential light profiles. It is thus applicable for describing bulges in spiral galaxies, dwarf elliptical galaxies, and both ordinary and giant elliptical galaxies. Moreover, the observed Sérsic quantities define the parameters of the density model. Here we provide simple equations, in terms of elementary and special functions, for the gravitational potential and force associated with this density profile. Furthermore, to match galaxies with partially depleted cores, and better explore the supermassive black hole/galaxy connection, we have added a power-law core to this density profile and derived similar expressions for the potential and force of this hybrid profile. Expressions for the mass and velocity dispersion, assuming isotropy, are also given. These spherical models may also prove appropriate for describing the dark matter distribution in haloes formed from ΛCDM cosmological simulations.     

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.     

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.     

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.     

GHASP: an Hα kinematic survey of spiral and irregular galaxies – V. Dark matter distribution in 36 nearby spiral galaxies

The results obtained from a study of the mass distribution of 36 spiral galaxies are presented. The galaxies were observed using Fabry–Perot interferometry as part of the GHASP survey. The main aim of obtaining high-resolution Hα 2D velocity fields is to define more accurately the rising part of the rotation curves which should allow to better constrain the parameters of the mass distribution. The Hα velocities were combined with low resolution H  i data from the literature, when available. Combining the kinematical data with photometric data, mass models were derived from these rotation curves using two different functional forms for the halo: an isothermal sphere (ISO) and a Navarro–Frenk–White (NFW) profile. For the galaxies already modelled by other authors, the results tend to agree. Our results point at the existence of a constant density core in the centre of the dark matter haloes rather than a cuspy core, whatever the type of the galaxy from Sab to Im. This extends to all types the result already obtained by other authors studying dwarf and low surface brightness galaxies but would necessitate a larger sample of galaxies to conclude more strongly. Whatever model is used (ISO or NFW), small core radius haloes have higher central densities, again for all morphological types. We confirm different halo scaling laws, such as the correlations between the core radius and the central density of the halo with the absolute magnitude of a galaxy: low-luminosity galaxies have small core radius and high central density. We find that the product of the central density with the core radius of the dark matter halo is nearly constant, whatever the model and whatever the absolute magnitude of the galaxy. This suggests that the halo surface density is independent from the galaxy type.     

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.     

Cuspy dark matter haloes and the Galaxy

The microlensing optical depth to Baade's Window constrains the minimum total mass in baryonic matter within the Solar circle to be greater than ∼     , assuming the inner Galaxy is barred with viewing angle ∼20°. From the kinematics of solar neighbourhood stars, the local surface density of dark matter is ∼     . We construct cuspy haloes normalized to the local dark matter density and calculate the circular-speed curve of the halo in the inner Galaxy. This is added in quadrature to the rotation curve provided by the stellar and ISM discs, together with a bar sufficiently massive so that the baryonic matter in the inner Galaxy reproduces the microlensing optical depth. Such models violate the observational constraint provided by the tangent-velocity data in the inner Galaxy (typically at radii     . The high baryonic contribution required by the microlensing is consistent with implications from hydrodynamical modelling and the pattern speed of the Galactic bar. We conclude that the cuspy haloes favoured by the cold dark matter cosmology (and its variants) are inconsistent with the observational data on the Galaxy.     

The extended rotation curve and the dark matter halo of M33

We present the 21-cm rotation curve of the nearby galaxy M33 out to a galactocentric distance of 16 kpc (13 disc scalelengths). The rotation curve keeps rising out to the last measured point and implies a dark halo mass ≳5×10¹⁰ M_⊙. The stellar and gaseous discs provide virtually equal contributions to the galaxy gravitational potential at large galactocentric radii, but no obvious correlation is found between the radial distribution of dark matter and the distribution of stars or gas.
Results of the best fit to the mass distribution in M33 picture a dark halo which controls the gravitational potential from 3 kpc outward, with a matter density which decreases radially as R ^−1.3. The density profile is consistent with the theoretical predictions for structure formation in hierarchical clustering cold dark matter (CDM) models, and favours lower mass concentrations than those expected in the standard cosmogony.     

Galaxy cluster mass profiles

Accurate measurements of the mass distribution in galaxy and cluster haloes are essential to test the cold dark matter (CDM) paradigm. The cosmological model predicts a universal shape for the density profile in all haloes, independent of halo mass. Its profile has a 'cuspy' centre, with no evidence for the constant density core. In this paper, we carry out a careful analysis of 12 galaxy clusters, using Chandra data to compute the mass distribution in each system under the assumption of hydrostatic equilibrium. Due to their low concentration, clusters provide ideal objects for studying the central cusps in dark matter haloes. The majority of the systems are consistent with the CDM model, but four objects exhibit flat inner density profiles. We suggest that the flat inner profile found for these clusters is due to an underestimation of the mass in the cluster centre (rather than any problem with the CDM model), since these objects also have a centrally peaked gas mass fraction. We discuss possible causes for erroneously low-mass measurements in the cores of some systems.     

Dark matter in dwarf spheroidals – II. Observations and modelling of Draco

We present stellar radial velocity data for the Draco dwarf spheroidal (dSph) galaxy obtained using the AF2/WYFFOS instrument combination on the William Herschel Telescope. Our data set consists of 186 member stars, 159 of which have good quality velocities, extending to a magnitude   V ≈19.5  with a mean velocity precision of ≈2 km s⁻¹. As this survey is based on a high-precision photometric target list, it contains many more Draco members at large radii. For the first time, this allows a robust determination of the radial behaviour of the velocity dispersion in a dSph.
We find statistically strong evidence of a rising velocity dispersion consistent with a dark matter halo that has a gently rising rotation curve. There is a <2 σ signature of rotation about the long axis, inconsistent with tidal disruption as the source of the rising dispersion. By comparing our data set with earlier velocities, we find that Draco probably has a binary distribution and fraction comparable to those in the solar neighbourhood.
We apply a novel maximum likelihood algorithm and fit the velocity data to a two parameter spherical model with an adjustable dark matter content and velocity anisotropy. Draco is best fit by a weakly tangentially anisotropic distribution of stellar orbits in a dark matter halo with a very slowly rising rotation law  ( v _circ∝ r ^0.17)  . We are able to rule out both a mass-follows-light distribution and an extended halo with a harmonic core at the 2.5 to 3 σ significance level, depending on the details of our assumptions about Draco's stellar binary population. Our modelling lends support to the idea that the dark matter in dwarf spheroidals is distributed in the form of massive, nearly isothermal haloes.     

The distribution function of dark matter in massive haloes

We study the distribution function (DF) of dark matter particles in haloes of mass range  10¹⁴–10¹⁵ M_⊙  . In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Λ cold dark matter (ΛCDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L , a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form     . This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (β₀ and  β_∞  ) and the scale of transition between them ( L ₀). The energy part   f _E ( E )  is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density   Q =ρ( r )/σ³( r )  .     

Are mergers responsible for universal halo properties?

N -body simulations of cold dark matter (CDM) have shown that, in this hierarchical structure formation model, dark matter halo properties, such as the density profile, the phase-space density profile, the distribution of axial ratio, the distribution of spin parameter and the distribution of internal specific angular momentum, follow 'universal' laws or distributions. Here, we study the properties of the first generation of haloes in a hot dark matter (HDM) dominated universe, as an example of halo formation through monolithic collapse. We find all these universalities to be present in this case also. Halo density profiles are very well fit by the Navarro, Frenk & White profile over two orders of magnitude in mass. The concentration parameter depends on mass as   c ∝ M ^0.2  , reversing the dependence found in a hierarchical CDM universe. However, the concentration–formation time relation is similar in the two cases: earlier forming haloes tend to be more concentrated than their later forming counterparts. Halo formation histories are also characterized by two phases in the HDM case: an early phase of rapid accretion followed by slower growth. Furthermore, there is no significant difference between the HDM and CDM cases concerning the statistics of other halo properties: the phase-space density profile; the velocity anisotropy profile; the distribution of shape parameters; the distribution of spin parameter and the distribution of internal specific angular momentum are all similar in the two cases. Only substructure content differs dramatically. These results indicate that mergers do not play a pivotal role in establishing the universalities, thus contradicting models which explain them as consequences of mergers.