Results on MeV-scale dark matter from a gram-scale cryogenic calorimeter operated above ground

Models for light dark matter particles with masses below 1 GeV/c$^2$ are a natural and well-motivated alternative to so-far unobserved weakly interacting massive particles. Gram-scale cryogenic calorimeters provide the required detector performance to detect these particles and extend the direct dark matter search program of CRESST. A prototype 0.5 g sapphire detector developed for the $\nu$-cleus experiment has achieved an energy threshold of $E_{th}=(19.7\pm 0.9)$ eV, which is one order of magnitude lower than previous results and independent of the type of particle interaction. The result presented here is obtained in a setup above ground without significant shielding against ambient and cosmogenic radiation. Although operated in a high-background environment, the detector probes a new range of light-mass dark matter particles previously not accessible by direct searches. We report the first limit on the spin-independent dark matter particle-nucleon cross section for masses between 140 MeV/c$^2$ and 500 MeV/c$^2$.


Introduction
The origin of dark matter (DM) is still unknown, despite plenty of compelling observational evidence [1]. The most popular model for DM in recent years suggests the existence of weakly interacting massive particles (WIMPs) which are produced in the early Universe [2]. Such particles are in particular appealing since the correct relic density of DM can be produced with a thermally-averaged annihilation cross section σ v which is amazingly close to the weak scale. This natural explanation is referred to as the "WIMP miracle" and has been the main driver for DM physics. Minimal supersymmetric (SUSY) models provide an attractive candidate for WIMPs: the neutralino often assumed as the lightest supersymmetric particle with a mass of O(1 GeV/c 2 ) to O(1 TeV/c 2 ) [3]. Lee and Weinberg constrain the mass of WIMPs to m χ ≥ 2 GeV [4]. Lighter particles would have a too small annihilation cross section σ v ∝ m 2 χ m 4 Z (m Z : mass of the Z-boson) which would lead to a too large relic abundance of DM and, thus, an overclosure of the Universe. Therefore, most experimental DM searches in the past decades focused on GeV-scale particles.
Despite enormous experimental effort, no clear and unambiguous signature for DM has been found by direct, indirect or accelerator searches so far.
Due to this lack of evidence, models which provide alternatives to WIMPs should be probed experimentally. In the past years, rising interest is drawn to light DM (lDM) which extends the allowed mass range of DM particles to O(keV). Asymmetric DM [5][6][7], scalar DM particles [8,9] and hidden sector DM [10] are examples for promising and naturally motivated lDM models. These theories are compatible with the observed relic density (Lee-Weinberg argumentation), the constraints from existing DM searches, in particular indirect searches which are sensitive to lDM, and limits from particle physics experiments.
While searches for axions probe particle masses of keVscale and below, the mass range between 1 MeV and 1000 MeV is only weakly constrained by lDM-electron scattering [11] or via Bremsstrahlung emission in nuclear recoils [12,13].
We report here results from the ν-cleus 0.5 g prototype detector [14] which enables for the first time to directly probe nuclear recoils induced by DM particles with masses below 500 MeV/c 2 .

Detector technology
Cryogenic detectors are sensitive to the temperature rise induced by a particle interaction. CRESST-type detectors use single crystals (e.g. CaWO 4 or Al 2 O 3 ) equipped with transition-edge-sensors (TES) made of thin tungsten films. The TES is usually directly coupled to the heat sink by a gold bond wire. Due to the operation at very low temperatures of ∼ 10 mK and the resulting electron-phonon decoupling in the thermometer film, such a cryogenic calorimeter is mainly sensitive to non-thermal phonons [15]. Very low energy thresholds, a high dynamic range (of up to 10 5 ) and a particle-type independent detector response are the main advantages of these devices. State-of-the-art CaWO 4 detectors of 300 g operated in CRESST-II reach thresholds of ∼ 300 eV [16].
The energy threshold of cryogenic detectors depends strongly on the mass M of the crystal which can be quantified by a scaling law discussed in [14]. This describes nicely the performance of existing CaWO 4 and Al 2 O 3 detectors and can be used to predict the sensitivity for detectors of various sizes, materials and geometries [17]. In case of cubic crystals the energy threshold scales as E th ∝ M 2/3 suggesting that thresholds of 10 eV and below are in reach for detectors with masses of O(1 g) [14].

