# Halo-independent comparison of direct dark matter detection data.

1. IntroductionThe presence of dark matter (DM) in the universe is now an established fact that has been confirmed once more by the recent precise measurements of the Planck satellite [1]. Many different particle candidates exist as possible explanations for the DM. A particular class of candidates, the WIMPs (for weakly interacting massive particles), is very actively searched for. WIMPs are particles with weakly interacting cross sections and masses in the 1 GeV/[c.sup.2]-10 TeV/[c.sup.2] range. Of particular interest are light WIMPs, with mass around 1-10 GeV/[c.sup.2].

At present, four direct dark matter search experiments (DAMA [2], CoGeNT [3-5], CRESST-II [6], and CDMS-II-Si [7]) have data that may be interpreted as signals from DM particles in the light wImPs range. DAMA [2] and CoGeNT [4] report annual modulations in their event rates, compatible with those expected for a DM signal [8, 9]. CoGeNT [3, 5], CRESST-II [6], and CDMS-II-Si [7] observe an excess of events above their expected backgrounds that may be interpreted as due to DM WIMPs.

However, other experiments do not observe significant excesses above their estimated background, thus setting upper limits on the interaction of WIMPs with nuclei. The most stringent limits on the average (unmodulated) rate for light WIMPs are set by the LUX [10], XENON10 [11], XENON100 [12], CDMS-II-Ge [13], and CDMSlite [14] experiments, with the addition of SIMPLE [15], PICASSO [16], and COUPP [17] for spin-dependent and isospin-violating interactions. CDMS-II-Ge [18] also constrains directly the amplitude of an annually modulated signal.

In order to compare a model for WIMPs with data from direct DM detection experiments, one needs to assume a value for the DM local density and velocity distribution in our galaxy. The Standard Halo Model (SHM) is usually assumed for the DM halo, corresponding to a truncated Maxwell-Boltzmann distribution for the DM velocity (see e.g., [19]). However, the parameters of this model are not known to great accuracy, and the model itself is not supported by data. Actually, quantitatively different velocity distributions are obtained from numerical simulations (see e.g., [20]). Various models and parameterizations for the DM velocity distribution in our galaxy have been proposed as alternatives to the SHM, either derived from astrophysical data or from N-body simulations (see e.g., [9] and references therein). Other authors have attempted to estimate the uncertainty in the determination of the properties of the DM halo and to quantify its effects on the interpretation of DM direct detection data (see e.g., [21-25]). Another approach is that of marginalizing over the parameters of the DM halo when computing bounds and allowed regions from the experimental data (see e.g., [26]). However, all these procedures maintain a certain degree of model dependence, for example, in the choice of the functional form of the parameterization of the halo. It is very important to notice here that the high velocity tail of the DM velocity distribution plays a crucial role in determining the number of DM particles that is above threshold for a given experiment, and therefore a way to analyze the data without the need to make any assumption on its shape is highly desirable.

The problem of comparing results from different direct detection experiments can indeed be formulated without the need to assume a velocity profile for the DM [27-34]. The basic idea is to factor out from the formulas used to compute the scattering rate all the astrophysical quantities such as the DM velocity distribution function. In this way the rate can be computed, for any model of particle interactions between the DM and the nuclei in the detector, with no need to assume a velocity profile for the DM, while rather allowing using the experimental data to constrain the unknown quantities. Such a "halo-independent" analysis was first proposed in [27], where many of the features of the method were presented, and was further developed in [28] which extended the analysis to annual modulations and [29] which showed how to include detector resolutions. The method was further generalized in [32] to more complicate particle interactions; that is, those where the scattering cross section has a nontrivial dependence on the DM velocity.

The halo-independent analysis is particularly useful to investigate the compatibility of the different experimental results in the light WIMP hypothesis, for which the details of the DM velocity distribution, especially at high velocities, are notably relevant. Here we will review this method, applying it to spin-independent interactions with both isospin-conserving and isospin-violating [35, 36] couplings, and to WIMPs with magnetic dipole moment. We will compare data from DAMA [2], CoGeNT [4], CRESST-II [6], CDMS-II-Si [7], CDMS-II-Ge low threshold analysis [13], CDMS-II-Ge annual modulation analysis [18], CDMSlite [14], XENON10 S2-only analysis [11], XENON100 [12], LUX [10], and SIMPLE [15], following the analysis described in [31-33]. This review summarizes the results presented in [31-33].

2. The Scattering Rate

What is observed at direct DM detection experiments is the WIMP-nucleus differential scattering rate, usually measured in units of counts/kg/day/keV. For a target nuclide T initially at rest, recoiling with energy ER after the scattering with a WIMP with mass m and initial velocity v, the differential rate is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)

Here [m.sub.T] is the target nuclide mass and [C.sub.T] is its mass fraction in the detector, and we denoted by v = [absolute value of v] the WIMP speed. d[[sigma].sub.T]/d[E.sub.R] is the differential scattering cross section. The dependence of the rate on the local characteristics of the DM halo is contained in the local DM density [rho] and the DM velocity distribution in the Earth's frame f(v, t), which is modulated in time due to Earth's rotation around the Sun [8, 9]. The distribution f(v, t) is normalized to [integral] [d.sup.3]vf(v, t) = 1. In the velocity integral, [v.sub.min]([E.sub.R]) is the minimum speed required for the incoming DM particle to cause a nuclear recoil with energy [E.sub.R]. For an elastic collision

[v.sub.min] = [square root of ([m.sub.T][E.sub.R]/2[[mu].sup.2.sub.T])], (2)

where [[mu].sub.T] = [mm.sub.T]/(m + [m.sub.T]) is the WIMP-nucleus reduced mass.

