First array of enriched Zn[Formula: see text]Se bolometers to search for double beta decay.

The R&D activity performed during the last years proved the potential of ZnSe scintillating bolometers to the search for neutrino-less double beta decay, motivating the realization of the first large-mass experiment based on this technology: CUPID-0. The isotopic enrichment in [Formula: see text]Se, the Zn[Formula: see text]Se crystals growth, as well as the light detectors production have been accomplished, and the experiment is now in construction at Laboratori Nazionali del Gran Sasso (Italy). In this paper we present the results obtained testing the first three Zn[Formula: see text]Se crystals operated as scintillating bolometers, and we prove that their performance in terms of energy resolution, background rejection capability and intrinsic radio-purity complies with the requirements of CUPID-0.

Introduction

Neutrino-less double beta decay (0DBD) is a hypothesized nuclear transition in which a nucleus decays emitting only two electrons. This process cannot be accommodated in the Standard Model, as the absence of emitted neutrinos would violate the lepton number conservation. For this reason, its prized observation would have several implications for particle physics, astrophysics and cosmology. According to the majority of theoretical frameworks, for 0DBD to happen neutrinos must be Majorana particles. This means that, in contrast to all the other known fermions, they must coincide with their own antiparticles. In this scenario, the half-life of the process T is determined by the effective Majorana neutrino mass m, a linear superposition of the three neutrinos mass eigenstates weighted by the elements of the first row of the PMNS neutrino mixing matrix and including two Majorana phases, providing the dependance: T.

Current lower limits on T exceed 10 years for several possible emitters, corresponding to upper bounds on m of a few hundreds of meV [1]. Moreover, recent measurements made by the KamLAND-Zen collaboration using a xenon-loaded liquid scintillator allowed to set a limit on T(Xe) of 1.110 years at 90 C.L., corresponding to m < 60–161 meV [2].

Among the several experimental approaches proposed for the search of 0DBD, cryogenic calorimeters [3] (historically also called bolometers) stand out for the possibility of achieving excellent energy resolution (of the order of 0.1), efficiency (>80) and intrinsic radio-purity. Moreover, the crystals that are operated as bolometers can be grown starting from most of the 0DBD emitters, enabling the test of different nuclei.

The CUORE-0 experiment proved the potential of the bolometric technique by reaching an unprecedented sensitivity on the half-life of  [4]. The source was embedded in 52 bolometers (for a total mass of 39 kg), operated between March 2013 and March 2015 at Laboratori Nazionali del Gran Sasso (LNGS, Italy). The combination of the 9.8 kg(Te)y exposure collected by CUORE-0 with the 19.75 kg(Te) years of its ancestor, Cuoricino, allowed to set a new limit on the T(Te): T4.010 years at 90 C.L., corresponding to m270–760 meV.

The evolution of CUORE-0, named CUORE [5], will start operation before the end of 2016 with 988 bolometers, pointing to a sensitivity of T(Te)9.510 years at 90 C.L. (m50–130 meV). This remarkable sensitivity will allow to touch, but not to fully explore, the m region of 10–50 meV, corresponding to the inverted mass hierarchy of neutrinos (see Refs. [1, 6–8] and references therein for a complete discussion on the relationship between the 0DBD and the neutrino mass hierarchy). The entire coverage of this region, indeed, requires a level of sensitivity of about 10–10 years, depending on the isotope [9]. The CUPID (CUORE Upgrade with Particle IDentification [10, 11]) interest group is defining the technological upgrades that will allow to reach this target with a CUORE-size bolometric detector. Such an ambitious goal poses a stringent requirement on the background, that must be close to zero at the tonneyear level of exposure.

