# New Directions in Direct Dark Matter Searches

Authors: Paolo Panci
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

I present the status of direct dark matter detection with specific attention to the experimental results and their phenomenological interpretation in terms of dark matter interactions. In particular I review a new and more general approach to study signals in this field based on nonrelativistic operators which parametrize more efficiently the dark matter-nucleus interactions in terms of a very limited number of relevant degrees of freedom. Then I list the major experimental results, pointing out the main uncertainties that affect the theoretical interpretation of the data. Finally, since the underlying theory that describes both the dark matter and the standard model fields is unknown, I address the uncertainties coming from the nature of the interaction. In particular, the phenomenology of a class of models in which the interaction between dark matter particles and target nuclei
is of a long range type is discussed.


---

## Body

## 1. Introduction

One of the most exciting open questions at the interface between particle physics and cosmology is the nature of Dark Matter (DM). The first person who provided evidence and inferred the presence of DM was a Swiss-American astrophysicist, Fritz Zwicky. He applied the virial theorem to the Coma cluster of galaxies and obtained evidence of unseen mass. Roughly 40 years following the discoveries of Zwicky and others, Vera Rubin and collaborators conducted an extensive study of the rotation curves of isolated spiral galaxies. They announced the discovery that the rotational curves of stars in spiral galaxies exhibit a characteristic flat behaviour at large distance in contrast with Kepler’s law. Many other evidence of unseen mass on distance scales of the size of galaxies and clusters of galaxies appeared throughout the years, but the most precise measurement of the total amount of DM comes from cosmological scales. In particular, the measurements of modern precision cosmology (the Cosmic Microwave Background (CMB) and the surveys of the Large Scale Structure (LSS) of the Universe) provide the current most relevant evidence. Apart from the qualitative agreement, it is the quantitative fitting of the wealth of available data that allows the amount of DM to be one of the cosmological parameters now most precisely measured (Ωχh2=0.1199±0.0027; see [1, Table 2]). Therefore, we have compelling evidence of unseen mass, but the microscopic features of this new kind of matter remain unknown yet. Direct and indirect searches may shed light on the nature of DM, and therefore a careful study of their phenomenology is fundamental. For a pedagogical review on this subject, see, for example, [2].Direct searches for DM aim at detecting the nuclear recoils arising from scattering between DM particles and target nuclei in underground detectors. DM direct detection experiments are providing exciting results in terms of measured features which have the right properties to be potentially ascribed to a DM signal. For example in addition to the long-standingDAMA results, nowadays there are other experiments, such as  CoGeNT,  CRESST-II, and  CDMS-Si   that start to see some anomalies in their counting rates. On the other hand, the situation in this field is extremely unclear and confusing, because, on top of these positive result experiments, the constraints coming from null results, like  XENON100,  COUPP,  PICASSO, and very recently  LUX, are very stringent and put the interpretation of the anomalies in terms of a DM interaction in serious trouble. Nevertheless, there are at least two main caveats when the results from the experiments commented upon above are interpreted. The first is that one has to treat with great care the fine experimental details associated with the results quoted by each experiment. The second caveat is instead associated with the interpretation of the data within a very simple-minded DM model. For instance, the DM-nucleus spin independent contact interaction is just a benchmark example. Upon relaxing some of these assumptions, the current complicated experimental puzzle can probably be solved.Thescope of this work is to present thestatus of direct DM detection with specific attention to the experimental results and their phenomenological interpretation in terms of DM interaction. In particular in Section 2, I review a new and more general approach to study signals in this field based on nonrelativistic operators which parameterize most efficiently the DM-nucleus interactions in terms of a very limited number of relevant degrees of freedom. In Section 3, I review the experimental results and their interpretation in terms of the “standard” spin independent (SI) interaction. I list then the main uncertainties that affect the theoretical interpretation of the data: this is a very promising area of research since only major advancements here can probably reconcile the complicated puzzle showed by the experiments up. Finally in Section 4, I pose my attention on the uncertainties coming from the nature of the interaction. In particular the phenomenology of a class of models in which the interaction between DM particles and target nuclei is of a long range type is discussed.


## 2. Basics of Direct Detection Computations

### 2.1. Kinematics

As already stated, when a DM particle scatters off a nucleus, depending on the DM properties, one can envision at least two distinct kinematics, the elastic and the inelastic one. The elastic scattering is represented by(1)χ+𝒩(A𝒩,Z𝒩)atrest⟶χ+𝒩(A𝒩,Z𝒩)recoil,
while the inelastic is
(2)χ+𝒩(A𝒩,Z𝒩)atrest⟶χ′+𝒩(A𝒩,Z𝒩)recoil.
In (1) and (2), χ and χ′ are two DM particle states, and A𝒩, Z𝒩 are, respectively, the mass and atomic numbers of the nucleus 𝒩.We know, thanks to rotational curves data, that the velocity ofχ in the vicinity of the Earth is of the order of 10-3c. Hence the scattering between a DM particle with velocity v and mass mχ and a nucleus at rest with mass m𝒩 occurs in deeply nonrelativistic regime. This is analogous with the collision between two billiard balls and therefore the recoil energy can be simply obtained by considering energy and momenta conservation separately. In the detector rest frame, such scattering would end up into a nuclear recoil of energy:
(3)ER=2μχ𝒩2m𝒩v2(1-(vt2/2v2)-1-(vt2/v2)cos⁡θ2),vt=0,elasticscattering,vt≠0,inelasticscattering,
where μχ𝒩=mχm𝒩/(mχ+m𝒩) is the DM-nucleus reduced mass, θ as usual is the scattering angle in the center-of-mass frame, and vt=2δ/μχ𝒩 is a threshold velocity. Here δ=mχ′-mχ is the mass splitting between χ and χ′, and the equation above holds for δ≪mχ′, mχ. Elastic scattering occurs for δ=0, while δ≠0 implies inelastic scattering. The minimal DM velocity providing the ER recoil energy can be then obtained by putting θ=-π in (3). One gets
(4)vmin⁡(ER)=m𝒩ER2μχ𝒩2(1+μχ𝒩δm𝒩ER).


### 2.2. Formalism of Nonrelativistic Operators

Having a disposal the basic quantityvmin⁡ that fully describes the kinematics of the DM-nucleus scatterings in (1) and (2) and the differential rate of nuclear recoil expected in a given detector can be achieved by weighting the differential cross section dσ𝒩/dER with the DM velocity distribution fE(v→) in the velocity range allowed by the kinematics. It reads as follows:
(5)dR𝒩dER=ξ𝒩m𝒩ρ⊙mχ∫vmin⁡(ER)vesc‍d3vvfE(v→)dσ𝒩dER(v,ER),
where ξ𝒩 are the mass fractions of different nuclides (ξ𝒩=103NAm𝒩ζ𝒩/kgA¯, where NA=6.022×1023 is Avogadro’s number, ζ𝒩 are the numeric abundances and A¯≡∑𝒩‍ζ𝒩A𝒩.), and ρ⊙≃0.3GeV/cm3 is the DM energy density at the Earth’s location. This is the canonical value routinely adopted in the literature (see, e.g., [2–4]). Recent computations have found a higher central value of it, still subject for some debate [5–7]. The integral in the right-hand side of the equation above is a key ingredient because it contains all the information related to the geometry of the DM halo, the nature of the DM-nucleus interaction, and in turn the nuclear response of the target. In the following we present the formalism of nonrelativistic (NR) operators and then we describe how to write down the main observables in terms of it.(i)
As already stated, the weight of the velocity integral is the DM velocity distributionfE(v→). In the rest frame of our Galaxy it can be roughly approximated with a Maxwell-Boltzmann distribution due to the fact, as pointed out in [8], that the violent relaxation of the gravitational potential at the formation of the Milky Way leads to fast mixing of the DM phase space elements. Therefore, DM particles were basically frozen in high-entropy configuration which are indeed Maxwell-Boltzmann-like. This has been roughly confirmed by some recent numeral simulation, but as one can see in [9, Figure 2], there are some deviations due to the DM assembly history of the Milky Way. Indeed one can see different features both at low and high velocities that of course are not present in the case of a Maxwell-Boltzmann distribution. Furthermore the geometry of the halo predicted by this kind of numerical simulation is not exactly spherical but tends to a triaxial configuration (see, e.g., [10]).Keeping in mind that this is not the truth because we have not directly measured the properties of the DM halo yet, what we know for sure is that, since we are setting in the Earth laboratory frame, one has to boost the DM-nucleus relative velocity with the drift velocity of the Sun around the galactic center and the Earth’s orbital velocity projected in the galactic plane. Therefore the DM velocity distributionfE(v→) in the Earth frame is related to the DM velocity distribution in the galactic one fG(w→) by the Galilean velocity transformation fE(v→)=fG(v→+v→E(t)), where
(6)v→E(t)=v→⊙+v→⊕(t).Herev→⊙=v→G+v→S is the sum of the galactic rotational velocity of the local system v→G=(0,v0,0) (here v0≃220±50km/s is the local DM velocity) and the Sun’s proper motion v→S≃(10,13,7) km/s with absolute value v⊙=233±50km/s, while v→⊕(t) is the time dependent Earth’s orbital velocity with period of one year, phase around June 2nd (when it is aligned to v→⊙), and size v⊕≃30km/s, and it is inclined of an angle γ≃60∘ with respect to the galactic plane. More details can be found in [11–13]. Since in (6) the Earth’s orbital velocity projected in the galactic plane v⊕cos⁡γ is relatively small compared to v⊙, we can approximate v→E(t) with its component directed toward the galactic center. We can then write
(7)vE(t)≃v⊙+Δvcos⁡[2π(t-ϕ)τ],whereΔv=v⊕cos⁡⁡γ≃15km/s. Here ϕ≃152.5 days (June 2nd) is the phase and τ≃365 days is the period of the Earth motion around the Sun.Therefore, in general, it is expected an annual modulation in the counting rate of direct detection experiment, as the incoming flux of DM particles toward us contains small oscillatory terms due to the Earth’s proper motion around the Sun during the years. More details on this model independent signature can be found, for instance, in [14, 15]. It is worth stressing that looking for an annual modulation of the counting rate in direct detection is very challenging from the experimental point of view, simply because the modulus of the time dependent component of the velocity is relatively small compared to the constant one; in better words, the size of the modulated signal that experiments like  DAMA and  CoGeNT are looking for is suppressed with respect to the unmodulated one, due basically to the collective motion of the Solar system around the galactic center.(ii)
In the nonrelativistic limit the differential cross section can be written in the usual form:(8)dσ𝒩dER(v,ER)=132π1mχ2m𝒩1v2|ℳ𝒩|2¯,where|ℳ𝒩|2¯ is the DM-nucleus matrix element that encodes all the information related to the nuclear physics and the nature of the interactions. Since, as already stated, the velocity of the DM particles in the vicinity of the Earth is much smaller than the speed of light, the framework of relativistic quantum field theory is not too much appropriate to study signals in this field, especially if we do not have any clue about the underlying theory that describes the DM-quark interactions. A more useful and general framework is actually the one based on nonrelativistic operators. Indeed, since, for the elastic scattering, the relevant degrees of freedom are the DM-nucleon relative velocity v→, the exchanged momentum q→, the nucleon spin s→N (N=p, n can be proton or neutron), and the DM one s→χ (if different from zero), the scattering amplitude will then be a rotationally invariant function of this variables. In this regard, with v→, q→, s→N, and s→χ [16], we can construct a basis for 16 operators, which includes all possible spin configurations (the complete list and the numbering of the nonrelativistic operators can be found in [17–19]). The DM-nucleon matrix element can then be written as a linear combination of these operators, with coefficients that may depend on the momenta only through the q2 or v2 scalars (DM models feature a long range interaction with the Standard Model fermions provide perhaps the most notable example. Indeed, in this case, the differential cross section is enhanced at smaller exchanged momenta, due to the negative power dependence of q. A specific realization is offered by models in which DM particles can carry small but nonzero electric charge, electric dipole moment, or magnetic one. Their phenomenology in the context of DM direct searches has been studied in [20–31] and references therein.). In particular it reads
(9)ℳN=∑i=116‍𝔠iN(λ,mχ)𝒪iNR,where the coefficients𝔠iN are functions of the free parameters of the underlying relativistic theory, such as mediator masses, mixing angles, and couplings constants (collectively denoted by λ), and the DM mass mχ. For instance, if the interaction between a fermionic DMχ and the nucleon N is described by the “standard” SI Lagrangian gN/Λ2χ¯χN¯N, the only nonrelativistic operator involved will be the identity (𝒪1NR=𝕝). Before moving on, it is worth stressing that the operators above, which entirely describe the nonrelativistic physics of the DM-nucleus scattering, are many and they can either depend on the exchanged momenta q→ or the relative velocity v→. As a consequence, focussing on just one operator (e.g., 𝒪1NR) is not the most model independent way to study signals in direct DM searches, since theories can predict several operators entering together with the possibility that some of them may even interfere.Since now the nucleus is of course made of protons and neutrons, one has to correct the DM-nucleon matrix element (9) with the nuclear responses, which are a sort of form factors that take account of the finite size of the target nuclei. According to [17, Equation (55)] we can then write the spin-averaged amplitude squared for scattering off a target nucleus 𝒩 with mass m𝒩 as
(10)|ℳ𝒩|2¯=m𝒩2mN2∑i,j=116‍∑N,N′=p,n‍𝔠iN𝔠jN′Fi,j(N,N′).TheFi,j(N,N′)(v,ER,𝒩) are the nuclear responses and they depend critically on the type of scattering nucleus 𝒩: they are also functions of mχ, v, and the nuclear recoil energy ER=q2/2m𝒩. In Figure 1, the nuclear responses for both light element (on the left panel) and heavy one (on the right panel), considering two completely different kinds of interactions, are depicted. On a more specific level, in thick red the total nuclear response for the standard SI interaction is shown. As we can see, in the energy range of interest in direct DM experiments (few keV), it is enhanced by the canonical A𝒩2 factor: this is due to the fact that, in this case, the incoming DM particles see the nucleus as a point-like object with A𝒩 scatter centers. On the other hand, we show in blue the nuclear responses of a completely different type of interaction which is both momenta and spin dependent (𝒪6NR=(s→χ·q→)(s→N·q→)). This type of interaction can give rise if the DM-nucleus interaction is mediated by a pseudoscalar particle. We can see that the nuclear responses are dominated by the scattering with protons, simply because neither the fluorine nor the iodine has unpaired neutrons and moreover the interaction goes to zero in the small momenta transferred limit.A complete set of these nuclear responses for each pair of nucleons(N,N′), each pair of nonrelativistic operators (i,j), and for several target nuclei 𝒩 has been provided in the appendices of [17]. This is extremely useful because, for the first time, different types of interactions can be studied in a more general ground. However, since this kind of computations is quite new, the uncertainties, especially for the spin, momenta, and velocity dependent interactions, are still quite large.Nuclear responses of the fluorine (a) and the iodine (b) targets considering two completely different types of interactions. On a more specific level, we show, in red, the nuclear responses for the “standard” SI contact interaction, while in blue those for another kind of interaction described in the nonrelativistic limit by the operator 𝒪6NR=(s→χ·q→)(s→N·q→). The different hatching refers instead to the possible choices of nucleon pairs in the nucleus. In particular, the solid, dashed, and dotted lines are for (p,p), (n,n), and (p,n) pairs, while the thick ones are obtained by summing over all of them (total nuclear responses). Notice that the total nuclear responses for the “standard” SI contact interaction (thick red lines) reduce to A𝒩2FHelm2(q), where FHelm2(q) is the customary Helm form factor [33].
(a)
(b)Having a disposal the general relation of the DM-nucleus matrix element (10), the differential cross section (8), and in turn the rate of nuclear recoil (5) can be rewritten, following [32], in a very general way. It reads
(11)dR𝒩dER=Xξ𝒩∑i,j=116∑N,N′=p,n𝔠iN(λ,mχ)00000000000000000000×𝔠jN′(λ,mχ)ℱi,j(N,N′)(ER,𝒩),
where we defined
(12)X≡ρ⊙mχ132π1mχ2mN2,(13)ℱi,j(N,N′)(ER,𝒩)≡∫vmin⁡⁡(ER)vescd3v1vfG(v→+v→E(t))Fi,j(N,N′)(v,ER,𝒩).


#### 2.2.1.  “Standard” SI Contact Interaction

