Fitting a Self-Interacting Dark Matter Model to Data Ranging From Satellite
Galaxies to Galaxy Clusters
Sudhakantha Girmohanta1,2and Robert Shrock1
1 C. N. Yang Institute for Theoretical Physics and Department of Physics and Astronomy,
Stony Brook University, Stony Brook, New York 11794, USA and
2 Tsung-Dao Lee Institute and School of Physics and Astronomy,
Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
We present a fit to observational data in an asymmetric self-interacting dark matter model using
our recently calculated cross sections that incorporate both t-channel and u-channel exchanges in
the scattering of identical particles. We find good fits to the data ranging from dwarf galaxies to
galaxy clusters, and equivalent relative velocities from ∼20 km/sec to >
∼103km/s. We compare
our results with previous fits that used only t-channel exchange contributions to the scattering.
I. INTRODUCTION
There is strong evidence for dark matter (DM), com-
prising about 85 % of the matter in the universe. Cold
dark matter (CDM) can account for structures on length
scales larger than ∼10 Mpc [1–6] (reviews include [7–
13].) However, problems have been noted with fits to
observational data on shorter length scales of ∼1−100
kpc using early CDM simulations without baryon feed-
back [14–16]. These problems included the prediction of
greater density in the central region of galaxies than was
observed (the core-cusp problem), a greater number of
dwarf satellite galaxies than were seen (the missing satel-
lite problem), and the so-called “too big to fail” problem
pertaining to star formation in dwarf satellite galaxies.
Models with self-interacting dark matter (SIDM) have
been shown to avoid these problems (some reviews in-
clude [17–19]). The extension of cold dark matter N-
body simulations to include baryon feedback can amelio-
rate these problems with pure CDM simulations [20–33].
Nevertheless, cosmological models with self-interacting
dark matter (SIDM) are of considerable interest in their
own right and have been the subject of intensive study
[17–19, 34–89].
In the framework of a particle theory of dark mat-
ter, the rate of DM-DM scatterings is given by Γ =
(σ/mDM)vrelρDM , where σ,mDM,vrel, and ρDM are the
DM-DM scattering cross section, DM particle mass, rela-
tive velocity of two colliding DM particles, and DM mass
density, respectively. Fits to observational data on the
scale of ∼1−10 kpc, with velocities vrel ∼20−200 km/s,
yield values σ/mDM ∼1 cm2/g, while fits to observations
of galaxy clusters on distance scales of several Mpc and
vrel ∼O(103) km/s yield smaller values of σ/mDM ∼0.1
cm2/g. This implies that viable SIDM models should
have cross sections that decrease as a function of vrel.
This property can be achieved in models in which DM
particles, denoted χhere, interact via exchange of a light
(Lorentz scalar or vector) mediator field, generically de-
noted ξ.
In models with asymmetric dark matter (ADM), after
the number asymmetry is established in the early uni-
verse, the DM self-interaction occurs via the reaction
χ+χ→χ+χ . (1.1)
Because of the identical particles in the final state, a
proper treatment necessarily includes both the t-channel
and the u-channel contributions to the scattering ampli-
tude. In [89], we presented differential and integrated
cross sections for the reaction (1.1) with both the t-
channel and u-channel terms included and discussed the
differences with respect to previous calculations that in-
cluded only the t-channel term. Identical-particle effects
have also been noted in [58, 88] in a field-theoretic con-
text and in [48, 79] in the context of solutions of the
Schr¨odinger equation for potential scattering. An inter-
esting question raised by our work in [89] is the following:
how do the fits to observational data change when one
uses the cross section with both t-channel and u-channel
contributions to the scattering, as contrasted with pre-
vious fits that used only the t-channel contributions? In
the present paper we address this question using the same
observational data set that was analyzed in [52].
II. CROSS SECTIONS
First, we review the basic properties of the SIDM
model with asymmetric dark matter that we used in [89].
In this model, the dark matter particle χis a spin-1/2
Dirac fermion, and the mediator, ξ, is a real scalar, ξ=φ,
or a vector, ξ=V. Both χand ξare singlets under the
Standard Model (SM). For the version of the model with
a real scalar mediator, we take the χ-φinteraction to be of
Yukawa form, as described by the interaction Lagrangian
LYuk =yχ[ ¯χχ]φ. In the version with a vector media-
tor, the DM fermion χis assumed to be charged under
a U(1)Vgauge symmetry with gauge field Vand gauge
coupling g. Since only the product of the U(1)Vcharge
of χtimes goccurs in the covariant derivative in this the-
ory, we may, without loss of generality, take this charge
to be unity and denote the product as gχ. The corre-
sponding interaction Lagrangian is L¯χχV =gχ[¯χγµχ]Vµ.
A Higgs-type mechanism is assumed to break the U(1)V
symmetry and give a mass mVto the gauge field V. For
arXiv:2210.01132v2 [hep-ph] 14 Feb 2023