Observational evidence for self-interacting cold dark matter David N. Spergel and Paul J. Steinhardt Princeton University, Princeton NJ 08544 USA ing cross-section may be due to strong, short-range in- Cosmological models with cold dark matter composed of teractions, similar to neutron-neutron scattering at low- weakly interacting particles predict overly dense cores in the energies, or weak interactions mediated by the exchange centers of galaxies and clusters and an overly large number of halos within the Local Group compared to actual obser- of light particles (although not so light as to produce a vations. We propose that the conflict can be resolved if the long-range force). Depending on the interaction and the cold dark matter particles are self-interacting with a large mean free path, the requisite mass for the dark matter is scattering cross-section but negligible annihilation or dissipa- in the range 1 MeV to 10 GeV. For the purposes of our tion. In this scenario, astronomical observations may enable proposal, only two-body scattering effects are important us to study dark matter properties that are inaccessible in the so either repulsive or attractive interactions are possible. laboratory. Exchanged particles should be massive enough that they are not radiated by the scattering of dark matter parti- cles in the halo. We are led to consider self-interactions because ordi- Flat cosmological models with a mixture of ordinary nary astrophysical processes are unlikely to resolve the baryonic matter, cold matter, and cosmological constant problems with standard, weakly interacting dark matter. Consider the dwarf galaxy problem. One might suppose (or quintessence) and a nearly scale-invariant, adiabatic 15 spectrum of density fluctuations are consistent with stan- that supernova explosions could cause the galactic core dard inflationary cosmology and provide an excellent fit density to be made smoother; but, while the explosions ≫ 1 suppress star formation in dwarf galaxies, numerical sim- to current observations on large scales ( 1 Mpc). How- 16 ever, an array of observations on galactic and subgalactic ulations find that starbursts in dwarfs are very ineffi- scales (≤ few Mpc) appears to conflict with the struc- cient at removing gas or matter from the core. One might ture formation predicted by analytical calculations and also consider whether the apparent overabundance of ha- numerical simulations. The predictions are based on the los found in simulations can be explained if the low ve- locity halos form primarily low surface brightness galax- standard view of cold dark matter as consisting of parti- 17 cles with weak self-interactions, as well as weak interac- ies, which are difficult to find. However, while low brightness galaxy surveys suggest a steeper luminosity tions with ordinary matter. 18 A generic prediction for weakly self-interacting dark function outside of groups, even these surveys do not matter, independent of other details of the cosmological find enough small galaxies to eliminate the discrepancy model, is that cold dark matter forms triaxial halos with between theory and observations. If star formation in dwarfs is sufficiently suppressed,19 then they should have dense cores and significant dense substructures within the 20 halo. Yet, lensing observations of clusters2 reveal central been detected as gas clouds in the local group or ex- ternal systems. HI surveys do not find large numbers of regions (roughly galactic scale) with nearly spherical low 21 density cores. Dwarf irregular galaxies appear to have small isolated gas clouds. Even if any of the processes arXiv:astro-ph/9909386v2 28 Feb 2000 low density cores3–6 with much shallower profiles than were successful in reducing the number of visible dwarfs, predicted in numerical simulations.7, 8 The persistence the dense small halos would still persist. When these ha- los fall onto galactic disks, they will heat the stellar disks of bars in high surface brightness galaxies imply that 12, 22, 23 galaxies like our own Milky Way also have low density and destroy them. These dense halos will also set- cores.9 Observations of the Local Group reveal less than tle to the centers of the central halo and produce a high one hundred galaxies,10 while numerical simulations11, 12 density core in galaxies and clusters. Since the baryon and analytical theory13, 14 predict that there should be fraction in the centers of low surface brightness galaxies is low,17 hydrodynamic processes are not likely to alter roughly one thousand discrete dark matter halos within 3, 4 the Local Group. their dark matter profiles. In this paper, we propose that the inconsistencies with The success of the cold dark matter model on large the standard picture may be alleviated if the cold dark scales suggests that a modification of the dark matter matter is self-interacting with a large scattering cross- properties may be the best approach for resolving the section but negligible annihilation or dissipation. The problems on small scales. If the dark matter is not cold, but warm (moderately relativistic), this alleviates some key feature is that the mean free path should be in the 24 range 1 kpc to 1 Mpc at the solar radius, where the dark of these discrepancies. However, the remarkably good matter density is about 0.4 GeV/cm3. The large scatter- agreement between standard cold dark matter (CDM) models and the observed power spectrum of Lyman α ab- 1 sorbers25 likely rules out warm dark matter candidates. boson must be large enough that there is no dissipation We propose that a better resolution is dark matter when dark matter particles scatter; this requires that 2 that is cold, non-dissipative, but self-interacting. There my > 450eV(mx/GeV )(v/200km/s) , where v is the are stringent constraints on the interactions between typical velocity of dark matter particles in the halo. This 26–28 dark matter and ordinary matter and on long range mass scale for my corresponds to a force that is short- forces between dark matter particles.29 However, as long range compared to the dark matter interparticle spacing as the dark matter annihilation cross-section is much (about 1 cm in the halo). Hence, we need only consider smaller than the scattering cross-section, there are rel- two-body interactions in our analysis. If my = gmx and atively few constraints on short-range dark matter self- αy = O(1), then the maximum dark matter mass is interactions. Carlson, Machacek & Hall30, 31 suggested 1/3 a self-interacting dark matter model in which the dark λ − m < 80 g 4/3 MeV. (3) matter particle is warm rather than cold. Their model x 1Mpc assumed that the dark matter plus ordinary matter sum to the critical density predicted by inflationary cosmol- Beyond what is expressed in the relations above, there ogy. Their purpose was to reduce the power on 10 Mpc is no significant constraint on how light the dark matter scales in the dark matter mass spectrum as required if particles can be. the normalization of the spectrum is to agree fluctua- The strong self-interaction might occur if the dark mat- tions measured by the COBE satellite. Subsequently, de ter consists of particles with a conserved global charge Laix et al.32 pointed out that the alteration cannot simul- (such as a hidden baryon number) interacting through taneously fit the IRAS power spectrum and the observed a hidden gauge group (e.g., hidden color). If the gauge properties of galaxies. Our proposal does not suffer from group is unbroken, then the particles experience strong this problem because we assume a cosmological model interactions which can be non-dissipative but the par- with a low matter density and a cosmological constant (or ticle number is conserved. M-theory and superstrings, quintessence) which satisfies the COBE constraint with- for example, suggest the possibility that dark matter out self-interactions. Our proposed self-interactions do fields reside on domain walls with gauge fields sepa- not change structure on the 10 Mpc scale but only on rated from ordinary matter by an extra (small) dimen- the 1 kpc scale. Consequently, our model satisfies the sion.35, 36 Similar scenarios can be constructed in purely constraints raised by de Laix et al. four-dimensional supergravity models. Note that, if the To be more specific, we suggest that the dark matter sum of the hidden baryon number and the ordinary sec- particles should have a mean free path between ∼ 1 kpc tor baryon number is zero, then Ωx = (mx/mproton)Ωb ≃ 2 to 1 Mpc at the solar radius in a typical galaxy (mean 0.19 (mx/4GeV) (using Ωbh =0.02 and h =0.65). The density 0.4 GeV/cm3), for reasons to be explained be- particles we suggest include light versions of WIMPZIL- 33 34 low. For a particle of mass mx, this implies an elastic LAS and Q-balls. scattering cross-section of How does the mean free path of the dark matter par- ticles affect astrophysics? Since interactions only alter −1 −25 2 mx λ the evolution of cold dark matter when the density inho- σXX =8.1 × 10 cm , (1) GeV 1Mpc mogeneities are large, the cosmic microwave background (CMB) and large-scale power spectrum measurements intriguingly similar to that of an ordinary hadron. (In are not sensitive to the self-interactions. So long as the this paper, we consider the case of dark matter particles dark matter is cold (Tx/mx < Φ, where Φ is the depth of scattering only from themselves but, in a forthcoming the gravitational potential), the dark matter will collapse paper, we consider the possibility that dark matter is to form a bound halo regardless of its collisional proper- a stable, neutral hadron.) If the dark matter particles ties.