The origin of the density distribution of disc galaxies: a new problem for the standard model of disc formation

We present new models for the formation of disc galaxies that improve upon previous models by following the detailed accretion and cooling of the baryonic mass, and by using realistic distributions of specific angular momentum. Under the assumption of detailed angular momentum conservation, the discs that form have density distributions that are more centrally concentrated than an exponential. We examine the influence of star formation, bulge formation, and feedback on the outcome of the surface brightness distributions of the stars. Low angular momentum haloes yield disc galaxies with a significant bulge component and with a stellar disc that is close to exponential, in good agreement with observations. High angular momentum haloes, on the other hand, produce stellar discs that are much more concentrated than an exponential, in clear conflict with observations. At large radii, the models reveal distinct truncation radii in both the stars and the cold gas. The stellar truncation radii result from our implementation of star formation threshold densities, and are in excellent agreement with observations. The truncation radii in the density distribution of the cold gas reflect the maximum specific angular momentum of the gas that has cooled. We find that these truncation radii occur at H  i surface densities of roughly 1 M_⊙ pc⁻², in conflict with observations. We examine various modifications to our models, including feedback, viscosity, and dark matter haloes with constant-density cores, but show that the models consistently fail to produce bulge less discs with exponential surface brightness profiles. This signals a new problem for the standard model of disc formation: if the baryonic component of the protogalaxies out of which disc galaxies form has the same angular momentum distribution as the dark matter, discs are too compact.     

Completely dark galaxies: their existence, properties and strategies for finding them

There are a number of theoretical and observational hints that large numbers of low-mass galaxies composed entirely of dark matter exist in the field. The theoretical considerations follow from the prediction of cold dark matter theory that there exist many low-mass galaxies for every massive one. The observational considerations follow from the observed paucity of these low-mass galaxies in the field but not in dense clusters of galaxies; this suggests that the lack of small galaxies in the field is due to the inhibition of star formation in the galaxies as opposed to the fact that their small dark matter haloes do not exist. In this work we outline the likely properties of low-mass dark galaxies, and describe observational strategies for finding them, and where in the sky to search. The results are presented as a function of the global properties of dark matter, in particular the presence or absence of a substantial baryonic dark matter component. If the dark matter is purely cold and has a Navarro, Frenk & White density profile, directly detecting dark galaxies will only be feasible with present technology if the galaxy has a maximum velocity dispersion in excess of 70 km s⁻¹, in which case the dark galaxies could strongly lens background objects. This is much higher than the maximum velocity dispersions in most dwarf galaxies. If the dark matter in galaxy haloes has a baryonic component close to the cosmic ratio, the possibility of directly detecting dark galaxies is much more realistic; the optimal method of detection will depend on the nature of the dark matter. A number of more indirect methods are also discussed.     

The role of dark matter in the galaxy mass-size relationship

The observed relationship between stellar mass and effective radius for early type galaxies, pointed out by many authors, is interpreted in the context of Clausius’ virial maximum theory. In this view, it is strongly underlined that the key of the above mentioned correlation is owing to the presence of a deep link between cosmology and the existence of the galaxy Fundamental Plane. Then the ultimate meaning is: understanding visible mass-size correlation and/or Fundamental Plane means understanding how galaxies form. The mass-size relationship involves baryon (mainly stellar) mass and its typical dimension related to the light, but it gets memory of the cosmological mass variance at the equivalence epoch. The reason is that the baryonic component virializes by sharing virial energy in about equal amount between baryons and dark matter, this sharing depending, in turn, on the steepness of the dark matter distribution. The general strategy consists in using the two-component tensor virial theorem for determining the virialized baryonic configurations. A King and a Zhao density profile are assumed for the inner baryonic and the outer dark matter component, respectively, at the end of the relaxation phase. All the considerations are restricted to spherical symmetry for simplicity. The effect of changing the dark-to-baryon mass ratio, m, is investigated inside a ΛCDM scenario. A theoretical mass-size relation is expressed for the baryonic component, which fits fairly well to the data from a recently studied galaxy sample. Finally, the play of intrinsic dispersion on the mass ratio, m, is discussed in the light of the cusp/core problem and some consequences are speculated about the existence of a limit m_l expected by the theory.     

The angular momentum content of dwarf galaxies: new challenges for the theory of galaxy formation

