Limits on dark matter effective field theory parameters with CRESST-II

CRESST is a direct dark matter search experiment, aiming for an observation of nuclear recoils induced by the interaction of dark matter particles with cryogenic scintillating calcium tungstate crystals. Instead of confining ourselves to standard spin-independent and spin-dependent searches, we re-analyze data from CRESST-II using a more general effective field theory (EFT) framework. On many of the EFT coupling constants, improved exclusion limits in the low-mass region (< 3–4 GeV/c2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c^2$$\end{document}) are presented.


Introduction
The elusive nature of dark matter remains one of the major unsolved mysteries in modern physics. One leading hypothesis is that dark matter consists of as yet undetected particles with interactions at the weak scale (or below) [1]. If the hypothesis is correct, the microscopic properties of dark matter might be revealed in the coming years using existing detection methods [2]. The direct detection technique will play a key role in this context [3]. It searches for nuclear recoils induced by the non-relativistic scattering of Milky Way dark matter particles in low-background detectors [4,5]. Detectors based on dual-phase time projection chambers have proven to be very effective in the search for dark matter particles heavier than about 10 GeV/c 2 [6]. WIMPs (for weakly interacting massive particles) are the leading dark matter candidate in this mass range [7]. On the other hand, when the dark matter particle mass is below few GeV/c 2 , cryogenic experiments provide the best sensitivity to dark matter-nucleon interactions because of the low energy threshold these detectors can achieve [6]. The experiment CRESST (for Cryogenic Rare Event Search with Superconducting Thermometers) has pioneered the search for sub-GeV/c 2 dark matter, and currently places the most stringent exclusion limits on the spin-independent dark matter-nucleon scattering cross-section for dark matter masses below 1.8 GeV/c 2 [8]. Dark matter in the GeV/c 2 mass range is expected in models where the present cosmological density of dark matter is explained in terms of freeze-in and asymmetric production (for a review, see [9]), or in Strongly Interacting Massive Particles models (see [10]).
From the theoretical side, the null result of present direct detection experiments is usually interpreted within a framework where dark matter either couples to the total nucleon content of the nucleus (spin-independent interaction) or to the nucleon spin content of the nucleus (spin-dependent interaction) [11]. This is a reasonable approach, since spinindependent and spin-dependent interactions are in general expected to give the leading contribution to the cross section for dark matter-nucleus scattering. However, when these standard interactions are forbidden or suppressed, such as in the case of dark matter-nucleon interactions mediated by a pseudo-scalar particle [12] or in the case of anapole interactions [13], the leading contribution to the dark matter-nucleus scattering cross section may have a different nature. The classification and characterisation of non-standard dark matter-nucleus interactions has driven the theoretical research in the field of dark matter direct detection in the past few years. In this context, the nonrelativistic effective theory of dark matter-nucleon interactions has played a key role [14][15][16]. Assuming that dark matter and nucleons are the only relevant degrees of freedom, this theory describes all possible dark matter-nucleon interactions which are compatible with the symmetries characterising the non-relativistic dark matter-nucleus scattering. Within this theoretical framework, data collected by the SuperCDMS and XENON100 experiments have been interpreted in [17,18], respectively. Furthermore, a likelihood analysis of different direct detection experiments and non-relativistic dark matter-nucleon interactions has been performed in [19][20][21][22]. The role of operator interference has extensively been discussed in [23].
The main goal of this work is to set exclusion limits on the coupling constants of the effective theory of dark matternucleon interactions using data collected by the CRESST-II Phase 2 experiment. This analysis generalises previous results found by the CRESST Collaboration focusing on standard interactions, and extends limits on non-standard interactions presented by other groups to an as yet unexplored mass range.
This article is organised as follows. We introduce the nonrelativistic effective theory of dark matter-nucleon interactions in Sect. 2 and the CRESST experiment in Sect. 3. Section 4 is devoted to the methods and data used in our analysis, while a summary of our results and conclusions is presented in Sect. 5.