To give a concrete example, in the following the coefficient and the nuclear responses for the “standard” SI contact interaction are showed explicitly. In particular, when experiments present results in terms of this interaction, they implicitly assume (probably inspired by supersymmetric neutralino scattering) the following DM-nucleon effective Lagrangian:(14)ℒSIN=λSINχ¯χN¯N.
From this the DM-nucleon matrix element can be obtained by contracting (14) with initial and final states of the scattering. Performing then the nonrelativistic limit one gets
(15)ℳSIN≡〈χN|ℒSIN|χN〉≃4λSINmχmN𝒪1NR,
where 𝒪1NR=𝕝 is the operator that describes the nonrelativistic limit of (14) and the coefficient is obviously given by 𝔠1N(λSIN,mχ)=4λSINmχmN.For the SI scattering, DM-proton and DM-neutron couplings are customarily assumed to be equal and thusλSIp=λSIn≡λSI (isospin violating process is not taken into account). Plugging back then the coefficient in the general equation (11) and defining σSI=λSI2/π·μN2 (μN is the DM-nucleon reduced mass) to make contact with the usual physical cross sections, the differential rate of nuclear recoil reduces to the very well-know compact form:
(16)dR𝒩dER=ξ𝒩m𝒩ρ⊙mχ[m𝒩2μN2σSIF𝒩tot(ER)ℐ(ER)],
where we define the total SI nuclear response as
(17)F𝒩tot(ER)=∑N,N′=p,n‍F1,1(N,N′)(ER,𝒩).
Here, since λSIp=λSIn, the total nuclear response reduces to the square of the customary Helm form factor [33] multiplied by the coherent enhancement of the cross section A𝒩2. In particular it reads F𝒩tot(ER)=A𝒩2FHelm2(ER). The function ℐ(ER)=∫vmin⁡vesc‍d3vfG(v→+v→E(t))/v is instead the velocity integral encountered many times in the literature. It is worth noticing that ℐ(ER) can be written in this way because in the “standard” SI interaction the operator that the describes the nonrelativistic limit of the effective Lagrangian does not carry any dependences on the relative velocity v→. In Section 3.2 the interpretation of the experimental results in terms of the “standard” SI interaction is briefly reviewed.


### 2.3.  Experimental Observables

Since, as already stated, the Earth’s orbital velocity projected in the galactic plane is relatively small compared to the drift velocity of the Sun, we can then expand the recoil rate (11), assuming that the velocity distribution is not strongly anisotropic. Then, by means then of (7), one gets
(18)dR𝒩dER(t)≃dR𝒩dER|vE=v⊙+∂∂vEdR𝒩dER|vE=v⊙Δvcos⁡[2π(t-ϕ)τ].In order now to properly reproduce the experimental rate of nuclear recoil and in turn the expected number of events in a certain energy bin, one should take into account the response of the detectors as well. It can be done by the following energy convolution and transformation:(19)dRdE′=∑𝒩‍ϵ(E′)∫0∞‍dER𝒦𝒩(q𝒩ER,E′)dR𝒩dER(ER),
where E′ is the detected energy and the functions 𝒦𝒩(q𝒩ER,E′) and ϵ(E′) are the energy resolution centered in q𝒩ER and the detector’s efficiency, respectively. Here the sum runs over the different species of the detectors (e.g.,  DAMA is a multiple target experiment composed by crystal of sodium and iodine) and q𝒩 is the so-called quenching factor that takes account of the partial recollection of the released energy in the detector. After convolving with all the experimental effects, (19) must be averaged over the energy bins of the detector. For each energy bin k, we then obtain the number of the unmodulated events N0kth and of the modulated ones Nmkth as
(20)N0kth=wk∫ΔEk‍dE′dRdE′|vE=v⊙,(21)Nmkth=wkΔv∫ΔEk‍dE′∂∂vEdRdE′|vE=v⊙,
where wk is the exposure (expressed kilograms per day) and ΔEk is the width of the kth-esime energy bin. N0kth and Nmkth are the relevant quantities that one can use for the analysis of the experiments which address the annual modulation effect, namely,  DAMA and  CoGeNT. For the other experiments, only N0kth is instead relevant.Collecting all the elements in the previous equations we expand (20) and (21) and write
(22)N0kth=X∑i,j=116∑N,N′=p,n𝔠iN(λ,mχ)𝔠jN′(λ,mχ)00000000000×ℱ~i,j(N,N′)(mχ,k)|vE=v⊙,Nmkth=XΔv∑i,j=116∑N,N′=p,n𝔠iN(λ,mχ)𝔠jN′(λ,mχ)∂∂vE0000000000000×ℱ~i,j(N,N′)(mχ,k)|vE=v⊙.Hereℱ~i,j(N,N′)(mχ,k) is a sort ofintegrated form factor that encapsulates all the information related to nuclear physics, astrophysics, and the detector dependency of the rate. There is one of these factors for each energy bin k of each experiment under consideration and for each choice of the operators pair (i,j) and nucleons pair (N,N′). It reads explicitly
(23)ℱ~i,j(N,N′)(mχ,k)=wk∑𝒩‍ξ𝒩∫ΔEk‍dE′ϵ(E′)000000000000000×∫0∞‍dER𝒦T(q𝒩ER,E′)0000000000000000000×ℱi,j(N,N′)(ER,𝒩).It is worth noticing that the linearity of the expected number of events in theintegrated form factors is a fundamental ingredient. This, indeed, lets us parametrize the model dependent part of the rate from the model independent one, encapsulated in the function ℱ~i,j(N,N′)(mχ,k). Therefore, if the experimental collaborations, instead of presenting the experimental results in terms of a specific DM model, release all of this finite number ofintegrated form factors, one will be easily able to obtain the expected number of events for any kind of interactions, whose particle physics is completely encapsulated in the coefficient 𝔠iN. A first example of the potentiality of such parametrization has been presented in [32]. In particular the authors provide a self-contained set of numerical tools based on theintegrated form factors above, to derive the bounds from some current direct detection experiments on virtually any arbitrary model of DM elastically scattering on target nuclei.


## 2.1. Kinematics

As already stated, when a DM particle scatters off a nucleus, depending on the DM properties, one can envision at least two distinct kinematics, the elastic and the inelastic one. The elastic scattering is represented by(1)χ+𝒩(A𝒩,Z𝒩)atrest⟶χ+𝒩(A𝒩,Z𝒩)recoil,
while the inelastic is
(2)χ+𝒩(A𝒩,Z𝒩)atrest⟶χ′+𝒩(A𝒩,Z𝒩)recoil.
In (1) and (2), χ and χ′ are two DM particle states, and A𝒩, Z𝒩 are, respectively, the mass and atomic numbers of the nucleus 𝒩.We know, thanks to rotational curves data, that the velocity ofχ in the vicinity of the Earth is of the order of 10-3c. Hence the scattering between a DM particle with velocity v and mass mχ and a nucleus at rest with mass m𝒩 occurs in deeply nonrelativistic regime. This is analogous with the collision between two billiard balls and therefore the recoil energy can be simply obtained by considering energy and momenta conservation separately. In the detector rest frame, such scattering would end up into a nuclear recoil of energy:
(3)ER=2μχ𝒩2m𝒩v2(1-(vt2/2v2)-1-(vt2/v2)cos⁡θ2),vt=0,elasticscattering,vt≠0,inelasticscattering,
where μχ𝒩=mχm𝒩/(mχ+m𝒩) is the DM-nucleus reduced mass, θ as usual is the scattering angle in the center-of-mass frame, and vt=2δ/μχ𝒩 is a threshold velocity. Here δ=mχ′-mχ is the mass splitting between χ and χ′, and the equation above holds for δ≪mχ′, mχ. Elastic scattering occurs for δ=0, while δ≠0 implies inelastic scattering. The minimal DM velocity providing the ER recoil energy can be then obtained by putting θ=-π in (3). One gets
(4)vmin⁡(ER)=m𝒩ER2μχ𝒩2(1+μχ𝒩δm𝒩ER).


## 2.2. Formalism of Nonrelativistic Operators

Having a disposal the basic quantityvmin⁡ that fully describes the kinematics of the DM-nucleus scatterings in (1) and (2) and the differential rate of nuclear recoil expected in a given detector can be achieved by weighting the differential cross section dσ𝒩/dER with the DM velocity distribution fE(v→) in the velocity range allowed by the kinematics. It reads as follows:
(5)dR𝒩dER=ξ𝒩m𝒩ρ⊙mχ∫vmin⁡(ER)vesc‍d3vvfE(v→)dσ𝒩dER(v,ER),
where ξ𝒩 are the mass fractions of different nuclides (ξ𝒩=103NAm𝒩ζ𝒩/kgA¯, where NA=6.022×1023 is Avogadro’s number, ζ𝒩 are the numeric abundances and A¯≡∑𝒩‍ζ𝒩A𝒩.), and ρ⊙≃0.3GeV/cm3 is the DM energy density at the Earth’s location. This is the canonical value routinely adopted in the literature (see, e.g., [2–4]). Recent computations have found a higher central value of it, still subject for some debate [5–7]. The integral in the right-hand side of the equation above is a key ingredient because it contains all the information related to the geometry of the DM halo, the nature of the DM-nucleus interaction, and in turn the nuclear response of the target. In the following we present the formalism of nonrelativistic (NR) operators and then we describe how to write down the main observables in terms of it.(i)
As already stated, the weight of the velocity integral is the DM velocity distributionfE(v→). In the rest frame of our Galaxy it can be roughly approximated with a Maxwell-Boltzmann distribution due to the fact, as pointed out in [8], that the violent relaxation of the gravitational potential at the formation of the Milky Way leads to fast mixing of the DM phase space elements. Therefore, DM particles were basically frozen in high-entropy configuration which are indeed Maxwell-Boltzmann-like. This has been roughly confirmed by some recent numeral simulation, but as one can see in [9, Figure 2], there are some deviations due to the DM assembly history of the Milky Way. Indeed one can see different features both at low and high velocities that of course are not present in the case of a Maxwell-Boltzmann distribution. Furthermore the geometry of the halo predicted by this kind of numerical simulation is not exactly spherical but tends to a triaxial configuration (see, e.g., [10]).Keeping in mind that this is not the truth because we have not directly measured the properties of the DM halo yet, what we know for sure is that, since we are setting in the Earth laboratory frame, one has to boost the DM-nucleus relative velocity with the drift velocity of the Sun around the galactic center and the Earth’s orbital velocity projected in the galactic plane. Therefore the DM velocity distributionfE(v→) in the Earth frame is related to the DM velocity distribution in the galactic one fG(w→) by the Galilean velocity transformation fE(v→)=fG(v→+v→E(t)), where
(6)v→E(t)=v→⊙+v→⊕(t).Herev→⊙=v→G+v→S is the sum of the galactic rotational velocity of the local system v→G=(0,v0,0) (here v0≃220±50km/s is the local DM velocity) and the Sun’s proper motion v→S≃(10,13,7) km/s with absolute value v⊙=233±50km/s, while v→⊕(t) is the time dependent Earth’s orbital velocity with period of one year, phase around June 2nd (when it is aligned to v→⊙), and size v⊕≃30km/s, and it is inclined of an angle γ≃60∘ with respect to the galactic plane. More details can be found in [11–13]. Since in (6) the Earth’s orbital velocity projected in the galactic plane v⊕cos⁡γ is relatively small compared to v⊙, we can approximate v→E(t) with its component directed toward the galactic center. We can then write
(7)vE(t)≃v⊙+Δvcos⁡[2π(t-ϕ)τ],whereΔv=v⊕cos⁡⁡γ≃15km/s. Here ϕ≃152.5 days (June 2nd) is the phase and τ≃365 days is the period of the Earth motion around the Sun.Therefore, in general, it is expected an annual modulation in the counting rate of direct detection experiment, as the incoming flux of DM particles toward us contains small oscillatory terms due to the Earth’s proper motion around the Sun during the years. More details on this model independent signature can be found, for instance, in [14, 15]. It is worth stressing that looking for an annual modulation of the counting rate in direct detection is very challenging from the experimental point of view, simply because the modulus of the time dependent component of the velocity is relatively small compared to the constant one; in better words, the size of the modulated signal that experiments like  DAMA and  CoGeNT are looking for is suppressed with respect to the unmodulated one, due basically to the collective motion of the Solar system around the galactic center.(ii)
In the nonrelativistic limit the differential cross section can be written in the usual form:(8)dσ𝒩dER(v,ER)=132π1mχ2m𝒩1v2|ℳ𝒩|2¯,where|ℳ𝒩|2¯ is the DM-nucleus matrix element that encodes all the information related to the nuclear physics and the nature of the interactions. Since, as already stated, the velocity of the DM particles in the vicinity of the Earth is much smaller than the speed of light, the framework of relativistic quantum field theory is not too much appropriate to study signals in this field, especially if we do not have any clue about the underlying theory that describes the DM-quark interactions. A more useful and general framework is actually the one based on nonrelativistic operators. Indeed, since, for the elastic scattering, the relevant degrees of freedom are the DM-nucleon relative velocity v→, the exchanged momentum q→, the nucleon spin s→N (N=p, n can be proton or neutron), and the DM one s→χ (if different from zero), the scattering amplitude will then be a rotationally invariant function of this variables. In this regard, with v→, q→, s→N, and s→χ [16], we can construct a basis for 16 operators, which includes all possible spin configurations (the complete list and the numbering of the nonrelativistic operators can be found in [17–19]). The DM-nucleon matrix element can then be written as a linear combination of these operators, with coefficients that may depend on the momenta only through the q2 or v2 scalars (DM models feature a long range interaction with the Standard Model fermions provide perhaps the most notable example. Indeed, in this case, the differential cross section is enhanced at smaller exchanged momenta, due to the negative power dependence of q. A specific realization is offered by models in which DM particles can carry small but nonzero electric charge, electric dipole moment, or magnetic one. Their phenomenology in the context of DM direct searches has been studied in [20–31] and references therein.). In particular it reads
(9)ℳN=∑i=116‍𝔠iN(λ,mχ)𝒪iNR,where the coefficients𝔠iN are functions of the free parameters of the underlying relativistic theory, such as mediator masses, mixing angles, and couplings constants (collectively denoted by λ), and the DM mass mχ. For instance, if the interaction between a fermionic DMχ and the nucleon N is described by the “standard” SI Lagrangian gN/Λ2χ¯χN¯N, the only nonrelativistic operator involved will be the identity (𝒪1NR=𝕝). Before moving on, it is worth stressing that the operators above, which entirely describe the nonrelativistic physics of the DM-nucleus scattering, are many and they can either depend on the exchanged momenta q→ or the relative velocity v→. As a consequence, focussing on just one operator (e.g., 𝒪1NR) is not the most model independent way to study signals in direct DM searches, since theories can predict several operators entering together with the possibility that some of them may even interfere.Since now the nucleus is of course made of protons and neutrons, one has to correct the DM-nucleon matrix element (9) with the nuclear responses, which are a sort of form factors that take account of the finite size of the target nuclei. According to [17, Equation (55)] we can then write the spin-averaged amplitude squared for scattering off a target nucleus 𝒩 with mass m𝒩 as
(10)|ℳ𝒩|2¯=m𝒩2mN2∑i,j=116‍∑N,N′=p,n‍𝔠iN𝔠jN′Fi,j(N,N′).TheFi,j(N,N′)(v,ER,𝒩) are the nuclear responses and they depend critically on the type of scattering nucleus 𝒩: they are also functions of mχ, v, and the nuclear recoil energy ER=q2/2m𝒩. In Figure 1, the nuclear responses for both light element (on the left panel) and heavy one (on the right panel), considering two completely different kinds of interactions, are depicted. On a more specific level, in thick red the total nuclear response for the standard SI interaction is shown. As we can see, in the energy range of interest in direct DM experiments (few keV), it is enhanced by the canonical A𝒩2 factor: this is due to the fact that, in this case, the incoming DM particles see the nucleus as a point-like object with A𝒩 scatter centers. On the other hand, we show in blue the nuclear responses of a completely different type of interaction which is both momenta and spin dependent (𝒪6NR=(s→χ·q→)(s→N·q→)). This type of interaction can give rise if the DM-nucleus interaction is mediated by a pseudoscalar particle. We can see that the nuclear responses are dominated by the scattering with protons, simply because neither the fluorine nor the iodine has unpaired neutrons and moreover the interaction goes to zero in the small momenta transferred limit.A complete set of these nuclear responses for each pair of nucleons(N,N′), each pair of nonrelativistic operators (i,j), and for several target nuclei 𝒩 has been provided in the appendices of [17]. This is extremely useful because, for the first time, different types of interactions can be studied in a more general ground. However, since this kind of computations is quite new, the uncertainties, especially for the spin, momenta, and velocity dependent interactions, are still quite large.Nuclear responses of the fluorine (a) and the iodine (b) targets considering two completely different types of interactions. On a more specific level, we show, in red, the nuclear responses for the “standard” SI contact interaction, while in blue those for another kind of interaction described in the nonrelativistic limit by the operator 𝒪6NR=(s→χ·q→)(s→N·q→). The different hatching refers instead to the possible choices of nucleon pairs in the nucleus. In particular, the solid, dashed, and dotted lines are for (p,p), (n,n), and (p,n) pairs, while the thick ones are obtained by summing over all of them (total nuclear responses). Notice that the total nuclear responses for the “standard” SI contact interaction (thick red lines) reduce to A𝒩2FHelm2(q), where FHelm2(q) is the customary Helm form factor [33].
(a)
(b)Having a disposal the general relation of the DM-nucleus matrix element (10), the differential cross section (8), and in turn the rate of nuclear recoil (5) can be rewritten, following [32], in a very general way. It reads
(11)dR𝒩dER=Xξ𝒩∑i,j=116∑N,N′=p,n𝔠iN(λ,mχ)00000000000000000000×𝔠jN′(λ,mχ)ℱi,j(N,N′)(ER,𝒩),
where we defined
(12)X≡ρ⊙mχ132π1mχ2mN2,(13)ℱi,j(N,N′)(ER,𝒩)≡∫vmin⁡⁡(ER)vescd3v1vfG(v→+v→E(t))Fi,j(N,N′)(v,ER,𝒩).


