The COSINUS project: perspectives of a NaI scintillating calorimeter for dark matter search

The R&D project COSINUS (Cryogenic Observatory for SIgnatures seen in Next-generation Underground Searches) aims to develop a cryogenic scintillating calorimeter using an undoped NaI-crystal as target for direct dark matter search. Dark matter particles interacting with the detector material generate both a phonon signal and scintillation light. While the phonon signal provides a precise determination of the deposited energy, the simultaneously measured scintillation light allows for particle identification on an event-by-event basis, a powerful tool to study material-dependent interactions, and to suppress backgrounds. Using the same target material as the DAMA/LIBRA collaboration, the COSINUS technique may offer a unique possibility to investigate and contribute information to the presently controversial situation in the dark matter sector. We report on the dedicated design planned for the NaI proof-of-principle detector and the objectives of using this detection technique in the light of direct dark matter detection.


Introduction
In the era of precision cosmology we know that dark matter [1] is five times more prevalent than baryonic matter in the Universe, and experimental evidence of dark matter so far solely relies on gravitational interaction.
Among a long list of hypothetical new particles WIMPs (Weakly Interacting Massive Particles) provide compelling a e-mail: florian.reindl@mpp.mpg.de b e-mail: karoline.schaeffner@lngs.infn.it arguments as this class of candidates, with a mass in the (GeV-TeV)/c 2 region and weak-scale interactions, may provide for a relic density that matches observation. Participation of such particles in weak-scale interactions would also allow for a direct detection in earth-bound detectors via the process of elastic scattering off atomic nuclei [2].
Since the expected recoil energies are in the sub-keV to keV regime, depending on the mass of the dark matter particle and the target material, the challenge of such direct searches is to combine an ultra-low background with a highly sensitive detection apparatus.
At present, the field of direct dark matter search is very active, with numerous experiments all over the world aiming for further increase in size and sensitivity within the next years.
Direct searches are designed to register signals in the detector induced by the interacting particles. The involved detection channels induced by such scattering processes include light, charge and/or phonon signals. Experiments are either based on a single or on a dual channel readout. Which of the three channels are considered depends on the respective target and the applied detection technique. From the first channel the energy of the interacting particle is determined. In most cases the second channel serves to identify the nature of the interacting particle, thereby providing a powerful tool to discriminate dark matter signals from background.
Dark matter searches may be further divided with respect to the approach they use for detection: a numerous class of experiments is optimized to look for nuclear recoil events in a certain energy window where dark matter particles can induce signals. The number of observed recoil events is con-sequently combined with model assumptions, concerning e.g. the density of dark matter, in order to derive properties of dark matter particles, as e.g. their mass and scattering crosssection. Somewhat different therefrom, dark matter may also be identified via an annual modulation signal caused by the seasonal variation of the Earth's velocity with respect to the dark matter halo [3].
At present the situation in the sector of direct dark matter detection is controversial and a definite discovery is still absent. Several experiments [4][5][6] detect events above the known background level which allow to infer properties of a dark matter particle. In particular, the DAMA/LIBRA collaboration observes in more than 13 annual cycles a statistically robust modulation signal using a total of up to 250 kg of radiopure sodium iodide (NaI(Tl)) crystals operated as room temperature scintillating detectors in the Laboratori Nazionali del Gran Sasso (LNGS), a deep underground site in central Italy.
At the same time, the DAMA/LIBRA signal is not consistent, in the standard elastic scattering scenario, with the null results of other direct detection experiments as listed in legend and caption of Fig. 3. However, in order to compare the results of these experiments with the DAMA/LIBRA modulation signal one has to take into account astrophysical assumptions as well as assumptions on the physics of the scattering process of the dark matter particle with the respective target materials. A pure modulation signal instead can be considered model-independent and, thus, has the potential to give a unique and robust signature of new particles. Nonetheless, future investigations are of pivotal importance to understand the origin/composition of the modulation signal observed by the DAMA/LIBRA collaboration either being due to dark matter or a different, so far not understood, physics phenomenon or background.
A series of R&D projects and experiments with alkali halide crystals as target are presently in construction phase or being upgraded [7][8][9][10][11]. In comparison to the before mentioned scintillation detector approaches, dual channel bolometric detectors using NaI-crystal have some advantages but also face challenges.