Effective theory of dark matter direct detection
In this section we briefly review the non-relativistic effective theory of dark matter-nucleon interactions as defined in [14]. The theory is based upon the following considerations: (1) in the non-relativistic limit, the amplitude for dark matter scattering off nucleons N in target nuclei, M χ N , can in general be expanded in powers of |q|/m N 1, where |q| is the momentum transferred in the scattering and m N is the nucleon mass. (2) Each term in this expansion must be invariant under Galilean transformations and Hermitian conjugation, and can be expressed in terms of basic invariants under the above symmetries [14]: S χ , S N , iq, and v ⊥ ≡ v+q/2μ N , where S χ (S N ) is the dark matter (nucleon) spin, and μ N and v are the dark matter-nucleon reduced mass and relative velocity, respectively. These considerations imply that in the one-body approximation, 1 the Hamiltonian for the interactions of dark matter with a nucleus T , H χ T , can be written as follows [14]: where the index j characterises the dark matter-nucleon interaction type and c 0 j (c 1 j ) is the associated isoscalar (isovector) coupling constant. The A nucleons in the target nucleus are labeled by the index i = 1, . . . , A, and 1 i 2×2 (τ i 3 ) is the identity (third Pauli matrix) in the i-th nucleon isospin space. Finally,Ô i j is a non-relativistic operator for interactions of type j between the dark matter particle and the i-th nucleon in the nucleus.
In the Hamiltonian (1), the operatorsÔ i j act on particle coordinates through the momentum transfer and transverse relative velocity operators,q andv ⊥ , respectively. At linear order inv ⊥ , and at second order inq, Eq. (1) includes 16 independent interaction operatorsÔ i j , listed in Table 1. Not all of them appear as leading operators in the non-relativistic limit of simplified models 2 [25][26][27]. In Table 1, the dark matter and nucleon spin operators are denoted byŜ χ and S N , respectively. Following [14], here we do not consider the operatorsÔ i  [25]. For simplicity, from here onwards we will omit the nucleon index i in the definitions. Table 1 Quantum mechanical operators defining the non-relativistic effective theory of dark matter-nucleon interactions [14]. Here we adopt the notation introduced in Sect. 2. Standard spin-independent and spin-dependent interactions correspond to the operatorsÔ 1 andÔ 4 , respectively. The operatorsÔ 17 andÔ 18 can only arise for spin 1 dark matter, and S is a symmetric combination of spin 1 polarisation vectors [25]. Following [14], here we do not consider the operatorsÔ i 2 and O i 16 (see text above Eq. (2) for further details). For simplicity, we omit the nucleon index in the operator definitionŝ In a dark matter direct detection experiment, the differential rate of nuclear recoil events per unit detector mass is given by: where is the minimum dark matter velocity required to deposit an energy E R in the detector, μ T and m T are the dark matter-nucleus reduced mass and target nucleus mass, respectively, and m χ is the dark matter mass. In Eq. (2), ρ χ is the local dark matter density, while f (v) is the dark matter velocity distribution in the detector rest frame. The sum in Eq. (2) runs over all elements in the detector. Each contribution is weighted by the corresponding mass fraction ξ T .
In Eq. (2), the differential cross section for dark matternucleus scattering, dσ T /dE R , depends on the isoscalar and isovector coupling constants, c 0 j and c 1 j , respectively, and on nuclear matrix elements of H χ T . For an explicit expression, see [28]. This very general description of the dark matternucleus scattering captures most of the particle physics scenarios that one can conceive. Important exceptions include models where the dark matter-nucleus scattering is inelastic [29], or scenarios where dark matter-nucleon interactions are mediated by particles with mass comparable or lighter than typical momentum transfers [30]. In addition, Eq. (1) cannot be used to describe effects related to meson exchange in nuclei, e.g. the "pion pole" [24,31]. However, such effects are known to be important only for momentum transfers comparable with the pion mass, and are therefore expected to be negligible in the mass range of interest for the present analysis [32]. Finally, Eq. (1) cannot account for operator mixing effects induced by the running of coupling constants. These can be predicted within ultraviolet complete models [33]. While CRESST data are here interpreted within a non-relativistic effective theory for dark matter-nucleon interactions, the use of relativistic effective theories for dark matter-quark and -gluon interactions has also been discussed in the literature. The relation between the relativistic and the non-relativistic approach is explained in Refs. [27,31,32]. Exclusion limits from CRESST data are separately computed for each of the non-relativistic operators that we consider. The impact of interference effects on direct detection exclusion limits has been discussed in, e.g. [34].
Let us now comment on some of the assumptions made while evaluating Eq. (2). Regarding the local dark matter density, we adopt the standard value of 0.3 GeV/c 2 /cm 3 , although slightly larger values are favoured by astronomical data, e.g., [35]. For the dark matter velocity distribution in the detector rest frame, we assume a Maxwellian velocity distribution with a circular speed of 220 km/s for the local standard of rest and a galactic escape velocity of 544 km/s (i.e. the so-called Standard Halo Model [36]). As far as the detector composition is concerned, here we consider the contribution of Oxygen and Calcium to the scattering cross section, but neglect Tungsten. For Tungsten, the nuclear response functions, or "form factors", associated with (most of) the operators in Table 1 15 contribute to the event rate in the present analysis. In the notation of [14], these operators generate the nuclear responses W τ,τ M , W τ,τ Φ and W τ,τ MΦ . In the zero momentum transfer limit, W τ,τ M measures the nucleon content of the nucleus (and is proportional to the standard spin-independent form factor), whereas W τ,τ Φ measures the nucleon spin orbit coupling content of the nucleus. Finally, the nuclear response W τ,τ MΦ arises from the interference of the nuclear currents underlying W τ,τ M and W τ,τ Φ . Conclusions based on Eq. (2) are affected by uncertainties in astrophysical and nuclear physics inputs. For kinematical reasons, only uncertainties on the nuclear response functions at zero momentum transfer are relevant for light dark matter. Whereas W τ τ M (0) is known exactly, being proportional to the square of the number of nucleons in the nucleus, uncertainties on W τ,τ Φ (0) and W τ,τ MΦ (0) must be assessed through nuclear structure calculations. In the case of Helium, the relative uncertainty on these response functions was found to be of a factor of 3 or so using an ab initio no core shell model approach [37]. Using large-scale nuclear structure calculations, similar results were found for the nuclear response functions of interest in the case of Xenon isotopes [38]. On the other hand, astrophysical uncertainties can play a crucial role in the search for light dark matter, especially those on the dark matter velocity distribution. For example, in the small mass limit, it has been found that exclusion limits can be modified by up to few orders of magnitude by variations in the astrophysical inputs that govern the dark matter and baryon mass profiles in the Milky Way [39]. In order to consistently compare our results with those in [17], we will present our exclusion limits focusing on the Standard Halo Model, and adopting the nuclear response functions for Oxygen and Calcium computed in [40] through a shell model calculation.

