ournal of C osmology and A stroparticle hysics J The unbearable lightness of being: CDMS versus XENON

. The CDMS-II collaboration has reported 3 events in a Si detector, which are consistent with being nuclear recoils due to scattering of Galactic dark matter particles with a mass of ∼ 8 . 6 GeV and a cross-section on neutrons of ∼ 2 × 10 − 41 cm 2 . While a previ-ous result from the XENON10 experiment has supposedly ruled out such particles as dark matter, we ﬁnd by reanalysing the XENON10 data that this is not the case. Some tension remains however with the upper limit placed by the XENON100 experiment, independently of astrophysical uncertainties concerning the Galactic dark matter distribution. We explore possible ways of ameliorating this tension by altering the properties of dark matter interactions. Nevertheless, even with standard couplings, light dark matter is consistent with both CDMS and XENON10/100.


Introduction
Identifying dark matter (DM) is a key challenge for modern astro-particle physics [1]. Its presence has been inferred solely through gravitational interactions but it is theoretically well-motivated to assume that DM consists of new particles arising in physics beyond the Standard Model (SM). It would then be natural for the strength of DM interactions to be less even than that of the weak interactions, consistent with the inferred collisionless behaviour of DM on cosmological scales. The relic abundance of weakly interacting massive particles (WIMPs) which decoupled from thermal equilibrium in the early universe is also naturally of order the observed DM abundance if they have a weak scale mass of O(100) GeV -the so-called 'WIMP miracle'. Another attractive possibility is that the DM particles inherited the same matter-antimatter asymmetry as baryons are observed to have -their mass would then have to be of O(5) GeV. In either case, the most direct way to identify DM is with well-shielded underground experiments which can detect the energy deposited by a recoiling nucleus when a passing DM particle from the Galactic halo scatters on it.
Many such searches have been conducted for over three decades but the overwhelming majority have focussed on DM particles of mass 100 m χ /GeV 1000, motivated simultaneously by the WIMP miracle and by the existence of a natural particle candidate viz. the neutralino in supersymmetric extensions of the SM. Less attention has been paid to searching for lighter particles, both in the absence of as clear a theoretical motivation for these, as well as the increased difficulty in distinguishing nuclear recoils from background at low energies. However in recent years results from several experiments have suggested m χ 10 GeV. This mass range is favoured by an excess of events over the expected background observed by CoGeNT [2] and CRESST-II [3] as well as the long-standing annual modulation signal seen by DAMA [4]. However, constraints set by other leading experiments, notably XENON10 [5], XENON100 [6] and CDMS-Ge [7,8], are in strong tension with these results and have cast doubt on the DM interpretation of these signals (see e.g. [9][10][11][12][13][14]).
The direct detection of 'light' DM is particularly challenging, because it is hard to reliably determine the detector response to a potential DM signal at such low recoil energies [15,16]. Moreover, the interpretation of the data is sensitive to the assumed properties (in particular, velocity distribution) of the DM halo [17][18][19][20]. To minimise these uncertainties, an ideal experiment to probe the low-mass region would have target nuclei with low mass JCAP07(2013)023 number and a low energy threshold. Si detectors satisfy both criteria, making them an obvious choice for this purpose. With increasing interest amongst theorists in developing particle physics models for light DM (and its co-genesis with baryons), the experimental community is becoming increasingly aware of the importance of looking in this mass region.
The CDMS-II collaboration has recently presented two analyses of data taken with Si detectors. While no candidate events were found in the small data set from the first 'Five-Tower Run' of CDMS-II [21], the longer exposure of 140.2 kg-days during 2007-08 [22] revealed 3 events in the DM search region, compared to a background expectation of 0.62 events. Considering only the number of events the signal is compatible with background at the ∼ 2σ level, but taking into account their distribution, a background+DM interpretation of the data is preferred over the background-only hypothesis with a probability of 99.8%. The highest likelihood is found for a DM particle with mass 8.6 GeV and cross-section 1.9 × 10 −41 cm 2 , which is significantly below the DAMA and CRESST-II regions, but still in tension with bounds from XENON10/100.
We investigate under what conditions the tension between CDMS and XENON10/100 can be ameliorated or even resolved. In particular, we find the bound from the 'S2-only analysis' of XENON10 to actually be weaker by nearly a factor of 10 than the published value [5]. However, the tension with XENON100 [6] remains and cannot easily be evaded by appealing to experimental or astrophysical uncertainties. We therefore consider various particle physics modifications of the interactions between DM and SM quarks. While a momentum-or velocity-dependence in the cross-section does not improve the situation, we find that isospin-dependent couplings as well as exothermic DM-scattering can bring XENON10/100 and CDMS into better agreement. We do not attempt, however, to find a common interpretation of the CDMS results and the signals found by DAMA, CoGeNT and CRESST-II.
In section 2, we introduce the standard framework for interpreting DM signals and review both the CDMS best-fit region and the bounds from XENON10/100. We pay particular attention to the XENON10 bound and its dependence on the form of the ionisation yield Q y , which is essential to understand the energy scale. In section 3, we briefly review v min -space, which enables comparison of results from different direct detection experiments without having to assume a specific DM velocity distribution. We show that CDMS and XENON10/100 probe much the same region of v min -space, hence the tension between them is independent both of astrophysical uncertainties concerning the DM halo and any momentumand velocity-dependence of the cross-section. We illustrate this by giving some examples for possible modifications. Then in section 4 we consider possible changes in the underlying particle physics, which can suppress the sensitivity of Xe targets relative to Si, and thus bring CDMS and XENON10/100 into better agreement. We present our conclusions in section 5.