The experience gained during the years of R&D activity for CUORE allowed to determine that the ultimate limit in the background suppression resides in the presence of -decaying isotopes located in the detector structure. Thus, the main breakthrough of CUPID with respect to CUORE will be the addition of independent devices to measure the light signals emitted either from Cherenkov radiation (in  [12]) or from scintillation (in scintillating crystals [13] like ZnSe [14, 15], ZnMoO [16-22], LiMoO [23, 24] and many others). The different light emission of electrons and particles will enable event-by-event rejection of interactions, suppressing the overall background in the region of interest for 0DBD of at least one order of magnitude. During the last years, several R&D activities demonstrated that the detection of Cherenkov light emitted by electrons crossing the crystals of CUORE can provide an effective background rejection [25]. Despite the encouraging results [26-33], none of these activities is yet mature for a next-generation experiment. A different approach was followed by the LUCIFER [34] and LUMINEU [35] collaborations, that tested different scintillating bolometers and proved that the read-out of the much more intense scintillation light (a few keV or tens of keV depending on the crystal, with respect to 100 eV produced by Cherenkov radiation) simplifies the suppression of the background. Moreover, choosing crystals containing high Q-value 0DBD emitters such as Se and Mo provides a natural suppression of the background produced by the environmental radioactivity, that drops above the 2615 keV -line of . The encouraging results obtained with large ZnSe bolometers motivated the endeavor in realizing a first large mass demonstrator of the scintillating bolometers technology: CUPID-0. A second phase of CUPID-0, not described in this paper, will include also enriched Li MoO and ZnMoO crystals, which showed excellent performance as scintillating bolometers. CUPID-0 aims at proving that the interactions due to particles can be efficiently suppressed, providing a path to ensure the achievement of the zero-background level required by CUPID.

CUPID-0

The CUPID-0 detector will consist of 30 ZnSe bolometers arranged in 5 towers, for a total mass of 13.2 kg (7.0 kg of ). The Q-value of (2997.9  0.3 keV [36]) is among the highest, and the half-life of the 2 mode is long enough (T=(9.6  0.1(stat)  1.0(syst))10 years [37]) to prevent 2 background contributions in the region of interest. To date, all the crystals (ZnSe and Ge light detectors) were grown and delivered to LNGS. The assembly of the CUPID-0 is expected to be completed before the end of Summer 2016.

The Se powder was enriched at URENCO Stable Isotope Group (Netherlands) to overcome the rather poor natural isotopic abundance of (8.73 [38]). We measured the radio-purity of the 96.3 enriched powder with HP-Ge spectroscopy and, since we found no evidences of contaminations, we set upper limits of 61, 110 and 74 Bq/kg at 90 C.L. for Th, U and U, respectively [39]. We measured with the same sensitivity also Zn metal, and found no signatures of Th, U and U.

The ZnSe synthesis, as well as the crystals growth were performed at the Institute for Scintillation Materials (ISMA, Kharkov, Ukraine). To prevent contaminations, the synthesis and the subsequent powder purification treatments were performed in Argon atmosphere using tools/containers selected according to their radio-purity. After the test run described in this paper, all the 95.4 enriched CUPID-0 crystals were cut, shaped and polished in order to obtain bolometers with 4.4 cm diameter and 5.5 cm height, for an average mass of 440.5 g. These operations were carried out in a clean room at LNGS following the protocols for radio-purity described in Ref. [40]. The procedures that we developed for the ZnSe synthesis and purification allowed to reach a yield of 98, with 2 irrecoverable loss of enriched Se. The yield of the crystals production was reduced by losses of material evaporated during the crystal growth. The crystals presented in this paper were grown using graphite crucible, with an irrecoverable loss of material of about 15. The growth procedure was then optimized in order to reduce the material loss down to 2 using glassy graphite crucible. The fraction of material discarded during the manufacturing process (about 40) can be reprocessed to obtain other ZnSe crystals. More details about the synthesis, crystal growth and processing can be found in Ref. [41].

