DARK MATTER 101

DARK MATTER 101

From production to detection

David G. Cerdeño
CONTENTS

1 Motivation for dark matter 1

1.1 Motivation for Dark Matter 1
1.1.1 Galactic scale 2
1.1.2 Galaxy Clusters 3
1.1.3 Cosmological scale 4
1.2 Dark Matter properties 5
1.2.1 Neutral 5
1.2.2 Nonrelativistic 5
1.2.3 NonBaryonic 6
1.2.4 Long-lived 6
1.2.5 Collisionless 7

2 Freeze Out of Massive Species 9

2.1 Cosmological Preliminaries 9
2.2 Time evolution of the number density 12
2.2.1 Freeze out of relativistic species 15
2.2.2 Freeze out of non-relativistic species 16
2.2.3 WIMPs 17
2.3 Computing the DM annihilation cross section 17
2.3.1 Special cases 18
v
vi CONTENTS

3 Direct DM detection 23
3.1 Preliminaries 23
3.1.1 DM flux 23
3.1.2 Kinematics 23
3.2 The master formula for direct DM detection 24
3.2.1 The scattering cross section 24
3.2.2 The importance of the threshold 25
3.2.3 Velocity distribution function 26
3.2.4 Energy resoultion 26
3.3 Exponential spectrum 26
3.4 Annual modulation 26
3.5 Directional detection 27
3.6 Coherent neutrino scattering 27
3.7 Inelastic 27

4 Neutrinos 29
4.1 Preliminaries - copied form internet 29
References 31
CHAPTER 1

MOTIVATION FOR DARK MATTER

The existence of a vast amount of dark matter (DM) in the Universe is supported by many
astrophysical and cosmological observations. The latest measurements indicate that ap-
proximately a 27% of the Universe energy density is in form of a new type of non-baryonic
cold DM. Given that the Standard Model (SM) of particle physics does not contain any vi-
able candidate to account for it, DM can be regarded as one of the clearest hints of new
physics.

1.1 Motivation for Dark Matter

Astrophysical and Cosmological observations have provided substantial evidence that point
towards the existence of vast amounts of a new type of matter, that does not emit or absorb
light. All astrophysical evidence for DM is solely based on gravitational effects (either
trough the observation of dynamical effects, deflection of light by gravitational lensing or
measurements of the gravitational potential of galaxy clusters), which cannot be accounted
for by just the observed luminous matter. The simplest way to solve these problems is the
inclusion of more matter (which does not emit light - and is therefore dark in the astro-
nomical sense1 ). Modifications in the Newtonian equation relating force and accelerations
have also been suggested to address the problem at galactic scales, but this hypothesis is

1 Since dark matter does not absorb light, a more adequate name would have been transparent matter.

Dark Stuff. 1
By D. G. Cerdeño, IPPP, University of Durham
2 MOTIVATION FOR DARK MATTER

Figure 1.1 Left) Vera Rubin. Right) Rotation curve of a spiral galaxy, where the contribution from
the luminous disc and dark matter halo is shown by means of solid lines.

insufficient to account for effects at other scales (e.g., cluster of galaxies) or reproduce the
anisotropies in the CMB.
No known particle can play the role of the DM (we will later argue that neutrinos con-
tribute to a small part of the DM). Thus, this is one of the clearest hints for Physics Beyond
the Standard Model and provides a window to new particle physics models. In the follow-
ing I summarise some of the main pieces of evidence for DM at different scales.
I recommend completing this section with the first chapters of Ref. [1] and the recent
article [2].

1.1.1 Galactic scale

Rotation curves of spiral galaxies Rotation curves of spiral galaxies are probably the
best-known examples of how the dynamical properties of astrophysical objects are affected
by DM. Applying Gauss Law to a spiral galaxy (one can safely ignore the contribution
from the spiral arms and assume a spherical distribution of matter in the bulge) leads to a
simple relation between the rotation velocity of objects which are gravitationally bound to
the galaxy and their distance to the galactic centre:
r
GM (r)
v= , (1.1)
r
where M (r) is the mass contained within the radius r. In the outskirts of the galaxy,
where we expect that M does not increase any more, we would therefore expect a decay
vrot ∝ r−1/2 .
Vera Rubin’s observations of rotation curves of spiral galaxies [3, 4] showed a very slow
decrease with the galactic radius. The careful work of Bosma [5], van Albada and Sancisi
[6] showed that this flatness could not be accounted for by simply modifying the relative
weight of the diverse galactic components (bulge, disc, gas), a new component was needed
with a different spatial distribution (see Fig. 1.1).
Notice that the flatness of rotation curves can be obtained if a new mass component is
introduced, whose mass distribution satisfies M (r) ∝ r in eq.(1.1). This is precisely the
 MOTIVATION FOR DARK MATTER 3

Figure 1.2 Left) Coma cluster and F. Zwicky, who carried out measurements of the peculiar
velocities of this object. Right) Modern techniques [7], based on gravitational lensing, allow for a
much more precise determination of the total mass of this object.

relation that one expects for a self-gravitational gas of non-interacting particles. This halo
of DM can extend up to ten times the size of the galactic disc and contains approximately
an 80% of the total mass of the galaxy.
Since then, flat rotation curves have been found in spiral galaxies, further strengthening
the DM hypothesis. Of course, our own galaxy, the Milky Way is no exception. N-body
simulations have proved to be very important tools in determining the properties of DM
haloes. These can be characterised in terms of their density profile ρ(r) and the velocity
distribution function f (v). Observations of the local dynamics provide a measurement of
the DM density at our position in the Galaxy. Up to substantial uncertainties, the local
DM density can vary in a range ρ0 = 0.2 − 1 GeV cm−3 . It is customary to describe
the DM halo in terms of a Spherical Isothermal Halo, in which the velocity distribution
follows a Maxwell-Boltzmann law, but deviations from this are also expected. Finally, due
to numerical limitations, current N-body simulations cannot predict the DM distribution at
the centre of the galaxy. Whereas some results suggest the existence of a cusp of DM in
the galactic centre, other simulations seem to favour a core. Finally, the effect of baryons
is not easy to simulate, although substantial improvements have been recently made.

Local probes

1.1.2 Galaxy Clusters

Peculiar motion of clusters. Fritz Zwicky studied the peculiar motions of galaxies in
the Coma cluster [8, 9]. Assuming that the galaxy cluster is an isolated system, the virial
theorem can be used to relate the average velocity of objects with the gravitational potential
(or the total mass of the system).
As in the case of galaxies, this determination of the mass is insensitive to whether ob-
jects emit any light or not. The results can then be contrasted with other determinations
that are based on the luminosity. This results in an extremely large mass-to-light ratio,
indicative of the existence of large amounts of missing mass, which can be attributed to a
DM component.
4 MOTIVATION FOR DARK MATTER

Figure 1.3 Left) Deep Chandra image of the Bullet cluster. Green lines represent mass contours
from weak lensing. Right) Dark filament in the system Abell 222/223, reconstructed using weak
lensing.

Modern determinations through weak lensing techniques provide a better gravitational

determination of the cluster masses [10, 7] (see Fig. 1.2). I recommend reading through
Ref.[9] for a derivation of the virial theorem in the context of Galaxy clusters.