The CRESST experiment
CRESST (Cryogenic Rare Event Search with Superconducting Thermometers) is a direct dark matter search experiment. The anticipated dark matter signals are nuclear recoils in a scintillating calcium tungstate (CaWO 4 ) target crystal. The target detectors are operated at a temperature of around 15mK. To shield the experiment from background signals, mainly induced by cosmic rays, the experiment is located at the underground laboratory of the LNGS (Laboratori Nazionali del Gran Sasso) in central Italy. This work uses data from phase 2 of CRESST-II, 3 which started in July 2013 and ended in August 2015. In total 18 detector modules with an overall mass of ∼ 5 kg were operated [41,42].
The module Lise was the detector module with the lowest trigger threshold for nuclear recoils (0.307 keV). Limits on the elastic spin-independent dark-matter-nucleon crosssection from this detector module were published in 2016 (see [42]). For dark matter masses below 2 GeV/c 2 these limits led the field at the time of the publication. 4 This work uses the same data as [42].
The detector module consists of a scintillating CaWO 4 crystal with a mass of 300 g (phonon detector) and an independent light detector. Most of the energy deposited by Fig. 1 All events from the detector module Lise after all cuts. The light yield, the fraction of light to phonon signal, is plotted against the (phonon) energy. Solid lines are 90% boundaries for electron recoils (blue) and nuclear recoils off oxygen (red) and tungsten (green). The red dashed line is the center of the oxygen band. The yellow area is the acceptance region Fig. 2 Energy spectrum of all events in the previously defined acceptance region from the detector module Lise. The range in the y-axis is chosen for reasons of clarity, although a few bins surpass the upper bound recoils leads to a phonon signal, which is thus used for the energy determination. The fraction of the energy that yields a scintillation light signal is called light yield and is used for discrimination between different types of recoils: while the phonon signal is fairly independent of the type of the scattered particle, the light signal depends strongly on the sort. The scattered particle can either be a calcium, tungsten or oxygen nucleus, or an electron. While most of the background events are electron recoils, the anticipated dark matter signals are nuclear recoils. Figure 1 shows the dataset used for this work. For more details on data preparation and cuts see [42]. The light yield is normalized to 1 for electron recoils with an energy of 122 keV. Since only nuclear recoils are considered as a possible dark matter signature in this analysis, the acceptance region is chosen accordingly: To avoid leakage from electron recoils, the upper bound is defined as the center of the oxygen band, which is the highest of the nuclear recoil bands in Fig. 1. The lower bound is chosen to be the 99.5% lower boundary of the lowest band, which is the tungsten band. The calcium band (not shown in Fig. 1) lies between the oxygen and the tungsten band. In terms of energy the acceptance region spans from threshold (307 eV) to 40 keV. The acceptance region, as well as all methods of data preparation and selection, have been defined and fixed before unblinding the data. Figure  2 shows the energy distribution of all events in the acceptance region. The double peak with measured energies of 6.0 keV and 6.6 keV is due to a 55 Fe X-ray source that was installed for the calibration of a detector module close to Lise.
Although this was an unintentional exposure, it doesn't significantly influence the sensitivity of the experiment because of the narrow width of the peaks. The smaller peaks at 2.7 keV and 8.1 keV are also understood. They originate from cosmogenic activation of tungsten and copper fluorescence respectively.