### 2.2.1.  “Standard” SI Contact Interaction

To give a concrete example, in the following the coefficient and the nuclear responses for the “standard” SI contact interaction are showed explicitly. In particular, when experiments present results in terms of this interaction, they implicitly assume (probably inspired by supersymmetric neutralino scattering) the following DM-nucleon effective Lagrangian:(14)ℒSIN=λSINχ¯χN¯N.
From this the DM-nucleon matrix element can be obtained by contracting (14) with initial and final states of the scattering. Performing then the nonrelativistic limit one gets
(15)ℳSIN≡〈χN|ℒSIN|χN〉≃4λSINmχmN𝒪1NR,
where 𝒪1NR=𝕝 is the operator that describes the nonrelativistic limit of (14) and the coefficient is obviously given by 𝔠1N(λSIN,mχ)=4λSINmχmN.For the SI scattering, DM-proton and DM-neutron couplings are customarily assumed to be equal and thusλSIp=λSIn≡λSI (isospin violating process is not taken into account). Plugging back then the coefficient in the general equation (11) and defining σSI=λSI2/π·μN2 (μN is the DM-nucleon reduced mass) to make contact with the usual physical cross sections, the differential rate of nuclear recoil reduces to the very well-know compact form:
(16)dR𝒩dER=ξ𝒩m𝒩ρ⊙mχ[m𝒩2μN2σSIF𝒩tot(ER)ℐ(ER)],
where we define the total SI nuclear response as
(17)F𝒩tot(ER)=∑N,N′=p,n‍F1,1(N,N′)(ER,𝒩).
Here, since λSIp=λSIn, the total nuclear response reduces to the square of the customary Helm form factor [33] multiplied by the coherent enhancement of the cross section A𝒩2. In particular it reads F𝒩tot(ER)=A𝒩2FHelm2(ER). The function ℐ(ER)=∫vmin⁡vesc‍d3vfG(v→+v→E(t))/v is instead the velocity integral encountered many times in the literature. It is worth noticing that ℐ(ER) can be written in this way because in the “standard” SI interaction the operator that the describes the nonrelativistic limit of the effective Lagrangian does not carry any dependences on the relative velocity v→. In Section 3.2 the interpretation of the experimental results in terms of the “standard” SI interaction is briefly reviewed.


## 2.2.1.  “Standard” SI Contact Interaction

To give a concrete example, in the following the coefficient and the nuclear responses for the “standard” SI contact interaction are showed explicitly. In particular, when experiments present results in terms of this interaction, they implicitly assume (probably inspired by supersymmetric neutralino scattering) the following DM-nucleon effective Lagrangian:(14)ℒSIN=λSINχ¯χN¯N.
From this the DM-nucleon matrix element can be obtained by contracting (14) with initial and final states of the scattering. Performing then the nonrelativistic limit one gets
(15)ℳSIN≡〈χN|ℒSIN|χN〉≃4λSINmχmN𝒪1NR,
where 𝒪1NR=𝕝 is the operator that describes the nonrelativistic limit of (14) and the coefficient is obviously given by 𝔠1N(λSIN,mχ)=4λSINmχmN.For the SI scattering, DM-proton and DM-neutron couplings are customarily assumed to be equal and thusλSIp=λSIn≡λSI (isospin violating process is not taken into account). Plugging back then the coefficient in the general equation (11) and defining σSI=λSI2/π·μN2 (μN is the DM-nucleon reduced mass) to make contact with the usual physical cross sections, the differential rate of nuclear recoil reduces to the very well-know compact form:
(16)dR𝒩dER=ξ𝒩m𝒩ρ⊙mχ[m𝒩2μN2σSIF𝒩tot(ER)ℐ(ER)],
where we define the total SI nuclear response as
(17)F𝒩tot(ER)=∑N,N′=p,n‍F1,1(N,N′)(ER,𝒩).
Here, since λSIp=λSIn, the total nuclear response reduces to the square of the customary Helm form factor [33] multiplied by the coherent enhancement of the cross section A𝒩2. In particular it reads F𝒩tot(ER)=A𝒩2FHelm2(ER). The function ℐ(ER)=∫vmin⁡vesc‍d3vfG(v→+v→E(t))/v is instead the velocity integral encountered many times in the literature. It is worth noticing that ℐ(ER) can be written in this way because in the “standard” SI interaction the operator that the describes the nonrelativistic limit of the effective Lagrangian does not carry any dependences on the relative velocity v→. In Section 3.2 the interpretation of the experimental results in terms of the “standard” SI interaction is briefly reviewed.


## 2.3.  Experimental Observables

Since, as already stated, the Earth’s orbital velocity projected in the galactic plane is relatively small compared to the drift velocity of the Sun, we can then expand the recoil rate (11), assuming that the velocity distribution is not strongly anisotropic. Then, by means then of (7), one gets
(18)dR𝒩dER(t)≃dR𝒩dER|vE=v⊙+∂∂vEdR𝒩dER|vE=v⊙Δvcos⁡[2π(t-ϕ)τ].In order now to properly reproduce the experimental rate of nuclear recoil and in turn the expected number of events in a certain energy bin, one should take into account the response of the detectors as well. It can be done by the following energy convolution and transformation:(19)dRdE′=∑𝒩‍ϵ(E′)∫0∞‍dER𝒦𝒩(q𝒩ER,E′)dR𝒩dER(ER),
where E′ is the detected energy and the functions 𝒦𝒩(q𝒩ER,E′) and ϵ(E′) are the energy resolution centered in q𝒩ER and the detector’s efficiency, respectively. Here the sum runs over the different species of the detectors (e.g.,  DAMA is a multiple target experiment composed by crystal of sodium and iodine) and q𝒩 is the so-called quenching factor that takes account of the partial recollection of the released energy in the detector. After convolving with all the experimental effects, (19) must be averaged over the energy bins of the detector. For each energy bin k, we then obtain the number of the unmodulated events N0kth and of the modulated ones Nmkth as
(20)N0kth=wk∫ΔEk‍dE′dRdE′|vE=v⊙,(21)Nmkth=wkΔv∫ΔEk‍dE′∂∂vEdRdE′|vE=v⊙,
where wk is the exposure (expressed kilograms per day) and ΔEk is the width of the kth-esime energy bin. N0kth and Nmkth are the relevant quantities that one can use for the analysis of the experiments which address the annual modulation effect, namely,  DAMA and  CoGeNT. For the other experiments, only N0kth is instead relevant.Collecting all the elements in the previous equations we expand (20) and (21) and write
(22)N0kth=X∑i,j=116∑N,N′=p,n𝔠iN(λ,mχ)𝔠jN′(λ,mχ)00000000000×ℱ~i,j(N,N′)(mχ,k)|vE=v⊙,Nmkth=XΔv∑i,j=116∑N,N′=p,n𝔠iN(λ,mχ)𝔠jN′(λ,mχ)∂∂vE0000000000000×ℱ~i,j(N,N′)(mχ,k)|vE=v⊙.Hereℱ~i,j(N,N′)(mχ,k) is a sort ofintegrated form factor that encapsulates all the information related to nuclear physics, astrophysics, and the detector dependency of the rate. There is one of these factors for each energy bin k of each experiment under consideration and for each choice of the operators pair (i,j) and nucleons pair (N,N′). It reads explicitly
(23)ℱ~i,j(N,N′)(mχ,k)=wk∑𝒩‍ξ𝒩∫ΔEk‍dE′ϵ(E′)000000000000000×∫0∞‍dER𝒦T(q𝒩ER,E′)0000000000000000000×ℱi,j(N,N′)(ER,𝒩).It is worth noticing that the linearity of the expected number of events in theintegrated form factors is a fundamental ingredient. This, indeed, lets us parametrize the model dependent part of the rate from the model independent one, encapsulated in the function ℱ~i,j(N,N′)(mχ,k). Therefore, if the experimental collaborations, instead of presenting the experimental results in terms of a specific DM model, release all of this finite number ofintegrated form factors, one will be easily able to obtain the expected number of events for any kind of interactions, whose particle physics is completely encapsulated in the coefficient 𝔠iN. A first example of the potentiality of such parametrization has been presented in [32]. In particular the authors provide a self-contained set of numerical tools based on theintegrated form factors above, to derive the bounds from some current direct detection experiments on virtually any arbitrary model of DM elastically scattering on target nuclei.


## 3. Direct Detection Status

### 3.1. Experimental Landscape

The elusive nature of DM particles makes their detection a challenge for experimentalist: in fact, considering typical atomic masses for the target nuclei ofm𝒩=100GeV and typical properties for the DM halo (Maxwell-Boltzmann distribution with dispersion velocity v0=220km/s), DM particles should induce tiny nuclear recoil in the range 1÷103 keV with a total rate lower than 1 cpd/kg/keV (here cpd refers to counts per days). Due to this rare phenomenon, the experimental priorities in this field are as follows:(i)
the detectors must work deeply underground in order to avoid the high rate interaction induced by cosmic rays scattering off on target nuclei,(ii)
they must use active shields and very clean materials against the residual radioactivity in the tunnel (mostlyα-particles, neutrons, and photons),(iii)
they must distinguish multiple scattering, simply because DM particles do not interact twice in the detector, being of course weakly interacting particles.A variety of different experimental techniques with the aim of measuring the tiny energy released by a DM interaction have been developed. In the following the main channels by which the scattered nucleus can deposit energy in the detectors are summarized as follows.(i)
Scintillation detection: a particle interacting within a scintillating target induces the emission of light produced by the deexcitation of exited atoms. This light can be detected by appropriate photomultipliers. Typically NaI(Tl) and xenon are used as scintillators.(ii)
Ionization detection: a particle interacting inside a target produces an amount of free electron-ion pairs that can be detected with a collecting drift field and a device sensitive to the electric charge.(iii)
Phonon detection: a particle interacting inside a detector releases a tiny energy deposition with a subsequently increase of the temperature. Cryogenic apparatus working at very low temperature (around few mK) are able to measure this small variation making this detection technique possible.Since various types of interacting particles release a different amount of energy in the channels commented upon above, for a better rejection of the background events, most of the experiments are designed to be sensitive to more than one of them. Indeed, thanks to this ability, the ratio between two channels can be used to distinguish between nuclear (due to a DM interaction) and electromagnetic recoils.The only two experiments that do not use this kind of experimental technique areDAMA and  CoGeNT that, perhaps not coincidentally, were the only two detectors (until the  CRESST-II and recently the  CDMS-Si results) observing an excess in their counting rate. In particular, since both the experiments are not able to disentangle the nuclear recoil signals from the electromagnetic ones, they infer the presence of DM in our halo by exploiting the model independent annual modulation signature in the counting rate. As follows, I review the most important experiments in this field.


#### 3.1.1. DAMA

TheDAMA experiment, located at the National Laboratory of Gran Sasso, is an observatory for rare processes made of highly radiopure scintillators (NaI(Tl) crystals). In particular the former  DAMA/NaI and the current  DAMA/Libra experiments have the main aim of investigating the presence of DM particles in the galactic halo by exploiting the model independent annual modulation induced by a DM interaction.TheDAMA detector is only able to measure the fraction of energy deposited in scintillation light, while the phonon excitation due to multiple nuclear interactions is not detected. This effect, as commented upon in Section 2.3, is taken into account by the so-called quenching factor. In particular, for NaI(Tl) crystals, one normally considers qNa≃0.3 and qI≃0.09 [34]. It is, however, known [34–36] that some scattered nuclei can travel long distances along the crystallographic axes and planes without colliding with other nuclei. This process is called channeling effect, and since no phonon excitation is produced, the scattered nuclei deposit all the energy in the detector electromagnetically (q𝒩≃1). The fraction of channeling depends on the nucleus itself, on the recoil energy and strongly on the temperature. It has been calculated, for instance, in [35] and, in particular, it has been found that, for low energy recoiling Na and I ions, the fraction of channeled events can be relatively large. Nevertheless, this result has not been confirmed by other theoretical [37–39] and experimental studies [40], which instead suggest that the channeling effect in NaI crystals is negligible.TheDAMA collaboration published results of the combined  DAMA/NaI and  DAMA/Libra experiments [41, 42], corresponding to an exposure of 1.17 ton·yr for a target. They observed a cosine-like modulation, present only in thesingle-hit events, with a measured period τ=(0.999±0.002) yr and a phase t′=(146±7) days well compatible with the roughly 152.30 days expected for DM signals (see the previous section). The modulation is present only in the low energy window (2–6) keVee (here, keVee refers to electron equivalent recoil energy. This must be converted with the quenching factor in order to get the total energy released to the nucleus by the scattering process (ER[keV]≡ER[keVee]/q𝒩)) and its amplitude Nm(2–6)exp⁡ is equal to (0.0116±0.0013) cpd/(kg·keV) at 8.9σ CL [42].


#### 3.1.2.  CoGeNT

TheCoGeNT experiment employs p-type point-contact Germanium detectors operating in the Soudan Underground Laboratory.LikeDAMA, the experiment performs just one of the techniques commented upon above; namely, the p-type point-contact germanium detectors only measure the fraction of energy deposited by incident particles in the form of ionization. The lack of energy measured is again taken into account by the quenching factor that can be extracted through the following empirical relation Edet⁡=0.2ER1.12 [43]. By virtue of its low energy threshold (0.4 keVee) and the ability to reject surface backgrounds, this type of detector is particularly sensitive to light DM candidates, although large background contamination may be present for these low energies.In 2010 the collaboration has reported a step rise of nuclear recoil spectrum at low energy which is not directly explainable in terms of known radioactive background [44]. The energy region probed by  CoGeNT partially overlaps with the one in which the  DAMA apparatus reported the annual modulation signal, and therefore it is natural to interpret the excess at low energy due to a DM interaction. In view of that, fifteen months of cumulative  CoGeNT data (442 live days) have been also examined in order to look for an annual modulation signature [45]. In particular the  CoGeNT data seem to favor a cosine-like modulation slightly shifted in phase with respect to the one measured by  DAMA. The modulation is again only present at low energy with a statistical significance of 2.8σ CL, limited by the short exposure.Very recently, the collaboration presented a 3.4-year data taking, with an improved analysis that allows a better discrimination between pure bulk and pure surface events [46, 47]. In the same energy region, where in 2010 they observed a step rise of the nuclear recoil spectrum, they still reported a preference for an annual modulation in the pure bulk counting rate. Although the statistical significance of the modulation is modest (2.2σ CL), the phase is compatible with the one observed by  DAMA.


#### 3.1.3.  CRESST-II