To properly reproduce the recoil rate measured by experiments, we need to take into account the characteristics of the detector. Most experiments do not measure the recoil energy directly but rather a detected energy E', often quoted in keVee (keV electron-equivalent) or in photoelectrons. The uncertainties and fluctuations in the detected energy corresponding to a particular recoil energy are expressed in a (target nuclide and detector dependent) resolution function [G.sub.T]([E.sub.R], E') that gives the probability that a recoil energy [E.sub.R] (usually quoted in keVnr for nuclear recoils) is measured as E'. The resolution function is often (but not always as the XENON and LUX experiments are a notable exception) approximated by a Gaussian distribution. It incorporates the mean value <E'> = [Q.sub.T][E.sub.R], which depends on the energy dependent quenching factor [Q.sub.T]([E.sub.R]), and the energy resolution [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Moreover, experiments have one or more counting efficiencies or cut acceptances, denoted here by [[epsilon].sub.1](E') and [[epsilon].sub.2]([E.sub.R]), which also affect the measured rate. Thus the nuclear recoil rate in (1) must be convolved with the function [[epsilon].sub.l](E')[[epsilon].sub.2]([E.sub.R])[G.sub.T]([E.sub.R], E'). The resulting differential rate as a function of the detected energy E' is

dR/dE' = [[epsilon].sub.1] (E') [summation over T] [[integral].sup.[infinity].sub.0] d[E.sub.R][[epsilon].sub.2] ([E.sub.R])[G.sub.T] ([E.sub.R], E') [d[R.sub.T]/d[E.sub.R]]. (3)

The rate within a detected energy interval [[E'.sub.1], [E'.sub.2]] follows as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

The time dependence of the rate (4) is generally well approximated by the first terms of a harmonic series,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (5)

where [t.sub.0] is the time of the maximum of the signal and [omega] = 2[pi]/yr. The coefficients [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are, respectively, the unmodulated and modulated components of the rate in the energy interval [[E'.sub.1], [E'.sub.2]].

3. Halo-Independent Method for Spin-Independent Interaction

The differential cross section for the usual spin-independent (SI) interaction is

d[[sigma.sub.T]/d[E.sub.R] = [[sigma].sup.SI.sub.T] ([E.sub.R]) [[m.sub.T]/2[[mu].sup.2.sub.T][v.sup.2]], (6)

with

[[sigma].sup.SI.sub.T] ([E.sub.R]) = [[sigma].sub.p] [[[mu].sup.2.sub.T]/[[mu].sup.2.sub.p]] [[[Z.sub.T] + [([A.sub.T] - [Z.sub.T]) [f.sub.n]/[f.sub.p]]].sup.2] [F.sup.2.sub.SI,T]([E.sub.R]). (7)

Here [Z.sub.T] and [A.sub.T] are, respectively, the atomic and mass number of the target nuclide T, [F.sub.SI,T]([E.sub.R]) is the nuclear spin- independent form factor (which we take to be the Helm form factor [37] normalized to [F.sub.SI,T](0) = 1), [f.sub.n] and [f.sub.p] are the effective WIMP couplings to neutron and proton, and [[mu].sub.p] is the WIMP-proton reduced mass. The WIMP-proton cross section [[sigma].sub.p] is the parameter customarily chosen to be constrained together with the WIMP mass m for SI interactions, as it does not depend on the detector, and thus bounds and allowed regions from different experiments can be compared on the same plot.

The isospin-conserving coupling [f.sub.n] = [f.sub.p] is usually assumed by the experimental collaborations. The isospin-violating coupling [f.sub.n]/[f.sub.p] = -0.7 [35, 36] produces the maximum cancellation in the expression inside the square bracket in (7) for xenon, thus highly suppressing the interaction cross section. This suppression is phenomenologically interesting because it weakens considerably the bounds from xenon-based detectors such as XENON and LUX which provide some of the most restrictive bounds.

Using this expression for the differential cross section and changing integration variable from [E.sub.R] to [v.sub.min] through (2), we can rewrite (4) as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (8)

where the velocity integral [??] is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9)

and we defined the response function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] ([v.sub.min]) for WIMPs with SI interactions as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

Introducing the speed distribution

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII (11)

we can rewrite the [??] function as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (12)

The velocity integral [??]([v.sub.min], t) has an annual modulation due to Earth's rotation around the Sun and can be separated into its unmodulated and modulated components as was done for the rate in 5)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)

Once the WIMP mass and interactions are fixed, the functions [[??].sup.0]([v.sub.min]) and [[??].sup.1]([v.sub.min]) are detector-independent quantities that must be common to all nondirectional direct DM experiments. Thus we can map the rate measurements and bounds of different experiments into measurements of and bounds on if [[??].sup.0]([v.sub.min]) and [[??].sup.1]([v.sub.min]) as functions of [v.sub.min].

For experiments with putative DM signals, in light of (8) we may interpret the measured rates [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in an energy interval [[E'.sub.1], [E'.sub.2]] as averages of the [[??].sup.i]([v.sub.min]) functions weighted by the response function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (14)

with i = 0,1 for the unmodulated and modulated component, respectively. Each such average corresponds to a point with error bars in the ([v.sub.min], [??]) plane. The vertical bars are given by [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] computed by replacing [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] with [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in (14). The [DELTA][R.sup.i] used here correspond to the 68% confidence interval. The horizontal bar shows the [v.sub.min] interval where the response function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] for the given experiment is sufficiently different from zero. Following [28, 29, 31, 33] the horizontal bar may be chosen to extend over the interval [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is the energy resolution and the function [v.sub.min](E') is obtained from [v.sub.min]([E.sub.R]) in (2) by using the recoil energy ER that produces the mean <E'> which is equal to the measured energy E'. When isotopes of the same element are present, like for Xe or Ge, the vmin intervals of the different isotopes almost completely overlap, and we take [v.sub.min,1] and [v.sub.min,2] to be the [C.sub.T]-weighted averages over the isotopes of the element. When there are nuclides belonging to very different elements, like Ca and O in CRESST-II, a more complicated procedure should be followed (see 28, 29] for details).

To determine the upper bounds on the unmodulated part of [??] set by experimental upper bounds on the unmodulated part of the rate, a procedure first outlined in [27, 28] maybe used. This procedure exploits the fact that, by definition, [[??].sup.0] is a nonincreasing function of [v.sub.min]. For this reason, the smallest possible [[??].sup.0](vmin) function passing by a fixed point ([v.sub.0], [[??].sub.0]) in the ([v.sub.min], [??]) plane is the downward step-function [[??].sub.0] [theta]([v.sub.0] - [v.sub.min]). In other words, among the functions passing by the point ([v.sub.0], [[??].sub.0]), the downward step is the function yielding the minimum predicted number of events. Imposing this functional form in (8) we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

The upper bound [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] on the unmodulated rate in an interval [[E'.sub.1], [E'.sub.2]] is translated into an upper bound [[eta].sub.lim]([v.sub.min]) on if [[??].sup.0] at [v.sub.0] by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

The upper bound so-obtained is conservative in the sense that any [[??].sup.0] function extending even partially above [[??].sup.lim] is excluded, but not every [[??].sup.0] function lying everywhere below [[??].sup.lim] is allowed [28].

The procedure just described does not assume any particular property of the DM halo. By making some assumptions, more stringent limits on the modulated part l can be derived from the limits on the unmodulated part of the rate (see [28, 38-40]), but we choose to proceed without making any assumption on the DM halo.

Figures 1 and 2 collect the results of the halo-independent analysis for a WIMP mass m = 7 GeV/[c.sup.2] and m = 9 GeV/[c.sup.2], respectively; the left and right columns correspond to isospin-conserving ([f.sub.n] = [f.sub.p]) and isospin-violating ([f.sub.n]/[f.sub.p] = -0.7) interactions, respectively; and the top, middle, and bottom rows show measurements and bounds for the unmodulated component [[??].sup.0] [c.sup.2], for the modulated component [[??].sup.1] [c.sup.2], and for both together, respectively, in units of [day.sup.-1]. The middle row also shows the upper bounds on [[??].sup.0] [c.sup.2] from the plots on the top row.

The bounds from CDMS-II-Ge, CDMS-II-Si, CDMSlite, XENON10, XENON100, and LUX are derived as 90% CL upper bounds using the maximum gap method [41]. The SIMPLE bound is derived as the 90% CL Poisson limit. The crosses show the DAMA modulation signal (green crosses), CoGeNT modulated (blue crosses) and unmodulated signal (plus an unknown flat background, dark red horizontal lines), CRESST-II, and CDMS-II-Si unmodulated signals (black and red crosses, resp.); the CDMS-II-Ge modulation bound is shown as a dark grey horizontal line with downward arrow. Only sodium is considered for DAMA (with quenching factor [Q.sub.Na] = 0.3), as for the DM masses considered here the WIMP scattering off iodine is supposed to be below threshold. For XENON10, limits produced by setting or not setting the electron yield [Q.sub.y] to zero below 1.4 keVnr (as in [12]) are obtained (solid and dashed orange line, resp.). For LUX, upper bounds considering 0, 1, 3, and 5 observed events are computed [33], corresponding (from bottom to top) to the magenta lines with different dashing styles in Figures 1 and 2.

The overlapping of the green and blue crosses in Figures 1 and 2 seems to indicate that the DAMA and CoGeNT modulation data are compatible one with the other. On the other hand, the three CDMS-II-Si points overlap or are below the CoGeNT and DAMA measurements of the modulated part of [??]. Thus, interpreted as a measurement of the unmodulated rate, the three CDMS-II-Si data points seem largely incompatible with the modulation of the signal observed by CoGeNT and DAMA, since a modulated signal is expected to be much smaller than the respective unmodulated component. For isospin-conserving interactions (left column of Figures 1 and 2), the experiments with a positive signal seem largely incompatible with the limits set by the other experiments, most notably by LUX, and by XENON10 at low [v.sub.min] values. In fact, only the DAMA, CoGeNT, and CDMS-II-Si data points at very low Vmin are below all the bounds. The compatibility of the DAmA, CoGeNT, and CRESST-II data with the exclusion bounds improves slightly for isospin-violating couplings with [f.sub.n]/[f.sub.p] = -0.7, for which the XENON and LUX limits are weakened (right column of Figures 1 and 2). However, still only the points at low [v.sub.min] are below the exclusion lines, while the DAMA and CoGeNT modulated data at high [v.sub.min] are now mostly excluded by the CDMS-II-Ge modulation bound. The improvement is better for the three CDMS-II-Si data that pass all the limits for DM with isospin-violating couplings.

In the top left panel of Figure 1 and the middle right panel of Figure 2, we show the predicted [[??].sup.0] [c.sup.2] (upper line) and [[??].sup.1] [c.sup.2] (lower line) from the SHM (assuming [v.sub.0] = 220 km/s for the DM velocity dispersion and [v.sub.esc] = 544 km/s for the galactic escape velocity), for particular values of the WIMP-proton cross section. We choose these cross sections so that (a) for the m = 7 GeV/[c.sup.2] isospin-conserving case in the top left panel of Figure 1, the CDMS-II-Si unmodulated data are well explained by the SHM; that is, the [[??].sup.0] curve passes through the red crosses in the figure (this happens for [[sigma].sub.p] = [10.sup.-40] [cm.sup.2]) and (b) for the m = 9 GeV/[c.sup.2] isospin-violating case in the middle right panel of Figure 2, the DAMA modulation data are well explained; that is, [[??].sup.1] passes through the green crosses (this occurs for [[sigma].sub.p] = 2 x [10.sup.-38] [cm.sup.2]). These plots show that the [[??].sup.0] ([v.sub.min]) of the SHM is a very steep function of [v.sub.min] and thus can be constrained at low [v.sub.min] values by the CDMSlite limit as well as by other upper limits on the unmodulated rate.

The procedure outlined in this section to compare data from different experiments in a halo-independent way can only be applied when the differential cross section can be factorized into a velocity dependent term, independent of the detector (e.g., it must be independent of [m.sub.T]), times a velocity independent term containing all the detector dependency. The differential cross section in (6) and (7) is of this form. In the case of a more general form of the differential cross section; instead, we can proceed as described in the following section.

4. Generalized Halo-Independent Method

Here we show how to define the response function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in (8) so that the halo- independent analysis can be extended to any type of interaction. Changing the order of the v and [E.sub.R] integrations in (4), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)

Here [E.sup.max.sub.R] = 2 [[mu].sup.2.sub.T] [v.sup.2]/[m.sub.T] is the maximum recoil energy a WIMP of speed v can impart in an elastic collision to a target nucleus T initially at rest. To make contact with the SI interaction method of the previous section, we have multiplied and divided by the factor [[sigma].sub.ref]/[v.sup.2], where [[sigma].sub.ref] is a target-independent reference cross section (i.e., a constant with the dimensions of a cross section) that coincides with [[sigma].sub.p] for SI interactions. In compact form, (17) reads

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (18)

where in analogy with (9) we defined

[??](v,t) [equivalent] [[rho][[sigma].sub.ref]/m] f(v,t), (19)

and we defined the "integrated response function" (the name stemming from (26))

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)

For simplicity, we only consider differential cross sections and thus integrated response functions that depend only on the speed v = [absolute value of v] and not on the whole velocity vector. This is true if the DM flux and the target nuclei are unpolarized and the detection efficiency is isotropic throughout the detector, which is the most common case. With this restriction,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (21)

where the function [??] is defined as in (11). We now define the function [??](v, t) by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (22)

with [??](v, t) going to zero in the limit of v going to infinity. This yields the usual definition of [??] (see (9) and 12))

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (23)

Using (22) in (21), the energy integrated rate becomes

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (24)

Integration by parts of (24) leads to an equation formally identical to (8) but which is now valid for any interaction

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)

The response function is now defined as the derivative of the "integrated response function" [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (26)

Notice that the boundary term in the integration by parts of (24) is zero because the definition of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in (20) imposes that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Similarly to what we did earlier for the SI interaction, we want again to compare average values of the if functions with upper limits. However, for a differential cross section with a general dependence on the DM velocity, it might not be possible to simply use (14) with [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] replaced by [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] to assign a weighted average of [[??].sup.0] or [[??].sup.1] to a finite [v.sub.min] range. This may happen because the width of the response function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in (26) at large [v.sub.min] is dictated by the high speed behavior of the differential cross section, and it might even be infinite. For example, if ([v.sup.2]d[[sigma].sub.T]/d[E.sub.R]) goes as [v.sup.n] for large v, with n a positive integer, then [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] also goes as [v.sup.n] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] goes as [v.sup.n-1.sub.m] for large [v.sub.min]. Thus, if n [greater than or equal to] l, the response function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] does not vanish for large [v.sub.min]. This implies that the denominator in (14) diverges.