The prototype calorimeter
The first prototype calorimeter made of Al 2 O 3 has a size of 5 × 5 × 5 mm 3 and a mass of 0.49 g. One side of the optically polished crystal is equipped with a specifically developed TES.
The design is similar to the one used for CRESST light detectors, but adjusted to the size of the crystal [17]. The sensor consists of a 200 nm thick tungsten film with an area of 0.0061 mm 2 and an aluminium phonon collector of thickness 1 μm and an area of 0.15 mm 2 . It is weakly coupled to the heat bath via a gold stripe and a bond wire (heat conductance: ∼ 10 pW/K at 10 mK). The dimension of the TES and the strength of the heat link are chosen such that the life time of non-thermal phonons in the detector τ n is much shorter than the decay time of the pulses. The former depends on size, geometry and material of the crystal and the TES, the latter is determined by the strength of the heat link to the thermal bath. A detector designed in this way operates in the calorimetric mode [15], i.e. the phonons impinging on the film are integrated over the duration of the non-thermal phonon signal. 1 A separate gold film sputtered on the crystal acts as an ohmic heater. A current applied to it heats the detector to the desired position in the superconducting phase transition. Artificial pulses of high energy (called control pulses [18]) which saturate the detector response are injected into the heater (here: every 10 s). Those allow to control and stabilize the operating point. In addition, artificial pulses of discrete, lower energies (called test pulses) are injected periodically to continuously monitor the detector response (here: every 60 s).
A fit of the pulse model [15] to the pulses of the prototype gram-scale detector confirms its calorimetric operation. The non-thermal phonon life-time in the crystal is τ n = (0.30 ± 0.01) ms and the dominant decay time is found to be τ dec = (3.64 ± 0.01) ms (for the non-thermal pulse component). Figure 1 (left) shows how the prototype detector is arranged in the experimental setup. The cubic crystal is installed on a copper plate and placed on three sapphire spheres (∅ 1 mm), pressed from the top with another sphere attached to a flexible bronze clamp. The sensor is placed on the top surface where it is electrically and thermally contacted to copper bond pads on the clamp via Al and Au wire bonds, respectively. At a distance of about 2 cm from the crystal, a 55 Fe calibration source (activity ∼ 0.2 Bq) is installed.

Experimental setup
The prototype is operated in a cryostat of the Max-Planck-Institute for Physics in Munich in a surface building. A commercial dilution refrigerator from Oxford Instruments (Kelvinox400) is used with a standard helium dewar. The experimental volume is surrounded by a 1 mm copper thermal screen, but no special care on radiopurity is taken in this setup. No dedicated shielding against ambient radiation is installed within or outside the cryostat. The ceiling of the room is made of ∼ 30 cm concrete, therefore a moderate reduction of the hadronic component of cosmogenic radiation is expected. In the present setup, no muon veto or other anti-coincidence detectors are installed.
The TES is biased with a constant current (1.0 µA) and read out with DC-SQUIDs (Jessy from Supracon). The output of the SQUID electronics is fed to a lownoise voltage amplifier (SR560 from Stanford Research Electronics). For data-taking, two systems have been used in parallel: (1) the standard CRESST data acquisition (DAQ), see e.g. [18], with a hardware trigger unit and 16 bit transient digitizer, and (2) a continuous data-taking system based on a 16 bit digitizer from National Instruments (NI USB-6218 BNC). The latter was used to record a continuous stream of the entire measurement with a sampling frequency of 20 kHz. Offline, a software trigger is applied to the data.

Measurement
The first measurement of the prototype gram-scale calorimeter in the MPI cryostat was primarily intended to de-monstrate an energy threshold in the 10 eV regime for the ν-cleus experiment [17]. Since a threshold of ∼ 20 eV was reached (see [14] and below), one order of magnitude lower than previous CRESST results [16], this first data can be used to search for new physics: a previously un-explored region of parameter space for DM-particle nucleus scattering can be explored (see Sect. 4). In this section we review the energy calibration and threshold analysis which is presented in detail in [14], and focus on dedicated aspects of the DM analysis.

Energy calibration
The measurement with the 0.49 g Al 2 O 3 prototype calorimeter had a total run-time of 5.31 h corresponding to an exposure of 0.11 g-days. The 55 Fe calibration source was placed in the cryogenic setup during the entire run. The dominant K α line (E lit = 5.895 keV) of 55 Mn is used to calibrate the energy spectrum. The linearisation of the pulse response is established by a truncated template fit which uses the pulse shape infomation only in the linear region of the detector response [14]. The method is well established and successfully used in previous CRESST ana-lyses [16,19]. The fit reproduces the K β line of 55 Mn at an energy of E obs = (6.485 ± 0.017) keV, well in agreement with the literature value of 6.490 keV [14]. The robustness of the linearisation at lower energies is studied and a moderate systematic error of 1.1 % on the energy calibration is derived [14]. Figure 2 (main frame) shows the final energy spectrum with the dominant X-ray lines of 55 Mn. Above 7 keV, a constant background rate of ∼ 1.2 × 10 5 counts/(kg keV day) is observed, while from ∼ 1 keV towards lower energies the spectrum is significantly rising to a rate of about 10 8 counts/[kg keV day] (see discussion below).