TheCRESST-II cryogenic dark matter experiment, located at the national laboratory of Gran Sasso, employs 300 g of scintillating CaWO4 crystals.In particular the eight moduli of the experiments measure the deposited energy in the form of phonons and scintillation light. The former provides a precise measure of the energy deposited, while the ratio between them gives an excellent rejection of the background events.In 2011 the collaboration completed a 730 kg·day data taking [48]. In particular, 67 events were found in the DM acceptance region with an expected background contributions in the same band which is not sufficient to explain all the observed events. The resulting statistical significance, at which the background-only hypothesis is rejected, is roughly 4.5σ CL.


#### 3.1.4.  CDMS

TheCDMS experiment, located at the Soudan Underground Laboratory, is composed by germanium and silicon cryogenic detectors.In particular 19 germanium and 11 silicon detectors measure the deposited energy in the form of phonons and ionization through superconducting technology. Thanks to the ability of detecting these two signals, the rejection of the electromagnetic recoils can be obtained by the ratio among them.In 2009 theCDMS collaboration [49] has reported no significant evidence for DM interaction. In particular, only germanium detectors were used with an exposure of 612 kg·days. The collaboration observed two events in the signal region with energy 12.3 keVnr and 15.5 keVnr (here, keVnr refers to nuclear recoil energy. Unlike the electron equivalent recoil energy, it already represents the total energy released to the nucleus by the scattering process (ER[keV]≡ER[keVnr])), against an expected background of radiogenic neutrons, cosmogenic neutrons, and misidentified surface events equal to roughly 0.9 events.More recently, the collaboration has also presented results based on the silicon analysis [50, 51]. In [51], a blind analysis of 140.2 kg·days revealed three DM candidate events in the signal region with an expected background of roughly 0.62 events. A profile likelihood ratio statistical test of the three events yields a 0.19% probability when the known-background-only hypothesis is tested, against the alternative DM + background hypothesis.


#### 3.1.5.  Xenon Experiments

Direct detection experiments based on liquid/gaseous xenon have done stunning progresses in the last decade. Indeed, thanks to the advantage in scale ability and to the good reconstruction of the three-dimensional coordinates, these experiments are among the largest in terms of fiducial mass.Discrimination between nuclear and electromagnetic recoils is achieved by the ratio between the scintillation signal due to a particle interaction in the liquid xenon and the subsequent ionization signal in the gas phase of the detector. Furthermore, the large mass number of the xenon nuclei makes them an excellent target for the detection of DM particles with SI interaction. However, the experiments based on liquid xenon are also sensitive to then-DM spin dependent interaction, by virtue of the unpaired neutron of the Xe129 and Xe131 isotopes.In the following the results of the main experiments using the double-phase xenon technology are summarized.(i)
XENON experiment: the  XENON detector is a two-phase time projection chamber located at the national laboratory of Gran Sasso. The first stage of the experiment was installed underground during March 2006, and it searched for DM interactions until October 2007. In 2008 the collaboration reported a blind analysis with an exposure of 5.4×58.6 kg·days, which yields no significance evidence for DM interactions [52].More recently theXENON10 experiment has been superseded by  XENON100 with more than one order of magnitude improvement in sensitivity. In [53], the collaboration reported the results of the last run, a blind analysis with an exposure of 34×224.6 kg·days which again yielded no evidence for DM interactions. In particular, the two observed events in the predefined nuclear recoil energy range of 6.6–43.3 keVnr were consistent with the background expectation of roughly 1 event.(iii)
LUX experiment: the  LUX experiment is a dual-phase xenon time-projection chamber operated at the Sanford Underground Research Facility in South Dakota. Like the  XENON100 experiment, it measures both the scintillation light and the ionization in the gas phase to disentangle the nuclear recoils from the electromagnetic ones. In [54], a nonblind analysis was performed on data, collected with an exposure of 118.3×85.3 kg·days. After cuts, 160 events were found within the signal energy region. The collaboration found that all of these events are compatible with the expected electron recoil background.


#### 3.1.6.  COUPP

TheCOUPP experiment operating at SNOLAB in Ontario is a 4 kg CF3I bubble chamber. By virtue of its unpaired proton, the fluorine nucleus gives an excellent sensitivity for p-DM spin-dependent interactions, while the iodine enhances the sensitivity for the SI ones.Particles entering the liquid in the superheated phase create an ionization tracks, around which the liquid vaporizes, forming microscopic bubbles. Bubbles grow in size as the chamber expands, making the detector able to record them both photographically and by pressure and acoustic measurements. The discrimination between nuclear and electromagnetic recoils can be achieved by choosing an appropriate chamber pressure and temperature. Indeed, under this condition, the abundant gamma-ray and beta-decay backgrounds do not nucleate bubbles.In [55], data obtained for an effective exposure to single recoil events of 437.4 kg·days (taking into account the 79.1% detection efficiency) were presented. Twenty single nuclear recoil events passing all the analysis cuts were observed over the three energy bins. Due to uncertainties in the neutron background estimation, the collaboration has not attempted any background subtraction and instead has treated all of them as DM events.


### 3.2. Interpretation in Terms of the “Standard” SI Interaction

There are many studies in literature that try to analyze the implications of the results ofDAMA and the anomalies in  CoGeNT and  CRESST-II and recently in  CDMS-Si together with the null results from other experiments in terms of specific DM models and interactions. A routine way to compare results from different DM experiments is by assuming a particular DM velocity distribution and a certain type of DM-nucleus interaction. The customary choices that can often be found in literature are (i) a truncated Maxwell-Boltzmann velocity distribution with isotropic velocity dispersion v0=220km/s and (ii) DM particles coupling through a contact SI interaction with equal strength to the protons and neutrons. In this case one customarily chooses the total DM-nucleon cross section σSI defined in Section 2.2.1 together with the DM mass mχ as free parameters, since the bounds and the allowed regions from different experiments can be compared on the same plot.In Figure2 a compilation of the allowed regions of the positive results experiments and the constraints coming from null results taken from different experimental collaborations are shown. Without entering in the details of the different analysis, we summarize the interpretation of the datasets of the positive results experiments listed above in terms of the “standard” SI contact interaction as follows.(i)
TheDAMA modulation signal effect has been shown to be compatible with two areas in the (mχ,σSI) parameter space, due to the different kinematics experienced by the two targets in the scintillator. In particular, a region pointing towards DM masses in the 10 GeV ballpark and towards a cross section σSI≃2×10-40 cm2 can be associated with a DM particle scattering off with sodium (see, e.g., the red contour in the top left-handed plot of Figure 2 denoted by  DAMA/Na). Another region around DM masses of roughly 60 GeV is instead due to the scattering with the iodine. In this case the cross section is enhanced by a larger coherent factor and indeed the favored value is σSI≃10-41 cm2 (red contour denoted by  DAMA/I).(ii)
The modulation inCoGeNT is mostly concentrated in the first two energy bins, close to the lower threshold of the detector. If this anomaly is interpreted in terms of SI contact interaction, it will be fitted by light DM candidates (around 9 GeV) with a total SI cross section σSI≃3×10-41 cm2 (see, e.g., the red contour in the top right-handed plot of Figure 2).(iii)
LikeDAMA, the  CRESST-II experiment is a multiple target detector and therefore more than one allowed region is in general expected. In particular, as one can see in Figure 2, the region pointing towards light DM masses (around 12 GeV) and σSI≃5×10-41 cm2 is due to the scattering with oxygen and calcium. On the other hand, the one pointing towards heavier DM mass (around 30 GeV) and cross section σSI≃10-42 cm2 can be associated with a DM particle scattering off with tungstenum. These regions are, for instance, denoted as gray contours in the top left-handed plot of Figure 2.(iv)
Very recently, also the 11 silicon detectors of theCDMS experiment reported an excess in their counting rate. Again, if the anomaly is interpreted in terms of SI contact interaction, it will be fitted by light DM particles (around 9 GeV) with a SI total cross section σSI≃2×10-41 cm2 (blue regions in the bottom left-handed plot of Figure 2).Figure 2
A compilation of the allowed regions and the constraints for the “standard” SI contact interaction in the (mχ,σSI) plane. The plots on the top row are taken from [53] (left) and from [46] (right), while those in the bottom row are taken from [51] (left) and from [54] (right). The galactic halo has been assumed in the form of an isothermal sphere with velocity dispersion v0=220km/s and local DM energy density ρ⊙=0.3GeV/cm3.The interesting feature is that all the experiments employing light target nuclei seem to be compatible with a DM interpretation, which, for the “standard” SI contact interaction, pin-points the properties of the DM particle quite precisely; namely, it leads to DM masses in the 10 GeV ballpark and total cross section in the range (2×10-41cm2–2×10-40 cm2).On the other hand, the constraints coming from null results are very stringent, and in particular theXENON100 and recently the  LUX results severely exclude the allowed parameter space of the positive results experiments. On a more specific level, before  LUX, the  XENON100 experiment provided the most stringent constraints for a DM mass above roughly 8 GeV, with a minimum of σSI≃1.8×10-44 cm2 (see, e.g., the thick blue line in the top left-handed plot of Figure 2). Nowadays this bound has been beat out by the one of  LUX by a factor ~2.5 (see the blue line in the bottom right-handed plot of Figure 2). Other relevant analyses can be found in [56–68].


### 3.3.  Uncertainties

Nevertheless, since the actual range of masses and cross section critically depends on several assumptions, one always bears in mind all the uncertainties that enter in this field when interpreting the data in terms of a specific DM interaction.(i)
A first class of uncertainties comes from the poor knowledge of the properties of the DM halo. The often used Maxwell-Boltzmann halo with isotropic velocity dispersion is only a benchmark choice, not the physical description of it. Other possible choices (spherical halo with nonisotropic velocity dispersion, axisymmetric halo, triaxial halo, and the possibility to have corotating halo with different number densities) must be taken into account, in order to have a much better idea of how the velocity distribution affects the final results. In particular, it has been recently demonstrated that the uncertainties on the astrophysical part of the rate can be removed by comparing different experimental results in thevmin⁡ space (see, e.g., [69–77]). In this way one can, for instance, relax a bit the tension between the positive results and the constraints showed in Figure 2 up. However, this is not still sufficient to reconcile the complicated puzzle the experimental data have left to us (see, e.g., [77–80]).(ii)
A second class of uncertainties comes from the experimental side. For instance, the direct measurements of quenching factor are performed with reference detectors, and thus the systematic uncertainties, important for all the detectors in direct DM searches, could be larger than what we normally assumed. These must be carefully investigated, because they produce a significant shift of the allowed regions in the total cross section/DM mass plane.(iii)
Together with these main uncertainties, the third and perhaps the most important class of them comes from the nature of the DM interaction and in turn from the nuclear responses of the target nuclei. Most of the models tested in recent years are based on the assumption that the mechanism of interaction is realized through a contact interaction. Deviations from this standard phenomenological approach are interesting to study, since, from more general conditions, the DM parameter space can be analyzed. In particular, there is a concrete possibility to remove the uncertainties coming from the nature of interaction thanks to the formalism of nonrelativistic operator reviewed in Section2.2. On a more specific level, if the experimental collaborations, in addition to presenting the results in terms of a specific DM-nucleus interaction, also provide all the model independentintegrated form factors, it will be extremely useful for the community. Indeed, in this way one can easily compute the expected number of events for any kind of interactions (e.g., including isospin violating interactions, momentum-dependent form factors, and velocity dependent form factors) and compare it directly with the experimental results.Therefore the data must be treated very prudently with a maximally conservative attitude, simply because slightly modification of the assumptions summarized above can affect the theoretical interpretation of the experimental results in a relevant way. In the next section, I will focus on the uncertainties coming from the nature of interaction and in particular I am going to show how the allowed regions and constraints are modified, if an exchanged momentum dependent DM-nucleus interaction is taken into account.


## 3.1. Experimental Landscape

The elusive nature of DM particles makes their detection a challenge for experimentalist: in fact, considering typical atomic masses for the target nuclei ofm𝒩=100GeV and typical properties for the DM halo (Maxwell-Boltzmann distribution with dispersion velocity v0=220km/s), DM particles should induce tiny nuclear recoil in the range 1÷103 keV with a total rate lower than 1 cpd/kg/keV (here cpd refers to counts per days). Due to this rare phenomenon, the experimental priorities in this field are as follows:(i)
the detectors must work deeply underground in order to avoid the high rate interaction induced by cosmic rays scattering off on target nuclei,(ii)
they must use active shields and very clean materials against the residual radioactivity in the tunnel (mostlyα-particles, neutrons, and photons),(iii)
they must distinguish multiple scattering, simply because DM particles do not interact twice in the detector, being of course weakly interacting particles.A variety of different experimental techniques with the aim of measuring the tiny energy released by a DM interaction have been developed. In the following the main channels by which the scattered nucleus can deposit energy in the detectors are summarized as follows.(i)
Scintillation detection: a particle interacting within a scintillating target induces the emission of light produced by the deexcitation of exited atoms. This light can be detected by appropriate photomultipliers. Typically NaI(Tl) and xenon are used as scintillators.(ii)
Ionization detection: a particle interacting inside a target produces an amount of free electron-ion pairs that can be detected with a collecting drift field and a device sensitive to the electric charge.(iii)
Phonon detection: a particle interacting inside a detector releases a tiny energy deposition with a subsequently increase of the temperature. Cryogenic apparatus working at very low temperature (around few mK) are able to measure this small variation making this detection technique possible.Since various types of interacting particles release a different amount of energy in the channels commented upon above, for a better rejection of the background events, most of the experiments are designed to be sensitive to more than one of them. Indeed, thanks to this ability, the ratio between two channels can be used to distinguish between nuclear (due to a DM interaction) and electromagnetic recoils.The only two experiments that do not use this kind of experimental technique areDAMA and  CoGeNT that, perhaps not coincidentally, were the only two detectors (until the  CRESST-II and recently the  CDMS-Si results) observing an excess in their counting rate. In particular, since both the experiments are not able to disentangle the nuclear recoil signals from the electromagnetic ones, they infer the presence of DM in our halo by exploiting the model independent annual modulation signature in the counting rate. As follows, I review the most important experiments in this field.


### 3.1.1. DAMA

TheDAMA experiment, located at the National Laboratory of Gran Sasso, is an observatory for rare processes made of highly radiopure scintillators (NaI(Tl) crystals). In particular the former  DAMA/NaI and the current  DAMA/Libra experiments have the main aim of investigating the presence of DM particles in the galactic halo by exploiting the model independent annual modulation induced by a DM interaction.TheDAMA detector is only able to measure the fraction of energy deposited in scintillation light, while the phonon excitation due to multiple nuclear interactions is not detected. This effect, as commented upon in Section 2.3, is taken into account by the so-called quenching factor. In particular, for NaI(Tl) crystals, one normally considers qNa≃0.3 and qI≃0.09 [34]. It is, however, known [34–36] that some scattered nuclei can travel long distances along the crystallographic axes and planes without colliding with other nuclei. This process is called channeling effect, and since no phonon excitation is produced, the scattered nuclei deposit all the energy in the detector electromagnetically (q𝒩≃1). The fraction of channeling depends on the nucleus itself, on the recoil energy and strongly on the temperature. It has been calculated, for instance, in [35] and, in particular, it has been found that, for low energy recoiling Na and I ions, the fraction of channeled events can be relatively large. Nevertheless, this result has not been confirmed by other theoretical [37–39] and experimental studies [40], which instead suggest that the channeling effect in NaI crystals is negligible.TheDAMA collaboration published results of the combined  DAMA/NaI and  DAMA/Libra experiments [41, 42], corresponding to an exposure of 1.17 ton·yr for a target. They observed a cosine-like modulation, present only in thesingle-hit events, with a measured period τ=(0.999±0.002) yr and a phase t′=(146±7) days well compatible with the roughly 152.30 days expected for DM signals (see the previous section). The modulation is present only in the low energy window (2–6) keVee (here, keVee refers to electron equivalent recoil energy. This must be converted with the quenching factor in order to get the total energy released to the nucleus by the scattering process (ER[keV]≡ER[keVee]/q𝒩)) and its amplitude Nm(2–6)exp⁡ is equal to (0.0116±0.0013) cpd/(kg·keV) at 8.9σ CL [42].


### 3.1.2.  CoGeNT