The scintillation light emitted by each ZnSe will be measured by two light detectors [42], facing the top and bottom surfaces of the crystal (see Fig. 1a). The light detectors are also bolometers made of disk-shaped Ge crystals (4.4 cm diameter) grown by UMICORE using the Czochralski technique. Since the energy emitted as scintillation light is only a few of the total, the Ge crystals must be very thin to ensure a low heat capacity and, thus, a higher signal. For this reason, the light detectors were cut in 170 m thick wafers, that were polished and etched on both sides. The face showing the best optical properties was coated with a thin layer (70 nm) of SiO in order to increase the light absorption [43, 44].

Fig. 1

a Lateral view of a single module: the ZnSe is placed between two light detectors (Ge-LD), not visible in this scheme because of their small thickness (about 170 m); all the detectors are held in the copper structure using PTFE elements. Top b and 3D c views of the CUPID-0 detector. d 3D view of the CUPID-0 detector hosted in the CUORE-0 cryostat

To convert temperature variations into readable voltage signals, each ZnSe and Ge crystal will be equipped with a Neutron Transmutation Doped (NTD) Ge thermistor using a semi-automatic gluing system. In addition, a Si Joule heater will be glued to the crystals for the offline correction of thermal drifts by heat pulses injection [45].

The ZnSe and Ge bolometers will be mounted in an Oxygen Free High Conductivity (OFHC) copper structure, serving as thermal bath to cool the detectors at 10 mK. To account for different thermal contractions of the detectors and copper, the crystals will be secured to the copper frames using PTFE elements, acting also as weak thermal coupling to the bath. The final layout of the CUPID-0 detector, that will be hosted in the same cryostat used for CUORE-0 with an upgrade for the reduction of the microphonic noise, is shown in Fig. 1.

To test the CUPID-0 assembly line, as well as to verify the compliance of the final detectors with the requirements on energy resolution, background rejection capability and intrinsic radio-purity, we mounted an array of three ZnSe crystals and performed a bolometric run. In this paper we present the results of the test run and discuss the perspectives for CUPID-0.

Detector

The detector consists of three ZnSe crystals (for a total ZnSe mass of about 1.32 kg) and four light detectors, arranged in a single tower following the layout depicted in Fig. 1.

Apart from the number of crystals, there are only few differences with the final CUPID-0 protocol for the detector assembly: (1) the cleaning procedure defined for the surface of the crystals in CUPID-0 was not yet applied to the present crystals, as we were interested mainly in the bolometric performance and crystal bulk contaminations, (2) the detector was mounted above-ground, shortly after the crystal growth, while the CUPID-0 towers will be assembled in a Radon-free underground clean room using crystals stored underground to reduce the cosmogenic activation, (3) during the measurement the ZnSe were permanently exposed to smeared sources of Sm (Q-value2.3 MeV) to assess the background rejection capability in a short measurement time.

We used two classes of NTD Ge thermistors, characterized by different dimensions: the large sensors (2.83.01.0 mm) were attached on ZnSe bolometers, while the smaller ones (2.82.00.5 mm) were glued to the light detectors that, because of their much smaller mass, require sensors with lower heat capacity. Two stages of the detector assembly are shown in Fig. 2.

Fig. 2

Two stages of the tower assembly. Left a light detector equipped with an NTD Ge thermistor and a Si Joule heater is mounted on top of a ZnSe bolometer, surrounded by a 3M VM2002 reflecting foil. Right the second ZnSe bolometer is placed on top of the previous light detector

The ZnSe array was anchored to the coldest point of a He/He dilution refrigerator located in the Hall C of LNGS. For technical reasons, the refrigerator could not be operated below 20 mK, which is however satisfactory for our purposes. The final CUPID-0 towers will reach a much lower base temperature (about 10 mK) exploiting the cryostat that was used for the CUORE-0 experiment.

