A low-threshold diamond cryogenic detector for sub-GeV Dark Matter searches

In this work we report the realization of the first low-threshold cryogenic detector that uses diamond as absorber for astroparticle physics applications. We tested two 0.175$\,$g CVD diamond samples, each instrumented with a W-TES. The sensors showed transitions at about 25 mK. We present the performance of the diamond detectors and we highlight the best performing one, where we obtained an energy threshold as low as 16.8 eV. This promising result lays the foundation for the use of diamond for different fields of applications where low threshold and excellent energy resolution are required, as i.e. light dark matter searches and BSM physics with coherent elastic neutrino nucleus scattering.


Introduction
It is undeniable that the identification of dark matter (DM) is one of the most urgent open questions in physics today. Several observational evidences based on gravitational effects (i.e. anisotropies in the cosmic microwave background, large-scale structure distributions, galaxy cluster velocity dispersions, gravitational lensing effect) support the assumption that 26% of the mass-energy density in the Universe is in the form of non-barionic and cold DM. For recent reviews see, e.g., [1,2]. However, despite the world-wide efforts and Corresponding author: canonica@mpp.mpg.de the tremendous experimental progresses, the nature of DM is still unresolved. A standard approach in the field of direct DM searches is to look for the scattering of hypothetical DM particles with the atomic nuclei of target materials, contained in low-background experiments typically located in underground laboratories. One of the crucial parameters in direct DM experiments is the energy threshold, since it drives the sensitivity to the detection of low-energy nuclear recoils induced by DM interactions. While the elastic scattering of DM particles with masses above a few GeVs can be easily accessed by experiments using large TPCs filled with noble-liquids (Argon [3,4] and Xenon [5][6][7]), light DM particles with masses in the sub-GeV range are not accessible to these experiments due to the kinematics of the process. Simply, light DM particles, that interact with Ar or Xe nuclei, induce nuclear recoils that are below the experimental threshold. Recently, the Migdal effect [8,9] has been proposed to extend the reach of these experiments to sub-GeV DM masses [10,11], but the existence of this effect has still to be experimentally demonstrated. An experimental technique complementary to TPCs is the one of cryogenic detectors. They enable the reach of low-energy thresholds by means of detecting of the phonons created following the particle interaction with the crystalline detector nuclei. As of today, Si [12], Ge [13] and CaWO 4 [14] crystals are the targets most widely used in state-of-the art experiments for direct DM searches. For the read-out of the tiny temperature variations induced by particle interactions, the crystal targets must be instrumented with highly sensitive temperature sensors. Among the available thermometers, Transition Edge Sensors using W arXiv:2203.11999v1 [physics.ins-det] 22 Mar 2022 films [15] and Neutron Transmutation Doped thermistors [16] have the properties that fulfill the requirements of a sub-GeV DM experiment, having obtained energy thresholds of O(10 eV). Thanks to their performance, cryogenic detector are currently leading the sensitivity of spin-independent DM-nucleons interactions for DM masses below few GeV [12][13][14] . However, to further extend the reach of cryogenic experiments and explore new region of parameters of lighter DM particles, experimental thresholds as low as few eV are required. Diamond targets have the potential to reach such a low energy threshold thanks to their superior thermal properties [17] and only recently it has been proposed as detection medium for sub-GeV DM searches [18]. Despite the outstanding thermal properties, diamond has not been operated as a cryogenic absorber for astroparticle physics applications until recently. This was mainly due to the limited availability of cheap, high-purity and large-size single crystal (SC) diamonds on the market. Our group has already reported the first operation of a diamond cryogenic absorber equipped with a W-TES, but due to the poor performance of the setup, it was not possible to obtain any information on the threshold achieved by the prototype [19]. Here we report the results on the cryogenic performance of a new prototype realized with a synthetic SC diamond equipped with a W-TES. We achieved a detector energy threshold of 16.8 eV. This result opens a wealth of opportunities for using high-purity lab-grown diamonds in field of light DM particle searches, but also for other sectors where low energy thresholds are needed (i.e. coherent elastic neutrino nucleus scattering). The low-threshold achieved will allow to study new properties of DM with the data presented in this work. We leave this topic for a future publication that is currently in preparation.