TheCoGeNT experiment employs p-type point-contact Germanium detectors operating in the Soudan Underground Laboratory.LikeDAMA, the experiment performs just one of the techniques commented upon above; namely, the p-type point-contact germanium detectors only measure the fraction of energy deposited by incident particles in the form of ionization. The lack of energy measured is again taken into account by the quenching factor that can be extracted through the following empirical relation Edet⁡=0.2ER1.12 [43]. By virtue of its low energy threshold (0.4 keVee) and the ability to reject surface backgrounds, this type of detector is particularly sensitive to light DM candidates, although large background contamination may be present for these low energies.In 2010 the collaboration has reported a step rise of nuclear recoil spectrum at low energy which is not directly explainable in terms of known radioactive background [44]. The energy region probed by  CoGeNT partially overlaps with the one in which the  DAMA apparatus reported the annual modulation signal, and therefore it is natural to interpret the excess at low energy due to a DM interaction. In view of that, fifteen months of cumulative  CoGeNT data (442 live days) have been also examined in order to look for an annual modulation signature [45]. In particular the  CoGeNT data seem to favor a cosine-like modulation slightly shifted in phase with respect to the one measured by  DAMA. The modulation is again only present at low energy with a statistical significance of 2.8σ CL, limited by the short exposure.Very recently, the collaboration presented a 3.4-year data taking, with an improved analysis that allows a better discrimination between pure bulk and pure surface events [46, 47]. In the same energy region, where in 2010 they observed a step rise of the nuclear recoil spectrum, they still reported a preference for an annual modulation in the pure bulk counting rate. Although the statistical significance of the modulation is modest (2.2σ CL), the phase is compatible with the one observed by  DAMA.


### 3.1.3.  CRESST-II

TheCRESST-II cryogenic dark matter experiment, located at the national laboratory of Gran Sasso, employs 300 g of scintillating CaWO4 crystals.In particular the eight moduli of the experiments measure the deposited energy in the form of phonons and scintillation light. The former provides a precise measure of the energy deposited, while the ratio between them gives an excellent rejection of the background events.In 2011 the collaboration completed a 730 kg·day data taking [48]. In particular, 67 events were found in the DM acceptance region with an expected background contributions in the same band which is not sufficient to explain all the observed events. The resulting statistical significance, at which the background-only hypothesis is rejected, is roughly 4.5σ CL.


### 3.1.4.  CDMS

TheCDMS experiment, located at the Soudan Underground Laboratory, is composed by germanium and silicon cryogenic detectors.In particular 19 germanium and 11 silicon detectors measure the deposited energy in the form of phonons and ionization through superconducting technology. Thanks to the ability of detecting these two signals, the rejection of the electromagnetic recoils can be obtained by the ratio among them.In 2009 theCDMS collaboration [49] has reported no significant evidence for DM interaction. In particular, only germanium detectors were used with an exposure of 612 kg·days. The collaboration observed two events in the signal region with energy 12.3 keVnr and 15.5 keVnr (here, keVnr refers to nuclear recoil energy. Unlike the electron equivalent recoil energy, it already represents the total energy released to the nucleus by the scattering process (ER[keV]≡ER[keVnr])), against an expected background of radiogenic neutrons, cosmogenic neutrons, and misidentified surface events equal to roughly 0.9 events.More recently, the collaboration has also presented results based on the silicon analysis [50, 51]. In [51], a blind analysis of 140.2 kg·days revealed three DM candidate events in the signal region with an expected background of roughly 0.62 events. A profile likelihood ratio statistical test of the three events yields a 0.19% probability when the known-background-only hypothesis is tested, against the alternative DM + background hypothesis.


### 3.1.5.  Xenon Experiments

Direct detection experiments based on liquid/gaseous xenon have done stunning progresses in the last decade. Indeed, thanks to the advantage in scale ability and to the good reconstruction of the three-dimensional coordinates, these experiments are among the largest in terms of fiducial mass.Discrimination between nuclear and electromagnetic recoils is achieved by the ratio between the scintillation signal due to a particle interaction in the liquid xenon and the subsequent ionization signal in the gas phase of the detector. Furthermore, the large mass number of the xenon nuclei makes them an excellent target for the detection of DM particles with SI interaction. However, the experiments based on liquid xenon are also sensitive to then-DM spin dependent interaction, by virtue of the unpaired neutron of the Xe129 and Xe131 isotopes.In the following the results of the main experiments using the double-phase xenon technology are summarized.(i)
XENON experiment: the  XENON detector is a two-phase time projection chamber located at the national laboratory of Gran Sasso. The first stage of the experiment was installed underground during March 2006, and it searched for DM interactions until October 2007. In 2008 the collaboration reported a blind analysis with an exposure of 5.4×58.6 kg·days, which yields no significance evidence for DM interactions [52].More recently theXENON10 experiment has been superseded by  XENON100 with more than one order of magnitude improvement in sensitivity. In [53], the collaboration reported the results of the last run, a blind analysis with an exposure of 34×224.6 kg·days which again yielded no evidence for DM interactions. In particular, the two observed events in the predefined nuclear recoil energy range of 6.6–43.3 keVnr were consistent with the background expectation of roughly 1 event.(iii)
LUX experiment: the  LUX experiment is a dual-phase xenon time-projection chamber operated at the Sanford Underground Research Facility in South Dakota. Like the  XENON100 experiment, it measures both the scintillation light and the ionization in the gas phase to disentangle the nuclear recoils from the electromagnetic ones. In [54], a nonblind analysis was performed on data, collected with an exposure of 118.3×85.3 kg·days. After cuts, 160 events were found within the signal energy region. The collaboration found that all of these events are compatible with the expected electron recoil background.


### 3.1.6.  COUPP

TheCOUPP experiment operating at SNOLAB in Ontario is a 4 kg CF3I bubble chamber. By virtue of its unpaired proton, the fluorine nucleus gives an excellent sensitivity for p-DM spin-dependent interactions, while the iodine enhances the sensitivity for the SI ones.Particles entering the liquid in the superheated phase create an ionization tracks, around which the liquid vaporizes, forming microscopic bubbles. Bubbles grow in size as the chamber expands, making the detector able to record them both photographically and by pressure and acoustic measurements. The discrimination between nuclear and electromagnetic recoils can be achieved by choosing an appropriate chamber pressure and temperature. Indeed, under this condition, the abundant gamma-ray and beta-decay backgrounds do not nucleate bubbles.In [55], data obtained for an effective exposure to single recoil events of 437.4 kg·days (taking into account the 79.1% detection efficiency) were presented. Twenty single nuclear recoil events passing all the analysis cuts were observed over the three energy bins. Due to uncertainties in the neutron background estimation, the collaboration has not attempted any background subtraction and instead has treated all of them as DM events.


## 3.1.1. DAMA

TheDAMA experiment, located at the National Laboratory of Gran Sasso, is an observatory for rare processes made of highly radiopure scintillators (NaI(Tl) crystals). In particular the former  DAMA/NaI and the current  DAMA/Libra experiments have the main aim of investigating the presence of DM particles in the galactic halo by exploiting the model independent annual modulation induced by a DM interaction.TheDAMA detector is only able to measure the fraction of energy deposited in scintillation light, while the phonon excitation due to multiple nuclear interactions is not detected. This effect, as commented upon in Section 2.3, is taken into account by the so-called quenching factor. In particular, for NaI(Tl) crystals, one normally considers qNa≃0.3 and qI≃0.09 [34]. It is, however, known [34–36] that some scattered nuclei can travel long distances along the crystallographic axes and planes without colliding with other nuclei. This process is called channeling effect, and since no phonon excitation is produced, the scattered nuclei deposit all the energy in the detector electromagnetically (q𝒩≃1). The fraction of channeling depends on the nucleus itself, on the recoil energy and strongly on the temperature. It has been calculated, for instance, in [35] and, in particular, it has been found that, for low energy recoiling Na and I ions, the fraction of channeled events can be relatively large. Nevertheless, this result has not been confirmed by other theoretical [37–39] and experimental studies [40], which instead suggest that the channeling effect in NaI crystals is negligible.TheDAMA collaboration published results of the combined  DAMA/NaI and  DAMA/Libra experiments [41, 42], corresponding to an exposure of 1.17 ton·yr for a target. They observed a cosine-like modulation, present only in thesingle-hit events, with a measured period τ=(0.999±0.002) yr and a phase t′=(146±7) days well compatible with the roughly 152.30 days expected for DM signals (see the previous section). The modulation is present only in the low energy window (2–6) keVee (here, keVee refers to electron equivalent recoil energy. This must be converted with the quenching factor in order to get the total energy released to the nucleus by the scattering process (ER[keV]≡ER[keVee]/q𝒩)) and its amplitude Nm(2–6)exp⁡ is equal to (0.0116±0.0013) cpd/(kg·keV) at 8.9σ CL [42].


## 3.1.2.  CoGeNT

TheCoGeNT experiment employs p-type point-contact Germanium detectors operating in the Soudan Underground Laboratory.LikeDAMA, the experiment performs just one of the techniques commented upon above; namely, the p-type point-contact germanium detectors only measure the fraction of energy deposited by incident particles in the form of ionization. The lack of energy measured is again taken into account by the quenching factor that can be extracted through the following empirical relation Edet⁡=0.2ER1.12 [43]. By virtue of its low energy threshold (0.4 keVee) and the ability to reject surface backgrounds, this type of detector is particularly sensitive to light DM candidates, although large background contamination may be present for these low energies.In 2010 the collaboration has reported a step rise of nuclear recoil spectrum at low energy which is not directly explainable in terms of known radioactive background [44]. The energy region probed by  CoGeNT partially overlaps with the one in which the  DAMA apparatus reported the annual modulation signal, and therefore it is natural to interpret the excess at low energy due to a DM interaction. In view of that, fifteen months of cumulative  CoGeNT data (442 live days) have been also examined in order to look for an annual modulation signature [45]. In particular the  CoGeNT data seem to favor a cosine-like modulation slightly shifted in phase with respect to the one measured by  DAMA. The modulation is again only present at low energy with a statistical significance of 2.8σ CL, limited by the short exposure.Very recently, the collaboration presented a 3.4-year data taking, with an improved analysis that allows a better discrimination between pure bulk and pure surface events [46, 47]. In the same energy region, where in 2010 they observed a step rise of the nuclear recoil spectrum, they still reported a preference for an annual modulation in the pure bulk counting rate. Although the statistical significance of the modulation is modest (2.2σ CL), the phase is compatible with the one observed by  DAMA.


## 3.1.3.  CRESST-II

TheCRESST-II cryogenic dark matter experiment, located at the national laboratory of Gran Sasso, employs 300 g of scintillating CaWO4 crystals.In particular the eight moduli of the experiments measure the deposited energy in the form of phonons and scintillation light. The former provides a precise measure of the energy deposited, while the ratio between them gives an excellent rejection of the background events.In 2011 the collaboration completed a 730 kg·day data taking [48]. In particular, 67 events were found in the DM acceptance region with an expected background contributions in the same band which is not sufficient to explain all the observed events. The resulting statistical significance, at which the background-only hypothesis is rejected, is roughly 4.5σ CL.


## 3.1.4.  CDMS

TheCDMS experiment, located at the Soudan Underground Laboratory, is composed by germanium and silicon cryogenic detectors.In particular 19 germanium and 11 silicon detectors measure the deposited energy in the form of phonons and ionization through superconducting technology. Thanks to the ability of detecting these two signals, the rejection of the electromagnetic recoils can be obtained by the ratio among them.In 2009 theCDMS collaboration [49] has reported no significant evidence for DM interaction. In particular, only germanium detectors were used with an exposure of 612 kg·days. The collaboration observed two events in the signal region with energy 12.3 keVnr and 15.5 keVnr (here, keVnr refers to nuclear recoil energy. Unlike the electron equivalent recoil energy, it already represents the total energy released to the nucleus by the scattering process (ER[keV]≡ER[keVnr])), against an expected background of radiogenic neutrons, cosmogenic neutrons, and misidentified surface events equal to roughly 0.9 events.More recently, the collaboration has also presented results based on the silicon analysis [50, 51]. In [51], a blind analysis of 140.2 kg·days revealed three DM candidate events in the signal region with an expected background of roughly 0.62 events. A profile likelihood ratio statistical test of the three events yields a 0.19% probability when the known-background-only hypothesis is tested, against the alternative DM + background hypothesis.


## 3.1.5.  Xenon Experiments

Direct detection experiments based on liquid/gaseous xenon have done stunning progresses in the last decade. Indeed, thanks to the advantage in scale ability and to the good reconstruction of the three-dimensional coordinates, these experiments are among the largest in terms of fiducial mass.Discrimination between nuclear and electromagnetic recoils is achieved by the ratio between the scintillation signal due to a particle interaction in the liquid xenon and the subsequent ionization signal in the gas phase of the detector. Furthermore, the large mass number of the xenon nuclei makes them an excellent target for the detection of DM particles with SI interaction. However, the experiments based on liquid xenon are also sensitive to then-DM spin dependent interaction, by virtue of the unpaired neutron of the Xe129 and Xe131 isotopes.In the following the results of the main experiments using the double-phase xenon technology are summarized.(i)
XENON experiment: the  XENON detector is a two-phase time projection chamber located at the national laboratory of Gran Sasso. The first stage of the experiment was installed underground during March 2006, and it searched for DM interactions until October 2007. In 2008 the collaboration reported a blind analysis with an exposure of 5.4×58.6 kg·days, which yields no significance evidence for DM interactions [52].More recently theXENON10 experiment has been superseded by  XENON100 with more than one order of magnitude improvement in sensitivity. In [53], the collaboration reported the results of the last run, a blind analysis with an exposure of 34×224.6 kg·days which again yielded no evidence for DM interactions. In particular, the two observed events in the predefined nuclear recoil energy range of 6.6–43.3 keVnr were consistent with the background expectation of roughly 1 event.(iii)
LUX experiment: the  LUX experiment is a dual-phase xenon time-projection chamber operated at the Sanford Underground Research Facility in South Dakota. Like the  XENON100 experiment, it measures both the scintillation light and the ionization in the gas phase to disentangle the nuclear recoils from the electromagnetic ones. In [54], a nonblind analysis was performed on data, collected with an exposure of 118.3×85.3 kg·days. After cuts, 160 events were found within the signal energy region. The collaboration found that all of these events are compatible with the expected electron recoil background.


## 3.1.6.  COUPP

TheCOUPP experiment operating at SNOLAB in Ontario is a 4 kg CF3I bubble chamber. By virtue of its unpaired proton, the fluorine nucleus gives an excellent sensitivity for p-DM spin-dependent interactions, while the iodine enhances the sensitivity for the SI ones.Particles entering the liquid in the superheated phase create an ionization tracks, around which the liquid vaporizes, forming microscopic bubbles. Bubbles grow in size as the chamber expands, making the detector able to record them both photographically and by pressure and acoustic measurements. The discrimination between nuclear and electromagnetic recoils can be achieved by choosing an appropriate chamber pressure and temperature. Indeed, under this condition, the abundant gamma-ray and beta-decay backgrounds do not nucleate bubbles.In [55], data obtained for an effective exposure to single recoil events of 437.4 kg·days (taking into account the 79.1% detection efficiency) were presented. Twenty single nuclear recoil events passing all the analysis cuts were observed over the three energy bins. Due to uncertainties in the neutron background estimation, the collaboration has not attempted any background subtraction and instead has treated all of them as DM events.


## 3.2. Interpretation in Terms of the “Standard” SI Interaction