We compute the specific angular momentum distributions for a sample of low-mass disc galaxies observed by Swaters. We compare these distributions to those of dark matter haloes obtained by Bullock et al. from high-resolution N -body simulations of structure formation in a ΛCDM universe. We find that although the disc mass fractions are significantly smaller than the universal baryon fraction, the total specific angular momenta of the discs are in good agreement with those of dark matter haloes. This suggests that discs form out of only a small fraction of the available baryons, but yet manage to draw most of the available angular momentum. In addition we find that the angular momentum distributions of discs are clearly distinct from those of the dark matter; discs lack predominantly both low and high specific angular momenta. Understanding these findings in terms of a coherent picture for disc formation is challenging. Cooling, feedback and stripping, which are the main mechanisms to explain the small disc mass fractions found, seem unable to simultaneously explain the angular momentum distributions of the discs. In fact, it seems that the baryons that make up the discs must have been born out of angular momentum distributions that are clearly distinct from those of ΛCDM haloes. However, the dark and baryonic mass components experience the same tidal forces, and it is therefore expected that they should have similar angular momentum distributions. Therefore, understanding the angular momentum content of disc galaxies remains an important challenge for our picture of galaxy formation.     

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

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 Effect of Central Baryonic Cores in Dark Halos on the Evaluation of Strong Lensing Probabilities

We present an estimate of the strong lensing probability by dark halos, with emphasis on the role of the baryonic matter arising purely from radiative cooling. We treat the contribution of the cooled baryons optimistically with all the cooled baryons confined within a central core, and including no feedback process from stellar evolution. Our two-component model provides a strong lensing probability that is in good agreement with the observed distribution of multiple images of quasars, provided that the cooled baryons are deposited within a spherical region of radius of 0.1 times the virial radius and follow an isothermal profile. It is pointed out that strong lensing may be used as an additional probe‘of baryon physics in dark halos though this may meanwhile complicate the test of the inner density profiles of dark matter in halos using the observed strong lensing probability.     

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.     

Properties of dark matter haloes

An overview is presented of the main properties of dark matter haloes, as we know them from observations, essentially from rotation curves around spiral and dwarf galaxies. Detailed rotation curves are now known for more than a thousand galaxies, revealing that they are not so flat in the outer parts, but rising for late-types, and falling for early-types. A well-established result now is that most bright galaxies are not dominated by dark matter inside their optical disks. Only for dwarfs and LSB (Low Surface Brightness galaxies) dark matter plays a dominant role in the visible regions. The 3D-shape of haloes are investigated through several methods that will be discussed: polar rings, flaring of HI planes, X-ray isophotes. It is not yet possible with rotation curves to know how far haloes extend, but tentatives have been made. It will be shown that the dark matter appears to be coupled to the gas in spirals and dwarfs, suggesting that dark baryons could play a major role in rotation curves. Theories proposing to replace the non-baryonic dark matter by a different dynamical or gravity law, such as MOND, have to take into account the dark baryons, especially since their spatial distribution is likely to be quite different from the visible matter.     

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.     

A dark matter disc in three cosmological simulations of Milky Way mass galaxies

Making robust predictions for the phase-space distribution of dark matter at the solar neighbourhood is vital for dark matter direct-detection experiments. To date, almost all such predictions have been based on simulations that model the dark matter alone. Here, we use three cosmological hydrodynamic simulations of bright, disc-dominated galaxies to include the effects of baryonic matter self-consistently for the first time. We find that the addition of baryonic physics drastically alters the dark matter profile in the vicinity of the solar neighbourhood. A stellar/gas disc, already in place at high redshift, causes merging satellites to be dragged preferentially towards the disc plane where they are torn apart by tides. This results in an accreted dark matter disc that contributes ∼0.25–1.5 times the non-rotating halo density at the solar position. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts Weakly Interacting Massive Particle (WIMP) capture in the Earth and Sun, boosts the annual modulation signal and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined.     

Constraining the Mass Distribution of Cluster Galaxies by Weak Lensing

Analysing the weak lensing distortions of the images of faint background galaxies provides a means to constrain the average mass distribution of cluster galaxies and potentially to test the extent of their dark matter haloes as a function of the density of their environment. The observable image distortions are a consequence of the interplay between the effects of a global cluster mass distribution and the perturbations resulting from individual cluster galaxies. Starting from a reconstruction of the cluster mass distribution with conventional techniques, we apply a maximum likelihood method to infer the average properties of an ensemble of cluster galaxies. From simulations this approach is found to be reliable as long as the galaxies including their dark matter haloes only contribute a small fraction to the total mass of the system. If their haloes are extended, the galaxies contain a substantial mass fraction. In this case our method is still applicable in the outer regions of clusters, where the surface mass density is low, but yields biased estimates of the parameters describing the mass profiles of the cluster galaxies in the central part of the cluster. In that case it will be necessary to resort to more sophisticated strategies by modelling cluster galaxies and an underlying global mass distribution simultaneously. We conclude that galaxy–galaxy lensing in clusters provides a unique means to probe the presence and extent of dark haloes of cluster galaxies.     