Dynamical systems. The Bullet Cluster (1E 0657-558) is a paradigmatic example of

the effect of dark matter in dynamical systems. It consists of two galaxy clusters which
underwent a collision. The visible components of the cluster, observed by the Chandra X-
ray satellite, display a characteristic shock wave (which gives name to the whole system).
On the other hand, weak-lensing analyses, which make use of data from the Hubble Space
Telescope, have revealed that most of the mass of the system is displaced from the visible
components. The accepted interpretation is that the dark matter components of the clusters
have crossed without interacting significantly (see e.g., Ref. [11, 12]).
The Bullet Cluster is considered one of the best arguments against MOND theories
(since the gravitational effects occur where there is no visible matter). It also sets an upper
bound on the self-interaction strength of dark matter particles.

DM filaments. Observations of the distribution of luminous matter at large scales have

shown that it follows a filamentary structure. Numerical simulations of structure formation
with cold DM have been able to reproduce this feature. To date, it is well understood
that DM plays a fundamental role in creating that filamentary network, gravitationally
trapping the luminous matter. Recently, the comparison of the distribution of luminous
matter in the Abell 222/223 supercluster with weak-lensing data has shown the existence
of a dark filament joining the two clusters of the system. That filament, having no visible
counterpart, is believed to be made of DM.

1.1.3 Cosmological scale

Finally, DM has also left its footprint in the anisotropies of the Cosmic Microwave Back-
ground (CMB). The analysis of the CMB constitutes a primary tool to determine the cos-
mological parameters of the Universe. The data obtained by dedicated satellites in the past
decades has confirmed that we live in a flat Universe (COBE), dominated by dark matter
 DARK MATTER PROPERTIES 5

Figure 1.4 Left) Contribution to the energy density for each of the components of the Universe.
Right) Planck temperature map.

and dark energy (WMAP), whose cosmological abundances have been determined with
great precision (Planck).
The abundance of DM is normally expressed in terms of the cosmological density pa-
rameter, defined as ΩDM h2 = ρDM /ρc where ρc is the critical density necessary to re-
cover a flat Universe and h = 0.7 is the normalised Hubble parameter. The most recent
measurements by the Planck satellite, combined with data obtained from Supernovae (that
trace the Universe expansion) yield
ΩCDM h2 = 0.1196 ± 0.0031 . (1.2)
Given that Ω ≈ 1, this means that dark matter is responsible for approximately a 26% of
the Universe energy density nowadays. Even more surprising is the fact that another exotic
component is needed, dark energy, which makes up approximately the 69% of the total
energy density (see Fig. 1.4).

1.2 Dark Matter properties

1.2.1 Neutral
It is generally argued that DM particles must be electrically neutral. Otherwise they would
scatter light and thus not be dark. Similarly, constrains on charged DM particles can be
extracted from unsuccessful searches for exotic atoms. Constraints on heavy millicharged
particles are inferred from cosmological and astrophysical observations as well as direct
laboratory tests [13, 14, 15]. Millicharged DM particles scatter off electrons and protons
at the recombination epoch via Rutherford-like interactions. If millicharged particles cou-
ple tightly to the baryonphoton plasma during the recombination epoch, they behave like
baryons thus affecting the CMB power spectrum in several ways [13, 14]. For particles
much heavier than the proton, this results in an upper bound of its charge  [14]
 ≤ 2.24 × 10−4 (M/1 TeV)1/2 . (1.3)
Similarly, direct detection places upper bounds on the charge of the DM particle [16]
 ≤ 7.6 × 10−4 (M/1 TeV)1/2 . (1.4)

1.2.2 Nonrelativistic

Numerical simulations of structure formation in the Early Universe have become a very
useful tool to understand some of the properties of dark matter. In particular, it was soon
6 MOTIVATION FOR DARK MATTER

found that dark matter has to be non-relativistic (cold) at the epoch of structure forma-
tion. Relativistic (hot) dark matter has a larger free-streaming length (the average distance
traveled by a dark matter particle before it falls into a potential well). This leads to incon-
sistencies with observations.
However, at the Galactic scale, cold dark matter simulations lead to the occurrence of
too much substructure in dark matter haloes. Apparently this could lead to a large number
of subhaloes (observable through the luminous matter that falls into their potential wells).
It was argued that if dark matter was warm (having a mass of approximately 2 − 3 keV)
this problem would be alleviated.
Modern simulations, where the effect of baryons is included, are fundamental in order
to fully understand structure formation in our Galaxy and determine whether dark matter
is cold or warm.

1.2.3 NonBaryonic
The results of the CMB, together with the predictions from Big Bang nucleosynthesis,
suggest that only 4 − 5% of the total energy budget of the universe is made out of ordi-
nary (baryonic) matter. Given the mismatch of this with the total matter content, we must
conclude that DM is non-baryonic.

Neutrinos. Neutrinos deserve special mention in this section, being the only viable non-
baryonic DM candidate within the SM. Neutrinos are very abundant particles in the Uni-
verse and they are known to have a (very small) mass. Given that they also interact very
feebly with ordinary matter (only through the electroweak force) they are in fact a com-
ponent of the DM. There are, however various arguments that show that they contribute in
fact to a very small part.
First, neutrinos are too light. Through the study of the decoupling of neutrinos in the
early universe we can compute their thermal relic abundance. Since neutrinos are relativis-
tic particles at the time of decoupling, this is in fact a very easy computation (we will come
back to this in Section 2.2.1), and yields
P
mi
Ων h2 ≈ i . (1.5)
91 eV

Using current upper bounds on the neutrino mass, we obtain Ων h2 < 0.003, a small
fraction of the total DM abundance.
Second, neutrinos are relativistic (hot) at the epoch of structure formation. As men-
tioned above, hot DM leads to a different hierarchy of structure formation at large scales,
with large objects forming first and small ones occurring only after fragmentation. This is
inconsistent with observations.

1.2.4 Long-lived
Possibly the most obvious observation is that DM is a long-lived (if not stable) particle.
The footprint of DM can be observed in the CMB anisotropies, its presence is essential
for structure formation and we can feel its gravitational effects in clusters of galaxies and
galaxies nowadays.
Stable DM candidates are common in models in which a new discrete symmetry is
imposed by ensuring that the DM particle is the lightest with an exotic charge (and there-
 DARK MATTER PROPERTIES 7

fore its decay is forbidden). This is the case, e.g., in Supersymmetry (when R-parity is
imposed), Kaluza-Klein scenarios (K-parity) or little Higgs models.
However, stability is not required by observation. DM particles can decay, as long
as their lifetime is longer than the age of the universe. Long-lived DM particles feature
very small couplings. Characteristic examples are gravitinos (whose decay channels are
gravitationally suppressed) or axinos (which decays through the axion coupling).

1.2.5 Collisionless
Dynamical systems, such as cluster collisions, set an upper bound to the self-interactions
of DM particles. Observations seem to suggest that the DM component in these objects is
mostly collision-less, thus behaving very differently than ordinary matter. Dark matter’s
lack of deceleration in the bullet cluster constrains its self-interaction cross-section σ/m <
1.25 cm2 g−1 ≈ 2 barn GeV−1 .
Notice however, that self-interacting dark matter with a cross section in the range 0.1 <
σ/m < 1 cm2 g−1 can be very beneficial in order to alleviate the problems with the amount
of substructure in numerical simulations of DM haloes.
[?]
CHAPTER 2