There are many studies in literature that try to analyze the implications of the results ofDAMA and the anomalies in  CoGeNT and  CRESST-II and recently in  CDMS-Si together with the null results from other experiments in terms of specific DM models and interactions. A routine way to compare results from different DM experiments is by assuming a particular DM velocity distribution and a certain type of DM-nucleus interaction. The customary choices that can often be found in literature are (i) a truncated Maxwell-Boltzmann velocity distribution with isotropic velocity dispersion v0=220km/s and (ii) DM particles coupling through a contact SI interaction with equal strength to the protons and neutrons. In this case one customarily chooses the total DM-nucleon cross section σSI defined in Section 2.2.1 together with the DM mass mχ as free parameters, since the bounds and the allowed regions from different experiments can be compared on the same plot.In Figure2 a compilation of the allowed regions of the positive results experiments and the constraints coming from null results taken from different experimental collaborations are shown. Without entering in the details of the different analysis, we summarize the interpretation of the datasets of the positive results experiments listed above in terms of the “standard” SI contact interaction as follows.(i)
TheDAMA modulation signal effect has been shown to be compatible with two areas in the (mχ,σSI) parameter space, due to the different kinematics experienced by the two targets in the scintillator. In particular, a region pointing towards DM masses in the 10 GeV ballpark and towards a cross section σSI≃2×10-40 cm2 can be associated with a DM particle scattering off with sodium (see, e.g., the red contour in the top left-handed plot of Figure 2 denoted by  DAMA/Na). Another region around DM masses of roughly 60 GeV is instead due to the scattering with the iodine. In this case the cross section is enhanced by a larger coherent factor and indeed the favored value is σSI≃10-41 cm2 (red contour denoted by  DAMA/I).(ii)
The modulation inCoGeNT is mostly concentrated in the first two energy bins, close to the lower threshold of the detector. If this anomaly is interpreted in terms of SI contact interaction, it will be fitted by light DM candidates (around 9 GeV) with a total SI cross section σSI≃3×10-41 cm2 (see, e.g., the red contour in the top right-handed plot of Figure 2).(iii)
LikeDAMA, the  CRESST-II experiment is a multiple target detector and therefore more than one allowed region is in general expected. In particular, as one can see in Figure 2, the region pointing towards light DM masses (around 12 GeV) and σSI≃5×10-41 cm2 is due to the scattering with oxygen and calcium. On the other hand, the one pointing towards heavier DM mass (around 30 GeV) and cross section σSI≃10-42 cm2 can be associated with a DM particle scattering off with tungstenum. These regions are, for instance, denoted as gray contours in the top left-handed plot of Figure 2.(iv)
Very recently, also the 11 silicon detectors of theCDMS experiment reported an excess in their counting rate. Again, if the anomaly is interpreted in terms of SI contact interaction, it will be fitted by light DM particles (around 9 GeV) with a SI total cross section σSI≃2×10-41 cm2 (blue regions in the bottom left-handed plot of Figure 2).Figure 2
A compilation of the allowed regions and the constraints for the “standard” SI contact interaction in the (mχ,σSI) plane. The plots on the top row are taken from [53] (left) and from [46] (right), while those in the bottom row are taken from [51] (left) and from [54] (right). The galactic halo has been assumed in the form of an isothermal sphere with velocity dispersion v0=220km/s and local DM energy density ρ⊙=0.3GeV/cm3.The interesting feature is that all the experiments employing light target nuclei seem to be compatible with a DM interpretation, which, for the “standard” SI contact interaction, pin-points the properties of the DM particle quite precisely; namely, it leads to DM masses in the 10 GeV ballpark and total cross section in the range (2×10-41cm2–2×10-40 cm2).On the other hand, the constraints coming from null results are very stringent, and in particular theXENON100 and recently the  LUX results severely exclude the allowed parameter space of the positive results experiments. On a more specific level, before  LUX, the  XENON100 experiment provided the most stringent constraints for a DM mass above roughly 8 GeV, with a minimum of σSI≃1.8×10-44 cm2 (see, e.g., the thick blue line in the top left-handed plot of Figure 2). Nowadays this bound has been beat out by the one of  LUX by a factor ~2.5 (see the blue line in the bottom right-handed plot of Figure 2). Other relevant analyses can be found in [56–68].


## 3.3.  Uncertainties

Nevertheless, since the actual range of masses and cross section critically depends on several assumptions, one always bears in mind all the uncertainties that enter in this field when interpreting the data in terms of a specific DM interaction.(i)
A first class of uncertainties comes from the poor knowledge of the properties of the DM halo. The often used Maxwell-Boltzmann halo with isotropic velocity dispersion is only a benchmark choice, not the physical description of it. Other possible choices (spherical halo with nonisotropic velocity dispersion, axisymmetric halo, triaxial halo, and the possibility to have corotating halo with different number densities) must be taken into account, in order to have a much better idea of how the velocity distribution affects the final results. In particular, it has been recently demonstrated that the uncertainties on the astrophysical part of the rate can be removed by comparing different experimental results in thevmin⁡ space (see, e.g., [69–77]). In this way one can, for instance, relax a bit the tension between the positive results and the constraints showed in Figure 2 up. However, this is not still sufficient to reconcile the complicated puzzle the experimental data have left to us (see, e.g., [77–80]).(ii)
A second class of uncertainties comes from the experimental side. For instance, the direct measurements of quenching factor are performed with reference detectors, and thus the systematic uncertainties, important for all the detectors in direct DM searches, could be larger than what we normally assumed. These must be carefully investigated, because they produce a significant shift of the allowed regions in the total cross section/DM mass plane.(iii)
Together with these main uncertainties, the third and perhaps the most important class of them comes from the nature of the DM interaction and in turn from the nuclear responses of the target nuclei. Most of the models tested in recent years are based on the assumption that the mechanism of interaction is realized through a contact interaction. Deviations from this standard phenomenological approach are interesting to study, since, from more general conditions, the DM parameter space can be analyzed. In particular, there is a concrete possibility to remove the uncertainties coming from the nature of interaction thanks to the formalism of nonrelativistic operator reviewed in Section2.2. On a more specific level, if the experimental collaborations, in addition to presenting the results in terms of a specific DM-nucleus interaction, also provide all the model independentintegrated form factors, it will be extremely useful for the community. Indeed, in this way one can easily compute the expected number of events for any kind of interactions (e.g., including isospin violating interactions, momentum-dependent form factors, and velocity dependent form factors) and compare it directly with the experimental results.Therefore the data must be treated very prudently with a maximally conservative attitude, simply because slightly modification of the assumptions summarized above can affect the theoretical interpretation of the experimental results in a relevant way. In the next section, I will focus on the uncertainties coming from the nature of interaction and in particular I am going to show how the allowed regions and constraints are modified, if an exchanged momentum dependent DM-nucleus interaction is taken into account.


## 4. Long Range Interaction

So far, theLHC has not reported any evidences of physics beyond the standard model. An optimistic point of view is that the new physics threshold is truly around the corner. A pessimistic and maybe more realistic point of view is that the  LHC results are instead telling us that the TeV scale is not a fundamental energy scale for new physics and therefore all the new theories beyond the standard model are, let us say, unnatural. Within this picture there is not any particular reason to expect new particles in the TeV mass range weakly interacting with ordinary matter. Therefore it should be time to ask ourselves whether there are other indications of the relevant energy scale. DM may play a central role in this picture. Indeed, the closeness between Ωχ and Ωb, usually referred to as a cosmic coincidence problem, might suggest a profound similarity between the dark and the ordinary sectors. Indeed, both ρχ and ρb scale as 1/a3 with the expansion of the Universe, and their ratio is independent of time. Why these two fractions are then comparable, if the ordinary and the dark sectors have a drastically different nature and different origin? A hidden parallel sector (mirror world (the idea of a mirror world was suggested before the advent of the standard model (see, e.g., [81, 82]). The idea that the mirror particles might constitute the DM of the Universe was instead discussed in [83, 84])) may shed light on this cosmic coincident problem [85]. Indeed, it is tantalizing to imagine that the dark world could be similarly complex (CP violating and asymmetric), full of forces (e.g., dark electromagnetism), and matter (dominant constituent with a mass around few GeV) that are invisible to us. For a review of mirror dark world; see, for example, [86, 87].To be more concrete the phenomenology of a complex dark sector, in which the matter fields are charged under an extraU(1) gauge group, is particularly interesting. In this case, the physics of the dark sector in itself can be as complicated as the one of our sector, providing at the same time also a feeble interaction between the two worlds. Indeed, thanks to the following renormalizable Lagrangian:
(24)ℒ=ϵϕ2FμνFμν′,
the new U(1) gauge boson ϕ may possess a small kinetic mixing ϵϕ with the ordinary photon. Here Fμν and Fμν′ are the field strength tensors of the ordinary and “dark” electromagnetism, respectively. One effect of this mixing is to cause DM particles to couple with ordinary particles with effective millicharge ϵϕe [88–90], and therefore a Rutherford-like interaction gives rise if the mass of the new U(1) gauge boson is smaller than the typical exchanged momentum in the interaction.


### 4.1. Derivation of the Main Equations

Apart from the theoretical motivations which are interesting and quite strong, what describes in the nonrelativistic limit the scattering between millicharged DM particles and the protons (notice that the interaction here maximally violates the isospin; e.g., the DM-neutron coupling is zero) is the following Yukawa potential:(25)V(r)=ϵϕ4πeeχre-mϕr,
where ϵϕe is the effective millicharge of the protons inside the nucleus felt by a DM particle with “dark” charge eχ. Here mϕ is the mass of the dark photon that acts like an electronic cloud which screens the charges of the particles involved in the scattering. The DM-proton matrix element can be then obtained by performing the Fourier transform of (25) which writes
(26)ℳp=4mpmχ∫0∞‍dr2πr2∫-1+1‍dcos⁡⁡θe-iqrcos⁡θV(r)=𝔠Yp(q2)𝒪1NR,
where q=(2m𝒩ER)1/2 is the exchanged momentum, θ is the scatter angle in the centre-of-mass frame, and 𝔠Yp(q2)=4ϵϕeeχmpmχ/(q2+mϕ2). Since 𝒪1NR=𝕝, the interaction, as already mentioned in Section 2.2, is SI like the usual case but with a coefficient which is instead dependent on the exchanged momentum in the scattering. Summing now over the total number of protons in the nucleus Z𝒩 and, if the dark sector is made of stable composite particles, over the total number of charge elementary constituents Zχ, the square of the DM-nucleus matrix element reads
(27)|ℳ𝒩|2¯=16m𝒩2mχ2(ϵϕZ𝒩eZχeχq2+mϕ2)2FHelm2(q2),
where FHelm2(q2) is the usual Helm form factor. Notice that in general this equation should be also multiplied for the form factor of the composite DM particle. For our computation we assume it to be equal to one. Plugging back the DM-nucleus matrix element in (8) the differential cross section will be
(28)dσ𝒩dER=8πm𝒩v2ϵϕ2α2kχ2(1q2+mϕ2)2Z𝒩2FHelm2(q2),
where α=e2/4π is the electromagnetic fine structure constant and kχ=Zχeχ/e is a factor that measures the strength of the DM-dark photon coupling.In order to have a rough comparison with the “standard” picture (SI contact interaction), the differential rate of nuclear recoil can be cast in the following form:(29)dR𝒩dER=ξ𝒩m𝒩ρ⊙mχ[m𝒩2μN2A𝒩2σϕγpℐ(ER)FHelm2(ER)]×𝒢(ER),
where
(30)σϕγp=16πα2ϵϕ2kχ2mN4μN2≃(ϵϕ10-5)2(kχ1)2(μN1GeV)210-40cm2
is a normalized total cross section that encapsulates all the dependences on the dark sector vertex and the kinetic mixing parameter, and the function
(31)𝒢(ER)=Z𝒩2A𝒩2[mN2q2+mϕ2]2=Z𝒩2A𝒩2[mN22m𝒩ER+mϕ2]2
measures the deviation of the allowed regions and constraints compared to the “standard” SI contact cross section. Here the function 𝒢(ER) is a sort of DM form factors which clearly exhibits two limits.(i)
Point-like limit (q≪mϕ): in this regime 𝒢(ER)=(Z𝒩/A𝒩)2mN4/mϕ4 is independent of ER and therefore the interaction is of a contact type. Indeed, the rate of nuclear recoil turns out to be proportional to α2ϵϕ2kχ2/mϕ4, which plays the same role of the Fermi’s constant in weak interaction. The factor (Z𝒩/A𝒩)2 shows up the fact that millicharged DM particles only couple with protons. Therefore, it is expected that the allowed regions and bounds in the plane (mχ,σϕγp) will be rigid shifted of a factor ~4 up with respect to the “standard” SI scenario, in which a DM particle symmetrically couples both with neutrons and protons.(ii)
Long range limit (q≫mϕ): in this regime 𝒢(ER)=(Z𝒩/A𝒩)2mN4/(4m𝒩2ER2) and therefore the differential cross section acquires an explicit dependence on the nuclear recoil energy, and a Rutherford-like cross section gives rise (dσ𝒩/dER∝1/q4). Experiments with low energy thresholds and light target nuclei (e.g.,  DAMA and  CoGeNT) are therefore more sensitive than the ones with high thresholds and heavy targets (e.g.,  XENON100). The compatibility among the experiments could therefore be improved.Considering typical nuclei (m𝒩~100GeV) and recoil energies (few keV) in the range of interest of the current experiments, the transition between the two limits is obtained for mϕ~𝒪(10) MeV. Since q2∝m𝒩, increasing m𝒩 the transition occurs at lower mϕ. Notice that once the long range regime is reached (mϕ≲10MeV), the differential cross section is independent of the mass of the mediator.


### 4.2. Results