Software trigger
For the evaluation of the pulse amplitudes in the linear region of the detector response (up to 600 eV), the optimum filter is used [20,21]. This method improves the reconstruction of a known signal in the presence of noise with a measured power spectrum and typically gives significantly better results than frequency [Hz] . The histogram shows the distribution of the filter output (maximum) of randomly sampled noise traces (right axis). The energy threshold E th is chosen such that a negligible amount of noise triggers is accepted. Data from [14] the template fit. The filter transfer function is calculated by the ratio of the power spectra of the template pulse and the noise. This function is used to weight the spectral components of the data sample (in the frequency domain) according to the respective signal-to-noise-ratio. Figure 3 shows the measured power spectra (left axis) and the resulting transfer function (right axis) of the optimum filter used for the analysis of the prototype calorimeter. After being applied in the frequency space, the filter output is transformed back to time domain and normalized to match the maximum of the original pulse. In case of this measurement, a baseline noise of σ b = (3.74± 0.21) eV [14] is reached with the optimum filter, a factor of 1.7 smaller compared to that obtained with the template fit. The trigger threshold is derived by a systematic study [14]. The optimum filter is applied to a set of noise traces, which are randomly recorded by the (standard) DAQ, and is evaluated at every position within a single trace. The maximum filter output for each trace is stored. Figure 4 shows a histogram of the maximum filter output for ∼ 400 baseline traces. The threshold of the software trigger has to be cho-sen such that noise triggers are negligible. In the case of this measurement, the trigger threshold is set to 13.0 mV [14]. To derive the trigger efficiency, artificial template pulses of discrete energies are added to the randomly chosen noise traces. After applying the optimum filter, the fraction of events above and below the trigger threshold are calculated for each simulated energy. The result is shown as dots in Fig. 4. These data points are well fit by an error function p trig which gives an energy threshold of E th = (19.7 ± 0.1(stat.)) eV. This value corresponds to about 5.3 σ b of the baseline noise. As expected, the width of the error function σ th = (3.83 ± 0.15) eV is in agreement with the baseline noise σ b evaluated at E = 0.
To fully exploit this improvement, a software-trigger algorithm for the continuously acquired data was developed. The trigger threshold was set accordingly to 19.7 eV. The entire data stream is processed with the optimum filter by an algorithm which avoids distortions in the finite Fourier transform (for details see [21]). An event is triggered whenever the filter output exceeds the chosen software threshold. Using the same filter for amplitude evaluation and triggering guarantees the consistency of the trigger definition and the energy calibration. The pulse data from hardware (standard DAQ) and software (continuous DAQ) trigger agree down to the hardware threshold of ∼ 40 eV. The pulse height evaluation of the optimum filter and the truncated standard event was studied on an event-by-event basis. In linear region of the pulse response, a maximum deviation of 2.8 % is observed [14] which adds to the systematic error of the energy calibration. The final value of the threshold is E th = (19.7±0.9) eV including statistical and systematic errors. For the following DM analysis we use exclusively the data recorded with the continuous DAQ.

Stability and data selection
Since the data was not blinded for the calibration measurement [14], we choose the most conservative approach for this DM analysis: we apply no data quality cuts which might introduce an energy-dependent efficiency or might bias the result.
The only cut applied is the so-called stability cut which rejects periods where the detector is not in its correct operating point. The cut is based on the pulse height of control pulses (saturated pulses) which are continuously injected into the ohmic heater every 10 s. This pulse height should remain constant if the operating point is stable. To be conservative all events within ± 5 min around every outlying control pulse are not used for the analysis. This removes 2.05 h from the total measuring time of 5.31 h. Figure 5 shows the stability of various detector parameters over the duration of the measurement. In the top frame the total event rate after the stability cut (grey crosses, left axis) The gaps in the data correspond to periods removed by the stability cut. The bottom frame shows the energy resolution as obtained from empty baseline samples (black error bars, left axis) and the energy of periodically injected test pulses corresponding to an energy of ∼ 120 eV (red dots, right axis) together with their average (1σ error bars). Over the entire measuring time the variance of the baseline is in agreement with σ b = (3.74 ± 0.21) eV. The test pulse response also remains constant over time demonstrating again the stability of the detector response, except for the period at 5.0 h when the base temperature of cryostat has changed. This period is however removed by the stability cut (compare to upper frame of Fig. 5).