Effective field theory data analysis
The goal of this analysis is to set limits on the coupling constants c i defined in Eq. 1. In principle the theory allows any linear combination of the operators, but we restrict ourselves to limits on single operators individually. The calculation of the expected spectra for each operator and for each isotope was executed based on a Matlab code 5 released by the CDMS Collaboration for a similar analysis [17]. For each operator, the spectra are then added up according to the abundance of each isotope in calcium tungstate (CaWO 4 ).
Only the most abundant isotopes of oxygen ( 16 O) and Calcium ( 40 Ca) are taken into account. The nuclear form factors from Sect. 2 have already been calculated for these two isotopes using shell model computations [40], but not for the much heavier tungsten isotopes. Therefore, tungsten is left out of this analysis.
In order to calculate limits, Yellin's optimum interval method is used [43,44]. It doesn't require a background model and the implementation of different spectrum shapes is unproblematic. Also, without taking into account all isotopes in the spectrum calculation (which leads to lower expected spectra), limit calculation is still possible and valid, but yields conservative results. 5 The code is cross-checked with [28].

Results and conclusions
The limits for the Wilson coefficients c 0 1 , c 0 3 , c 0 5 , c 0 8 , c 0 11 , c 0 12 and c 0 15 are shown in Fig. 3. The results are compared to limits from the CDMS [17] and the Xenon100 experiment [18].
The isotopes that are taken into account for this analysis, 16 O and 40 Ca, both contain the same number of protons and neutrons. Consequently, we are not sensitive to isovector operators in this analysis. Also, the overall spin of both of the two isotopes is zero. For this reason, the strongest limits, compared to the CDMS experiment, are obtained for the coefficients c 0 i for operators that contain no dependence on the nucleon spin, namely the operators O 1 , O 5 , O 8 and O 11 . However, limits on the coeficcients c 0 i for the operators O 3 , O 12 and O 15 are also provided. Figure 3 shows that this present analysis sets leading limits on all of the above-mentioned Wilson coefficients for dark matter masses below 3-4 GeV/c 2 . Similar to the standard spin-independent analysis [42], the low energy threshold of the CRESST detectors leads to particularly high sensitivity for low particle masses. The upgrade to CRESST-III, with a threshold that is even lower and improvements on the detector design, holds considerable potential to further improve these limits, especially for low dark matter particle masses. It should also be noted that in the effective field theory, where differences in nuclear properties of the target nuclei play a decisive role, the comparison of multiple different complementary experiments becomes even more important than in the standard spin-independent analysis.