arXiv:astro-ph/9909386v2 28 Feb 2000 ocretosrain nlresae ( fit scales large excellent on an observations current provide to and cosmology stan- inflationary with consistent dard are adiabatic fluctuations scale-invariant, density of nearly spectrum a and constant quintessence) cosmological and (or , cold matter, baryonic t in inaccessible are that enab properties laboratory. may matter observations dark study astronomical larg to scenario, us a this d In with the or self-interacting annihilation if tion. negligible resolved are but be cross-section obser- particles can scattering actual matter conflict to dark the compared number that cold Group large propose We Local overly an the vations. and within clusters halos and of in cores dense of overly centers predict particles interacting weakly ag p o1Mca h oa ais hr h dark the where the radius, GeV/cm in solar 0.4 be the about at should is Mpc density path 1 matter free to The mean kpc 1 the range dissipation. that cross- is or feature scattering annihilation key dark large negligible cold a but the with section if self-interacting alleviated is be matter may picture standard the within halos Group. matter Local dark the discrete thousand one roughly theory analytical and ao e,lnigosrain fclusters of observations lensing Yet, with the within halos halo. substructures triaxial dense significant forms and matter cores cosmological dark dense cold the that of is details model, other of independent matter, interac- weak matter. as ordinary well parti- with as of tions self-interactions, consisting weak as matter with dark cles the cold on of based view are and standard predictions calculations The analytical simulations. by numerical predicted formation ture aaislk u w ik a lohv o density low that have also imply Way galaxies Milky brightness own cores. surface our like high galaxies in bars of n ude galaxies, hundred one vr naryo bevtoso aatcadsubgalactic and ( galactic on scales observations of array an ever, rdce nnmrclsimulations. have numerical to in appear predicted cores galaxies density irregular low Dwarf low spherical cores. nearly with density scale) galactic (roughly regions ltcsooia oeswt itr fordinary of mixture a with models cosmological Flat nti ae,w rps htteicnitniswith inconsistencies the that propose we paper, this In eei rdcinfrwal efitrcigdark self-interacting weakly for prediction generic A omlgclmdl ihcl akmte opsdof composed matter dark cold with models Cosmological 9 ≤ bevtoso h oa ru eells than less reveal Group Local the of Observations e p)apast oflc ihtestruc- the with conflict to appears Mpc) few 3–6 bevtoa vdnefrsl-neatn oddr mat dark cold self-interacting for evidence Observational ihmc hloe rfie than profiles shallower much with 10 3 14 13, hl ueia simulations numerical while rdc htteesol be should there that predict rneo nvriy rneo J054USA 08544 NJ Princeton University, Princeton ,8 7, 3 h ag scatter- large The . ai .SegladPu .Steinhardt J. Paul and Spergel N. David ≫ h persistence The 2 Mpc). 1 eelcentral reveal 1 issipa- How- 1 12 11, the he le e 1 endtce sgscod ntelclgroup local the in clouds gas as detected been enlsses Isresd o n ag ubr of numbers large clouds. find gas not isolated do small surveys HI systems. ternal in formation star If suppressed, sufficiently discrepancy is observations. dwarfs the and eliminate theory to between galaxies small enough find o alot aatcdss hywl etteselrdisk stellar the heat them. will ha- destroy they these and disks, When galactic onto persist. fall still dwarfs, los would visible halos of small number dense the the reducing in successful were suppose might One explosions supernova problem. that dwarf the matter. the dark resolve Consider interacting to weakly standard, unlikely with are problems processes astrophysical nary parti- matter dark halo. of the scattering in the cles by they radiated that not enough massive are be possibl our should are of particles interactions Exchanged attractive purposes or important the repulsive For are either effects GeV. so scattering 10 two-body is to only matter MeV dark proposal, 1 the range for the mass a the requisite in produce and the to interaction path, free the as mean on light Depending so force). not long-range (although exchange particles the by light mediated of interactions in- low- weak at short-range or scattering strong, energies, neutron-neutron to to due similar be teractions, may cross-section ing oesadteosre oe pcrmo Lyman (CDM) of spectrum matter power observed dark the and cold models standard between agreement some cold, discrepancies. alleviates not these is this the matter of relativistic), dark resolving (moderately the for warm If matter but approach scales. dark best small the on the problems of be modification may a properties that suggests scales rgtesglx uvy ugs tee luminosity groups, steeper of a outside suggest function surveys galaxy ve- brightness galax- low brightness the surface if low ies, explained primarily be form halos can locity simulations in ha- found of overabundance might los apparent One the core. whether the consider from also matter or gas removing at cient hi akmte profiles. matter dark galaxies their brightness surface the low low, Since of high is centers a clusters. the produce and in and galaxies fraction halo in central the core of density centers the to tle explosions sim- numerical the galaxies, ulations while dwarf in but, formation smoother; star suppress made be to density eaeldt osdrsl-neatosbcueordi- because self-interactions consider to led are We h ucs ftecl akmte oe nlarge on model matter dark cold the of success The 17 hc r icl ofid oee,wielow while However, find. to difficult are which 17 16 yrdnmcpoessaentlkl oalter to likely not are processes hydrodynamic n htsabrt ndaf r eyineffi- very are dwarfs in starbursts that find 2 2 23 22, 12, 24 oee,termral good remarkably the However, 21 15 hs es ao ilas set- also will halos dense These ,4 3, 18 vni n fteprocesses the of any if Even ol as h aatccore galactic the cause could vnteesresd not do surveys these even 19 ter hnte hudhave should they then 20 rex- or α ab- e. s sorbers25 likely rules out warm 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 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. If the dark matter mean path were much longer than scatter through strong interactions similar to low-energy ∼ 1 Mpc, the typical dark matter particle would not ex- neutron-neutron scattering, then the cross-section is σ = perience any interactions as it moves through a halo. In 4πa2, where a is the scattering length. For neutrons, the this regime, the usual, triaxial cold with scattering length is more than 100 times its Compton dense core forms through gravitational collapse. On the −1 wavelength. Using the estimate a ≈ 100fmx , we obtain other hand, if the dark matter mean free path is much smaller than 1 kpc, then the dark matter behaves as a 1/3 λ collisional gas and this alters the halo evolution signifi- m =4 f 2/3 GeV. (2) x 1Mpc cantly. The dark matter will shock: this will heat up the low entropy material that would usually collapse to form Alternatively, the self-interaction may be weak but a core and produce a shallower density profile. Since col- longer-range, as in the case of the exchange of a light lisions tend to make the dark matter velocity distribution intermediate vector boson of mass my, in which case the isotropic, the halo can not be triaxial and will only be el- 2 4 cross-section is σ ≈ αymx/my. The mass of the vector liptical if flattened by significant rotation. Since dark

2 halos form with little angular momentum, if the dark locity dispersion in the dwarf core and v is the typical matter is not dissipative, then all halos will be nearly particle velocity in the halo.) This could produce a ram spherical. X-ray observations of clusters37 reveal that pressure drag that can slow the dwarf halo and cause it most halos are moderately ellipsoidal. This implies that to spiral into the cores of larger halos. Numerical stud- the collision time scale for dark matter near the half-mass ies of this process are required to determine whether the radius of clusters must be longer than the Hubble time: observed fraction of denser dwarfs are likely to survive one of the strongest constraints on this model. Studies after a Hubble time. suggest that polar ring galaxies38 are only mildly triax- The halos of large (e.g. L∗) galaxies moving through ial and oblate with the equatorial plane of the dark halo groups and clusters are less prone to destruction. When a nearly coinciding with that of the stellar body. If the cluster dark matter particle strikes a galaxy dark matter dark matter has an isotropic distribution function and particle, the probability that the recoiling particles will the form a disk, then the dark matter will form escape from the galaxy is significantly less than unity. a slightly flattened halo. For cross-sections near the smaller end of our suggested In our scenario, we consider a mean free path in the range, most collisions take place within the half mass ra- intermediate regime, larger than 1 kpc, but smaller than dius, where the escape velocity from a galaxy halo is com- ∼ 1 Mpc. Particles in this range have 1 − 103 interac- parable to the characteristic recoil velocity. Thus, many tions per Hubble time in the local neighborhood, which collisions will not lead to recoil energies large enough to is overdense by 106 relative to the mean density of the escape. At the smaller end of our cross-section range, . At the virial radius of a typical , the probability of a typical particle experiencing a col- which is overdense by ∼ 200 relative to the mean density lision during a Hubble time approaches unity only for of the universe, the typical particle has less than 1 colli- galaxies that fall deep into the cluster core. Towards the sion per Hubble time. Thus, near the virial radius, halos larger end of our suggested range, the galaxies are opaque can have anisotropic velocity ellipsoids and will be triax- to dark matter and the collision products are also likely ial. However, in the inner halo of galaxies, dark matter to experience multiple collisions within a massive galaxy is collisional. Infalling dark matter is scattered before halo. Because of these multiple collisions, the collision reaching the center of the galaxy so that the orbit dis- products are unlikely to escape from the galaxy. This tribution is isotropic rather than radial. These collisions effect is more important in massive halos as the optical increase the entropy of the dark matter phase space dis- depth through a galaxy scales as the one-third power of tribution and lead to a dark matter halo profile with a its mass, Thus, for large galaxies, the background cluster shallower density profile. The characteristic scale for the or group will primarily heat the entire dark halo rather core would correspond to an “optical depth” of 1, the than evaporate dark matter from it. “photosphere” of the dark matter. The presence of collisions will lead to energy transport When a dwarf halo with a low velocity dispersion falls within the dark matter halo, which eventually leads to into a larger, high velocity dispersion halo, the high ve- core collapse.39 We can obtain an estimate of the core col- locity particles will scatter off the low velocity particles. lapse time from Quinlan’s40 Fokker-Planck simulations of After the collision, neither particle is likely to be bound to the evolution of an isolated cluster of interacting particles the dwarf. As dark matter is slowly removed, the dwarf with a central density profile of r−1 and an outer profile halo expands, this also makes the halo more vulnerable of r−3. These models have a temperature inversion in to tidal stripping and shock heating. This process will the core and undergo two stages of core collapse. During slowly evaporate all substructure in the larger halo, par- the first stage, the inner region expands as heat is trans- ticularly, in the centers of galaxies, groups and clusters. ported inwards. After 0.1 half-mass relaxation times, the However, near the half-mass radius, the collision time is inner 1% of the mass has moved out in radius by factor 2. longer than the Hubble time, so substructure will only be After roughly 3 half-mass relaxation times, the whole sys- destroyed in the inner portions of halos. This dark mat- tem collapses as heat is transported from the virial radius ter evaporation will protect galaxy disks from dynamical outwards. We suspect that this second stage is delayed heating by collisions with dwarfs in the halo. Dwarfs with in our cosmological context as the infall of new mate- high central densities evaporate more slowly as particles rial is constantly adding heat near the dark halo photo- at optical depth greater than ∼ 1 are shielded by other sphere. If we model our Galaxy as starting with a den- particles from collisions. Intriguingly, all of our Galaxy’s sity profile corresponding to the Navarro-Frenk-White fit 10 7 −1 dwarf companions have very high phase space densities. to CDM models with V200 = 225 km s , Ω0 = 0.3, −1 −1 As the central density and particle velocities increase, not H0 =65kms Mpc and dimensionless concentration all collisions will lead to particles being deflected out of parameter c = 8, then its half mass radius is 220 kpc. the dwarf halo. A small fraction of the particles will be Since the density at the half mass radius is 300 times rescattered within the spheroidal core of the dwarf and smaller than the local dark matter density, the particle a fraction of the momentum absorbed. (The phase space mean free path at that radius is between 0.3 - 300 Mpc for capture is suppressed by (σ/v)3, where σ is the ve- and at that radius, the particle is in the weakly interact-

3 ing regime. Thus, for our Galaxy, the core collapse time [2] J.A. Tyson, G.P. Kochanski, & I.P. Dell’antonio, ApJ, (roughly 3 half-mass relaxation times) is between 4.5 - 498, L107 (1998). 6000 Gyr, or perhaps significantly longer if the collapse [3] B. Moore, Nature, 370, 629 (1994). stage is delayed by the infall of new material. Hence, for [4] R.A. Flores & J.A. Primack, ApJ,427, L1 (1994). most of our range of parameters, the collapse time for our [5] W.J. G. De Blok & S.S. McGaugh, Mon. Not. Roy. As- 290 Galaxy exceeds the lifetime of the universe, yet there is tron. Soc. , 533 (1997). [6] J. Dalcanton & R. Bernstein, in Galaxy Dynamics: from sufficient number of interactions to lower the dark matter the Early Universe to the Present, eds.: F. Combes, G. A. density in the inner 5 kpc of our Galaxy. As the particle Mamon, and V. Charmandaris, ASP Conference Series, mean free path approaches the lower bound (0.3 Mpc) or 23 (1999). the upper bound (300 Mpc) at the half-mass radius, our [7] J.F. Navarro, C.S. Frenk, & S.D.M. White, ApJ, 462, estimates suggest that one or the other condition is not 563 (1996). satisfied, but more accurate methods are needed to deter- [8] B. Moore, F. Governato, T. Quinn, J. Stadel, & G. Lake, mine the precise range. [Note added in proof: A. Burkert ApJ Letters, 499, L5 (1998). (astro-ph/0002409) has constructed an N-body code that [9] V.P. Debattista, & J.A. Sellwood, .ApJ (Letters), 493, simulates self-interactions and obtains results consistent L5 (1998). with these estimates based on the Fokker-Planck approx- [10] M. Mateo, Ann. Rev. Astron. Astrophys., 36, 435 (1998). imation.] [11] A.A. Klypin, A.V. Kravtsov, O. Valenzuela, & F. Prada, ApJ, 522, 82 (1999). To summarize, our estimated range of σ/m for the 2 [12] B. Moore, S. Ghigna, F. Governato, G. Lake, G., T. dark matter is between 0.45-450 cm /g or, equivalently, Quinn, J. Stadel, & P. Tozzi, ApJ Letters, 524, L19 −(25−22) 2 8 × 10 cm /GeV. Numerical calculations are es- (1999). sential for checking our approximations and refining our [13] W.H. Press, & P. Schechter, ApJ, 187, 425 (1974). estimates. Even without numerical simulations, we can [14] G. Kauffmann, S.D.M. White, & B. Guiderdoni, Mon. already make a number of predictions for the properties Not. Royal Astron. Soc. 264, 201 (1993). of galaxies in a self-interacting dark matter cosmology: [15] A. Dekel & J. Silk, ApJ 303, 39 (1986). (1) the centers of halos are spherical; (2) dark matter ha- [16] M.M. MacLow & A. Ferrara, ApJ 513, 142 (1999). los will have cores; and (3) there are few dwarf galaxies in [17] J. Dalcanton, D.N. Spergel, & F.J. Summers, ApJ, 482, groups but dwarfs persist in lower density environments; 659 (1997). and, (4) the halos of dwarf galaxies and galaxy halos [18] J. Dalcanton, D.N. Spergel, J.E. Gunn, J.E., M. Schmidt, & D.P. Schnedider, . Astron. J. 114, 635 (1997). in clusters will have radii smaller than the gravitational [19] J.V. Kepner, A. Babul, & D.N. Spergel, . ApJ 487, 61 tidal radius (due to collisional stripping). Intriguingly, (1997). current observations appear to be consistent with all of [20] L. Blitz, D.N. Spergel, P.J. Teuben, D. Hartmann, & these predictions. W.B. Burton, . ApJ 514, 818 ( 1999). [21] M. Zwann, F. Briggs, D. Sprayberry, E. Sorar, ApJ 490, 173 (1997). ACKNOWLEDGEMENTS [22] G. Toth & J.P. Ostriker, ApJ, 389, 5 (1992). [23] M. Weinberg, Mon. Not. Royal Astron. Soc. 299, 499 We would like to thank R. Dave, J. Dalcanton, G. (1998). Dvali, J. Goodman, E. Kolb, J. March-Russell, J. [24] R. Schaeffer, and J. Silk, ApJ 332, 1 (1988). Miralda-Escude, J. Ostriker, J. Peebles, J. Silk, S. [25] R.A.C. Croft, D.H. Weinberg, M. Pettini, L. Hernquist, & N. Katz, ApJ 520, 1 (1999). Tremaine, M. Turner and N. Turok for discussions. We [26] R. Bernabei, et al, Phys. Rev. Lett. 83, 4918 (1999). are grateful to the West Anglia Great Northern Railway [27] C.W. Stubbs, Phys. Rev. Letters, 70, 119 (1993). whose train delay provided the opportunity for initial [28] G. Starkman, A. Gould, R. Esmailzadeh, and S. Di- discussions. This work was initiated at the Newton In- mopoulos, Phys. Rev. D41, 3594 (1990). stitute for Mathematical Sciences. DNS acknowledges [29] B.-A. Gradwohl & J.A. Frieman, ApJ, 398, 407 (1992). the NASA Theory program and the NASA MAP Satel- [30] E.D. Carlson, M.E. Machacek, & L.J. Hall, ApJ 398, 43 lite program for support and PJS was supported by the (1992). US Department of Energy grant DE-FG02-91ER40671. [31] M.E. Machacek, ApJ, 431, 41 (1994). [32] A.A. de Laix, R.J. Scherrer, and R.K. Schaeffer, ApJ 452, 495 (1995). [33] E.W. Kolb, D.J.H. Chung,. & A. Riotto, hep-ph/9810361 (1998) [34] A. Kusenko & M. Shaposhnikov, Phys.Lett. B418, 46 (1998). [1] N., Bahcall, J.P. Ostriker, S. Perlmutter, S., & P.J. Stein- [35] G. Dvali, Phys. Lett. B459, 489 (1999). hardt, Science, 284, 1481 (1999). [36] G. Dvali & G. Gabadadze, Phys. Lett. B460 47 (1999). [37] J.J. Mohr, A.E. Evrard, D.G. Fabricant, & M.J. Geller, ApJ 447, 8 (1995).

4 [38] P.D. Sackett, “The Shape of Dark Matter Halos”, to ap- pear in Galaxy Dynamics, eds. D. Merritt, J.A. Sellwood, and M. Valluri (1999) [39] J. Goodman, Ap. J. 313, 576 (1987). [40] G. Quinlan, New Astronomy, 1, 255 (1996).

5