FREEZE OUT OF MASSIVE SPECIES

In this chapter we will address the computation of the relic abundance of dark matter
particles, making special emphasis in the case of thermal production in the Early Universe.

2.1 Cosmological Preliminaries

This section does not intend to be a comprehensive review on Cosmology, but only an
introduction to some of the elements that we will need for the calculation of Dark Matter
freeze-out.
We can describe our isotropic and homogeneous Universe in terms of the Friedman-
Lemaı̂tre-Robertson-Walker (FLRW) metric, which is exact solution of Einstein’s field
equations of general relativity

dr2
 
ds2 = dt2 − a2 (t) + r 2
(dθ 2
+ sin θdφ2
) = gµν dxµ dxν . (2.1)
1 − kr2
The constant k = {−1, 0, +1} corresponds to the spatial curvature, with k = 0 corre-
sponding to a flat Universe (the choice we will be making in these notes). Remember that
the affine connection, defined as
1 µσ
Γµνλ = g (gσ ν, λ + gσλ,ν − gνλ,σ ) , (2.2)
2
is greatly simplified, since most of the derivatives vanish.
Dark Stuff. 9
By D. G. Cerdeño, IPPP, University of Durham
10 FREEZE OUT OF MASSIVE SPECIES

In the following we are going to work with a radiation-dominated Universe. Notice that
matter-radiation equality only occurs very late (when the Universe is approximately 60 kyr)
and dark matter freeze-out occurs before BBN. The Hubble parameter for a radiation-
dominated Universe reads
2
1/2 T
H = 1.66 g∗ , (2.3)
MP
where MP = 1.22 × 1019 GeV.
It is customary to define the dimensionless parameter x = m/T (where m is a mass pa-
rameter that we will later associate to the DM mass) and extract the explicit x dependence
from the Hubble parameter to define H(m) as follows
1/2 m2
H(m) = 1.66 g∗ = Hx2 . (2.4)
MP
In this section we will try to compute the time evolution of the number density of dark
matter particles, in order to be able to compute their relic abundance today and what this
implies in the interaction strength of dark matter particles. The phase space distribution
function f describes the occupancy number in phase space for a given particle in kinetic
equilibrium, and distinguishes between fermions and bosons.
1
f= , (2.5)
e(E−µ)/T ± 1
where the (−) sign corresponds to bosons and the (+) sign to fermions. E is the energy
and µ the chemical potential. For species in chemical equilibrium, the chemical potential
is conserved in the interactions. Thus, for processes such as i + j ↔ c + d we have
µi +µj = µc +µd . Notice then that all chemical potentials can be expressed in terms of the
chemical potentials of conserved quantities, such as the baryon chemical potential µB . The
number of independent chemical potentials corresponds to conserved particle numbers.
This implies, for example, that given a particle with µi , the corresponding antiparticle
would have the opposite chemical potential −µi . For the same reason, since the number of
photons is not conserved in interactions, µγ = 0
Using the expression of the phase space distribution function (2.5), and integrating in
phase space, we can compute a series of observables in the Universe. In particular, the
number density of particles, n, the energy density, ρ, and pressure, p, for a dilute and
weakly-interacting gas of particles with g internal degrees of freedom read
Z
g
n = f (p) d3 p, (2.6)
(2π)3
Z
g
ρ = E(p) f (p) d3 p, (2.7)
(2π)3
|p|2
Z
g
p = 3
f (p) d3 p. (2.8)
(2π) 3E(p)
It is customary (and very convenient) to define densities normalised by the time de-
pendent volume a(t)−3 . The reason for this is that in the absence of number changing
processes, the density remains constant with time evolution (or redshift). Notice that since
the evolution of the Universe is isoentropic, the entropy density s = S/a3 has precisely
that dependence. Applying this prescription to the number density of particles, we define
the yield as a fraction of the number density and the entropy density as
n
Y = . (2.9)
s
 COSMOLOGICAL PRELIMINARIES 11

Notice that, in the absence of number-changing processes, the yield remains constant.
The evolution of the entropy density as a function of the temperature is given by 1

2π 2
s= g∗s T 3 , (2.10)
45
where the effective number of relativistic degrees of freedom for entropy is
 3  3
X Ti 7 X Ti
g∗s = g + g . (2.11)
T 8 T
bosons fermions

Remember also that we can express the energy density as

π2
ρ= g∗ T 4 , (2.12)
30
in terms of the relativistic number of degrees of freedom
 4  4
X Ti 7 X Ti
g∗ = g + g . (2.13)
T 8 T
bosons fermions

In these two equations, T is the temperature of the plasma (in equilibrium) and Ti is the
effective temperature of each species.
Solving the integral in eq. (2.6) explicitly for relativistic and non-relativistic particles,
and expressing the results in terms of the Yield results in the following expressions.

relativistic species
gef f
n= ζ(3)T 3 , (2.14)
π2
where gef f = g for bosons and gef f = 43 g for fermions2 . Then, using eq. (2.9), the
Yield at equilibrium reads
45 gef f gef f
Yeq = 4
ζ(3) ≈ 0.278 . (2.15)
2π g∗s g∗s

non-relativistic species
 3/2
mT
n = gef f e−m/T . (2.16)
2π

Then the Yield at equilibrium reads

45  π 1/2 gef f  m 3/2 −m/T

Yeq = e . (2.17)
2π 4 8 g∗s T

1 To arrive at this equation, one can calculate s = (p + ρ)/T for fermions and bosons, using the corresponding

expression for the phase space distribution function.

2 We are using here the approximation E ≈ |~ p| in the relativistic limit, and the integrals 0∞ p2 /(ep − 1)dp =
R
R∞ 2 p
2ζ(3), and 0 p /(e + 1)dp = 3ζ(3)/2, in terms or Riemann’s Zeta function. Remember also that ζ(3) ≈
1.202.
 12 FREEZE OUT OF MASSIVE SPECIES

EXAMPLE 2.1

It is easy to estimate the value of the Yield that we need in order to reproduce the
correct DM relic abundance, Ωh2 ≈ 0.1, since
ρχ 2 mχ nχh2 mχ Y∞ s0 h2
Ωh2 = h = = , (2.18)
ρc ρc ρc
where Y∞ corresponds to the DM Yield today and s0 is todays entropy density. We
can assume that the Yield did not change since DM freeze-out and therefore
mχ Yf s0 h2
Ωh2 = . (2.19)
ρc
Using the measured value s0 = 2970 cm−3 , and the value of the critical density
ρc = 1.054×10−5 h2 GeV cm−3 , as well as Plancks result on the DM relic abundance,
Ωh2 ≈ 0.1, we arrive at
 
−10 1 GeV
Yf ≈ 3.55 × 10 . (2.20)
mχ
In Figure 2.1 represent the yield as a function of x for non-relativistic particles, us-
ing expression (2.17). As we can observe, the above range of viable values for Yf
correspond to xf ≈ 20.
Notice that this is a crude approximation and we will soon be making a more careful
quantitative treatment.