The signals were amplified by high gain voltage amplifiers with input stages based on Silicon JFETs, filtered by a 6-pole low-pass Bessel–Thomson filter [46-48] with roll-off rate of , and fed into an 18 bit NI-6284 PXI ADC unit. Due to the different time development of heat and light pulses (shown in Fig. 3) we set a cut-off frequency of 63 Hz for ZnSe and 200 Hz for light detectors. When the software trigger of each ZnSe fired, we saved on disk 5 s long time-windows sampled at 1 kHz. The light detectors were acquired with a shorter window of 250 ms and a higher sampling frequency (2 kHz).

Fig. 3

Response of ZnSe bolometers obtained by averaging tens of pulses with energy of 2615 keV. Inset corresponding scintillation light detected by the top (blue, continuous line) and bottom (orange, dotted line) light detectors. The side of the top light detector facing the ZnSe bolometer was coated with SiO

In addition to the derivative trigger, we used a second trigger that forced the acquisition of the light detectors every time the trigger of a ZnSe fired.

To maximize the signal-to-noise ratio, the waveforms were processed offline with the optimum filter algorithm [49, 50]. Unlike the classical nuclear pulse shape processing (gaussian, semi-gaussian, CR-(RC) etc.) the optimum filter estimates the signal at each frequency with weights equal the signal-to-noise ratio at the same frequency, embedding in its shape the proper timing parameters. In our case it behaves as a pass-band filter with the lower frequency optimized for 1/f noise rejection, while the higher frequency is limited by signal roll-off, minimizing the contribution of white noise. Since light pulses feature a worse signal-to-noise ratio, we used a dedicated algorithm, exploiting the known time-delay between heat and light signals, to compute their amplitude [51].

The general features of the ZnSe bolometers are reported in Table 1. The working resistance of the NTD Ge thermistor R depends on the detector temperature: the lower the temperature, the larger the resistance. The electronics channels are optimized for R of the order of tens/hundreds of M (corresponding to a crystal temperature lower than 10 mK) but the high base temperature of the cryostat prevented the achievement of such a resistance, somehow spoiling the detector performance. Given the low values of R (reported in Table 1), we estimated a noise contribution from the electronics of 7.0 keV FWHM for ZnSe-1, 14.3 keV for ZnSe-2 and 18.3 keV for ZnSe-3. These values are in full agreement with the measured energy resolution of the detector baseline (Table 1), leading to the consideration that a lower cryostat temperature would have provided a much better baseline resolution, as observed in other ZnSe prototypes [15]. The analysis of the noise power spectra of the detectors allowed to infer that the degradation of the baseline resolution could be mainly ascribed to the series noise (both 1/f and white) of the very front-end system, operated in unmatched conditions. Nevertheless, our primary interest was the study of the detector performance at high energy where, as shown in the next section, the electronics noise is negligible. Furthermore, also the response of the bolometer, defined as the voltage signal produced by an energy deposit of 1 MeV, was affected by the higher than usual detector temperature. A lower temperature of operation, to be expected in the next runs, will allow to increase the energy conversion gain of the detectors, providing an enhancement of the baseline resolution. To face both the situations, low temperature with large detector impedance and warm temperature with small detector impedance, we developed a new set of preamplifiers able to match the latter condition with a noise smaller by a factor of 2.5 and input capacitance larger by a similar factor with respect to the actual setup.

Table 1

Resistance of the NTD Ge thermistor in working conditions R, FWHM noise energy resolution after the optimum filter, voltage signal (in V) produced by an energy release of 1 MeV

	R_work (MΩ)	Baseline noise (keV FWHM)	Response (μV/MeV)
Zn^82Se-1	0.20	7.0	7.1
Zn^82Se-2	0.22	14.1	3.7
Zn^82Se-3	0.17	18.6	2.9

Energy resolution

The ZnSe bolometers were energy calibrated using Th and K sources, emitting rays up to 2615 keV. We fitted the most intense peaks to check the linearity of the detector response, as well as to determine the energy resolution (see Fig. 4 for ZnSe-1).

Fig. 4