In our analysis we consider as free parameters the normalized total cross sectionσϕγp, the mass of the dark photon mϕ, and of course one of the DM particle mχ. The velocity distribution has been assumed to be Maxwell-Boltzmann-like with a velocity dispersion v0=220km/s. For the local DM density we have chosen ρ⊙=0.3GeV/cm3. For all the positive results experiments we use a different technique for analyzing the datasets with respect to the one often found in literature. Specifically, we adopt the same approach of [91], in which the null hypothesis is tested. From this kind of statistical test we can then extract the level at which the absence of signal on top of the background in  CRESST-II and the absence of modulation in  DAMA and  CoGeNT is rejected.In Figure3 the allowed regions of the positive results experiments and constraints coming from null results in the plane (mχ,σϕγp) are shown. For a better visualization of the transition from the contact to the long range regime, we have fixed three values of the dark photon mass; namely, mϕ=1GeV (left plot corresponding to the contact limit), mϕ=30MeV (central plot), and mϕ=0 (right plot corresponding to the pure long range regime). In all plots the solid green contours individuate the regions compatible with the  DAMA annual modulation [41, 42], without considering the channeling effect. The short-dashed blue contour refers to the region derived from the  CoGeNT annual modulation signal data published in [45] (in [46, 47] the collaboration reanalyzed the data finding that the background of surface events was underestimated. In view of that we expect that the  CoGeNT allowed regions in the (mχ,σϕγp) plane will shrink around its best value, like the one reported for the “standard” SI contact interaction (see, e.g., the red contour in the top right-handed plot of Figure 2)). The dashed brown contours individuate the allowed regions compatible with the  CRESST-II excess [48]. For each experiment the contour lines cover at least one area in the (mχ,σϕγp) plane. Specifically, they refer to regions where the absence of annual modulation can be rejected at 7σ CL (outer region) and 8σ CL (inner region) for  DAMA, and 1σ CL for  CoGeNT. For  CRESST-II the absence of an excess is instead excluded at 3σ CL (outer region) and 4σ CL (inner region).Normalized total cross section as a function of the dark matter mass. For a better visualization of the transition from the contact to the long range regime, we have fixed three values of the dark photon mass; namely, mϕ=1GeV (a), mϕ=30MeV (b), and mϕ=0 (c). In all plots we show the regions where the absence of annual modulation can be rejected at 7σ CL (outer region) and 8σ CL (inner region) forDAMA (solid green contours) and 1σ CL forCoGeNT (solid green contours). ForCRESST-II (solid green contours) the absence of an excess is instead excluded at 3σ CL (outer region) and 4σ CL (inner region). Constraints, derived at 5σ CL, are shown as dashed gray lines forCDMS-Ge and magenta lines forXENON100. These results have been presented in [24, 25].
(a)
(b)
(c)Very recently, also the 11 silicon detectors of theCDMS experiment have reported an excess in their counting rate. Since the results presented in Figure 3 are based on [24] published in 2011, we do not include a full statistical analysis of the  CDMS-Si datasets. We attempt however an analytical comparison between the results arising from the interaction studied here and the standard SI contact picture (blue regions in the bottom left-handed plot of Figure 2). More precisely, defining 〈𝒢〉≡𝒢(〈ER〉), where 〈ER〉≃10keV is the average recoil energy of the three events observed [51], the  CDMS-Si allowed region in the (mχ,σSI) plane can be projected in the parameter space of Figure 3, through the relation (mχ,σϕγp)≡(mχ,σSI/〈𝒢〉). Recalling now that the best fit value of the SI cross section is σSIbest≃2×10-41 cm2, we expect that the  CDMS-Si allowed region in the (mχ,σϕγp) plane will still point towards light DM candidate, but with a normalized total cross section σϕγp=ASi2/ZSi2mϕ4/mN4σSIbest≃10-40 cm2 in the point-like regime (mϕ=1GeV), and σϕγp=ASi2/ZSi24mSi2〈ER〉2/mN4σSIbest≃3×10-47 cm2 in the long range one (mϕ=0).Constraints, derived at 5σ CL, are shown as dashed gray lines for  CDMS-Ge and magenta lines for  XENON100. In particular, in order to bracket as much as possible the uncertainties coming from the poor knowledge of the detection efficiency close to the lower threshold; for the  XENON100 bound, we adopt two approaches: (i) a case (dotted lines) in which the constraints are computed with the nominal value of the lower threshold of 4 PHE and of the effective luminosity [92], which is a sort of quenching factors in liquid xenon and (ii) another case (dashed lines) in which the constraints are computed by considering a higher threshold of 8 PHE: this is in order to determine a situation which is nearly independent of the poor knowledge of the effective luminosity (especially at very low nuclear recoil) and more important on the statistical distribution of the few PHE collected close to the threshold of the detector. For similar discussions; see [93, 94].In October 2013, also theLUX collaboration announced their first DM search results. With the data collected in just 85 live-days, they were able to set the most stringent bound on the SI contact cross section, with respect to the preexisting limits. We do not perform a full statistical analysis of their data, since, as already pointed out, the results in Figure 3 are taken from [24]. Nevertheless, since both  XENON and  LUX are based on the same double-phase xenon technology, we can estimate the bound on σϕγp, just by rescaling the magenta lines in Figure 3 with the exposures of the two experiments. Naively, we then expect that the  LUX bound will be a factor wLUX/wXENON100≃1.5 more stringent with respect to  XENON100.We can see, as expected from the discussion above, that the agreement among the positive results experiments increases moving towards long range interaction without being excluded by both theXENON100 and  LUX bounds, if our conservative choice for the lower threshold is taken into account. In particular, a common region pointing towards a DM mass around 15 GeV with a normalized total cross section σϕγp~10-46 cm2 has been found. Since the normalized total cross section σϕγp depends collectively on the kinetic mixing parameter ϵϕ and on the DM-dark photon coupling encoded in the parameter kχ (see (30)), keeping fixed, for instance, two values of kχ, namely, kχ=(1,10); the best fit for the kinetic mixing parameter would then be ϵϕ~(10-8,10-9).On the other hand, the significance of theDAMA region alone gets lower, due to the fact that for pure long range interaction (q≫mϕ) the 1/q4 drop-off of the unmodulated signal rapidly overshoots the measured total rate (see [41, Figure 1]), that of course we treat as an intrinsic constraint in our analysis. This can be appreciated in more details in Figure 4, where again the  DAMA allowed regions are shown, but in the (mχ,mϕ) plane. As one can see a 99% confidence level lower bound on the mass of the dark photon around 10 MeV is possible to be placed. This is due to the fact that the unmodulated signal, unlike the modulated one, is particularly sensitive to the energy drop-off of the cross section, which is typical for Rutherford-like scattering. The situation does not dramatically vary if we change the properties of the dark matter halo. Indeed, if we consider a more physical triaxial halo with the same velocity dispersion in the major axis (right plot of Figure 4), the allowed regions are only slightly enhanced.The bounds on the mediator mass as a function of the DM one obtained by comparing the  DAMA modulated signal with the total rate. The allowed regions are reported for CLs of 90%, 99%, and 99.99% in shaded magenta. In (a) an isothermal sphere with velocity dispersion v0=220km/s and a local DM energy density ρ⊙=0.3GeV/cm3 has been assumed. In (b) we instead show the allowed regions by considering a more physical triaxial halo with the same velocity dispersion in the major axis and with ρ⊙=0.84GeV/cm3. These results have been presented in [24, 25].
(a)
(b)Nevertheless, it is worth noticing that, in a dark world filled by more than one species in thermal equilibrium, the DM particles responsible for the annual modulation may possess smaller velocity dispersion. For instance, in the context of mirror DM, the dominant species can be the mirror hydrogen, like our world, while the mass of the DM particles favored by direct detection experiments is around 16 GeV (mirror oxigen). If they are in thermal equilibrium, the relationv0=1/2(AH/AO)1/2·220km/s≲100km/s then will hold [21]. In this case the constraints in Figure 4 are less important, because the raise of the total rate for low nuclear recoil energy is under threshold (see, e.g., [21, Figures 5(a) and 5(b)]).Since, as already stated, the normalized total cross section and in turn the rate depend collectively on the kinetic mixing parameterϵϕ and on the DM-dark photon coupling kχ; in the next section, a complementary class of constraints which are relevant for DM models featured by a long range interaction will be presented.


### 4.3.  Complementary Constraints

A first class of complementary constraints that solely depends on the properties of the dark photon (kinetic mixing parameterϵϕ and its mass mϕ) is the one coming from beam dump neutrino experiments and supernovæ observations. In Figure 5 a compilation of the bounds in the (mϕ,ϵϕ) plane is shown. Specifically, on the left panel the shaded blue regions are taken from [95, 96], while those in the right one represent the latest constraints on the dark photon properties presented in [97] (notice that the labels ϵ and mA′ in the plots of Figure 5 stand for the kinetic mixing parameter ϵϕ and the dark photon mass mϕ, resp.). As is apparent, the bounds coming from supernovæ observations (namely, the energy loss observed from SN1987a) are the most stringent. They exclude small kinetic mixing parameter and light mass for the dark photon. Indeed, in this case, if dark photons were produced in relevant amount in the centre of a supernova via mixing with the ordinary photon, the subsequent emission of them would shorten so that the predicted neutrinos boost that it would become inconsistent with the measurements made by KamiokaNDE (see, e.g., [95, 98, 99]).Allowed regions and a compilation of constraints in the (mϕ,ϵϕ) plane. (a) The upper (orange) and the lower (red) allowed regions are derived by combining the  DAMA and  CoGeNT annual modulation dataset. Specifically, the absence of modulation is here rejected at 8σ CL. The two regions differ in the choice of parameter kχ; namely kχ=1 (upper) and kχ=10 (lower). In blue are instead represented the constraints coming from beam dumped neutrino experiments and supernovæ observations taken from [95, 96]. (b) A compilation of the latest constraints in the (mϕ,ϵϕ) plane. This plot has been taken from the following review [97]. It is worth pointing out that for mϕ≲10-6 eV the parameter space of the positive results experiments in direct DM searches is basically unbounded (see, e.g., [97, Figure 7]).
(a)
(b)The allowed regions of theDAMA and  CoGeNT experiments projected in the (mϕ,ϵϕ) plane are instead shown in the left plot of Figure 5. Since, as already stated, the normalization of the rate depends not only on the properties of the dark photon but also on its coupling with the DM particles, we present the results for two benchmark values of kχ, for example, kχ=1 (upper orange region) and kχ=10 (lower orange region). We can see that, in the “symmetric” case kχ=1, only the dark photons with mϕ≳100MeV can simultaneously satisfy the constraints and give a reasonable fit of the positive results of experiments. On the other hand, for kχ≳50, the whole range is basically unbounded. It is worth stressing that such large values of kχ can be easily obtained either in models of composite DM particles with large Zχ (e.g., mirror models) or in those with a strongly coupled dark sector. Furthermore, it is also relevant to point out that, for pure long range interaction (mϕ=0), the majority of the bounds, and in particular those coming from supernovæ observations, do not apply, being the direct production of dark photon forbidden by kinematical reason (see, e.g., [97, Figure 7]).A second class of constraints which instead solely depends on the properties on the dark sector itself is the one coming from the self-interactions. Indeed, since for this class of models the DM-DM scattering, given by(32)V(r)=αkχre-mϕr,
is not suppressed by ϵϕ, the self-interactions can easily reach high values, making the dynamics of virialized astrophysical objects affected. This is particularly relevant in the limit mϕ→0; indeed, in this case the self-interaction (32) could be of the same order of the electromagnetic scattering.A first example of such bounds comes from the observations of colliding galaxy clusters, like the bullet cluster [100], which points towards collision-less DM. In particular a quite robust bound 〈σT〉/mχ≲1.25cm2/g on the size of the self-interaction has been placed by [101]. In order to translate such bound in a constraint on the dark photon mass, one has to compute the weighted energy transfer cross section, which measures the rate at which energy is transferred in the system. It reads
(33)〈σT〉=∬d3v1d3v2f(v1)f(v2)σT(vrel⁡),
where σT(vrel⁡)=∫dΩdσ/dΩ(1-cos⁡θ)‍ is the two-body energy transfer cross section. Here f(v) is the DM velocity distribution, assumed to be Maxwellian, and vrel⁡=|v→1-v→2| is the relative velocity of the DM particles involved in the scattering. Considering now the typical velocity of collision in the bullet cluster of 4700 km/s, a DM mass of 10 GeV, and two values of kχ=(1,10), the bound of 1.25 cm2/g is exceeded if the mass of the dark photon is smaller than (1,20) MeV.A second example of bounds, which is in principle relevant for this kind of models, comes from the fact that quite large self-interaction rapidly drives the DM halo into a spherical configuration, due to the prompt equipartition of the energy in the system. In particular, the efficiency of such process can be roughly estimated by(34)τrel⁡=1〈ΓT〉,
where
(35)〈ΓT〉=∬d3v1d3v2f(v1)f(v2)nχvrel⁡σT(vrel⁡)vrel⁡2v02,
which is telling us the typical time at which the self-interactions affect the dynamics of a virialized astrophysical object with number density nχ and dispersion velocity v0. Indeed, if it is smaller than the age of the object, a spherical configuration will tend to form, making such scenario excluded by few elliptical DM halos observations [102–106]. In particular, in the case of galaxy clusters we get that the relaxation time is always bigger than the age of the object (τob~1010 years), making therefore the mass of the dark photon basically unbounded. On the other hand, if we consider a smaller halo, like the one of dwarf galaxies, which are characterized by larger number density and smaller velocity dispersion, one can in principle put a very stringent lower bound on the mass of the dark photon of the order of 100 MeV (see also [107, 108]).Nevertheless, it is worth stressing that the derivations of the constraints coming from self-interaction are affected by several uncertainties both from the theoretical and experimental side. Indeed, since the phenomenology of this class of models is completely different compared to the standard one, an N-body/hydrodynamical simulation is needed. This is especially true in the case of multicomponent dark sector, in which the different equipartitions of the energy among the DM particles can generate a sort of dark electric and magnetic fields in the long range regime. Furthermore, since the self-interaction needed to change their dynamics is in general of the order of the Thomson scattering (σem~10-24 cm2), from trivial analysis dimension of such large cross section, the following rough estimation yields: the self-interaction is of a long range type in most of the virialized astrophysical objects under the assumption that the DM-dark photon coupling is of the order of α. In this case, the probability to radiate a dark photon in the scattering process is then different from zero and therefore it might well be possible that the DM sector is dissipative like ours. The time at which energy is transferred in the system (τrel⁡) is no longer a good indicator, since the relevant quantity that describes the dynamical evolution of the system is now the cooling time: in particular, for a DM sector composed by heavy and light species, the dissipation time due to the soft emission of the dark photon (dark bremsstrahlung) can be smaller than the age of the virialized astrophysical objects (see, e.g., [109–111]). In this scenario the system is no longer stable, and in general it starts to collapse. By virtue of this fact, we do not consider this last class of constraints, since dedicated and careful analysis involving also numerical simulations is needed.


## 4.1. Derivation of the Main Equations

Apart from the theoretical motivations which are interesting and quite strong, what describes in the nonrelativistic limit the scattering between millicharged DM particles and the protons (notice that the interaction here maximally violates the isospin; e.g., the DM-neutron coupling is zero) is the following Yukawa potential:(25)V(r)=ϵϕ4πeeχre-mϕr,
where ϵϕe is the effective millicharge of the protons inside the nucleus felt by a DM particle with “dark” charge eχ. Here mϕ is the mass of the dark photon that acts like an electronic cloud which screens the charges of the particles involved in the scattering. The DM-proton matrix element can be then obtained by performing the Fourier transform of (25) which writes
(26)ℳp=4mpmχ∫0∞‍dr2πr2∫-1+1‍dcos⁡⁡θe-iqrcos⁡θV(r)=𝔠Yp(q2)𝒪1NR,
where q=(2m𝒩ER)1/2 is the exchanged momentum, θ is the scatter angle in the centre-of-mass frame, and 𝔠Yp(q2)=4ϵϕeeχmpmχ/(q2+mϕ2). Since 𝒪1NR=𝕝, the interaction, as already mentioned in Section 2.2, is SI like the usual case but with a coefficient which is instead dependent on the exchanged momentum in the scattering. Summing now over the total number of protons in the nucleus Z𝒩 and, if the dark sector is made of stable composite particles, over the total number of charge elementary constituents Zχ, the square of the DM-nucleus matrix element reads
(27)|ℳ𝒩|2¯=16m𝒩2mχ2(ϵϕZ𝒩eZχeχq2+mϕ2)2FHelm2(q2),
where FHelm2(q2) is the usual Helm form factor. Notice that in general this equation should be also multiplied for the form factor of the composite DM particle. For our computation we assume it to be equal to one. Plugging back the DM-nucleus matrix element in (8) the differential cross section will be
(28)dσ𝒩dER=8πm𝒩v2ϵϕ2α2kχ2(1q2+mϕ2)2Z𝒩2FHelm2(q2),
where α=e2/4π is the electromagnetic fine structure constant and kχ=Zχeχ/e is a factor that measures the strength of the DM-dark photon coupling.In order to have a rough comparison with the “standard” picture (SI contact interaction), the differential rate of nuclear recoil can be cast in the following form:(29)dR𝒩dER=ξ𝒩m𝒩ρ⊙mχ[m𝒩2μN2A𝒩2σϕγpℐ(ER)FHelm2(ER)]×𝒢(ER),
where
(30)σϕγp=16πα2ϵϕ2kχ2mN4μN2≃(ϵϕ10-5)2(kχ1)2(μN1GeV)210-40cm2
is a normalized total cross section that encapsulates all the dependences on the dark sector vertex and the kinetic mixing parameter, and the function
(31)𝒢(ER)=Z𝒩2A𝒩2[mN2q2+mϕ2]2=Z𝒩2A𝒩2[mN22m𝒩ER+mϕ2]2
measures the deviation of the allowed regions and constraints compared to the “standard” SI contact cross section. Here the function 𝒢(ER) is a sort of DM form factors which clearly exhibits two limits.(i)
Point-like limit (q≪mϕ): in this regime 𝒢(ER)=(Z𝒩/A𝒩)2mN4/mϕ4 is independent of ER and therefore the interaction is of a contact type. Indeed, the rate of nuclear recoil turns out to be proportional to α2ϵϕ2kχ2/mϕ4, which plays the same role of the Fermi’s constant in weak interaction. The factor (Z𝒩/A𝒩)2 shows up the fact that millicharged DM particles only couple with protons. Therefore, it is expected that the allowed regions and bounds in the plane (mχ,σϕγp) will be rigid shifted of a factor ~4 up with respect to the “standard” SI scenario, in which a DM particle symmetrically couples both with neutrons and protons.(ii)
Long range limit (q≫mϕ): in this regime 𝒢(ER)=(Z𝒩/A𝒩)2mN4/(4m𝒩2ER2) and therefore the differential cross section acquires an explicit dependence on the nuclear recoil energy, and a Rutherford-like cross section gives rise (dσ𝒩/dER∝1/q4). Experiments with low energy thresholds and light target nuclei (e.g.,  DAMA and  CoGeNT) are therefore more sensitive than the ones with high thresholds and heavy targets (e.g.,  XENON100). The compatibility among the experiments could therefore be improved.Considering typical nuclei (m𝒩~100GeV) and recoil energies (few keV) in the range of interest of the current experiments, the transition between the two limits is obtained for mϕ~𝒪(10) MeV. Since q2∝m𝒩, increasing m𝒩 the transition occurs at lower mϕ. Notice that once the long range regime is reached (mϕ≲10MeV), the differential cross section is independent of the mass of the mediator.


## 4.2. Results