Detector design
The objective of COSINUS (Cryogeni Observatory for SIgnatures seen in Next-generation Underground Searches) is the development of a cryogenic scintillating calorimeter using a NaI-crystal as target -thereby providing a NaI detector with the possibility for active particle identification on an event-by-event basis thanks to the dual channel detection approach. The use of an alkali halide crystal for a scintillating calorimeter (with dual-channel readout) has already been studied e.g. in [12,13]. At first glance, the low Debye temperature of NaI is not favorable for a detector relying on good phonon propagation properties. However, a thermometer with low heat capacity and high sensitivity as e.g. a transition edge sensor (TES) has a realistic chance to provide a sufficient level of sensitivity [12,13]. A major challenge to overcome is the hygroscopic nature of NaI which asks for a handling of the detector -from the production to the assembling phase -in dry atmosphere using glove box systems and, hence, an adapted production procedure of the detector. Another approach is to coat the fragile NaI-crystal with a protective layer, as was studied by the authors of [14].
While COSINUS will use undoped NaI-crystals, experiments which are only measuring scintillation light typically apply 0.1 % of thallium (Tl) as dopant (e.g. [6,[9][10][11]15]). However, we want to stress that such small amount of Tldopant does not affect the potential interaction with dark matter, but solely serves to enhance scintillation light output (at room temperature).
Particles scattering in NaI-crystals do not only create phonons, but also scintillation photons. At temperatures of few milli-Kelvin, both the phonon signal in the crystal and the signal arising from the absorption of scintillation photons in a suitable light absorber can be detected by the use of superconducting thin film thermometers: as thin film technology is difficult to be directly applied to the hygroscopic NaI-crystal, a carrier disk made of a more robust scintillating material (e.g. CdWO 4 ) is planned to be used instead. The carrier disk (diameter of 40-50 mm and about 1 mm thickness) is designed to exceed in diameter the size of the NaI-crystal and carries the thermometer: a transition edge sensor (TES) of CRESST-type consisting of a thin tungsten film (200 nm, W-TES) directly evaporated onto the carrier crystal. 1 The target crystals of undoped NaI can weigh up to ∼200 g, depending on the achievable performance in function of the crystal mass. Using a carrier crystal allows to exclude the fragile NaI-crystal from the whole TES production procedures (thin film production by evaporation and sputtering processes, photo-lithography, chemical etching and intermediate cleaning). Once the TES on the carrier crystal is ready, the NaI is attached to it by a small amount of e.g. silicon oil, as the last step of detector mounting in a dedicated glove box. The delicate NaI-crystal is pressed to the carrier only by its own weight, in this way avoiding the need for metal clamping of the crystal. The phonon propagation from the NaI to the TES on the carrier occurs via the silicon oil interface.
In order to efficiently reject any α-related background, e.g. a recoiling nucleus, which could mimic a dark matter signal, there must not be any non-scintillating surfaces in the line-of-sight of the target crystal. An elegant way to tag and reject such backgrounds is a completely active surround-ing of the crystal [17]. Thus in COSINUS, the light detector is planned to consist of a beaker-shaped work-piece made from high purity silicon serving two purposes: scintillation light detection and a fully active surrounding of the target crystal. The beaker will exhibit dimensions of about 40-50 mm in diameter and height, the wall-thickness will be about 600 µm. The front face of the polished silicon beaker will be equipped with a W-TES, optimized in size for the purpose of light detection. Despite being quite macroscopic devices, such kind of light detectors have shown to be highly performing, achieving a baseline noise of below σ 10 eV 2 [18].
With the help of a dedicated copper structure the target crystal will be held inside the beaker-shaped light detector. Since the carrier crystal is chosen to be slightly larger in diameter than the target crystal and the light detector, an almost 100 % coverage of the target crystal should be achievable. 3 A schematic drawing thereof can be seen in Fig. 1.
A particle interaction in the target crystal mainly induces a thermal signal detected by the W-TES. This so-called phonon-signal provides a precise measurement of the deposited energy, independent of the type of particle [19,20], often referred to as unquenched channel. 4 Interactions which take place in the carrier crystal itself can be discriminated from energy deposits in the NaI-crystal by pulse-shape analysis [18,21].
Simultaneously to the phonon signal a small fraction of the energy deposited in the NaI-crystal is emitted in form of scintillation light and detected by the high-purity silicon light-absorbing beaker. The light signal allows to identify the type of interacting particle, as the amount of scintillation light strongly depends on it. Therefore, the light to phonon ratio, referred to as light yield, is characteristic for each type of event. Betas and gammas produce the most light and get assigned a light yield value of one by definition. Other types of particles (e.g. α-particles, neutrons inducing nuclear recoils off iodine and sodium) exhibit a lower light yield quantified by the so-called quenching factor [22,23]. The QF is defined as the ratio of scintillation light produced by an interacting particle of type X to the scintillation light produced by a gamma of the same deposited energy.
The dual channel detection approach of COSINUS aims to allow for a highly efficient discrimination and rejection of background events (e − /γ s, α-particles) from potential signal events (nuclear recoils) in the NaI target. The low temperatures required for detector operation will be supplied by a 3 He/ 4 He-dilution refrigerator. The test cryostat of the Max-Planck-Institute is the favored choice for all the measurements necessary for the first basic prototype testing. The hall C CUORE/CUPID R&D dilution refrigerator, since equipped with passive shielding for means of background reduction, is considered for final prototype testing and a longer performance run. Both refrigerators are located in the LNGS underground site for obvious reasons of shielding against cosmic radiation.
The results from CsI measurements give a first reference point on the already achieved performance of an alkali halide crystal operated as cryogenic calorimeter. With a sequence of tests we commissioned the carrier crystal-absorber crystal concept and we achieved, with two different CsI crystals, an energy resolution of 0.95 keV/0.70 keV, corresponding to a trigger threshold of 4.7 keV/3.5 keV [21]. 5