FWHM energy resolution as a function of the energy (ZnSe-1) for the most intense peaks produced by Th and  K sources. The point at zero energy is the baseline energy resolution reported in Table 1. The black line is the fit function: FWHM(E) = FWHM+E. The green dotted lines indicate the Q-value

Since the natural radioactivity drops above the 2615 keV line, we could not rely on calibration peaks at the Q-value of . For this reason, we estimated the energy resolution in the ROI by fitting the energy resolution of the calibration peaks as a function of the energy with the following function:where accounts for the contribution of the electronic noise (see Table 1), and and are arbitrary coefficients. At the Q-value we derived a FWHM energy resolution of 30.1±1.7 keV (ZnSe-1), 29.7±1.4 keV (ZnSe-2) and 30.2±1.7 keV (ZnSe-3), proving the reproducibility of the ZnSe bolometers.

These values are about a factor 3 larger than the initial target of the experiment (10 keV FWHM [34]). Even if a further cooling of the crystals is expected to improve the resolution, we will demonstrate in the next sections that ZnSe bolometers with 30 keV FWHM would anyway comply with the requirements of CUPID-0.

Performance of light detectors

The sensitivity of light detectors is crucial for CUPID-0, as the background rejection capability relies on the different light emission of electrons and particles. In this run, all the light detectors were permanently exposed to Fe X-rays sources (characterized by peaks at 5.9 and 6.4 keV) for energy-calibration. To further investigate the reproducibility of light detectors, we assembled a second array, with three ZnSe interleaved by four Ge disks. Two of these light detectors were equipped with large NTD Ge thermistors, similar to those attached on ZnSe crystals, to investigate the effect of the sensor heat capacity on the bolometric performance. The features of the light detectors are reported in Table 2.

Table 2

Features of light detectors equipped with small NTD Ge sensors: and are the rise and decay times, computed as the time difference between the 90 and 10 of the leading edge and as the time difference between the 30 and 90  of the trailing edge respectively. The other parameters were defined in Table 1. In the last two lines we report the features of light detectors equipped with large NTD Ge sensors

R_work (MΩ)	Response (μV/MeV)	Baseline noise (eV FWHM)	τ_r (ms)	τ_d (ms)
0.63	1.0 × 10^3	134	1.7	3.4
1.45	1.6 × 10^3	92	1.9	5.7
0.71	1.0 × 10^3	103	1.8	5.4
0.89	1.7 × 10^3	76	1.7	5.1
0.85	1.7 × 10^3	89	1.7	5.2
0.72	1.1 × 10^3	108	1.8	5.4
0.65	1.2 × 10^3	260	1.8	9.2
0.37	0.6 × 10^3	160	1.9	11.0

As explained in Sect. 3, the high cryostat temperature could produce an electronic noise limiting the detector energy resolution. Nevertheless, it is worth highlighting that, in absence of any signal amplification, energy resolutions as good as those reported in Table 2 have never been achieved with wide area Ge light detectors equipped with conventional NTD Ge thermistors. The detector response in terms of signal height, as well as the energy resolution, showed variations lower than 40 across the 6 light detectors equipped with small NTD Ge sensors, proving the reproducibility of the light detectors performance in view of CUPID-0. These results can mainly be ascribed to the improvement of the interface between the Ge disks and the NTD Ge thermistors provided by dedicated surface treatments (etching and polishing), as well as by the semi-automatic gluing system that enhanced the reproducibility of the detectors features

Background study

Cosmogenic activation

To investigate the cosmogenic activation of ZnSe crystals we exploited the comparison between data extracted from the cosmogenic activation calculation software Activia [52] and the measurements performed with this array. The isotope of (enriched) Selenium with the highest activation rate is Se (Q-value864 keV, T118 d) which decays via electron capture with a rather complex combination of de-excitation ’s and X-rays. None of the enriched crystals showed evidences of Se, and we set a 90  C.L. upper limit on its activity of 14 Bq/kg. All the other isotopes are expected to have a much lower activation rate. Among them, the only emitters that could produce background events for the 0DBD of are V (Q-value 4012 keV, T16 d) and Co (Q-value 4566 keV, T79 days). Both these isotopes have activation rates about 3 orders of magnitude lower than Se. Furthermore, due to their short half-lives, they will not participate to the background of CUPID-0, as the crystals will be permanently stored underground until operations.