2.2 Time evolution of the number density

The evolution of the number density operator can be computed by applying the covariant
form of Liuvilles operator to the corresponding phase space distribution function. Formally
speaking, we have
L̂[f ] = C[f ], (2.21)
where L̂ is the Liouville operator, defined as
∂ ∂
L̂ = pµ µ
− Γµσρ pσ pρ µ , (2.22)
∂x ∂p
and C[f ] is the collisional operator, which takes into account processes which change the
number of particles (e.g., annihilations or decays). In the expression above, gravity enters
through the affine connection, Γµσρ .
One can show that in the case of a FRW Universe, for which f (xµ , pµ ) = f (t, E), we
have
∂ ∂
L̂ = E − Γ0σρ pσ pρ
∂t ∂E
∂ 2 ∂
= E − H|p| . (2.23)
∂t ∂E
Integrating over the phase space we can relate this to the time evolution of the number
density
Z Z
g L̂[f ] 3 g C[f ] 3
d p = d p, (2.24)
(2π)3 E (2π)3 E
 TIME EVOLUTION OF THE NUMBER DENSITY 13

Figure 2.1 Equilibrium yield as a function of the dimensionless variable, x, for non-relativistic
particles. The green band represents the freeze-out value, Yf , for which the correct thermal relic
abundance is achieved (for masses of order 1-1000 GeV.

EXAMPLE 2.2

We can show that

Z
g L̂[f ] 3 dn
d p= + 3Hn . (2.25)
(2π)3 E dt

Regarding the collisional operator, it encodes the microphysical description in terms of

Particle Physics, and incorporates all number-changing processes that create or deplete
particles in the thermal bath. For simplicity, let us concentrate in annihilation processes,
where SM particles (A, B) can annihilate to form a pair of DM particles (labelled 1, 2), or
vice-versa (A, B ↔ 1, 2). The phase space corresponding to each particle is defined as
gi d3 pi
dΠi = , (2.26)
(2π)3 2Ei
from where
Z Z
g C[f ] 3
d p = − dΠA dΠB dΠ1 dΠ2 (2π)4 δ(pA + pB − p1 − p2 )
(2π)3 E
|M12→AB |2 f1 f2 (1 ± fA )(1 ± fB ) − |MAB→12 |2 fA fB (1 ± f1 )(a ± f2 )
 
Z
= − dΠA dΠB dΠ1 dΠ2 (2π)4 δ(pA + pB − p1 − p2 )

|M12→AB |2 f1 f2 − |MAB→12 |2 fA fB .
 
(2.27)
14 FREEZE OUT OF MASSIVE SPECIES

The terms (1 ± fi ) account for the viable phase space of the produced particles, taking
into account whether they are fermions (−) or bosons (+). Assuming no CP violation
in the DM sector (T invariance) |M12→AB |2 = |MAB→12 |2 ≡ |M|2 . Also, energy
conservation in the annihilation process allows us to write EA + EB = E1 + E2 , thus,
EA +EB E1 +E1
fA fB = fAeq fBeq = e− T = e− T = f1eq f2eq . (2.28)

In the first equality we have just used the fact that SM particles are in equilibrium. This
eventually leads to
Z
g C[f ] 3
d p = −hσvi n2 − n2eq ,


3
(2.29)
(2π) E
where we have defined the thermally-averaged cross-section as
Z
1
hσvi ≡ 2 dΠA dΠB dΠ1 dΠ2 (2π)4 δ(pA + pB − p1 − p2 )|M|2 f1eq f2eq . (2.30)
neq

Collider enthusiasts would realise that this expression is similar to that of a cross-section,
but we have to consider that the “initial conditions” do not correspond to a well-defined
energy, but rather we have to integrate to the possible energies that the particles in the
thermal bath may have. This explains the extra integrals in the phase space of incident
particles with a distribution function given by f1eq f2eq . We are thus left with the familiar
form of Boltzmann equation,
dn
+ 3Hn = −hσvi n2 − n2eq .


(2.31)
dt
Notice that this is an equilibrium-restoring equation. If the right-hand-side of the equation
dominates, then n traces its equilibrium value n ≈ neq . However, when Hn > hσvin2 ,
then the right-hand-side can be neglected and the resulting differential equation dn/n =
−3da/a implies that n ∝ a−3 . This is equivalent to saying that DM particles do not
annihilate anymore and their number density decreases only because the scale factor of the
Universe increases.
It is also customary to define the dimensionless variable 3
m
x= . (2.32)
T

EXAMPLE 2.3

Using the yield defined in equation (2.9) we can simplify Boltzmann equation. Notice
that
d a3 n
     
dY d n 1 2 3 dn 1 dn
= = = 3a ȧn + a = 3Hn + . (2.33)
dt dt s dt a3 s a3 s dt s dt

Here we have used that the expansion of the Universe is iso-entropic and thus a3 s
remains constant. Also we use the definition of the Hubble parameter H = aȧ . This

3 It
is important to point that this definition of x is not universal; some authors use T /m and care should be taken
when comparing results from different sources in the literature.
 TIME EVOLUTION OF THE NUMBER DENSITY 15

allows us to rewrite Boltzmann equation as follows

dY
= −shσvi Y 2 − Yeq
2


. (2.34)
dt
Now, since a ∝ T −1 and s ∝ T 3 ,
d 3 d d a
(a s) = 0 → (aT ) = 0 → =0, (2.35)
dt dt dt x
which in turns leads to
dx
= Hx , (2.36)
dt
and thus
dY dY dx dY
= = Hx . (2.37)
dt dx dt dx
Using the results of Example (2.3) we can express Boltzmann equation (2.31) as
dY −sxhσvi
Y 2 − Yeq2


=
dx H(m)
−λhσvi
Y 2 − Yeq
2


= 2
, (2.38)
x
where we have used the expression of the entropy density (2.10) in the last line and defined
2π 2 MP g∗s
λ ≡ m
45 1.66 g∗1/2
g∗s
≈ 0.26 1/2 MP m . (2.39)
g∗
Eq. (2.38) is a Riccati equation, without closed analytical form. Thus, to calculate its
solutions we have to rely on numerical methods. However, it is possible to solve it approx-
imately.

2.2.1 Freeze out of relativistic species

The freeze-out of relativistic species is easy to compute, since the yield (2.15) has no
dependence on xf . Neutrinos are a paradigmatic example of relativistic particles and one
must in principle consider their contribution to the total amount of dark matter (after all,
they are dark).
Since neutrinos decouple while they are still relativistic, their yield reads
gef f
Yeq ≈ 0.278 . (2.40)
g∗s
Neutrinos decouple at a few MeV, when the species that were still relativistic are e± , γ, ν
and ν̄. Thus, the number of relativistic degrees of freedom is g∗ = g∗s = 10.75. For one
neutrino family, the effective number of degrees of freedom is gef f = 3g/4 = 3/2. Using
these values, the relic density today an be written as
2
P
i mνi Y∞ s0 h
Ωh2 =
ρc
P
i mνi
≈ . (2.41)
91 eV
16 FREEZE OUT OF MASSIVE SPECIES

P
Notice that in order for neutrinos to be the bulk of dark matter, we would need i mνi ≈
9 eV , which is much bigger than current upper limits (for example, obtainedP from cos-
mological observations). Notice, indeed, that if we consider the current bound i mνi ≤
0.3 eV we can quantify the contribution of neutrinos to the total amount of dark matter,
resulting in Ωh2 ≤ 0.003. This is less than a 3% of the total dark matter density.

2.2.2 Freeze out of non-relativistic species

We can define the quantity

∆Y ≡ Y − Yeq . (2.42)

Boltzmann equation (2.38) is now easier to solve, at least approximately, as follows

For early times, 1 < x  xf , the yield follows closely its equilibrium value, Y ≈
Yeq , and we can assume that d∆Y /dx = 0. We then find
dYeq
dx x2
∆Y = − . (2.43)
Yeq 2λhσvi
Thus, at freeze-out we obtain
x2f
∆ Yf ≈ , (2.44)
2λhσvi
where in the last line we have used that for large enough x, using eq. (2.17) implies
dYeq
dx ≈ −Yeq .

For late times, x  xf , we can assume that Y  Yeq , and thus ∆Y∞ ≈ Y∞ , leading
to the following expression,
d∆Y λhσvi
≈ − 2 ∆2Y , (2.45)
dx x
This is a separable equation that we integrate from the freeze-out time up to nowa-
days. In doing so, it is customary to expand the thermally averaged annihilation cross
section in powers of x−1 as hσvi = a + xb .
Z ∆Y ∞ Z x∞
d∆Y λhσvi
=− dx . (2.46)
∆Yf ∆2Y xf x2

Taking into account that x∞  xf , this leads to

 
1 1 λ b
= + a+ . (2.47)
∆ Y∞ ∆ Yf xf 2xf

The term 1/∆Yf is generally ignored (if we are only aiming at a precision up to a few
per cent [17]) . We can check that this is a good approximation using the previously
derived (2.44) for xf ≈ 20 (which, as we saw in Fig. 2.1 is the value for which the
equilibrium Yield has the right value). This leads to
xf
∆Y∞ = Y∞ =  . (2.48)
b
λ a+ 2xf
 COMPUTING THE DM ANNIHILATION CROSS SECTION 17

The relic density can now be expressed in terms of this result as follows

mχ Y∞ s0 h2
Ωh2 =
ρc
10−10 GeV−2
≈ b
a + 40
3 × 10−27 cm3 s−1
≈ b
. (2.49)
a + 40

This expression explicitly shows that for larger values of the annihilation cross sec-
tion, smaller values of the relic density are obtained.

2.2.3 WIMPs
Equation (2.49) implies that in order to reproduce the correct relic abundance, dark matter
particles must have a thermally averaged annihilation cross section (from now on we will
shorten this to simply annihilation cross section when referring to hσvi) of the order of
hσvi ≈ 3 × 10−26 cm3 s−1 .
We can now consider a simple case in which dark matter particles self-annihilate into
Standard Model ones through the exchange (e.g., in an s-channel) of a gauge boson. It is
easy to see that if the annihilation cross section is of order hσvi ∼ G2F m2W IM P , where
GF = 1.16 × 10−5 GeV−2 , then the correct relic density is obtained for masses of the
order of ∼ GeV.

2.3 Computing the DM annihilation cross section

In the previous sections we have derived a relation between the thermally averaged annihi-
lation cross section and the corresponding dark matter relic abundance. This is very useful,
since it provides an explicit link with particle physics. A central point in that calculation
was the expansion in velocities of the thermally averaged annihilation cross section.

3 b0 15 c
hσvi = ha + bv 2 + cv 4 + . . .i = a + + + ... . (2.50)
2x 8 x2
Notice that in the expressions of the previous section we have defined b ≡ 3b0 /2. As we
also mentioned before, DM candidates tend to decouple when xf ≈ 20. For this value, the
rms velocity of the particles is about c/4, thus corrections of order x−1 can in general not
be ignored (they can be of order 5 − 10%). Moreover, some selection rules can actually
lead to a = 0 for some particular annihilation channels and in that case hσvi is purely
velocity-dependent.
It is important to define correctly the relative velocity that enters the above equation. In
Ref. [17] an explicitly Lorentz-invariant formalism is introduced where

d3 p1
Z Z
g1 C[f1 ] 3 = − hσviMøl (dn1 dn2 − dneq eq
1 dn2 ) , (2.51)
2π E1

where hσviMøl n1 n2 is invariant under Lorentz transformations and equals vlab n1,lab n2,lab
in the rest frame of one of the incoming particles. In our case the densities and Møller
18 FREEZE OUT OF MASSIVE SPECIES

velocity refer to the cosmic comoving frame. In terms of the particle velocities ~vi = p~i /Ei ,
1/2
vMøl = |~v1 − ~v2 |2 + |~v1 × ~v2 |2

. (2.52)
The thermally-averaged product of the dark matter pair-annihilation cross section and their
relative velocity hσvMøl i is most properly defined in terms of separate thermal baths for
both annihilating particles [17, 18],

d p1 d3 p2 σvMøl e−E1 /T e−E2 /T

R 3
hσvMøl i(T ) = R , (2.53)
d3 p1 d3 p2 e−E1 /T e−E2 /T

where p1 = (E1 , p1 ) and p2 = (E2 , p2 ) are the 4-momenta of the two colliding particles,
and T is the temperature of the bath. The above expression can be reduced to a one-
dimensional integral which can be written in a Lorentz-invariant form as [17]
Z ∞ √ 
1 2
√ s
hσvMøl i(T ) = 4 2 ds σ(s)(s − 4mχ ) sK1 , (2.54)
8mχ T K2 (mχ /T ) 4m2χ T

where s = (p1 + p2 )2 and Ki denote the modified Bessel function of order i. In comput-
ing the relic abundance [19] one first evaluates eq. (2.54) and then uses this to solve the
Boltzmann equation. The freeze out temperature can be computed by solving iteratively
the equation  r 
mχ 45 −1/2
xf = ln hσv Møl i(x f ) xf (2.55)
2π 3 2g∗ GN
√
where g∗ represents the effective number of degrees of freedom at freeze-out ( g∗ ≈ 9).
As explained in the previous section, one finds that the freeze-out point xf ≡ mχ /Tf is
approximately xf ∼ 20.
The procedure can be simplified if we consider that the annihilation cross section can
be expanded in plane waves. For example, consider the dark matter annihilation process
χχ → ij and assume that the thermally averaged annihilation cross section can be ex-
pressed as hσviij ≈ aij + bij x. It can then be shown that the coefficients aij and bij can
be computed from the corresponding matrix element. For example,

1 1
 Z 
1 Nc 2
aij = 2 β(s, mi , mj ) d cos θCM |Mχχ→ij | , (2.56)
mχ 32π 2 −1 s=4m2 χ

where θCM denotes the scattering angle in the CM frame, N c = 3 for q̄q final states and 1
otherwise, and
1/2  1/2
(mi + mj )2 (mi − mj )2

β(s, mi , mj ) = 1− 1− (2.57)
s s
The contribution for each final state is calculated separately.

2.3.1 Special cases

The derivation of equation (2.49) relied on the expansion of hσvi in terms of plane waves.
This expansion can be done when hσvi varies slowly with the energy (we can express this
in terms of the centre of mass energy s). However, there are some special cases in which
this does not happen and which deserve further attention.
 COMPUTING THE DM ANNIHILATION CROSS SECTION 19

Figure 2.2 Relativistic thermal average near a threshold (thick solid line) compared to the result
fro the expansion in powers of x−1 (thin line). Figure from Ref. [17].

Annihilation thresholds
A new annihilation channel χ + χ → A + B opens up when 2mχ ≈ mA + mB . In
this case the expansion in velocities of hσvi diverges (at the threshold energy) and it
is no longer a good approximation [17]. Notice in particular that below the threshold,
the expression of aij in Equation (2.56) is equal to zero (as it is only evaluated for
s > 4m2χ ). A qualitative way of understanding this is of course that DM particles have
a small velocity, which is here approximated to zero. In the limit of zero velocity, the
total energy available is determined by the DM mass.
However, we are here ignoring that a fraction of DM particles (given by their thermal
distribution in the Early Universe) have a kinetic energy sufficient to annihilate into
heavier particles (above the threshold). In other words, hσvi is different from zero
below the corresponding thresholds. A very good illustration of this effect is shown
in Ref. [17] and is here reproduced in Fig. 2.2.
The thin solid line corresponds to the approximate expansion in velocities and shows
that not only hσvi is zero below the threshold, but also diverges at the threshold,
thereby not leading to a good solution. Expression (2.54), represented by a thick solid
line, still provides a good solution .

Resonances
The annihilation cross section is not a smooth function of s in the vicinity of an s-
channel resonance. Thus, the velocity expansion of hσvi will fail (although once
more, expression (2.54) still provides a good solution). For a Breit-Wigner resonance
(due to a particle φ) we have

4πw m2φ Γ2φ

σ= B i B f , (2.58)
p2 (s − m2φ )2 + m2φ Γ2φ
20 FREEZE OUT OF MASSIVE SPECIES

in terms of the centre of mass momentum p = 1/2(s − 4m2 )1/2 and the statistical
factor w = (2J + 1)/(2S + 1)2 . The quantities Bi,f correspond to the branching
fractions of the resonance into the initial and final channel.
We can define the kinetic energy per unit mass in the lab frame, , as

(E1,lab − m) + (E2,lab − m) 2 − 4m2

= = , (2.59)
2m 4m2
and rewrite the expression for σ in the lab frame (we want to use Equation (3.21)
in Ref. [17] to compute hσvMøl i). Summing to all final states, and using vlab =
21/2 (1 + )1/2 /(1 + 2), we obtain

8πw γφ2
σvlab = b φ () , (2.60)
m2 ( − 2φ )2 + γφ2

with the definitions b() = Bi (1 − Bi )(1 + )1/2 /(1/2 (1 + 2), γφ = mφ Γφ /4m2 ,

and φ = (m2φ − 4m2 )/4m2 .
It can be shown that in the case of a very narrow resonance, γφ  1, the expression
above can be approximated as

8πw
σvlab = bφ ()πγφ δ( − φ ) , (2.61)
m2
the relativistic formula for the thermal average then reads [17]

16πw x 1/2 p
hσvMøl i = 2 2 πγφ φ (1 + 2eφ )K1 (2x 1 + φ )bφ (eφ )θ(φ ) . (2.62)
m K2 (x)

Notice that φ > 0 when m < 2mφ , i.e., when the mass of the DM is not enough
to enter the resonance. The reason is easy to understand. Only through the extra
kinetic energy provided by the thermal bath, the resonance condition can be satisfied.
However, when the mass of the DM exceeds the resonance condition, the kinetic
energy only takes us further away from the resonant condition and the thermalised
cross section tends to vanish. In other words, the centre of mass rest energy exceeds
mφ /2. This can be seen in Figure 2.3.
For a large width the expression has to be computed numerically and can be found in
[17].

Coannihilations
When deriving Boltzmann equation (2.31) we have only considered one exotic species,
but this needs not be the case. In fact, in most particle models for DM, there are more
exotic species that we need to take into account. Notice that, in principle, this would
lead to a system of coupled Boltzmann equations. If we label exotic species as χi ,
with i = 0, 1 . . . k, and SM particles as A, B, we have to consider all number chang-
ing processes for each species,

(i) χi + χj → A + B
(ii) χi + A → χj + B
(iii) χj → χi + A
 COMPUTING THE DM ANNIHILATION CROSS SECTION 21

Figure 2.3 Relativistic thermal average in a resonance (thick solid line) compared to the result fro
the expansion in powers of x−1 (thin line). Figure from Ref. [17].

If we consider the (usual) case in which the DM is protected by a symmetry (e.g.,

in the case of Supersymmetric theories) and that the exotic particles all must decay
eventually into the lightest one χ0 , then, we must only trace the evolution of the total
Pk
number density of exotic species, n = i=0 ni . Under this assumption, processes
(ii) and (iii) do not need to be considered, as they do not change the number of exotics.
This is correct as long as the rate of these is faster than the expansion of the Universe.
Regarding process (i) we have to be aware that the cross section σij is going to appear
multiplied by the corresponding number densities, ni nj . Now, we are considering
the case in which both particles i and j are non-relativistic and as a consequence, ni,j
are Boltzmann suppressed, ni,j /e−mi,j /T . Thus, unless mj ≈ mi , the abundance of
χj is negligible and only the process χi + χj → A + B is important (and we are back
to the case of a single exotic).
However, when mj ≈ mi , there can be coannihilation effects and particle j may serve
as a channel through which particles i can be more effectively depleted. This is the
case, e.g., of the stau and the neutralino in supersymmetric theories.
CHAPTER 3

DIRECT DM DETECTION

3.1 Preliminaries

3.1.1 DM flux
We can easily estimate the flux of DM particles through the Earth. The DM typical velocity
is of the order of 300 km s−1 ∼ 10−3 c. Also, the local DM density is ρ0 = 0.3 GeV cm−3 ,
thus, the DM number density is n = ρ/m.
vρ 107
φ= ≈ cm−2 s−1 (3.1)
m m
These particles interact very weakly with SM particles.
Assuming a typical WIMP cross section σ

3.1.2 Kinematics

Direct DM detection is based on the search of the scattering between DM particles and
nuclei in a detector. This process is obly observable through the recoiling nucleus, with an
energy ER . DM particles move at non-relativistic speeds in the DM halo. Thus, the dy-
namics of their elastic scattering off nuclei are easily calculated. In particular, the recoiling
energy of the nucleus is given by
1 4mχ mN 1 + cos θ
ER = mχ v 2 (3.2)
2 (mχ + mN )2 2
Dark Stuff. 23
By D. G. Cerdeño, IPPP, University of Durham
24 DIRECT DM DETECTION

It can be checked that for DM particles with a mass of the order of 100 GeV, this leads
to recoil energies of approximately ER ∼ 100 keV. Notice also that the maximal energy
transfer occurs on a head-on-collision and when the DM mass is equal to the target mass.
In such a case
1 1 1  mχ 
ERmax
= mχ v 2 = mχ × 10−6 = keV (3.3)
2 2 2 1 GeV
where we have used that in a DM halo the typical velocity is v ∼ 10−3 c.
Experiments must therefore be very sensitive and be able to remove an overwhelming
background of ordinary processes which lead to nuclear recoils of the same energies.

3.2 The master formula for direct DM detection

The total number of detected DM particles, N , can be understood as the product of the DM
flux (which is equal to the DM number density, n, times its speed, v), times the effective
area of the target (i.e., the number of targets NT times the scattering cross-section, σ), all
of this multiplied by the observation time, t,

N = t n v NT σ . (3.4)

We will be interested in determining the spectrum of DM recoils, i.e., the energy depen-
dence of the number of detected DM particles. Thus,
dN dσ
= t n v NT . (3.5)
dER dER
Now, the DM velocity is not unique, and in fact DM particles are described by a local
velocity distribution, f (~v ), where ~v is the DM velocity in the reference frame of the detec-
tor. We therefore have to integrate to all possible DM velocities, with their corresponding
probability density, Z
dN dσ
= t n NT vf (~v ) d~v , (3.6)
dER vmin dE R

where q
vmin = mχ ER /2µ2χN (3.7)
is the minimum speed necessary to produce a DM recoil of energy ER , in terms of the
WIMP-nucleus reduced mass, µχN . Using n = ρ/mχ and NT = MT /mN (where MT
is the total detector mass and mN is the mass of the target nuclei), and defining the exper-
imental exposure  = t MT , we arrive at the usual expression for the DM detection rate
Z
dN ρ dσ
= vf (~v ) d~v . (3.8)
dER mχ mN vmin dER

3.2.1 The scattering cross section

The scattering takes place in the non-relativistic limit. The cross section is therefore ap-
proximately isotropic (angular terms being suppressed by v 2 /c2 ∼ 10−6 . This implies
that
dσ σ
∗
= constant = (3.9)
d cos θ 2
 THE MASTER FORMULA FOR DIRECT DM DETECTION 25

On the other hand,

max 1 + cos θ∗ dER max

ER
ER = ER → = (3.10)
2 d cos θ∗ 2
From this, we can see that
dσ dσ d cos θ∗ σ mN σ
= ∗
= max = (3.11)
dER d cos θ dER ER 2µχN v 2
Notice finally that the momentum transfer from WIMP interactions reads (remember
that we are considering non-relativistic processes and thus we neglect the kinetic energy of
the nucleus) p
q = 2 mN ER (3.12)
and is typically of the order of the MeV. The equivalent de Broglie length would be λ ∼
2π~/p ∼ 10 − 100 fm. For light nuclei, the DM particle sees the nucleus as a whole,
without substructure, only for heavier nuclei we have to take into account a suppression
form factor. The nuclear form factor, F 2 (ER ), accounts for the loss of coherence
dσ mN σ0 2
= F (ER ) (3.13)
dER 2µχN v 2
Finally, the scattering cross section receives different contributions, depending on the mi-
croscopic description of the interaction.
In the end, we can
Z
dN ρ f (~v )
= σ0 F 2 (ER ) d~v . (3.14)
dER 2 mχ µχN vmin v
The inverse mean velocity
Z
f (~v )
η(vmin ) = d~v . (3.15)
vmin v
is the main Astrophysical input.

3.2.2 The importance of the threshold

From the kinematics of the DM-nucleus interaction, we see that, for a given recoil energy
ER , we require a minimal velocity of DM particles, given by expression (3.7).
Thus, given that experiments are only sensitive to DM interactions above a certain en-
ergy threshold, ET , this means that we are only probing a part of the WIMP velocity
distribution function (for a given DM mass). Conversely, given that DM particles have a
maximum velocity in the halo (otherwise they become unbound and escape the galaxy),
the experimental energy threshold is a limitation to explore low-mass WIMPs.

EXAMPLE 3.1

Consider a germanium experiment and a xenon experiment with a threshold of 2 keV.

Given the escape velocity in a typical isothermal halo, vesc = 554 km s−1 , determine
the minimum DM mass that these experiments can probe.
This is the reason that experiments loose sensitivity for small masses.
26 DIRECT DM DETECTION

3.2.3 Velocity distribution function

It is customary to consider the Isothermal Spherical Halo, which assumes that the Milky
Way (MW) halo is an isotropic, isothermal sphere with density profile ρ ∝ r−2 . The
velocity distribution, in the galactic rest frame, for such a halo reads
1 v |2
|~
fgal (~v ) = 3/2
e− 2σ2 , (3.16)
(2πσ)

where the√one-dimensional velocity dispersion, σ, is related to the circular speed, vc , as

σ = vc / 2. The canonical values are vc = 220km s−1 , with a statistical error of order
10% (see references in [20])
Now, in order to use it for direct detection experiments we need to carry out a Galilean
transformation ~v → ~v + ~vE , such that

f (~v ) = fgal (~v + ~vE (t)). (3.17)

where ~vE (t) is the velocity of the Earth with respect to the Galactocentric rest frame.

~vE (t) = ~vLRF + ~v + ~vorbit (t) (3.18)

Notice that vE includes contributions from the speed of the Local Standard of Rest vLSR ,
the peculiar velocity of the Sun with respect to vLSR , and the Earths velocity around the
Sun, which has an explicit time dependence.
Notice that if we work with the SHM, the angular integration in the computation of
direct detection rates can be easily done as follows
Z Z Z Z
f (~v ) 3 1 v |2 +|~
|~ v E |2 |~
v | |~
vE | cos θ
d v = dφ d cos θ dv v 2 3/2
e− 2σ2 e σ2
v (2πσ )
2σ 2
 
|~v | |~vE |
Z
v |2 +|~
|~ v E |2
−
= 2π dv v e 2σ 2 sinh
|v||~vE |(2πσ)3/2 σ2
√  
|~v | |~vE |
Z
2 v |2 +|~
|~ v E |2
= dv √ e− 2σ2 sinh (3.19)
πσ|~vE | σ2

3.2.4 Energy resoultion

3.3 Exponential spectrum

η(vmin ) ∝ e−αER (3.20)

3.4 Annual modulation

vE (t) = v0 [1.05 + 0.07 cos(2π(t − tp )/1yr)] (3.21)

The amplitude of the modulation in the velocity is very small, approximately a 7% (vorbit cos θorbit /(vLRF +
v ) ≈ 15/220
This eventually implies a modulation in η(vmin ), which means that the amount of DM
particles in the tail of the distribution can change substantially. This, in turn, leads to a
modulation in the detected rate.
 DIRECTIONAL DETECTION 27

3.5 Directional detection

The current experimental situation regarding directional DM detection has been sum-
marised in the recent review article, Ref. [20].

3.6 Coherent neutrino scattering

Solar neutrinos might leave a signal in DD experiments, either through their coherent scat-
tering with the target nuclei or through scattering with the atomic electrons.
In general, the number of recoils per unit energy can be written
Z
dR  dφν dσν
= dEν , (3.22)
dER mT dEν dER
where  is the exposure and mT is the mass of the target electron or nucleus. If several iso-
topes are present, a weighted average must be performed over their respective abundances.
The SM neutrino-electron scattering cross section is
G2F me

dσνe
= (gv + ga )2 + (3.23)
dER 2π
 2 
ER me ER
(gv − ga )2 1 − + (ga2 − gv2 ) ,
Eν Eν2
where GF is the Fermi constant, and
1 1
gv;µ,τ = 2 sin2 θW − ; ga;µ,τ = − , (3.24)
2 2
for muon and tau neutrinos. In the case νe + e → νe + e, the interference between neutral
and charged current interaction leads to a significant enhancement:
1 1
gv;e = 2 sin2 θW + ; ga;e = + . (3.25)
2 2
The neutrino-nucleus cross section in the SM reads
G2
 
dσνN mN ER
= F Q2v mN 1 − F 2 (ER ), (3.26)
dER 4π 2Eν2
where F 2 (ER ) is the nuclear form factor, for which we have taken the parametrisation
given by Helm [21]. Qv parametrises the coherent interaction with protons (Z) and neu-
trons (N = A − Z) in the nucleus:
Qv = N − (1 − 4 sin θW )Z. (3.27)

3.7 Inelastic

[22]
The WIMP needs to have sufficient speed to interact with the nucleus and promote to
an excited state (with energy separation δ)
1
µχN v 2 > δ (3.28)
2
28 DIRECT DM DETECTION

This leads to the condition

r  
1 mN ER
vmin = +δ (3.29)
2mN ER µχN

Therefore, the main effect at a given experiment is to limit the sensitivity only to a part
of the phase space of the halo. This favours heavy nuclei (since they can transfer more
energy to the outgoing WIMP) and can account for observation in targets such as iodine
(DAMA/LIBRA) while avoiding observation in lighter ones such as Ge (CDMS)
CHAPTER 4

NEUTRINOS

4.1 Preliminaries - copied form internet

The neutrino is a state that is produced in a weak interaction. It is, by definition, a flavour
eigenstate, in the sense that a neutrino is always produced with, or absorbed to give, a
charged lepton of electron, muon or tau flavour. The neutrino that is generated with the
charged electron is the electron neutrino, and so on. However, as with the quarks and the
CKM matrix, it is possible that the flavour eigenstates (states with definite flavour) are not
identical to the mass eigenstates (states which have definite mass).
Mass eigenstates are therefore superpositions of the flavour eigenstates. Denoting α =
e, µ, τ and i = 1, 2, 3, we can write
X
|νii = Uiα |νiα (4.1)
α

If this is the case, the following can happen : suppose one generates a neutrino at a source.
This neutrino will have definite flavour, but will be produced as a linear combination of
states of definite mass. The states of definite mass will propagate out of the source towards
the detector. If the states have different masses, then the phase between the states will
change with distance from the source. At the detector, the mass states will have different
relative phases to those the mass states had at the source, and when we go to detect them,
it is possible that we will detect a flavour state which was not present in the beam to begin
with.

Dark Stuff. 29
By D. G. Cerdeño, IPPP, University of Durham
REFERENCES

[1] J. Silk and Others, Particle Dark Matter: Observations, Models and Searches. 2010.
[2] G. Bertone and D. Hooper, A History of Dark Matter, .
[3] V. C. Rubin and J. Ford, W. Kent, Rotation of the Andromeda Nebula from a Spectroscopic
Survey of Emission Regions, Astrophys. J. 159 (1970) 379.
[4] V. C. Rubin, N. Thonnard, and J. Ford, W. K., Rotational properties of 21 SC galaxies with a
large range of luminosities and radii, from NGC 4605 /R = 4kpc/ to UGC 2885 /R = 122 kpc/,
Astrophys. J. 238 (1980) 471.
[5] A. Bosma, 21-cm line studies of spiral galaxies. II. The distribution and kinematics of neutral
hydrogen in spiral galaxies of various morphological types., Astron. J. 86 (1981) 1825.
[6] T. S. van Albada, J. N. Bahcall, K. Begeman, and R. Sancisi, Distribution of dark matter in the
spiral galaxy NGC 3198, Astrophys. J. 295 (1985) 305.
[7] R. Gavazzi, C. Adami, F. Durret, et al., A weak lensing study of the Coma cluster, Astron.
Astrophys. 498 (2009) L33–L36, [arXiv:0904.0220].
[8] F. Zwicky, Die Rotverschiebung von extragalaktischen Nebeln, Helv. Phys. Acta 6 (1933)
110–127.
[9] F. Zwicky, On the Masses of Nebulae and of Clusters of Nebulae, Astrophys. J. 86 (1937) 217.
[10] J. M. Kubo, A. Stebbins, J. Annis, et al., The Mass Of The Coma Cluster From Weak Lensing
In The Sloan Digital Sky Survey, Astrophys. J. 671 (2007) 1466–1470,
[arXiv:0709.0506].
[11] M. Markevitch, A. H. Gonzalez, D. Clowe, et al., Direct constraints on the dark matter
self-interaction cross-section from the merging galaxy cluster 1E0657-56, Astrophys. J. 606
(2003) 819–824, [astro-ph/0309303].
[12] D. Clowe, M. Bradac, A. H. Gonzalez, et al., A direct empirical proof of the existence of dark
matter, Astrophys. J. 648 (2006) L109–L113, [astro-ph/0608407].
Dark Stuff. 31
By D. G. Cerdeño, IPPP, University of Durham
32 REFERENCES

[13] S. D. McDermott, H.-B. Yu, and K. M. Zurek, Turning off the Lights: How Dark is Dark
Matter?, Phys. Rev. D 83 (2010) 063509, [arXiv:1011.2907].
[14] C. Kouvaris, Composite Millicharged Dark Matter, Phys. Rev. D 88 (2013) 015001,
[arXiv:1304.7476].
[15] A. D. Dolgov, S. L. Dubovsky, G. I. Rubtsov, and I. I. Tkachev, Constraints on millicharged
particles from Planck, Phys. Rev. D 88 (2013) 117701, [arXiv:1310.2376].
[16] E. Del Nobile, M. Nardecchia, and P. Panci, Millicharge or Decay: A Critical Take on
Minimal Dark Matter, arXiv:1512.05353.
[17] P. Gondolo and G. Gelmini, Cosmic abundances of stable particles: Improved analysis, Nucl.
Phys. B 360 (1991) 145–179.
[18] M. Srednicki, R. Watkins, and K. A. Olive, Calculations of relic densities in the early
universe, Nucl. Phys. B 310 (1988) 693–713.
[19] T. Nihei, L. Roszkowski, and R. R. de Austri, Exact Cross Sections for the Neutralino-Slepton
Coannihilation, JHEP 07 (2002) 46, [hep-ph/0206266].
[20] F. Mayet, A. M. Green, J. B. R. Battat, et al., A review of the discovery reach of directional
Dark Matter detection, Phys. Rep. 627 (2016) 1–49, [arXiv:1602.03781].
[21] R. H. Helm, Inelastic and Elastic Scattering of 187-Mev Electrons from Selected Even-Even
Nuclei, Phys. Rev. 104 (1956) 1466–1475.
[22] D. Smith and N. Weiner, Inelastic Dark Matter, Phys. Rev. D 64 (2001) 043502,
[hep-ph/0101138].