Performance from simulated data
We base the performance estimate of a COSINUS detector module on the following assumptions: -light detector energy resolution from standard performance of silicon beaker light detector operated in CRESST-II phase 2 (σ = 0.11 keV ee ) -phonon detector resolution of σ = 0.2 keV (corresponding to a threshold of 5σ = 1 keV), -4 % of deposited energy in the NaI detected in the light detector, -a gross exposure (before cuts) of 100 kg-days. Based on results from CRESST-II [19,24], we conservatively estimate an efficiency of 50 % down to an energy of 2 keV and a linear decrease to 20 % efficiency at the threshold of 1 keV, -a flat background contribution down to energy threshold of 1 count/(keV kg days) and an additional contribution from 40 K at a level of 600 µBq.
In the following the above assumptions will be motivated in more detail.
In COSINUS we aim to bring the performance of a NaIbased detector in line with other scintillating bolometers, rendering an energy resolution of σ = 0.2 keV, corresponding to an energy threshold of 1 keV, feasible. From our understanding the key points to be addressed are the optimization of the sensitivity of the TES itself and of the interface between NaI-crystal and carrier crystal.
In [21] we measured two CsI-crystals; in average around 7 % of the energy deposited in the crystal were determined as scintillation light (for an e − /γ -event). From measurements in [13] it is known that CsI roughly emits 3.5 times more energy in form of scintillation light as NaI. However, for COSINUS beaker-shaped light detectors will be used which are found to collect more than twice as much light the light detector of conventional CRESST design used for the measurement in [21]. Thus, in total we expect 7 %/3.5 × 2 = 4 % of the deposited energy in a NaI-crystal to be detected in the light detector.
Using the above assumptions a simulation is carried out -the outcome is depicted in the light yield-energy plane in Fig. 2 with data points arising from simulated contributions of e − /γ -events in black and a potential signal in red. To stick to the DAMA/LIBRA background budget, an overall flat background contribution of 1 count/(keV kg days) is assumed.
NaI-crystals typically show a contamination with 40 K which undergoes an electron capture, either directly to the ground state of 40 Ar (0.2 % branching ratio) or to an excited state of 40 Ar de-exciting by the emission of an 1.46 MeV gamma (10.55 % branching ratio). Both processes are accompanied by the energy releases of 3.2 and 3 keV. This contamination makes up for the line at about 3 keV in the e − /γ -band in Fig. 2 corresponding to a 40 K-activity of 600 µBq as measured by the DAMA/LIBRA collaboration [25].
The bands in Fig. 2 are calculated on the basis of the bands determined in the CsI measurement (a detailed description on the underlying model is e.g. given in [22]), the experience of beaker-shaped light detectors operated in CRESST and the assumptions listed above. Thereby, the lines depict the lower and upper 90 %-boundaries for the bands corresponding to different event classes to be observed in NaI (electron recoils in black, nuclear recoils off Na and I in Fig. 2 Simulated data for an exposure before cuts of 100 kg-days applying certain performance assumptions as discussed in detail in the text. The bands correspond to the following event classes: electron recoils in black, scatterings off Na and I in blue and green. Depicted in black dots are events originating from a flat background contribution of 1 c/(keV kg days) and a 40 K contamination with an activity of 600 µBq. The contribution from 40 K is visible as a line at around 3 keV in the e − /γ -band. A hypothetical WIMP scattering off I and Na nuclei yields the events marked in red. In addition, the magenta boxes indicate the regions which may contribute, apart from the e − /γ -band, to the DAMA/LIBRA energy range of positive modulation detection corresponding to 2-6 keV ee blue and green, respectively). Thus, in between two corresponding lines 80 % of the events of the respective event class are expected. The widths of the bands is dominated by two effects: firstly by the finite baseline resolution of phonon and light detector and secondly by Poissonian fluctuations in the number of scintillation photons produced. The mean of a band obviously depends on the type of particle and the respective quenching factor. For scatterings off Na and I we use the energy-dependent values reported by Tretyak [23].
In summary, the overall e − /γ -background (black dots) is given by the 40 K contamination on top of the constant background level of 1 count/(keV kg days). With the performance and backgrounds assumed, we expect five counts below the mean of the Na-band for a gross exposure of 100 kg-days and an energy threshold of 1 keV. Thereby, one half of the counts originates from the constant background and the other half from the 40 K contamination, which roughly corresponds to 0.9 % of the activity in the full double-peak. The expected leakage quickly drops with increasing energy, due to the enhanced separation between e − /γ -band and Na-band. Additionally, the above numbers show that the leakage is to a large extent caused by the 40 K contamination yielding a considerably diminishing leakage for energies above 3 keV, where we expect less than one leakage event (for 100 kg-days before cuts), which is commonly denoted background-free. Considering leakage to the full 80 % Na-band (as depicted in Fig. 2) the values change to 26 events expected above the Fig. 3 Projected sensitivity for a NaI-based experiment using COSI-NUS detector technology [solid blue band, 1σ confidence level (C.L.)] for spin-independent elastic dark matter nucleus scattering. Recent results from experiments using silicon and germanium targets are drawn in green [4,[28][29][30][31][32][33], results from CRESST-II (CaWO 4 ) are depicted in red [19,24,34,35] and exclusion limits using liquid noble gases in magenta [36][37][38][39]. Limits drawn in cyan correspond to bubble chamber technology and experiments with CsI target [7,8,40]. The light blue shaded regions correspond to the interpretation of the DAMA/LIBRA (NaI(Tl)) modulation signal by Savage et al. [27]. The benchmark point (blue cross) indicates the mass and cross-section chosen for the simulated WIMP contribution presented in Fig. 2. Gray-shaded regions in parameter space will be affected by coherent neutrino nucleus scattering on NaI mainly originating from solar neutrinos [41] threshold of 1 keV, out of which 13 events are attributed to 40 K-origin. For the complete Na-band the anticipated leakage drops below one event at 3.9 keV.
Because of the high discrimination power we only anticipate a few ten leakage events in the full nuclear recoil bands for the radiopurity levels reported by the DAMA/LIBRA experiment (1 count/(keV kg days) + 40 K: 600 µ Bq), while at the same time expecting O(10 3 ) signal events. Hence, we can tolerate a higher level of e − /γ -background (as e.g. the two times higher background recently reported for NaIcrystals by the KIMS collaboration [26]), in particular originating from 40 K, while still maintaining high sensitivity for a nuclear recoil signal -a distinct feature of COSINUS compared to other NaI-based dark matter searches [6,[9][10][11].
The red events in Fig. 2 result from a contribution of a hypothetical WIMP of 10 GeV/c 2 and 0.0002 pb which is consistent with the interpretation of the DAMA/LIBRA modulation signal by Savage et al. [27]. The blue-colored islands in Fig. 3 correspond to the DAMA/LIBRA signal regions for scatterings off Na and I. For clarity we indicated the assumed dark matter particle mass and nucleon cross-section for the simulated data presented here in form of a blue benchmark point (see Fig. 3).
What concerns the distribution of the hypothetical WIMP events we find in the simulation a total of 2386 events in the energy window 1-6 keV. About 45 % of these events are present in the energy interval from 1 to 2 keV and about 53 % within 2-6 keV. These numbers clearly show the benefit of a low threshold. Despite the conservative assumption on the cut efficiency being 20 % at threshold energy, the number of expected events roughly doubles when lowering the threshold from 2 to 1 keV, which is attributed to the anticipated exponential rise of the dark matter recoil spectrum. Above an energy of 6 keV (up to 100 keV) we find a total of 46 events.
Since WIMPs are expected to scatter coherently off the nucleus as a whole, the cross-section scales quadratically with the atomic mass number (A 2 ), thus preferring the heavy I over the light Na. As a consequence, the WIMP events in Fig. 2 appear almost symmetric to the I recoil band. However, the energy transferred in a scattering for the rather light WIMP considered here is enhanced for light nuclei due to kinematic reasons. Consequently, the majority of the 46 events above 6 keV is found in the Na recoil band. For this reason, cryogenic detectors providing low thresholds in combination with light target nuclei currently lead the field of direct dark matter searches in the low mass regime [24,28].
At last we want to discuss and underline the prominent features of such COSINUS detector in contrast to a purely scintillating detector considering the DAMA/LIBRA experiment as representative. In Fig. 2 we display magenta colored boxes indicating the different regions that may contribute to the positive modulation signal observed by DAMA/LIBRA in an energy window from 2 to 6 keV ee . Due to the lack of particle discrimination it remains unknown if the positive modulation signal is made from particles scattering off the electrons or purely off the nuclei in NaI(Tl).
For interactions with electrons the box extends from 2 to 6 keV (e − /γ -band), the same energy range as for DAMA/LIBRA. Instead, interactions on nuclei are quenched in the light channel, hence the energy interval for Na and I recoil events has to be corrected for by the respective quenching factor (energy-dependent QF are taken from [23]). For Na-recoils (QF ≈ 0.3) the region is restricted to about 6-20 keV as qualitatively indicated by the magenta box in the Na-recoil band. For I-recoils, due to the even higher light quenching effect (QF ≈ 0.1), the box is confined in an energy region of about 20-60 keV.
As already mentioned, a distinctive feature of the COSI-NUS technology is the unquenched phonon channel directly measuring the deposited energy quasi independent of the type of interacting particle. Taking into account the aimed for energy threshold of 1 keV, this would result in an improvement in detection threshold by a factor of about six for Narecoil events and a factor of about 20 for I-recoils in comparison to the sensitivity demonstrated by DAMA/LIBRA.
It deserves mentioning that we already achieved, by operating two scintillating calorimeters based on CsI an energy threshold as low as 4.7 and 3.5 keV, respectively [21]. Thus, such detectors already indicate an increased sensitivity for nuclear recoil events in comparison to the DAMA/LIBRA experiment with further refinements anticipated for the future [21].

Perspective and conclusion
As can be seen, competitive sensitivity also to other direct dark matter searches may be reached, in particular in the mass range of 1-10 GeV/c 2 (see Fig. 3). This makes NaI an attractive target for direct dark matter search by itself, beyond offering a material-independent cross-check of the DAMA/LIBRA signal. Of special benefit is the presence of light (Na) and heavy (I) target nuclei which is, however, also given for CaWO 4 . A CaWO 4 -target (CRESST-II/III) offers similar discrimination capabilities and background level, but presumably better phonon properties. The latter are decisive for the achievable energy threshold and resolution rendering CaWO 4 better suited for the low-mass range. Currently, NaIcrystals do not reach radiopurity levels of Ge-based experiments (SuperCDMS). This is also true for liquid noble gas experiments (LUX/XENON) which, in addition, are easier to scale to large target masses, but fall behind cryogenic experiments regarding energy resolution and threshold. For masses of O(≥10 GeV/c 2 ), however, exposure outweighs the achievable threshold, making liquid noble gases the leading technology in this mass range.
The comparison of results from different experiments, as depicted in Fig. 3, only holds under certain assumptions concerning the dark matter halo and the interaction of dark matter with standard model particles. This consideration is further augmented by the use of different target materials, as the impact of the mentioned uncertainties significantly depends on the target material.
Obviously, material dependencies will be ruled out in the evaluation of COSINUS and DAMA/LIBRA data as both are using NaI as target. 6 As Fig. 3 clearly shows, the anticipated sensitivity of a COSINUS detector is about two orders of magnitude below the interpretation of the DAMA/LIBRA claim, assuming a standard dark matter halo and elastic WIMP-nucleus scattering (light-blue regions in Fig. 3 corresponding to recoils off Na and I, respectively [27]).
The enhanced sensitivity, thereby, is driven by two key factors. Firstly, the particle discrimination via the simultaneous measurement of phonon signal and scintillation. Secondly, the better energy resolution going along with a lower threshold (of 1 keV, independent of the type of particle). The latter is of special benefit due to the expected exponential rise of the dark matter recoil spectrum towards low energies.
Thus, even with a moderate exposure, COSINUS technology has the potential to add knowledge on the underlying nature of the DAMA/LIBRA signal, in particular on the question whether the signal originates from nuclear recoils or not.
For a hypothetical dark matter particle interacting with the electrons of the target material the advantage of particle discrimination vanishes, as the dominant background will also be found in the e − /γ -band. This renders radiopurity the key issue, for the experimental set-up and the NaI-crystals. The benefit of an excellent energy resolution and a low threshold still persist. In this scenario a potential future experiment (based on COSINUS technology) would require a higher target mass (O(10-100 kg)) measuring for multiple annual cycles. Large-scale bolometric experiments are feasible, as e.g. proven by CUORE [42].
Thus, a potential future dark matter experiment, increased in target mass and based on COSINUS technology, has promising prospects to give new insight on the long-standing DAMA/LIBRA claim.