We performed a similar analysis to investigate the activation of Zn. In this case, the isotope with highest activation rate is Zn (Q-value1352 keV and T244.26 d). The signature of Zn was observed in all the ZnSe crystals with similar activities of 3.06±0.44 mBq/kg. Also in this case, V and Co would be the only dangerous emitters produced by cosmogenic activations, with activation rates much lower than Zn. As said before, due to the low activation and the short half-lives, these isotopes will not be of concern for CUPID-0.

Nevertheless, since this test run was performed shortly after the crystal production and delivery, we expect Co to contribute about 1 /day in the total energy spectrum of the array above 2615 keV.

Alpha background rejection

The large light output of ZnSe crystals and the excellent performance of the cryogenic light detectors provided a very efficient background rejection. The three crystals featured a similar light yield (LY), defined as the amount of light (in keV) measured when an interaction of 1 MeV occurs in the ZnSe bolometer. ZnSe crystals operated at cryogenic temperatures show a rather peculiar (and not yet fully understood) feature: the LY of electrons is lower with respect to the one of particles. With this array we measured LY ranging from 3.3 to 5.2 keV/MeV, to be compared with LY ranging from 9.1 to 14.1 keV/MeV.

Even if the different LY allows to identify and separate electrons and events, we used a more powerful estimator for background rejection: the shape of the light pulses. As shown in the inset of Fig. 5, the light pulses produced by interactions are slower than light pulses produced by ’s of the same energy.

Fig. 5

Shape parameter of a light detector as a function of the energy released in ZnSe-1 (530 h). The red and blue lines indicate the 2 (continuous) and 3 (dotted) and bands respectively. events produced by the smeared Sm source (below 3 MeV) and by contaminations of the crystal bulk (peaks above 5 MeV) can be easily rejected, in particular in the region of interest for the 0DBD (green lines). The other ZnSe showed similar results. Inset time development of light pulses produced by (blue) and (red) interactions with energy of about 2.6 MeV

To exploit this difference, we used a shape-sensitive parameter (SP) computed on the optimum filtered pulse and defined as:where y is the pulse, A and i its amplitude and maximum position, s the ideal signal pulse scaled to unitary amplitude and aligned to y, the right width at half maximum of s. The Shape Parameter of light pulses computed for the best light detector is shown for ZnSe-1 in Fig. 5.

Because of the worse signal-to-noise ratio, the and bands become wider at lower energies. To compute the discrimination capability at the 0DBD energy, we divided the and bands in intervals and made Gaussian fits to derive the mean value () and the standard deviation () of the Shape Parameter in each interval. We fitted the energy dependence of (E) and (E) with polynomial functions and defined the Discrimination Power (DP) as a function of the energy as:and found at the Q-value. The same analysis was made on the other two crystals, obtaining a DP of 11 and 10 for ZnSe-2 and ZnSe-3 respectively.

We report in Fig. 6 the zoom at high energy of the spectrum collected with the ZnSe array in about 530 h, before and after the rejection of the background and of the events occurring simultaneously in more than one crystal.

Fig. 6

High-energy spectrum of the ZnSe array collected in 530 h before (gray) and after (red) the rejection of the background. The green lines indicate a FWHM region centered on the Q-value

The number of events above the Tl peak is very small if we consider that the cryogenic set-up used for this test did not provide an efficient shield against Bi produced in Rn-contaminated air, and that the crystals were measured immediately after production (when the activation could still give rise to an important background). For CUPID-0 we foresee an even better result, as the cryostat used by CUORE-0 will be equipped with a more effective Rn-free system and a more efficient shielding. Moreover, storing the crystals underground will suppress the activity of isotopes like Co (T79 d) that, as explained in the previous section, was expected to produce a non-negligible background in this run.