General framework and discussion of experimental results
We first review the standard method of analysing DM direct detection experiments and then apply this to the silicon data from CDMS-II (CDMS-Si in short) and the XENON10/100 experiments. The event rate in a given detector depends on the flux of Galactic DM particles passing through the detector and the DM-nucleon scattering cross-section. In the laboratory frame the differential event rate with respect to recoil energy E R is given by

JCAP07(2013)023
where ρ is the local halo DM density, m χ and m N are the DM and target nucleus mass respectively, f (v) is the local DM velocity distribution evaluated in the Galactic rest frame, v = |v| and v E (t) is the velocity of the Earth relative to the Galactic rest frame [23,24]. We have separated the differential scattering cross-section into the piece for scattering off a point like nucleus, dσ/dE R , and the nuclear form factor F (E R ), which accounts for the finite size of the nucleus. For elastic scattering the minimum velocity required for a DM particle to transfer the energy E R to a given nucleus N is: where µ is the reduced mass of the DM-nucleus system. For a detector with energy resolution ∆E R , the number of events N (E 1 , E 2 ) expected in an energy range [E 1 , E 2 ] is: where Ex is the exposure of the detector, (E R ) is the detector acceptance and is the detector response function [25]. We now focus on whether CDMS-Si is compatible with XENON10 and XENON100. We consider only spin-independent scattering in this work since only a very small fraction of naturally occurring Si contains an isotope with nuclear spin. The usual parameterisation for a spin-independent cross-section is where µ nχ is the reduced DM-nucleon mass, C T (A, Z) ≡ (f p /f n Z + (A − Z)), A and Z are the mass and charge number of the target nucleus and f n,p denote the effective DM coupling to neutrons and protons, respectively. Experiments typically quote results in terms of σ n , the DM-neutron cross-section at zero momentum transfer. The above parameterisation of the cross-section holds e.g. for DM-nucleus scattering mediated by a heavy CP -even scalar boson or by a heavy vector mediator. However, other interactions can give rise to a cross-section that cannot be parameterised in this way. We consider two such examples in section 3. Unless otherwise stated, we adopt for the astrophysical parameters the Standard Halo Model (SHM), i.e. a Maxwell-Boltzmann distribution for f (v) with v 0 = 220 km s −1 and v esc = 544 km s −1 , and take ρ = 0.3 GeV cm −3 . We use the Helm form factor [26] and take f n /f p = 1. The results of our calculations for this choice of parameters are shown in figure 1.
For XENON100 we take the data from an exposure of 20.9 kg-years [6], combined with the recent measurements of L eff [27]. We calculate the energy resolution under the assumption that it is dominated by Poisson fluctuations in the number of photoelectrons and obtain exclusion limits using the 'maximum gap method' [28]. As noted already, CDMS-Si found 3 events after background rejection in data from an exposure of 140.2 kg-days [22]. For our analysis we assume an energy resolution of ∆E R = 0.3 keV [29] and use the detector  . CDMS-Si confidence region (68% and 90% CL) together with the 90% exclusion curves from XENON10 and XENON100 from our analysis, assuming standard elastic spin-independent scattering and equal coupling to protons and neutrons. Our CDMS-Si confidence region and XENON100 bound agree well with the results from the respective collaborations [6,22], however, our XENON10 bound is significantly weaker than the published one [5]. We consider three choices of the ionisation yield Q y at low energy to illustrate the corresponding variation of the extracted bound.
acceptance from [22]. We take the normalised background distributions from [30] and rescale the individual contributions in such a way that 0.41, 0.13 and 0.08 events are expected from surface events, neutrons and 206 Pb, respectively.
To calculate the allowed parameter region, we employ the extended maximum likelihood method [31]. The contours we show correspond to 68% and 90% confidence level (CL) -these agree well with the ones shown by the CDMS-II collaboration [22]. In excellent agreement with their analysis, we find that the best-fit point is m χ = 8.55 GeV with σ n = 2 × 10 −41 cm 2 , and that the background+DM hypothesis is favoured over the background-only hypothesis with a probability of 99.6%. We find that varying the normalisation of the background within the range suggested in [22] while keeping its shape fixed does not sensibly affect either the significance of the background+DM hypothesis, nor the shape of the best-fit region.
For XENON10, we use the S2-only analysis [5]. The sensitivity of this analysis depends strongly on the behaviour of the ionisation yield Q y at low energies [16]. To calculate the central bound (thick dashed blue line) in figure 1 we adopt the original choice of Q y made in [5], assuming that Q y vanishes for E R < 1.4 keV. We assume that the energy resolution is given by ∆E R = E R / E R Q y (E R ) and use the maximum gap method to set a bound. We have checked that the 'p max method' [28] and the 'binned Poisson method' [25,32] give comparable bounds. In every case, we find bounds that are significantly weaker than the published one [5] -even when we assume exactly the same values of Q y -but in good agreement with other independent analyses [10,16,33,34]. 1 It has been suggested that the JCAP07(2013)023 XENON10 data be analysed by converting recoil spectra to photoelectrons and using Poisson statistics to model the detector response [35]. We find that this procedure gives even weaker bounds at low DM mass. Therefore we assume Gaussian fluctuations in order to enable fair comparison with other results in the literature.
We also consider values of Q y that have recently been determined by the XENON100 collaboration by comparing calibration data to Monte Carlo simulations. We adopt the range of values given in figure 3 of [27], which we have reproduced together with the choice of Q y from [5] in figure 1 (inset). Most of the sensitivity for detecting low-mass DM relies on recoil energies below 3 keV, where Q y cannot be reliably extracted from the data. We therefore extrapolate Q y in two different ways that bracket the range of possibilities: the pessimistic choice (light dashed blue line in figure 1) assumes that Q y follows the lower bound of the range given in [27] and then drops to zero at 2 keV. The optimistic choice (light dotted blue line in figure 1) assumes that Q y follows the upper bound of the range, continuing to rise up to 8 e − /keV below 3 keV, and then drops to 0 at 1 keV. Even for the optimistic case we do not extrapolate Q y below 1 keV, i.e. we follow the approach of [5] and neglect upward fluctuations. Changing Q y also changes the energy scale shown in figure 2 of [5] so we convert the energies of the observed events appropriately when setting the bound. It is clear from figure 1 that although the XENON10 S2-only analysis does constrain the CDMS-Si parameter region, an unambiguous bound cannot be placed at low mass because of uncertainty in the value of the ionisation yield Q y at low energy.
Henceforth we present only one XENON10 bound, adopting the original choice of Q y from [5]. As this is significantly weaker than the published bound [5], CDMS-Si and XENON10/100 are actually consistent at about 90% CL. Nevertheless there is some tension between these results. We will now explore possible ways to ameliorate this tension.