We can regularize the behavior of the response function at large [v.sub.min] by using, for example, the function [v.sup.r.sub.min] [??]([v.sub.min], t) with integer r [greater than or equal to] n, instead of just [??] ([v.sub.min], t). Since this new function is common to all experiments, we can use it to compare the data in [v.sub.min] space. (While any other function that goes to zero fast enough would be equally good to regularize [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], as, for instance, an exponentially decreasing function, the power law [v.sup.-r.sub.min] does not require the introduction of an arbitrary [v.sub.min] scale in the problem.) In fact, by multiplying and dividing the integrand in (25) by [v.sup.r.sub.min],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (27)

we can define the average of the functions [v.sup.r.sub.min] [[??].sup.i] ([v.sub.min])

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (28)

(i = 0 for the unmodulated and i = 1 for the modulated component, see (13)). Notice that exploiting the definition of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in (26), we can write this relation in terms of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] instead of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (29)

where in the integration by parts the finite term [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] vanishes since by assumption r has been appropriately chosen to regularize the integral of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], that is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Equations (28) or (29) allow translating rate measurements in a detected energy interval [[E'.sub.1], [E'.sub.2]] into averaged values of [V.sup.r.sub.min] [[??].sup.i] ([v.sub.min]) in a finite [v.sub.min interval [[v.sub.min,1], [v.sub.min,2]]. This is now the interval outside which the integral of the new response function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (and not of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is negligible. We choose to use 90% central quantile intervals; that is, we determine vmin 1 and [v.sub.min,2] such that the area under the function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] to the left of [v.sub.min,1] is 5% of the total area, and the area to the right of [v.sub.min,2] is also 5% of the total area. In practice, the larger the value of r, the smaller is the width of the [[v.sub.min,1], [v.sub.min,2]] interval, designated by the horizontal error bar of the crosses in the ([v.sub.min,[eta]]) plane. However, r cannot be chosen arbitrarily large, because large values of r give a large weight to the low velocity tail of the [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] function, and this tail depends on the low energy tail of the resolution function [G.sub.T]([E.sub.R], E') in (20), which is never well known. Therefore too large values of r make the procedure very sensitive to the way in which the tails of the [G.sub.T]([E.sub.R], E') function are modeled. This is shown more explicitly in Section 5 (see also Figure 3), where we apply this procedure to a particular WIMP-nuclei interaction. In the figures, the horizontal placement of the vertical bar in the crosses corresponds to the maximum of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The extension of the vertical bar, unless otherwise indicated, shows the 1a interval around the central value of the measured rate.

The upper limit on the unmodulated part of [v.sup.r.sub.min][??] is simply [v.sup.r.sub.min][[??].sup.lim]([v.sub.min]), where [[??].sup.lim]([v.sub.min]) is computed as described in Section 3 by using a downward step-function [[??].sub.0] [theta]([v.sub.0] - [v.sub.min]) for [[??].sup.0] ([v.sub.0]) to determine the maximum value of the step [[??].sub.0]. Given the definition of the response function R in the general case in terms of H, (26), the downward step-function choice for [[??].sub.0] yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (30)

From this equation we find the maximum value of [[??].sub.0] at [v.sub.0] allowed by the experimental upper limit on the unmodulated rate [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (31)

In the figures, rather than drawing the new averages [bar.[v.sup.r.sub.min][[??].sup.i]] [c.sup.2] and the limits [v.sup.r.sub.min][[??].sup.lim]([v.sub.min])[c.sup.2], we may draw [v.sup.r.sub.min] [bar.[v.sup.r.sub.min][[??].sup.i]] [c.sup.2] and [[??].sup.lim] ([v.sub.min])[c.sup.2] (in units of [day.sup.-1]), so that a comparison can be easily made with the results obtained for the SI interaction shown in the previous section.

5. Application to Magnetic Dipole Dark Matter

We apply here the generalized halo-independent method to a Dirac fermion DM candidate that interacts only through an anomalous magnetic dipole moment [[lambda].sub.[chi]] (see e.g., [42-61])

[L.sub.int] = [[[lambda].sub.[chi]]/2] [bar.[chi]] [[sigma].sub.[mu]v] [chi] [F.sup.[mu]v]. (32)

The differential cross section for scattering of a magnetic dipole dark matter (MDM) with a target nucleus is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (33)

Here [alpha] = [e.sup.2]/4[pi] is the electromagnetic fine structure constant, [m.sub.p] is the proton mass, [S.sub.T] is the spin of the target nucleus, and [[??].sub.T] is the magnetic moment of the target nucleus in units of the nuclear magneton e/(2[m.sub.p]) = 0.16 [(GeV/[c.sup.2]).sup.-1]. The first term corresponds to the dipole-nuclear charge coupling, and [F.sub.E,T]([E.sub.R]) is the corresponding nuclear charge form factor. We take it to be the Helm form factor [37] normalized to [F.sub.E,T](0) = 1. The second term, which we call "magnetic", corresponds to the coupling of the DM magnetic dipole to the magnetic field of the nucleus, and the corresponding nuclear form factor is the nuclear magnetic form factor [F.sub.M,T]([E.sub.R]). This magnetic form factor is not identical to the spin form factor that accompanies spin-dependent interactions, in that the magnetic form factor includes the magnetic currents due to the orbital motion of the nucleons in addition to the intrinsic nucleon magnetic moments (spins). For the light WIMPs considered here, the magnetic term is negligible for all the target nuclei we consider except Na. This term is more important for lighter nuclei, such as Na and Si, but only [sup.23]Na has a nonnegligible magnetic dipole moment, [[??].sub.Na] = 2.218. The magnetic form factor for this nuclide is taken from [62] as explained in [32].

The spin-independent part of the differential cross section has two terms, with different v dependences. Therefore, had we computed the rate with the method used to get to (8), we would have obtained two terms in the rate each containing a different function of [v.sub.min] multiplied by detector dependent coefficients. It would have been impossible in this way to translate a rate measurement or bound into only one of the two [v.sub.min] functions. In such a situation, the approach presented in Section 3 cannot be applied and one needs to resort to the generalized method of Section 4. The function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] has in this case a [v.sup.2] dependence for large values of v, with [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] scaling as [v.sub.min]. More precisely we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (34)

where we defined [[sigma].sub.ref] = [alpha][[lambda].sup.2.sub.[chi]]. As a consequence,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (35)

The denominator of (28) is therefore

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (36)

where r can be any number larger than 2. To obtain this result, we first integrated the [theta] term in (35) with respect to [v.sub.min] and then used (2) to change integration variable from [E.sub.R] to [v.sub.min] again.

In Figure 3, we illustrate the effect of various choices of r on the response function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] for MDM for several energy bins and experiments: the first energy bin of DAMA/LIBRA [2], 2 to 2.5keVee, the 7 to 9keVCDMS-II used for the Si data [7], and the first, 0.43 to 1.11keVee, and last, 2.49 to 3.18keVee, of CoGeNT [3, 4]. We also include [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] from (10) for the standard SI interaction (gray dashed line) for a comparison. The normalization of each curve is arbitrary. For r = 0, the MDM response function is divergent and goes like v at large velocities, given the [v.sup.2] behavior of ([v.sup.2]d[[sigma].sub.T]/d[E.sub.R]) (see discussion after (26)). The divergent behavior is much more pronounced in the low energy bins. The choice r = 3 is already enough to regularize the divergent behavior but still yields too large [v.sub.min] intervals. For growing values of r, the peak of the response function shifts towards low velocities (mostly in the low energy bins), due to the [v.sup.-r.sub.min] factor. This peak, when far from the [v.sub.min] interval where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is nonnegligible, is unreliable as it is due to the low energy tail of the detector energy resolution function [G.sub.T]([E.sub.R], E'), which determines the low velocity tail of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (see (20)) and is never well known. We found the optimum r value by trial and error and for MDM we find that r = 10 is an adequate choice (see Figure 3) to get a localized response function in [v.sub.min] space without relying on how the low energy tail of the energy resolution function is modeled. The choice of r is dictated by the lowest energy bins, where the function [v.sup.-r.sub.m] is largest. Higher energy bins are less sensitive to the choice of r.

Let us remark that this way of comparing data is not an inherent part of the halo-independent method but is only due to the choice of finding averages over measured energy bins to translate putative measurements of a DM signal. So far a better way of presenting the data has not been found, and more work is necessary to make progress in this respect.

Figures 4,5, and 6 show the measurements of and bounds on [v.sup.10.sub.min] [[??].sup.0]([v.sub.min]) and [v.sup.10.sub.min] [[??].sup.1] ([v.sub.min]) for a WIMP with magnetic dipole interactions and mass m = 6 GeV/[c.sup.2], 9 GeV/[c.sup.2], and 15 GeV/[c.sup.2], respectively. These masses are motivated by previous studies on MDM as a potential explanation for the putative DM signal found by DAMA, CoGeNT, and CRESST-II (see e.g., [57]). The averages (indicated by the crosses) and upper bounds are multiplied by [v.sup.-10.sub.min] so that the vertical axis has the usual [??][c.sup.2] units of [day.sup.-1] and the bounds show [[??].sup.lim] ([v.sub.min]) (as usual for SI interactions). Figures 4, 5, and 6 include the DAMA modulation signal (green crosses), CoGeNT modulated (blue crosses), and unmodulated signal (plus an unknown flat background, darkred horizontal lines), CDMS-II-Si unmodulated rate signal (red crosses and limit line), CDMS-II-Ge unmodulated rate limit (light blue line) and modulation bound (dark grey horizontal line with downward arrow), XENON100 limit (purple line), and XENON10 S2-only limit without [Q.sub.y] suppression below 1.4 keVnr (orange line). The crosses represent the averages [bar.[v.sup.10.sub.min][[??].sup.i]] (i = 0 for the unmodulated and i = 1 for the modulated parts of the velocity integral) over the [v.sub.min] intervals indicated by the horizontal bar of each cross multiplied by [v.sup.-10.sub.min]. The lines represent upper limits on [[??].sup.0]([v.sub.min]). The CDMS-II-Ge modulation limit is instead an upper limit on [bar.[v.sup.10.sub.min][[??].sup.i]] multiplied by [v.sup.-10.sub.min].

The measurements and limits in Figures 4, 5, and 6 for MDM move to larger [v.sub.min] values as the WIMP mass increases, as expected from the relation (2) between [v.sub.min] and the recoil energy. As shown in Figure 4, for a WIMP of mass m = 6 GeV/[c.sup.2] the three CDMS-II-Si points are largely below the XENON10 and XENON100 upper limits, but they move progressively above them as m increases to 9 GeV/[c.sup.2], see Figure 5, and are almost entirely excluded by them for m = 15 GeV/[c.sup.2] in Figure 6. The addition of the recent LUX bound, which is not shown here, poses however a greater threat to the DM interpretation of the CDMS-II-Si excess even for very low DM masses. Moreover, the three CDMS-II-Si points overlap with or are below the CoGeNT and DAMA measurements of [[??].sup.1], and therefore these appear incompatible with the interpretation of the CDMS-II-Si data as a measurement of [[??].sup.0], as one usually expects that if [[??].sup.0] [much greater than] [[??].sup.1]. For all three WIMP masses shown in the figures, the DAMA and CoGeNT modulation measurements seem compatible with each other. However, the upper limits on the unmodulated part of the rate imposed by XeNoN10 and XENON100 (plus CDMSlite and LUX, which are not shown here) reject the MDM interpretation of the DAMA and CoGeNT modulation signal, except for the lowest energy bins.

6. Conclusions

In order to interpret and compare data from different experiments, a model is often (if not always) needed. Concerning direct DM detection experiments, one needs to assume a model of particle interactions between DM particles and nuclei in the detectors, as well as a model for the dark halo, most notably the DM velocity distribution. For light WIMPs with ~10 GeV/[c.sup.2] mass, as those pointed by direct detection experiments with positive signals in the assumption of SI interactions and the SHM, the details of the high velocity region of the DM velocity distribution are crucial in comparing positive and negative results. For this reason, a framework to analyze the direct detection data independently on the properties of the DM halo is an important tool to address the compatibility of the different experiments.

In this work we have reviewed the halo-independent method to compare data from direct DM detection experiments, as introduced and developed in [27-34]. We followed closely the treatment in [31-33], which present the most updated analyses of DM direct detection experiment data in the context of the halo-independent method. We applied the halo-independent analysis to SI interactions with both isospin-conserving and isospin-violating couplings. In both cases the situation seems to be of disagreement between most of the experiments with positive signals (DAMA, CoGeNT, and CRESST-II) and those with negative results (most notably LUX, XENON, and CDMS-II). The three CDMS-II-Si events seem however compatible with all the limits for DM with isospin-violating couplings. DAMA and CoGeNT modulation data sets seem to agree with one another, but they appear to be incompatible with the CDMS-II-Si events when these are interpreted as measurements of the unmodulated rate.

We have also shown the results of the halo-independent analysis in the assumption of WIMPs interacting with nuclei only through a magnetic dipole moment. In this case the scattering cross section has a more complicated dependence on the DM velocity, thus requiring a generalized version of the method as it was originally devised. The conclusions for this DM candidate are similar as above, with only the lowest energy data points of DAMA, CoGeNT, and CDMS-II-Si lying below the exclusion bounds. The situation is somewhat better for very light WIMPs (m ~ 6 GeV/[c.sup.2]), especially for the CDMS-II-Si events.

The halo-independent analysis is a promising framework to compare different direct detection experiments without making assumptions on the DM halo. This feature is highly desirable given the crucial role played by the DM velocity distribution in the galaxy in determining the total scattering rate at direct detection experiments. This analysis allows to directly compare the recoil spectra measured by different experiments in [v.sub.min] space, together with bounds from null experiments. These spectra indicate the integrated DM velocity distribution [??] favored by the experiments as a function of [v.sub.min] (see (23)).

At present this framework presents some drawbacks, which could be addressed and improved in future work. For instance, the relation between the [??] function that one wants to fit and the observed rates is an integral equation, (25). So far it has been assumed that [??]([v.sub.min], t) is approximately constant in any [v.sub.min] interval where the response functions [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are significantly different from zero, so that it could be extracted from the integral in the form of the average in (14) or (28) and (29). However, this is not necessarily a good assumption. Secondly, the degree of agreement or disagreement between two data sets cannot be readily quantified (in a statistical sense) in the current halo-independent analysis. Finally, the method provides no information on the consistency of modulated and unmodulated signals, even when these are measured by the same experiment as for CoGeNT. By making some (mild) assumptions, more stringent limits on the modulated part [[??].sup.1] can be derived from the limits on the unmodulated part of the rate [28, 38-40]. However, with no additional assumptions the only way one can bound the modulated rate with the unmodulated rate is by imposing the most general inequality [[??].sup.0] > [[??].sup.1].

http://dx.doi.org/10.1155/2014/604914

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The author thanks L. Marcus and L. Prosperi for useful discussions. Partial support from Department of Energy under Award no. DE-SC0009937 is also acknowledged.

References

[1] P. A. R. Ade, N. Aghanim, C. Armitage-Caplan et al., "Planck 2013 results--XVI. Cosmological parameters," http://arxiv.org/abs/1303.5076.

[2] R. Bernabei, P. Belli, F. Cappella et al., "New results from DAMA/LIBRA," European Physical Journal C, vol. 67, no. 1, pp. 39-49, 2010.

[3] C. E. Aalseth, P. S. Barbeau, N. S. Bowden et al., "Results from a search for light-mass dark matter with a p-type point contact germanium detector," Physical Review Letters, vol. 106, no. 13, Article ID 131301, 4 pages, 2011.

[4] C. E. Aalseth, P. S. Barbeau, J. Colaresi et al., "Search for an annual modulation in a p-type point contact germanium dark matter detector," Physical Review Letters, vol. 107, no. 14, Article ID 141301, 5 pages, 2011.

[5] C. E. Aalseth, P. S. Barbeau, J. Colaresi et al., "CoGeNT: a search for low-mass dark matter using p-type point contact germanium detectors," Physical Review D, vol. 88, no. 1, Article ID 012002, 20 pages, 2013.

[6] G. Angloher, M. Bauer, I. Bavykina et al., "Results from 730 kg days of the CRESST-II dark matter search," The European Physical Journal C, vol. 72, article 1971, 2012.

[7] R. Agnese, Z. Ahmed, A. J. Anderson et al., "Dark matter search results using the silicon detectors of CDMS II," Physical Review Letters, vol. 111, no. 25, Article ID 251301, 6 pages, 2013.

[8] A. K. Drukier, K. Freese, and D. N. Spergel, "Detecting cold dark-matter candidates," Physical Review D, vol. 33, no. 12, pp. 3495-3508, 1986.

[9] K. Freese, M. Lisanti, and C. Savage, "Annual modulation of dark matter: a review," Reviews of Modern Physics, vol. 85, no. 4, pp. 1561-1581, 2013.

[10] D. S. Akerib, H. M. Araujo, X. Bai et al., "First results from the LUX dark matter experiment at the Sanford Underground Research Facility," Physical Review Letters, vol. 112, no. 9, Article ID 091303, 7 pages, 2014.

[11] J. Angle, E. Aprile, F. Arneodo et al., "A search for light dark matter in XENON10 data," Physical Review Letters, vol. 107, no. 5, Article ID 051301, 5 pages, 2011.

[12] E. Aprile, M. Alfonsi, K. Arisaka et al., "Darkmatter results from 225 live days of XENON100 data," Physical Review Letters, vol. 109, no. 18, Article ID 181301, 6 pages, 2012.

[13] Z. Ahmed, D. S. Akerib, S. Arrenberg et al., "Results from a low-energy analysis of the CDMS II germanium data," Physical Review Letters, vol. 106, no. 13, Article ID 131302, 5 pages, 2011.

[14] R. Agnese, A. J. Anderson, M. Asai et al., "CDMSlite: a search for low-mass WIMPs using voltage-assisted calorimetric ionization detection in the SuperCDMS experiment," Physical Review Letters, vol. 112, no. 4, Article ID 041302, 6 pages, 2014.

[15] M. Felizardo, T. A. Girard, T. Morlat et al., "Final analysis and results of the phase II SIMPLE dark matter search," Physical Review Letters, vol. 108, no. 20, Article ID 201302, 5 pages, 2012.

[16] S. Archambault, E. Behnke, P. Bhattacharjee et al., "Constraints on low-mass WIMP interactions on [sup.19]F from PICASSO," Physics Letters B, vol. 711, no. 2, pp. 153-161, 2012.

[17] E. Behnke, J. Behnke, S. J. Brice et al., "First dark matter search results from a 4-kg CF3 I bubble chamber operated in a deep underground site," Physical Review D, vol. 86, no. 5, Article ID 052001, 9 pages, 2012.

[18] Z. Ahmed, D. S. Akerib, A. J. Anderson et al., "Search for annual modulation in low-energy CDMS-II data," http://arxiv.org/abs/1203.1309.

[19] C. Savage, G. Gelmini, P. Gondolo, and K. Freese, "Compatibility of DAMA/LIBRA dark matter detection with other searches," Journal of Cosmology and Astroparticle Physics, vol. 2009, no. 4, article 010, 2009.

[20] Y.-Y. Mao, L. E. Strigari, R. H. Wechsler, H.-Y. Wu, and O. Hahn, "Halo-to-Halo similarity and scatter in the velocity distribution of darkmatter," The Astrophysical Journal, vol. 764, no. 1, article 35, 2013.

[21] A. M. Green, "Effect of halo modeling on WIMP exclusion limits," Physical Review D, vol. 66, no. 8, Article ID 083003, 9 pages, 2002.

[22] C. McCabe, "The astrophysical uncertainties of dark matter direct detection experiments," Physical Review D, vol. 82, no. 2, Article ID 023530, 13 pages, 2010.

[23] A. M. Green, "Dependence of direct detection signals on the WIMP velocity distribution," Journal of Cosmology and Astroparticle Physics, vol. 1010, article 034, 2010.

[24] A. M. Green, "Astrophysical uncertainties on direct detection experiments," Modern Physics Letters A, vol. 27, no. 3, Article ID 1230004, 20 pages, 2012.

[25] M. Fairbairn, T. Douce, and J. Swift, "Quantifying astrophysical uncertainties on dark matter direct detection results," Astroparticle Physics, vol. 47, pp. 45-53, 2013.

[26] C. Arina, "Bayesian analysis of multiple direct detection experiments," http://arxiv.org/abs/1310.5718.

[27] P. J. Fox, J. Liu, and N. Weiner, "Integrating out astrophysical uncertainties," Physical Review D, vol. 83, no. 10, Article ID 103514, 12 pages, 2011.

[28] M. T. Frandsen, F. Kahlhoefer, C. McCabe, S. Sarkar, and K. Schmidt-Hoberg, "Resolving astrophysical uncertainties in dark matter direct detection," Journal of Cosmology and Astroparticle Physics, vol. 2012, no. 1, article 024, 2012.

[29] P. Gondolo and G. B. Gelmini, "Halo independent comparison of direct dark matter detection data," Journal of Cosmology and Astroparticle Physics, vol. 2012, no. 12, article 015, 2012.

[30] M. T. Frandsen, F. Kahlhoefer, C. McCabe, S. Sarkar, and K. Schmidt-Hoberg, "The unbearable lightness of being: CDMS versus XENON," Journal of Cosmology and Astroparticle Physics, vol. 2013, no. 7, article 023, 2013.

[31] E. Del Nobile, G. B. Gelmini, P. Gondolo, and J.-H. Huh, "Halo-independent analysis of direct detection data for light WIMPs," Journal of Cosmology and Astroparticle Physics, vol. 2013, no. 10, article 026, 2013.

[32] E. Del Nobile, G. Gelmini, P. Gondolo, and J.-H. Huh, "Generalized halo independent comparison of direct dark matter detection data," Journal of Cosmology and Astroparticle Physics, vol. 2013, no. 10, article 048, 2013.

[33] E. Del Nobile, G. B. Gelmini, P. Gondolo, and J.-H. Huh, "Update on lightWIMP limits: LUX, lite and light," Journal of Cosmology and Astroparticle Physics, vol. 2014, no. 3, article 014, 2014.

[34] E. Del Nobile, G. B. Gelmini, P. Gondolo, and J.-H. Huh, "Direct detection of light anapole and magnetic dipole DM," http://arxiv.org/abs/1401.4508.

[35] A. Kurylov and M. Kamionkowski, "Generalized analysis of the direct weakly interacting massive particle searches," Physical Review D, vol. 69, no. 6, Article ID 063503, 12 pages, 2004.

[36] J. L. Feng, J. Kumar, D. Marfatia, and D. Sanford, "Isospin-violating dark matter," Physics Letters B, vol. 703, no. 2, pp. 124-127, 2011.

[37] R. H. Helm, "Inelastic and elastic scattering of 187-Mev electrons from selected even-even nuclei," Physical Review, vol. 104, no. 5, pp. 1466-1475, 1956.

[38] J. Herrero-Garcia, T. Schwetz, and J. Zupan, "On the annual modulation signal in dark matter direct detection," Journal of Cosmology andAstroparticlePhysics, vol. 2012, no. 3, article 005, 2012.

[39] J. Herrero-Garcia, T. Schwetz, and J. Zupan, "Astrophysics independent bounds on the annual modulation of dark matter signals," Physical Review Letters, vol. 109, no. 14, Article ID 141301, 5 pages, 2012.

[40] N. Bozorgnia, J. Herrero-Garcia, T. Schwetz, and J. Zupan, "Halo-independent methods for inelastic dark matter scattering," Journal of Cosmology and Astroparticle Physics, vol. 2013, no. 7, article 049, 2013.

[41] S. Yellin, "Finding an upper limit in the presence of an unknown background," Physical Review D, vol. 66, no. 3, Article ID 032005, 7 pages, 2002.

[42] M. Pospelov and T. ter Veldhuis, "Direct and indirect limits on the electro-magnetic form factors of WIMPs," Physics Letters B, vol. 480, no. 1-2, pp. 181-186, 2000.

[43] H. An, S. L. Chen, R. N. Mohapatra, S. Nussinov, and Y. Zhang, "Energy dependence of direct detection cross section for asymmetric mirror dark matter," Physical Review D, vol. 82, no. 2, Article ID 023533, 8 pages, 2010.

[44] K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, and M. Kamionkowski, "Dark-matter electric and magnetic dipole moments," Physical Review D, vol. 70, no. 8, Article ID 083501, 15 pages, 2004.

[45] K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, and M. Kamionkowski, "Dark-matter electric and magnetic dipole moments," Physical Review D, vol. 73, no. 8, Article ID 089903, 2 pages, 2006.

[46] V Barger, W.-Y. Keung, and D. Marfatia, "Electromagnetic properties of dark matter: dipole moments and charge form factor," Physics Letters B, vol. 696, no. 1-2, pp. 74-78, 2011.

[47] S. Chang, N. Weiner, and I. Yavin, "Magnetic inelastic dark matter," Physical Review D, vol. 82, no. 12, Article ID 125011, 8 pages, 2010.

[48] W. S. Cho, J.-H. Huh, I.-W. Kim, J. E. Kim, and B. Kyae, "Constraining WIMP magnetic moment from CDMS II experiment," Physics Letters B, vol. 687, no. 1, pp. 6-10, 2010.

[49] W. S. Cho, J.-H. Huh, I.-W. Kim, J. E. Kim, and B. Kyae, "Constraining WIMP magnetic moment from CDMS II experiment," Physics Letters B, vol. 694, no. 4-5, pp. 496-497, 2011.

[50] J. H. Heo, "Minimal Dirac fermionic dark matter with nonzero magnetic dipole moment," Physics Letters B, vol. 693, no. 3, pp. 255-258, 2010.

[51] S. Gardner, "Shedding light on dark matter: a Faraday rotation experiment to limit a dark magnetic moment," Physical Review D, vol. 79, no. 5, Article ID 055007, 14 pages, 2009.

[52] E. Masso, S. Mohanty, and S. Rao, "Dipolar dark matter," Physical Review D, vol. 80, no. 3, Article ID 036009, 8 pages, 2009.

[53] T. Banks, J. F. Fortin, and S. Thomas, "Direct detection of dark matter electromagnetic dipole moments," http://arxiv.org/abs/1007.5515.

[54] J. F Fortin and T. M. P. Tait, "Collider constraints on dipole-interacting dark matter," Physical Review D, vol. 85, no. 6, Article ID 063506, 6 pages, 2012.

[55] K. Kumar, A. Menon, and T. M. P. Tait, "Magnetic fluffy dark matter," Journal of High Energy Physics, vol. 2012, no. 2, article 131, 2012.

[56] V Barger, W.-Y. Keung, D. Marfatia, and P.-Y. Tseng, "Dipole moment dark matter at the LHC," Physics Letters B, vol. 717, no. 1-3, pp. 219-223, 2012.

[57] E. Del Nobile, C. Kouvaris, P. Panci, F. Sannino, and J. Virkajarvi, "Light magnetic dark matter in direct detection searches," Journal of Cosmology and Astroparticle Physics, vol. 2012, no. 8, article 010, 2012.

[58] J. M. Cline, Z. Liu, and W. Xue, "Millicharged atomic dark matter," Physical Review D, vol. 85, no. 10, Article ID 101302, 6 pages, 2012.

[59] N. Weiner and I. Yavin, "How dark are majorana WIMPs? Signals from magnetic inelastic dark matter and Rayleigh dark matter," Physical Review D, vol. 86, no. 7, Article ID 075021, 16 pages, 2012.

[60] S. Tulin, H. B. Yu, and K. M. Zurek, "Three exceptions for thermal dark matter with enhanced annihilation to [gamma][gamma]," Physical Review D, vol. 87, no. 3, Article ID 036011, 12 pages, 2013.

[61] J. M. Cline, G. D. Moore, and A. R. Frey, "Composite magnetic dark matter and the 130 GeV line," Physical Review D, vol. 86, no. 11, Article ID 115013, 9 pages, 2012.

[62] T W. Donnelly and I. Sick, "Elastic magnetic electron scattering from nuclei," Reviews of Modern Physics, vol. 56, no. 3, pp. 461-566, 1984.

Eugenio Del Nobile

Department of Physics and Astronomy, UCLA, 475 Portola Plaza, Los Angeles, CA

90095, USA

Correspondence should be addressed to Eugenio Del Nobile; delnobile@physics.ucla.edu

Received 26 December 2013; Accepted 3 March 2014; Published 1 April 2014

Academic Editor: Anselmo Meregaglia

Printer friendly Cite/link Email Feedback | |

Author: | Del Nobile, Eugenio |
---|---|

Publication: | Advances in High Energy Physics |

Date: | Jan 1, 2014 |

Words: | 9000 |

Previous Article: | Quantum-spacetime scenarios and soft spectral lags of the remarkable GRB130427A. |

Next Article: | The annual modulation signature for dark matter: DAMA/LIBRA-phase1 results and perspectives. |

Topics: |