In our analysis we consider as free parameters the normalized total cross sectionσϕγp, the mass of the dark photon mϕ, and of course one of the DM particle mχ. The velocity distribution has been assumed to be Maxwell-Boltzmann-like with a velocity dispersion v0=220km/s. For the local DM density we have chosen ρ⊙=0.3GeV/cm3. For all the positive results experiments we use a different technique for analyzing the datasets with respect to the one often found in literature. Specifically, we adopt the same approach of [91], in which the null hypothesis is tested. From this kind of statistical test we can then extract the level at which the absence of signal on top of the background in  CRESST-II and the absence of modulation in  DAMA and  CoGeNT is rejected.In Figure3 the allowed regions of the positive results experiments and constraints coming from null results in the plane (mχ,σϕγp) are shown. For a better visualization of the transition from the contact to the long range regime, we have fixed three values of the dark photon mass; namely, mϕ=1GeV (left plot corresponding to the contact limit), mϕ=30MeV (central plot), and mϕ=0 (right plot corresponding to the pure long range regime). In all plots the solid green contours individuate the regions compatible with the  DAMA annual modulation [41, 42], without considering the channeling effect. The short-dashed blue contour refers to the region derived from the  CoGeNT annual modulation signal data published in [45] (in [46, 47] the collaboration reanalyzed the data finding that the background of surface events was underestimated. In view of that we expect that the  CoGeNT allowed regions in the (mχ,σϕγp) plane will shrink around its best value, like the one reported for the “standard” SI contact interaction (see, e.g., the red contour in the top right-handed plot of Figure 2)). The dashed brown contours individuate the allowed regions compatible with the  CRESST-II excess [48]. For each experiment the contour lines cover at least one area in the (mχ,σϕγp) plane. Specifically, they refer to regions where the absence of annual modulation can be rejected at 7σ CL (outer region) and 8σ CL (inner region) for  DAMA, and 1σ CL for  CoGeNT. For  CRESST-II the absence of an excess is instead excluded at 3σ CL (outer region) and 4σ CL (inner region).Normalized total cross section as a function of the dark matter mass. For a better visualization of the transition from the contact to the long range regime, we have fixed three values of the dark photon mass; namely, mϕ=1GeV (a), mϕ=30MeV (b), and mϕ=0 (c). In all plots we show the regions where the absence of annual modulation can be rejected at 7σ CL (outer region) and 8σ CL (inner region) forDAMA (solid green contours) and 1σ CL forCoGeNT (solid green contours). ForCRESST-II (solid green contours) the absence of an excess is instead excluded at 3σ CL (outer region) and 4σ CL (inner region). Constraints, derived at 5σ CL, are shown as dashed gray lines forCDMS-Ge and magenta lines forXENON100. These results have been presented in [24, 25].
(a)
(b)
(c)Very recently, also the 11 silicon detectors of theCDMS experiment have reported an excess in their counting rate. Since the results presented in Figure 3 are based on [24] published in 2011, we do not include a full statistical analysis of the  CDMS-Si datasets. We attempt however an analytical comparison between the results arising from the interaction studied here and the standard SI contact picture (blue regions in the bottom left-handed plot of Figure 2). More precisely, defining 〈𝒢〉≡𝒢(〈ER〉), where 〈ER〉≃10keV is the average recoil energy of the three events observed [51], the  CDMS-Si allowed region in the (mχ,σSI) plane can be projected in the parameter space of Figure 3, through the relation (mχ,σϕγp)≡(mχ,σSI/〈𝒢〉). Recalling now that the best fit value of the SI cross section is σSIbest≃2×10-41 cm2, we expect that the  CDMS-Si allowed region in the (mχ,σϕγp) plane will still point towards light DM candidate, but with a normalized total cross section σϕγp=ASi2/ZSi2mϕ4/mN4σSIbest≃10-40 cm2 in the point-like regime (mϕ=1GeV), and σϕγp=ASi2/ZSi24mSi2〈ER〉2/mN4σSIbest≃3×10-47 cm2 in the long range one (mϕ=0).Constraints, derived at 5σ CL, are shown as dashed gray lines for  CDMS-Ge and magenta lines for  XENON100. In particular, in order to bracket as much as possible the uncertainties coming from the poor knowledge of the detection efficiency close to the lower threshold; for the  XENON100 bound, we adopt two approaches: (i) a case (dotted lines) in which the constraints are computed with the nominal value of the lower threshold of 4 PHE and of the effective luminosity [92], which is a sort of quenching factors in liquid xenon and (ii) another case (dashed lines) in which the constraints are computed by considering a higher threshold of 8 PHE: this is in order to determine a situation which is nearly independent of the poor knowledge of the effective luminosity (especially at very low nuclear recoil) and more important on the statistical distribution of the few PHE collected close to the threshold of the detector. For similar discussions; see [93, 94].In October 2013, also theLUX collaboration announced their first DM search results. With the data collected in just 85 live-days, they were able to set the most stringent bound on the SI contact cross section, with respect to the preexisting limits. We do not perform a full statistical analysis of their data, since, as already pointed out, the results in Figure 3 are taken from [24]. Nevertheless, since both  XENON and  LUX are based on the same double-phase xenon technology, we can estimate the bound on σϕγp, just by rescaling the magenta lines in Figure 3 with the exposures of the two experiments. Naively, we then expect that the  LUX bound will be a factor wLUX/wXENON100≃1.5 more stringent with respect to  XENON100.We can see, as expected from the discussion above, that the agreement among the positive results experiments increases moving towards long range interaction without being excluded by both theXENON100 and  LUX bounds, if our conservative choice for the lower threshold is taken into account. In particular, a common region pointing towards a DM mass around 15 GeV with a normalized total cross section σϕγp~10-46 cm2 has been found. Since the normalized total cross section σϕγp depends collectively on the kinetic mixing parameter ϵϕ and on the DM-dark photon coupling encoded in the parameter kχ (see (30)), keeping fixed, for instance, two values of kχ, namely, kχ=(1,10); the best fit for the kinetic mixing parameter would then be ϵϕ~(10-8,10-9).On the other hand, the significance of theDAMA region alone gets lower, due to the fact that for pure long range interaction (q≫mϕ) the 1/q4 drop-off of the unmodulated signal rapidly overshoots the measured total rate (see [41, Figure 1]), that of course we treat as an intrinsic constraint in our analysis. This can be appreciated in more details in Figure 4, where again the  DAMA allowed regions are shown, but in the (mχ,mϕ) plane. As one can see a 99% confidence level lower bound on the mass of the dark photon around 10 MeV is possible to be placed. This is due to the fact that the unmodulated signal, unlike the modulated one, is particularly sensitive to the energy drop-off of the cross section, which is typical for Rutherford-like scattering. The situation does not dramatically vary if we change the properties of the dark matter halo. Indeed, if we consider a more physical triaxial halo with the same velocity dispersion in the major axis (right plot of Figure 4), the allowed regions are only slightly enhanced.The bounds on the mediator mass as a function of the DM one obtained by comparing the  DAMA modulated signal with the total rate. The allowed regions are reported for CLs of 90%, 99%, and 99.99% in shaded magenta. In (a) an isothermal sphere with velocity dispersion v0=220km/s and a local DM energy density ρ⊙=0.3GeV/cm3 has been assumed. In (b) we instead show the allowed regions by considering a more physical triaxial halo with the same velocity dispersion in the major axis and with ρ⊙=0.84GeV/cm3. These results have been presented in [24, 25].
(a)
(b)Nevertheless, it is worth noticing that, in a dark world filled by more than one species in thermal equilibrium, the DM particles responsible for the annual modulation may possess smaller velocity dispersion. For instance, in the context of mirror DM, the dominant species can be the mirror hydrogen, like our world, while the mass of the DM particles favored by direct detection experiments is around 16 GeV (mirror oxigen). If they are in thermal equilibrium, the relationv0=1/2(AH/AO)1/2·220km/s≲100km/s then will hold [21]. In this case the constraints in Figure 4 are less important, because the raise of the total rate for low nuclear recoil energy is under threshold (see, e.g., [21, Figures 5(a) and 5(b)]).Since, as already stated, the normalized total cross section and in turn the rate depend collectively on the kinetic mixing parameterϵϕ and on the DM-dark photon coupling kχ; in the next section, a complementary class of constraints which are relevant for DM models featured by a long range interaction will be presented.


## 4.3.  Complementary Constraints

A first class of complementary constraints that solely depends on the properties of the dark photon (kinetic mixing parameterϵϕ and its mass mϕ) is the one coming from beam dump neutrino experiments and supernovæ observations. In Figure 5 a compilation of the bounds in the (mϕ,ϵϕ) plane is shown. Specifically, on the left panel the shaded blue regions are taken from [95, 96], while those in the right one represent the latest constraints on the dark photon properties presented in [97] (notice that the labels ϵ and mA′ in the plots of Figure 5 stand for the kinetic mixing parameter ϵϕ and the dark photon mass mϕ, resp.). As is apparent, the bounds coming from supernovæ observations (namely, the energy loss observed from SN1987a) are the most stringent. They exclude small kinetic mixing parameter and light mass for the dark photon. Indeed, in this case, if dark photons were produced in relevant amount in the centre of a supernova via mixing with the ordinary photon, the subsequent emission of them would shorten so that the predicted neutrinos boost that it would become inconsistent with the measurements made by KamiokaNDE (see, e.g., [95, 98, 99]).Allowed regions and a compilation of constraints in the (mϕ,ϵϕ) plane. (a) The upper (orange) and the lower (red) allowed regions are derived by combining the  DAMA and  CoGeNT annual modulation dataset. Specifically, the absence of modulation is here rejected at 8σ CL. The two regions differ in the choice of parameter kχ; namely kχ=1 (upper) and kχ=10 (lower). In blue are instead represented the constraints coming from beam dumped neutrino experiments and supernovæ observations taken from [95, 96]. (b) A compilation of the latest constraints in the (mϕ,ϵϕ) plane. This plot has been taken from the following review [97]. It is worth pointing out that for mϕ≲10-6 eV the parameter space of the positive results experiments in direct DM searches is basically unbounded (see, e.g., [97, Figure 7]).
(a)
(b)The allowed regions of theDAMA and  CoGeNT experiments projected in the (mϕ,ϵϕ) plane are instead shown in the left plot of Figure 5. Since, as already stated, the normalization of the rate depends not only on the properties of the dark photon but also on its coupling with the DM particles, we present the results for two benchmark values of kχ, for example, kχ=1 (upper orange region) and kχ=10 (lower orange region). We can see that, in the “symmetric” case kχ=1, only the dark photons with mϕ≳100MeV can simultaneously satisfy the constraints and give a reasonable fit of the positive results of experiments. On the other hand, for kχ≳50, the whole range is basically unbounded. It is worth stressing that such large values of kχ can be easily obtained either in models of composite DM particles with large Zχ (e.g., mirror models) or in those with a strongly coupled dark sector. Furthermore, it is also relevant to point out that, for pure long range interaction (mϕ=0), the majority of the bounds, and in particular those coming from supernovæ observations, do not apply, being the direct production of dark photon forbidden by kinematical reason (see, e.g., [97, Figure 7]).A second class of constraints which instead solely depends on the properties on the dark sector itself is the one coming from the self-interactions. Indeed, since for this class of models the DM-DM scattering, given by(32)V(r)=αkχre-mϕr,
is not suppressed by ϵϕ, the self-interactions can easily reach high values, making the dynamics of virialized astrophysical objects affected. This is particularly relevant in the limit mϕ→0; indeed, in this case the self-interaction (32) could be of the same order of the electromagnetic scattering.A first example of such bounds comes from the observations of colliding galaxy clusters, like the bullet cluster [100], which points towards collision-less DM. In particular a quite robust bound 〈σT〉/mχ≲1.25cm2/g on the size of the self-interaction has been placed by [101]. In order to translate such bound in a constraint on the dark photon mass, one has to compute the weighted energy transfer cross section, which measures the rate at which energy is transferred in the system. It reads
(33)〈σT〉=∬d3v1d3v2f(v1)f(v2)σT(vrel⁡),
where σT(vrel⁡)=∫dΩdσ/dΩ(1-cos⁡θ)‍ is the two-body energy transfer cross section. Here f(v) is the DM velocity distribution, assumed to be Maxwellian, and vrel⁡=|v→1-v→2| is the relative velocity of the DM particles involved in the scattering. Considering now the typical velocity of collision in the bullet cluster of 4700 km/s, a DM mass of 10 GeV, and two values of kχ=(1,10), the bound of 1.25 cm2/g is exceeded if the mass of the dark photon is smaller than (1,20) MeV.A second example of bounds, which is in principle relevant for this kind of models, comes from the fact that quite large self-interaction rapidly drives the DM halo into a spherical configuration, due to the prompt equipartition of the energy in the system. In particular, the efficiency of such process can be roughly estimated by(34)τrel⁡=1〈ΓT〉,
where
(35)〈ΓT〉=∬d3v1d3v2f(v1)f(v2)nχvrel⁡σT(vrel⁡)vrel⁡2v02,
which is telling us the typical time at which the self-interactions affect the dynamics of a virialized astrophysical object with number density nχ and dispersion velocity v0. Indeed, if it is smaller than the age of the object, a spherical configuration will tend to form, making such scenario excluded by few elliptical DM halos observations [102–106]. In particular, in the case of galaxy clusters we get that the relaxation time is always bigger than the age of the object (τob~1010 years), making therefore the mass of the dark photon basically unbounded. On the other hand, if we consider a smaller halo, like the one of dwarf galaxies, which are characterized by larger number density and smaller velocity dispersion, one can in principle put a very stringent lower bound on the mass of the dark photon of the order of 100 MeV (see also [107, 108]).Nevertheless, it is worth stressing that the derivations of the constraints coming from self-interaction are affected by several uncertainties both from the theoretical and experimental side. Indeed, since the phenomenology of this class of models is completely different compared to the standard one, an N-body/hydrodynamical simulation is needed. This is especially true in the case of multicomponent dark sector, in which the different equipartitions of the energy among the DM particles can generate a sort of dark electric and magnetic fields in the long range regime. Furthermore, since the self-interaction needed to change their dynamics is in general of the order of the Thomson scattering (σem~10-24 cm2), from trivial analysis dimension of such large cross section, the following rough estimation yields: the self-interaction is of a long range type in most of the virialized astrophysical objects under the assumption that the DM-dark photon coupling is of the order of α. In this case, the probability to radiate a dark photon in the scattering process is then different from zero and therefore it might well be possible that the DM sector is dissipative like ours. The time at which energy is transferred in the system (τrel⁡) is no longer a good indicator, since the relevant quantity that describes the dynamical evolution of the system is now the cooling time: in particular, for a DM sector composed by heavy and light species, the dissipation time due to the soft emission of the dark photon (dark bremsstrahlung) can be smaller than the age of the virialized astrophysical objects (see, e.g., [109–111]). In this scenario the system is no longer stable, and in general it starts to collapse. By virtue of this fact, we do not consider this last class of constraints, since dedicated and careful analysis involving also numerical simulations is needed.


## 5. Conclusion

Direct DM searches are now characterized by tantalizing results and hints that make this field very active both from the theoretical and experimental side. In particular, in addition to the long standingDAMA results, nowadays there are other experiments, like  CoGeNT,  CRESST-II, and  CDMS-Si that are starting to observe anomalies in their counting rates. On the other hand, the increasingly stringent constraints coming from null result experiments put in serious trouble the theoretical interpretation of the data, at least in terms of the simple-minded SI contact interaction. In this work I discussed thestatus of direct DM detection with specific attention to the experimental results and their phenomenological interpretation in terms of DM interaction. In particular, in the first part I presented a new and more general approach to study signals in this field based on nonrelativistic operators. Then I reviewed the experimental results and their interpretation in terms of the “standard” SI interaction pointing out all the uncertainties which enter in this field. In the last part of this work, I investigated a fermionic dark matter particle carrying a small millicharged and analyzed its impact on direct detection experiments. I showed that this kind of long range interaction can accommodate the positive experimental results. By assuming a conservative choice for the lower threshold of the  XENON100 and  LUX experiments I have demonstrated that this candidate is not ruled out. I also determined the complementary class of constraints which are relevant for millicharged DM particles with long range forces.Finally, I would like to propose a possible direction to pursue in order to make sense of the current exciting experimental panorama based on the formalism of nonrelativistic operators. Indeed, as we have seen in the first part of this work, it allows us to describe the DM-nucleus interactions in terms of a very limited number of relevant degrees of freedom. In this way, it is possible to parametrize the model-dependent part of the rate from the model-independent one encapsulated in a sort ofintegrated form factors that encode all the dependences on the astrophysics, nuclear physics, and experimental details. Therefore, since one is ignorant or agnostic about the underlying theory, I would like to encourage a synergy between nuclear physicists and experimentalists in order to provide a complete set ofintegrated form factors defined in (23). It would be extremely useful for the community, because, in this way, one can compute the expected number of events for any kind of interaction (e.g., including isospin-violating interactions, momentum-dependent form factors, and velocity-dependent form factors) and compare it directly with the experimental results. Providing theintegrated form factors will thus be the first step towards a model-independent analysis in direct DM searches.



---

*Source: 681312-2014-04-16.xml*