Analysing the experiments in v min -space
We have seen that many parameters need to be specified before a theoretical prediction for the number of scattering events in a direct detection experiment can be compared with the observed number. While a parameter such as the local DM density affects all experiments in the same way, other parameters can change the number of events in one experiment while having no impact on another experiment. A useful technique to gain insight into this involves mapping the experimental result into v min -space [36]. If experiments probe different regions of this space, they will be affected differently by varying parameters such as the local escape velocity v esc ; conversely, if experiments probe the same region of v min -space, then modifying such parameters cannot improve agreement between the experiments. We first apply this technique in the usual way to astrophysical parameters, before applying it also to momentum-dependent interactions. Our discussion and notation closely follow [37].

Varying astrophysical parameters
After substituting the usual parameterisation of the cross-section for spin-independent scattering from eq. (2.5) into eq. (2.1), we see that direct detection experiments do not directly probe the local velocity distribution f (v), but rather the velocity integral g( where v min is defined by eq. (2.2). For our purposes, it will be convenient to absorb the DM mass, cross-section and density into this definition and consider the 'rescaled velocity integral':g (v min ) = ρ σ n m χ g(v min ) . Since g(v min ) is a monotonically decreasing function of v min [38], we can use experimental null results to place upper bounds on g(v min ) for v min =v min using the standard techniques for setting an exclusion limit, but with [36,37] The resulting bounds from XENON10 and XENON100 are shown as the blue (dashed) and purple (dot-dashed) lines in figure 2, respectively. In order to deduce information ong(v min ) from CDMS-Si, we need to infer the differential event rate and use the formula: The differential event rate is found most easily by binning both the observed events and the background expectation. We subtract the background using the Feldman-Cousins technique [39], which allows us to calculate an upper and lower bound on the signal event rate at 1σ CL. The orange data points in figure 2 show the result of this procedure for m χ = 9 GeV  (upper panels) and m χ = 7 GeV and m χ = 11 GeV (lower left and lower right panels respectively). Binning the data introduces a certain arbitrariness so we check the robustness of our results by considering two choices of the bin width: 2 keV and 3 keV for the upper left and right panels of figure 2 respectively. The inferred values forg(v min ) agree well, implying that our conclusions are largely independent of the choice of bin width. In all cases, the highest bin is in significant tension with the XENON100 bound except for the case m χ = 7 GeV, corresponding to the least constrained mass in figure 1.
We observe from figure 2 that all three experiments probe essentially the same region of v min -space. This suggests that it will not be possible to significantly improve the consistency of CDMS-Si and XENON10/100 by varying astrophysical parameters. To explicitly demonstrate that this is so, we consider two variations in astrophysical parameters. In the left panel of figure 3 we keep the usual Maxwell-Boltzmann velocity distribution but choose v 0 = 250 km/s and v esc = 650 km/s, which are at the upper end of the allowed range for these parameters (see e.g. [17] and references within). Although we see that the CDMS-Si region and XENON10/100 bounds move towards lower masses by ∼ 1 GeV, the tension between the experiments remains essentially unchanged. As a more radical modification we consider the effect of a 'debris flow' in the presence of which g(v min ) is given by [40,41] , (3.4) where v flow is the velocity of the flow in the Galactic rest frame and α is the contribution of the debris to the local DM density. In the right panel of figure 3 we have chosen α = 0.2 and v flow = 400 km/s. We see again that although the preferred mass changes, the overall agreement does not.

Non-standard momentum and velocity dependence
We have seen that CDMS-Si and XENON10/100 probe the same region of v min -space so the tension between these experiments is independent of astrophysical uncertainties. We now consider the effect of modifying the particle physics to introduce an additional momentumor velocity-dependence of the cross-section which can change the recoil energy spectrum. Non-trivial momentum dependence can be studied using the simple parameterisation [42] dσ where (dσ/dE R ) 0 is the differential cross-section introduced in eq. (2.5), q = √ 2 m N E R is the momentum transfer and q ref is an arbitrary normalisation. The crucial observation is that the momentum transfer q is related to the minimal velocity v min by v min = q 2µ . (3.6) For light DM, m χ m N , hence µ ≈ m χ independent of m N . Consequently, experiments that probe the same range of v min will also probe the same range of q (see also [43]). We cannot therefore significantly shift the parameter region favoured by CDMS-Si relative to the XENON10/100 bounds by introducing a momentum dependence. This is illustrated in the left panel of figure 4 for q ref = 10 MeV and n = −2, which corresponds to longrange interactions arising from a (nearly) massless mediator. As expected, no significant improvement is found.
In many realistic models, the cross-section depends not only on the momentum transfer but also on the velocity v of the incoming DM particle. This is the case e.g. for dipole interactions [44][45][46] and anapole interactions [47,48]. Since we know already that both CDMS-Si and XENON10/100 probe the same part of the DM velocity distribution, such an additional velocity dependence is not expected to change the picture significantly. We confirm this expectation for the case of anapole interactions, which lead to a differential scattering cross-section [48]: Our results are shown in the right panel of figure 4. As expected, because of the momentumand velocity-suppression, the CDMS-Si favoured parameter region and the XENON10/100 bounds are moved to much larger cross-sections, but their relative position remains unchanged. For the anapole operator we do in fact observe a slight shift of the CDMS-Si region compared to the XENON10/100 bounds. This shift can be traced back to the fact that for anapole interactions DM particles couple to protons only leading to a factor of Z 2 rather than A 2 in the cross-section. We explore the effect of different DM couplings to protons and neutrons in more detail in the next section.

Reducing the tension between CDMS-Si and XENON10/100
As we have seen CDMS-Si and XENON10/100 cannot be brought into better agreement by modifying either the DM velocity distribution or the velocity/momentum dependence of the cross-section. To weaken the constraints from XENON10/100, we need to reduce the enhancement of the cross-section for heavy nuclei. In this section, we discuss two possible modifications of DM interactions that can increase the sensitivity of light targets compared to heavy ones: inelastic DM [49] and isospin-dependent couplings [50][51][52][53].
In the former case, DM-nucleon interactions require the transition between two DM states of slightly different mass. The minimum velocity required for a recoil of energy E R is then: where δ is the mass splitting between incoming and outgoing DM state. If the incoming DM particle is lighter (δ > 0), scattering will be enhanced for heavy targets. only the heavier state is populated initially (δ < 0), lighter targets will be favoured. This second case, referred to as exothermic DM [54,55], thus seems a promising option to fully reconcile CDMS-Si and XENON10/100. We present our results in figure 5. As expected, the XENON10/100 bounds are strengthened for inelastic DM and weakened for exothermic DM. In fact, a relatively small splitting of δ = −50 keV is sufficient to bring CDMS-Si and XENON10/100 into good agreement.
Another possible route to concordance is based on the observation that the strength of the XENON10/100 bounds results partially from the factor A 2 in the cross-section, following from the assumption that DM couples equally to protons and neutrons. If DM couples to protons only (as for anapole interactions), the cross-section will be proportional to Z 2 , thus favouring targets with a larger ratio of protons to neutrons. If the coupling to neutrons is slightly negative, f n /f p < 0, targets with a large fraction of neutrons will suffer more strongly from destructive interference between protons and neutrons. Such negative values of f n /f p can arise e.g. in theories with a light Z that mixes with the SM gauge bosons [56].
To study this possibility, we scan simultaneously over f n /f p , σ n and m χ and calculate the likelihood for each set of parameters, using the maximum likelihood method described in [10]. In particular, we assume that the likelihood for XENON10 and XENON100 for a given set of parameters are given by L Xe = exp(−N max ), where N max is the number of events expected in the largest interval determined using the maximum gap method. Since we use an extended maximum likelihood method to fit CDMS-Si, the minimum value of the likelihood depends on an arbitrary normalisation constant and does not carry physical significance. We cannot therefore perform a goodness-of-fit analysis for our model, but we can infer confidence regions for parameter estimation.
The results of this analysis are shown in the left panel of figure 6. As expected [56], the best-fit point corresponds to f n /f p −0.7, which strongly suppresses the bounds from XENON10/100, as shown in the right panel of figure 6. Note that in this particular case, the strongest constraints on CDMS-Si arise from SIMPLE [57]  run [58] (not shown). For f n /f p = −0.7 these experiments require σ n 10 −39 cm 2 at m χ 10 GeV [37] and therefore do not significantly constrain the CDMS-Si preferred region.
In spite of the preference for f n /f p −0.7, we observe that much larger values of f n /f p still give a good fit to the data. At 1σ confidence level, we find −0.76 < f n /f p < −0.58 and the 2σ confidence region extends up to f n /f p −0.2. To illustrate this point, we show the cases f n /f p = −0.5 and f n /f p = −0.2 in figure 7. We conclude that little fine-tuning is required to suppress the bounds from XENON10/100, in particular we do not require a precise cancellation of proton and neutron contributions and are therefore not sensitive to potential differences in the form factors as discussed in [59]. We would also emphasise that even though f n = f p is strongly disfavoured compared to f n /f p = −0.7 when performing a parameter estimation, this does not imply that such a value cannot give a good fit to the data. In fact, as we have shown above, CDMS-Si and XENON10/100 are compatible at the 90% CL even for f n = f p .

Conclusions
The report by the CDMS-II collaboration [22] of 3 events consistent with nuclear recoils from scattering of Galactic DM particles with a mass of ∼ 8.6 GeV is both exciting and intriguing. This low-mass region has been attracting increasing theoretical interest as DM particles with such a mass can have a natural connection to baryons by virtue of sharing a primordial asymmetry. However there appears, prima facie, to be a discordance with the upper limits placed by the XENON10/100 experiments.
We have demonstrated that in fact the XENON10/100 experiments do not exclude the entire CDMS parameter region; in particular, the bound from XENON10 is significantly weaker than reported by the collaboration [5]. Moreover, this bound is rather sensitive to the (unknown) behaviour of the ionisation yield at the relevant low recoil energies.
Nevertheless there is tension between the CDMS results and the XENON10/100 bounds and we have shown explicitly that this is independent of astrophysical uncertainties concerning the velocity distribution of halo DM e.g. the presence of a 'debris flow'. Allowing a JCAP07(2013)023 possible momentum or velocity dependence in the DM scattering cross-section will also not alleviate the tension as we have illustrated for both long-range and anapole interactions.
However the tension is reduced if scattering of DM on heavy targets like Xe is suppressed compared to light ones like Si, e.g. because DM couplings are isospin-dependent, leading to interference between scattering rates on protons and neutrons. As another example for enhancing scattering off light targets, we have briefly discussed exothermic DM.
In conclusion, the DM interpretation of the 3 events reported by the CDMS-II collaboration is not entirely excluded by XENON10/100. This mass region is very interesting theoretically and can be accessed with relatively small experiments employing light target nuclei. In view of our ignorance concerning the possible interactions of DM, it is essential that experimental searches employ as many different target nuclei as is possible.
Noted added in proof. The XENON10 Collaboration has issued an erratum [60] acknowledging the error in their published analysis [5] which we have pointed out in the present paper (section 2). The corrected limit is stated to be in good agreement with our result shown in figure 1.