Dark matter halo structure in CDM hydrodynamical simulations

We have carried out a comparative analysis of the properties of dark matter haloes in N -body and hydrodynamical simulations. We analyse their density profiles, shapes and kinematical properties with the aim of assessing the effects that hydrodynamical processes might produce on the evolution of the dark matter component. The simulations performed allow us to reproduce dark matter haloes with high resolution, although the range of circular velocities is limited. We find that for haloes with circular velocities of [150–200] km s⁻¹ at the virial radius, the presence of baryons affects the evolution of the dark matter component in the central region, modifying the density profiles, shapes and velocity dispersions. We also analyse the rotation velocity curves of disc-like structures and compare them with observational results.     

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.     

Marriage à-la-MOND: Baryonic dark matter in galaxy clusters and the cooling flow puzzle

I start with a brief introduction to MOND phenomenology and its possible roots in cosmology—a notion that may turn out to be the most far reaching aspect of MOND. Next I discuss the implications of MOND for the dark matter (DM) doctrine: MOND’s successes imply that baryons determine everything. For DM this would mean that the puny tail of leftover baryons in galaxies wags the hefty DM dog. This has to occur in many intricate ways, and despite the haphazard construction history of galaxies—a very tall order. I then concentrate on galaxy clusters in light of MOND, which still requires some yet undetected cluster dark matter, presumably in some baryonic form (CBDM). This CBDM might contribute to the heating of the X-ray emitting gas and thus alleviate the cooling flow puzzle. MOND, qua theory of dynamics, does not directly enter the microphysics of the gas; however, it does force a new outlook on the role of DM in shaping the cluster gas dynamics: MOND tells us that the cluster DM is not cold dark matter, is not so abundant, and is not expected in galaxies; it is thus not subject to constraints on baryonic DM in galaxies. The mass in CBDM required in a whole cluster is, typically, similar to that in hot gas, but is rather more centrally concentrated, totally dominating the core. The CBDM contribution to the baryon budget in the universe is thus small. Its properties, deduced for isolated clusters, are consistent with the observations of the “bullet cluster”. Its kinetic energy reservoir is much larger than that of the hot gas in the core, and would suffice to keep the gas hot for many cooling times. Heating can be effected in various ways depending on the exact nature of the CBDM, from very massive black holes to cool, compact gas clouds.     

Dark haloes of spiral galaxies: ISO photometry

We exclude hydrogen-burning stars, of any mass above the hydrogen-burning limit and any metallicity, as significant contributors to the massive haloes deduced from rotation curves to dominate the outer parts of spiral galaxies. We present and analyse images of four nearly edge-on bulgeless spiral galaxies (UGC 711, NGC 2915, UGC 12426, UGC 1459) obtained with ISOCAM (The CAMera instrument on board the Infrared Space Observatory ) at 14.5 and 6.75 μm. Our sensitivity limit for detection of any diffuse infrared emission associated with the dark haloes in these galaxies is a few tens of μJy per 6 × 6 arcsec² pixel, with this limit currently set by remaining difficulties in modelling the non-linear behaviour of the detectors. All four galaxies show zero detected signal from extended non-disc emission, consistent with zero halo-like luminosity density distribution. The 95 per cent upper limit on any emission, for NGC 2915 in particular, allows us to exclude very low mass main-sequence stars ( M  > 0.08 M⊙) and young brown dwarfs (≲1 Gyr) as significant contributors to dark matter in galactic haloes. Combining our results with those of the Galactic microlensing surveys, which exclude objects with M  < 0.01 M⊙, excludes almost the entire possible mass range of compact baryonic objects from contributing to Galactic dark matter.     

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.     

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

Halo masses for optically selected and for radio-loud AGN from clustering and galaxy–galaxy lensing

We compute two-point correlation functions and measure the shear signal due to galaxy–galaxy lensing for 80 000 optically identified and 5700 radio-loud active galactic nuclei (AGN) from Data Release 4 of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed   M _*  . Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive haloes than optical AGN: in our samples, their mean halo masses are  1.6 × 10¹³  and  8 × 10¹¹  h ⁻¹ M_⊙  , respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter haloes of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in haloes where gas cools quasi-statically.