Diamond as cryogenic target for DM detection
Cryogenic calorimeters have been developed to explore different sectors of particle physics and astrophysics, from direct DM searches [12][13][14], to neutrinoless double beta decay [20][21][22] and neutrino mass measurements [23,24]. One of their characteristics is that they can be realized with many different materials, with the only requirement to have crystalline structure with low enough heat capacity at the operating temperature O(10 mK).
The Debye temperature of diamond, that can be used as a figure of merit of the thermal properties of the material [25], is among the highest available in nature (2220 K) and it is significantly higher than other established target materials used in DM direct detection experiments (e.g. 645 K for Si, 371 K for Ge and 250 K for CaWO 4 ). This parameter ensures high phonon density of states, that is key for achieving the ultimate sensitivity in the detection of small energy deposits. Moreover, diamond is made of light nuclei (Z=6) and this ensures a good match in the kinematics of the interaction of the light DM particle. A complete review of the properties of diamond as cryogenic sub-GeV DM detector can be found in [18]. Thanks to the developments in the synthetic grown diamond technologies of the last decades, the barrier of the availability of high-quality crystals has been overcome. We can nowadays easily find on the market SC crystals grown using the Chemical Vapour Deposition (CVD) technology [26]. In this technique the crystal is grown starting from a source gas (typically methane) that is deposited over a proper seed substrate. The samples used in this work have been realized using the heteroepitaxial technique on a foreign multilayer substrate made of Ir/YSZ/Si [27]. This, in contrast to the homoepitaxial technique, does not require the use of a diamond seed to start the growth process. However, it requires only a lattice-matched substrate, enabling the growth of large volume diamonds more easily.

Experimental setup
In this work, the crystals used as absorbers are made by two identical SC diamond samples of 0.175 g each and (2 × 5 × 5) mm 3 , grown at AuDiaTec [28]. To detect the phonons produced in the absorber following a particle interaction, both crystals are instrumented with a W-TES. The TES design is very similar to the ones used for the CRESST-III experiment [29]. They are made of a thin strip of W with two large Al pads partially overlapping the W layer. These Al pads have two different features: they serve as phonon collectors [15] and as ohmic contacts. They are connected via a pair of 25 µm Al bond wires through which the bias current is injected. The W film is also connected by a long and thin strip of Au to a thicker Au bond pad on which a 25 µm Au wire is bonded. This connection serves as weak thermal link between the sensor and the heat bath at ∼10 mK. On the same surface, but separated from the TES, we also evaporated a heater. The heater is made of a thin strip of Au with two Al pads deposited on top. These pads are also bonded with a pair of 25 µm Al bond wires through which a tunable current can be injected to maintain the TES at the desired operating temperature. The heater is also used to inject artificial pulses to monitor the detector response over time and to refine the energy calibration during data analysis. In Figure 1 (top) we show one of the two diamond crystals after the TES fabrication. Each of the two diamond crystals is housed in a Cu structure, held in position by a pair of bronze clamps (see Figure 1 (bottom)). Four sapphire balls sitting on the Cu holder are used as spacers and provide a pointlike contact between the crystal and the holder itself. At a distance of about 0.5 cm from the crystal, a 55 Fe source (activity ∼0.3 Bq) is installed for calibration purposes. The samples were cooled down in a dilution refrigerator at the Max Planck institute for Physics in Munich, Germany, in a surface building, using a Kelvinox-400HA dilution refrigerator from Oxford Instruments (see [30] and references therein for details of the cryogenic infrastructure). Electrical connections to the TES were made by superconducting NbTi twisted pairs wires. The TES readout circuit uses a reference resistor at mK and a commercial SQUID system sensor as front end amplifier, combined with a CRESST-like detector control system [31]. The two detectors were readout using an Applied Physics System model 581 DC SQUID.