Results and discussion
We define the linear region of the pulse response which extends from the energy threshold of 19.7-600 eV as our region-of-interest for DM search (ROI). The final spectrum is depicted in Fig. 2 (inset) with a binning of 5 eV. All 511 events in the ROI are considered conservatively as candidate events for DM particle-nucleus scattering, although it is assumed that these pulses originate from backgrounds. Since we do not apply data-quality cuts, 2 the remaining events (higher energetic particle pulses, test pulses, control pulses and possible artifacts) are unspecified in the analysis. Again conservatively, we remove the sampling time (614.28 ms per pulse) of the 5779 events which lie outside the ROI from the live-time. This results in a net live-time of 2.27 h which corresponds to a net exposure of 0.046 g-days.
The constant background level of ∼ 1.2 × 10 5 counts/(kg keV day) is not unexpected due to lack of shielding against ambient radiation and due to the operation of the detector above ground. In addition, an X-ray calibration source was present during the entire measurement. The steep rise of the event rate towards threshold can be caused by source-related Auger electrons. Also, backgrounds such as beta/gamma or alpha decays on surfaces surrounding the detector could contribute to the exponentially increasing spectrum (see [14] for a more detailed discussion). The presence of exponentially rising, unknown backgrounds prohibits an unambiguous discovery of DM, however promising techniques are on their way to significantly reduce and understand remaining backgrounds [17].
Since the threshold of the prototype detector is lowered by more than one order of magnitude compared to previous macroscopic devices, the mass range of DM particles below 500 MeV/c 2 can be probed for the first time despite a relatively high background level. An upper limit on the elastic spin-independent DM particle-nucleon cross section is derived. For the limit calculation we use the Yellin optimal interval method [22] which was developed to derive DM limits in the presence of backgrounds with unknown energy spectra. The sensitivity is based on a comparison of the observed spectrum [see Fig. 2 (inset)] and the expected recoil spectrum from DM particles of a certain mass m χ in Al 2 O 3 .
For the calculation of the recoil spectrum we use the standard astrophysical parameters for the DM halo of the Milky Way: a Maxwellian velocity distribution, an asymptotic velocity of 220 km/s and a galactic escape velocity of 544 km/s. The local DM density at Earth position is assumed to be 0.3 GeV/cm 3 . The Helm form factor [23] is used in the calculation of the scattering cross section. Both elements of Al 2 O 3 are considered as targets. To take into account the finite energy resolution of the calorimeter the calculated spectrum is convolved with the Gaussian resolution (σ b = (3.74 ± 0.21) eV at threshold).
The measurement presented here extends the reach of direct DM search experiments for lDM masses to below 500 MeV/c 2 . The achieved threshold of E th = (19.7 ± 0.9) eV allows to probe elastic scattering of DM particles down to masses of 140 MeV/c 2 . Figure 6 shows the upper Fig. 6 Parameter space of DM mass vs. spin-independent DM particlenucleon cross-section. The result of this work (solid red) explores a new mass region between 140 and 500 MeV/c 2 and extends the reach of DM direct searches. In comparison selected experimental results are shown: CRESST-II Phase 2 (dashed red) [16], CDMS-lite (full green) [24], EDELWEISS (dashed green) [25], XENON100 low mass analyis (dotted blue) [26], PandaX (full blue) [27], LUX (dashed blue) [28] and DAMIC (black) [29]. The shaded grey area indicates the neutrino-floor calculated for Al 2 O 3 limit with 90 % confidence level achieved in this work (red line). Below 500 MeV/c 2 we explore a new region of parameter space via the direct detection of nuclear recoils and improve existing limits from the detection of bremstrahlung emission which accompanies nuclear recoils [12,13]. For masses below ∼ 1 GeV/c 2 O-recoils dominate the recoil spectrum, while for higher masses Al-recoils provide the highest sensitivity. The transition is visible as a weak inflection in the exclusion curve.
In comparison, selected results of direct DM search experiments on the elastic spin-independent DM particle-nucleon cross-section are shown.

Outlook
Gram-scale cryogenic calorimeters have demonstrated their high potential for DM search. An Al 2 O 3 prototype has explored a new range of DM particles below 500 MeV/c 2 . This first result can be improved in two directions: (1) detector performance. A recently developed scaling law [14] predicts energy thresholds in the 1-10 eV regime for cubic detectors of 1 g. Modifications of the TES sensor are foreseen in order to reach this goal [14]. (2) Background level. The shielding of the setup used for the measurement can be improved significantly. We expect a reduction of the background rate by two orders of magnitude when installing a low-radioactivity Pb shielding around the experimental volume. A 4π active cryogenic veto around the calorimeter would significantly reduce surface-related backgrounds (see MC results in [17]) which are suspected to cause backgrounds at lowest energies.
The next generation of gram-scale calorimeters which will be developed for the ν-cleus experiment [17] may extend the reach of direct search experiments to DM particle masses of O(10 MeV/c 2 ) and probe new models for DM. Due to hints that backgrounds at lowest energies are not dominated by cosmogenic radiation, large improvements in sensitivity may be achieved in above-ground or shallow laboratories.