Crystal contaminations

The radio-purity of ZnSe crystals is of primary concern for the achievement of the target sensitivity, as the long-living isotopes of the natural U and Th chains may give rise to an irreducible background in the region of interest around 3 MeV.

Due to their very short range, particles can give a clear indication of bulk or surface contaminations. If the -decaying isotope is located in the crystal bulk, both the particle and the nuclear recoil are absorbed by the bolometer, giving rise to a peak at the Q-value of the decay. On the contrary, if the contaminant is located on the crystal surface, ’s (and with much lower probability nuclear recoils) may leave the crystal before being stopped and the event is recorded with an energy lower than the Q-value. Only for very shallow surface contamination the line (about 100 keV below the Q-value) becomes visible.

To extract the crystal contaminations in the isotopes of the U and Th chains, we analyzed the region of the total energy spectrum, reported in Fig. 7. The same analysis was performed on the energy spectra of each ZnSe to highlight possible differences, and the results are shown in Table 3. In some cases we had no evidences of contaminations, thus we reported the 90  C.L. upper limit. For this purpose, we defined the signal as the number of events falling in the energy region [Q-value , Q-value ] and the background as the average number of events falling in the 3 side-bands of this interval. Following the Feldman–Cousins approach, we computed the 90  C.L. upper limit on the number of events and, correcting for the branching ratio of the nuclide, we inferred the upper limit on the activity. The lines not quoted in the table are the two 6 MeV lines of Th because they produce a pile-up event rejected as a deformed signal by pulse shape cuts, and the two lines that in both chains are summed to the following beta emission (Bi-Po events) producing a continuum above 8 MeV.

Fig. 7

region of the ZnSe array collected in 530 h. The spectrum was energy-calibrated using the nominal energy of the most prominent peaks

Table 3

Crystal contaminations in the isotopes of U and Th chains derived from the analysis of the region (see Fig. 7). We made the same analysis on the energy spectrum of each ZnSe bolometer to highlight differences among crystals. Isotopes with long half-lives, that can break the secular equilibrium, are highlighted in bold type

	Zn^82Se-1 (μBq/kg)	Zn^82Se-2 (μBq/kg)	Zn^82Se-3 (μBq/kg)	Array (μBq/kg)
^232 Th	13 ± 4	13 ± 4	<5	7 ± 2
^228 Th	32 ± 7	30 ± 6	22 ± 4	26 ± 2
^224Ra	29 ± 6	26 ± 5	23 ± 5	27 ± 3
^212Bi	31 ± 6	31 ± 6	23 ± 5	29 ± 3
^238 U	17 ± 4	20 ± 5	<10	10 ± 2
^234 U+^226 Ra	42 ± 7	30 ± 6	23 ± 5	33 ± 4
^230 Th	18 ± 5	19 ± 5	17 ± 4	18 ± 3
^218Po	20 ± 5	24 ± 5	21 ± 5	21 ± 2
^210 Pb	100 ± 11	250 ± 17	100 ± 12	150 ± 8

The energy spectrum reported in Fig. 7 was reconstructed using a Monte Carlo simulation (not shown), obtaining activities compatible with those reported in Table 3.

The position and shape of the peaks suggested that most of the contaminants are uniformly diffused in the crystal bulk, but we could not completely exclude the hypothesis that part of the contaminations is located deep in the crystal surface. For simplicity, however, the activities reported in Table 3 were normalized to the mass of the crystals.

The only isotope showing a clear evidence for a surface contamination was Po that, because of its rather short half-life (138 days) was likely produced by Pb. As shown in Fig. 7, Pb produced two structures at 5.3 MeV () and 5.4 MeV ( + nuclear recoil). The presence of the first peak could be explained by the simulation only assuming a shallow contamination of the crystal surface or of the inert materials surrounding the detector. The activity of Pb shown in Table 3 was computed attributing the entire event rate (surface + bulk) to a contamination of the crystal bulk. The presence of this contaminant, however, is not of concern for CUPID-0, as none of the daughters of Pb/Po produces dangerous events.