Measurement and data processing
During the cryogenic measurement, the diamond samples were measured simultaneously in a stable working Fig. 2 Spectrum of events coming from the iron source (black) and events coming from the artificial pulses sent by a heater (red) in data set 1. point for 58.4 hours. We collected statistics for a total exposure of 4.3×10 −4 kg·d. The TESs showed a superconducting transition at about 25 mK. In the following, the data sets acquired with the diamond samples #1 and #2 are referred to as data set 1 and 2 respectively. Given the small size of the samples and the limited dynamic range of the TESs, a sub-set of measurements was carried out in under-performing conditions. The goal was to extend the sensor dynamic range such that events from the 55 Fe at 5.9 keV were clearly visible and could be used for calibration purposes. The artificial pulses sent through the heater during this measurement were used to calculate the conversion factor between the deposited energy and the injected heater voltage in the entire energy range from threshold up to the iron energy region. Thanks to this calibration of the heater response, we were able to reconstruct the energy scale also in different operational conditions, namely during the full data set acquired with the optimized settings. As an example, events from the 55 Fe source and from the artificial pulses injected through the heater are shown in Fig. 2 for detector 1. After the data set for calibration was acquired, the detectors were operated at their optimal operating point, focusing on achieving a low energy threshold. Data were acquired in a continuously recorded stream and processed with a software trigger after data taking [30]. Triggered events were filtered with an optimum filter [32] algorithm which improves the reconstruction of the amplitude of a pulse with a given pulse shape in the presence of noise. The noise conditions of each detector were inferred from empty traces which were acquired with the continuous recording process. For the representation of the pulse shape, we created a standard event (SEV) by averaging a selection of clean particle events. Following the directions described in [33], the SEV could be used to reproduce the shape of the ac-quired pulses at different energies. With the aim of estimating the noise level of the baselines of each detector, we superimposed the SEV with a fixed amplitude on a set of empty traces. The distribution of this amplitude evaluated with the filter was then fitted with a gaussian function. The sigma of this distribution was used to determine the resolution of the baseline. The obtained values in Volt and their respective calibrated values in eV can be found in Tab. 1 for the two diamond samples. Following the approach described in [34], the threshold was computed accepting a fixed number of noise triggers for exposure units which is negligible compared to the particle event rate. For this analysis in both modules the number of accepted noise trigger events is around 1% of the total number of triggered events. The final threshold values obtained with this method can be inferred from Fig. 3, where the noise trigger curves are shown as a function of the trigger threshold. Taking into account the calibration factors, these threshold values correspond to 19.7 eV and 16.8 eV for data set 1 and data set 2, respectively. Fig. 3 Threshold evaluation of the two sets of data. The threshold is estimated with a function that gives the number of noise triggers for a given threshold (see [34] for details). The thresholds that have been calculated are 14.9 mV for data set 1 and 18.7 mV for data set 2. Taking into account the calibration functions, they correspond to 19.7 eV and 16.8 eV respectively. This result, although obtained with a prototype, shows for the first time the excellent performance of diamond when operated as cryogenic detector. We believe that future work on the optimization of the W-TES design to better match the properties of thermal signals in diamonds will allow for further reduction of the energy threshold, making the goal of reaching O(eV) values within reach.  Table 1 Values for threshold and baseline resolution of the two data sets. Each value is shown in mV and in eV.

Conclusions
This work shows the advanced performance achieved with two detectors made of diamonds operated as cryogenic detectors. The results presented in this work represent a major milestone for future sub-GeV DM investigations, since the low-energy thresholds achieved with these prototypes already compete with the most advanced detector currently operating in the field of light DM searches. A dedicated future work on the exploration of new region of the parameter space for DM interactions with regular matter obtained with the data presented in this work is currently under preparation.