In contrast to Po, the signature of isotopes belonging to U and Th chains could be explained both by crystal bulk contaminations, and by deep contaminations of the crystal surface (> 0.1m). If the contaminations were located in the crystal surface, we would expect a background reduction after the surface processing. The ZnSe crystals, indeed, will be treated with a cleaning procedure similar to the one developed for TeO crystals by the CUORE collaboration [40], that proved to be very effective.

Finally, the GEANT-4 simulation allowed to infer the expected background contribution to the region of interest for CUPID-0. We implemented the geometry described in Fig. 1, and assumed a discrimination power of 12 (see Sect. 6.2), and an energy resolution of 30 keV at the Q-value of (see Sect. 4).

If all the contaminants were diffused in the crystal bulk, we would expect a background at the Q-value of 410 counts/keV/kg/year. Since 0DBD electrons are expected to be contained in a single crystal (with a containment efficiency of 80 ) while the events produced by Tl decay mainly occur in more than one detector, we can exploit coincidences among nearby crystals to suppress this background to the level of 2.310 counts/keV/kg/year. Taking advantage of the small time delay between Tl and its mother (Bi), we can reach a level of 110 counts/keV/kg/y with a negligible dead-time (2 ) [15, 19], in compliance with the requirements of CUPID-0.

If all the contaminants were located only on the crystal surface, and in the unlikely event that crystal polishing does not change crystal surface activity, the background in the region of interest would be of 610 counts/keV/kg/y (coincidence suppression included), dominated by Tl interactions. In this case, the background suppression obtained exploiting the time delay between Tl and Bi would be slightly less effective, and would allow to reach the level of 310 counts/keV/kg/year. However, following what observed in CUORE [53, 54], crystal surface treatment allows a reduction of surface activity by more than a factor 6. In this case the expected counting rate could be even lower than 10 counts/keV/kg/year.

Conclusions and perspectives

The results presented in this paper allowed to assess the performance of the first ZnSe array in view of CUPID-0. We demonstrated that the assembly line guarantees the reproducibility both of light detectors, characterized by an unprecedented sensitivity, and of ZnSe bolometers. We derived the energy resolution of ZnSe detectors, obtaining 30 keV FWHM at the Q-value of 0DBD, and we proved that the excellent performance of light detectors allows to completely disentangle and reject the background due to interactions.

CUPID-0 is expected to run at least for 1 y of live-time to prove the potential of the scintillating bolometers technology, and the high number of emitters will allow to reach a remarkable sensitivity on the 0DBD. Thanks to the low crystal contaminations, the expected background should be lower than 1.510 counts/keV/kg/year, mainly produced by the cryogenic setup [34]. As explained in Ref [55], when the background is low enough to produce a number of events in the ROI of the order of unity along the experiment life (“zero background” approximation), the sensitivity scales as:where C.L. is the confidence level, the Avogadro constant, a the isotopic abundance, the stoichiometric coefficient of the 0DBD candidate, W the atomic weight. is the efficiency (80  for these ZnSe), M the detector active mass, t the live-time and the fraction of signal event that fall into the considered energy region (0.76 when E is the detector FWHM energy resolution).

Given a background index of 1.510 counts/keV/kg/year, we expect 0.6 events in 1 year in a FWHM energy interval centered around the Q-value. With a lower cryostat temperature, we expect CUPID-0 to feature a better energy resolution and, as a consequence, an even lower number of background events in the region of interest. Thus, the sensitivity can be computed in the “zero background” approximation and results 9.310 y at 90  C.L. in 1 year of data taking. The sensitivity will increase linearly with the live-time of the experiment as long as the conditions will be compatible with the zero background approach.