Development of 100 Mo -containing scintillating bolometers for a high-sensitivity neutrinoless double-beta decay search.

This paper reports on the development of a technology involving 100 Mo -enriched scintillating bolometers, compatible with the goals of CUPID, a proposed next-generation bolometric experiment to search for neutrinoless double-beta decay. Large mass ( ∼ 1 kg ), high optical quality, radiopure 100 Mo -containing zinc and lithium molybdate crystals have been produced and used to develop high performance single detector modules based on 0.2-0.4 kg scintillating bolometers. In particular, the energy resolution of the lithium molybdate detectors near the Q-value of the double-beta transition of 100 Mo (3034 keV) is 4-6 keV FWHM. The rejection of the α -induced dominant background above 2.6 MeV is better than 8 σ . Less than 10 μ Bq/kg activity of 232 Th ( 228 Th ) and 226 Ra in the crystals is ensured by boule recrystallization. The potential of 100 Mo -enriched scintillating bolometers to perform high sensitivity double-beta decay searches has been demonstrated with only 10 kg × d exposure: the two neutrino double-beta decay half-life of 100 Mo has been measured with the up-to-date highest accuracy as T 1 / 2 = [6.90 ± 0.15(stat.) ± 0.37(syst.)] × 10 18 years . Both crystallization and detector technologies favor lithium molybdate, which has been selected for the ongoing construction of the CUPID-0/Mo demonstrator, containing several kg of 100 Mo .

Introduction

Neutrinoless double-beta () decay, a yet-to-be-observed nuclear transition, consists in the transformation of an even-even nucleus into a lighter isobar containing two more protons with emission of two electrons and no other particles, resulting in a violation of the total lepton number by two units: (e.g. see Ref. [1]). This hypothetical transition is energetically allowed for 35 nuclei [2]. The detection of decay would have profound implications for our understanding of nature, proving that neutrinos are their own antiparticles (Majorana fermions), fixing the absolute neutrino mass scale and offering also a clue for the creation of matter abundance in the primordial universe (see recent reviews [1, 3] and references therein). It is to remark that this process is much more than a neutrino-physics experiment, because decay is a powerful, inclusive test of lepton number violation. Non-conservation of the total lepton number is as important as baryon number violation and is naturally incorporated by many theories beyond the Standard Model (SM). The current most stringent lower limits on the decay half-lives are in the range of 10– [1, 4]. The SM allowed process two neutrino double-beta () decay is the rarest observed nuclear transition and it has been measured in 11 nuclides with the half-lives in the range of – [5].

There are a number of proposed next-generation decay experiments, based on upgrades of the most promising current technologies (e.g. see Refs. [1, 6–9]). The goal of these future searches is to improve by up to two orders of magnitude the present best limits on the half-life with a sensitivity to the effective Majorana neutrino mass (a measure of the absolute neutrino mass scale) at the level of 10–20 meV, covering the so-called inverted hierarchy region of the neutrino mass pattern. The bolometric approach is amongst the most powerful methods to investigate decay. Particularly, one of the most stringent constrains on the effective Majorana neutrino mass [1] have been set by the results of CUORICINO and CUORE-0, precursors of the Cryogenic Underground Observatory for Rare Events (CUORE) [10], which studies the candidate isotope with the help of bolometers. CUORE, a ton-scale decay experiment, is now taking data in the Gran Sasso National Laboratories (Italy) and will be in operation for several years. A large group of interest is proposing a next-generation bolometric experiment, CUORE Upgrade with Particle ID (CUPID) [11, 12], to follow CUORE after the completion of its physics program. The nuclei , , and are the candidates considered for CUPID. A selection of the CUPID technologies, isotopes and materials is foreseen in 2018/2019.

Scintillating bolometers, the devices used in the present work, are favorable nuclear detectors for the conduction of sensitive decay searches, as they offer high detection efficiency, excellent energy resolution (at the level of ), efficient /() particle separation and potentially low intrinsic background [9, 13–18]. The isotope is one of the most promising candidates, since its signal is expected at [19] (Q-value of the transition), while the environmental background mainly ends at 2615 keV. The candidate is embedded in zinc and lithium molybdate crystals ( and ), working both as low-temperature bolometers and scintillators. An auxiliary bolometer, consisting of a thin Ge wafer, faces each -containing crystal in order to detect the scintillation light. The energy region above is dominated by events produced by radioactive contamination of surfaces, especially particles (e.g. as shown by the CUORICINO [20] and CUORE-0 [21]). Scintillation light yield from alpha interactions is usually quenched when compared to the () interactions of the same energy [22]. Combined with the fact that the thermal response for and () interactions are nearly equivalent, this allows for dual channel scintillating bolometer readouts to perform an effective event-by-event active background rejection [7–9, 17].

The development of a reproducible crystallization and detector technologies is needed for the scintillating bolometer technique to be applicable to a large-scale experiment, like CUPID. The specific requirements to be fulfilled by a crystallization technology of -containing scintillators are [23]: large enough crystal boule size; limited losses of the high-cost enriched isotope in the purification-crystallization chain; good optical properties; high scintillation yield; exceptionally low radioactive contamination. The size of a boule should be enough to produce at least one – scintillation element. The volume of the -containing crystal is bounded to the aformentioned value in order to avoid a significant impact on background from random coincidences of the decay events of [24-26]. Irrecoverable losses of the enriched material are acceptable at the level of a few % taking into account that the price of the enriched isotope is 80 $/g [27]. High transmittance (no less than 30 cm absorption length at the emission maximum) is welcome to reduce the amount of the trapped light and therefore to improve the scintillation light yield [28]. and crystals have a reasonable scintillation yield, at the level of 1 keV/MeV, which do not require ultra-low-noise bolometric light detectors. (Baseline noise at the level of a few hundreds of eV are sufficient to provide efficient light-assisted particle identification) According to Monte Carlo simulations of experiments based on -containing scintillating bolometers [7–9, 29, 30], a crystal bulk contamination of the order of 0.01 mBq/kg of would result to a minor contribution to the background in the region of interest (ROI; e.g. FWHM wide centered at ), at the level of [9]. As far as is concerned, a specific activity of even an order of magnitude higher would provide the significantly lower contribution to the background (e.g. see [29, 30]). The total activity of other radionuclides from the U/Th chains should not be higher than few mBq/kg to avoid pile-up effects. The main demands concerning the detector performance at the ROI are [8, 9, 29]: better than 10 keV FWHM energy resolution (5 keV FWHM is the CUPID goal [11]); at least 99.9% rejection of -induced events (with ()s acceptance larger than 90%) to suppress this background component to less than .

Preliminary results have been achieved in the past with bolometers containing molybdenum with natural isotopic composition in the Gran Sasso underground laboratory in Italy [8, 31–34], in the Modane underground laboratory in France [29, 35–37] and in an aboveground cryogenic laboratory located at CSNSM (Orsay, France) [23, 28, 38–41]. In the latter set-up, the first small -enriched two-detector array has been tested recently [42]. Most of these R&D activities were conducted in the framework of the scintillating-bolometer research programs of LUCIFER [43] – focused on ZnSe for the decay candidate but involving also -containing scintillators – and of LUMINEU [44], dedicated to the investigation of .

The present work represents a crucial step forward in the development of radiopure scintillating bolometers based on and crystals grown from -enriched molybdenum. A protocol for crystal growth was developed, and several prototypes were tested showing excellent energy resolution, efficient background rejection power and remarkable radiopurity. The results described here prove in particular that the technology is mature enough to carry out a pilot experiment on a several-kilogram scale. This technology demonstrator will provide essential information for the choice of the CUPID technique by clarifying the merits and the drawbacks of the option.

R&D on natural and -enriched zinc and lithium molybdates

Important milestones were achieved by LUMINEU in the R&D on zinc molybdate scintillators: the development of a molybdenum purification procedure [23]; the growth of large () [35] and small (0.17 kg) [42] crystals with the help of the low-temperature-gradient Czochralski (LTG Cz) method [45, 46]; further optimization of the growth process [40]. The R&D goal has been accomplished by the successful development of a large-mass crystal boule ( in weight, enrichment is ) shown in Fig. 1 (top left). Even though there is still room to improve the crystal quality – the boule exhibits a faceted structure and contains inclusions, mainly in the bottom part – the developed crystallization technology ensures the growth of reasonably good quality scintillators with a mass of about – which represents more than 80% yield from the initial charge of the powder in the crucible – and below 4% irrecoverable losses of the enriched material.

Fig. 1

Photographs of the first large-mass -enriched scintillators: the crystal boule with the cut scintillation element enrZMO-t (top panels), and the boule of crystal with the produced sample enrLMO-t (bottom panels). Both scintillation elements were cut from the top part of the boules. Color and transparency of the enrZMO-t crystal are different from the ones of the boule due to artificial light source and grinded side surface. The photo on the top left panel is reprinted from [37]

In the present work, we report about the study of four massive (0.3–) crystals operated as scintillating bolometers at . Two scintillation elements have been cut from a boule containing molybdenum of natural isotopic composition [35], while the other two are obtained from the boule (Fig. 1, top left). The information about the applied molybdenum purification, the size and the mass of the produced samples are listed in Table 1. The size of crystals is chosen to minimize the material losses and to produce similar-size samples from each boule. According to [47], a hexagonal shape of elements (e.g. see Fig. 1, top right) should provide a higher light output than the cylindrical one.

Table 1

Zinc and lithium molybdate crystal scintillators grown by the LTG Cz method from molybdenum with natural isotopic composition and enriched in . The molybdenum compound has been purified by single or double sublimation with subsequent double recrystallization in aqueous solutions. A compound supplied by NRMP (see text) was used to produce all Li-containing scintillators, except LMO-3 (produced from Alfa Aesar ). The position in the crystal boule is given for those samples cut from the same boule. The crystal ID is represented by the abbreviation of the chemical compound with an extra “enr” to mark enriched samples and/or “t” or “b” to indicate the position in the boule and a number to distinguish boules of the same material

Scintillator	Molybdenum sublimation	Boule crystallization	Crystal ID	Position in boule	Size (⊘×h mm)	Mass (g)
ZnMoO_4	Single	Double	ZMO-t	Top	50×40	336
			ZMO-b	Bottom	50×40	334
Zn^100MoO_4	Double	Single	enrZMO-t	Top	60×40	379
			enrZMO-b	Bottom	60×40	382
Li_2MoO_4	Single	Single	LMO-1	–	40×40	151
	Single	Double	LMO-2	–	50×40	241
	Single	Single	LMO-3	–	50×40	242
Li2100MoO_4	Double	Triple	enrLMO-t	Top	44×40	186
			enrLMO-b	Bottom	44×44	204
	Double	Double	enrLMO-2t	Top	44×46	213
			enrLMO-2b	Bottom	44×44	207

Because of some experienced difficulties with the crystallization process,1 which prevented us from obtaining top quality large-mass crystals, the LUMINEU collaboration initiated an R&D on the production of large-mass radiopure lithium molybdate scintillators [41, 48]. Thanks to the low and congruent melting point of , the growth process is expected to be comparatively easier than that of . However, the chemical affinity of lithium and potassium results in a considerably high contamination of () in crystal scintillators, as it was observed in early studies of this material [49]. Despite the low of , random coincidences of and events of Mo can produce background in the ROI [25]. In particular, a contamination level around of in a detector with dimension provides the same background counting rate in the ROI as the random coincidences of the events. So, in addition to the LUMINEU specifications on U/Th contamination, the acceptable activity in crystals is of the order of a few mBq/kg. Therefore, the R&D on scintillators included the radioactive screening and selection of commercial lithium carbonate samples, the optimization of the LTG Cz crystal growth and the investigation of the segregation of radioactive elements in the crystallization process.

Three samples of high purity lithium carbonate were measured by HPGe spectrometry at the STELLA facility of the Gran Sasso underground laboratory (Italy): (1) 99.99% purity grade powder produced by Novosibirsk Rare Metal Plant (NRMP, Novosibirsk, Russia) [50]; (2) 99.998% lithium carbonate by Puratronic (Alfa Aesar GmbH & Co KG, France) [51]; (3) 99.99% raw material by Sigma-Aldrich (USA) [52]. The results are given in Table 2. The lithium carbonate produced by NRMP, the material of highest radiopurity, was selected for crystals production. Due to the high contamination, the Sigma-Aldrich material was rejected for further investigation.

Table 2

Radioactive contamination of commercial lithium carbonate compounds measured by low background HPGe spectrometry. Errors are given in parenthesis at 68% C.L., upper limits – at 95% C.L.

Chain	Nuclide	Activity in Li_2CO_3 powder (mBq/kg)
		NRMP	Alfa Aesar	Sigma-Aldrich
^232Th	^228Ra	≤2.9	≤14	16(8)
	^228Th	≤3.7	12(4)	13(4)
^238U	^226Ra	≤3.3	705(30)	53(6)
	^40K	≤42	≤42	210(70)

Even first attempts of the growth by the LTG Cz technique were successful providing high quality crystal boules with masses of 0.1–0.4 kg [41]. The growing conditions have been optimized to extend the crystal size up to 100 mm in length and 55 mm in diameter (0.5–0.6 kg mass) [53] allowing us to produce two large scintillating elements of about 0.2 kg each from one boule. For the present study, we developed three large scintillators by using highly purified molybdenum oxide and high purity grade lithium carbonate. Two of them have been grown from the NRMP compound by applying a single (LMO-1 sample in Table 1) and a double (LMO-2) crystallization, while the last one (LMO-3) was grown by the single crystallization from the Alfa Aesar Li-containing powder.

Once the LTG Cz growth of crystals containing molybdenum of natural isotopic composition was established, we started to process molybdenum enriched in . Figure 1 (bottom left) shows a first large-mass () -enriched crystal boule grown at the beginning of 2016. The crystal was produced by a triple crystallization due to an accident that happened during the second crystal growth process. The second massive crystal boule (; see Fig. 4 in [53]) was grown by double crystallization at the end of May 2016. Both enriched crystals demonstrated high optical quality and have the size required for the production of two similar transparent scintillation elements with masses of 0.2 kg (see Table 1 and Fig. 1, bottom right). Two cylindrical samples produced from the first crystal boule were used for the bolometric tests described in the present work.2

Underground tests of -containing scintillating bolometers

-containing scintillating bolometers

The bolometers were fabricated from the crystal scintillators listed in Table 1. Each scintillating crystal was equipped with one or two epoxy-glued Neutron Transmutation Doped (NTD) Ge temperature sensors [55], whose resistance exponentially depends on temperature as . and are two parameters depending on the doping, the compensation level and on the geometry in the case of . In our samples, is derived to be 0.5. In the present work we used high resistance (HR) and low resistance (LR) sensors with typical parameter values , and , , respectively. Therefore, HR NTDs have a resistance of at working temperature, while an order of magnitude lower resistance is typical for LR NTDs. The NTD Ge thermistors, biased with a constant current, act as temperature-voltage transducers. The thermal link to the bath was provided by Au bonding wires which give also the electrical connection with the NTD Ge sensors. In addition, each crystal was supplied with a small heater made of a heavily-doped Si [56], through which a constant Joule power can be periodically injected by a pulser system to stabilize the bolometer response over temperature fluctuations [56, 57].

The detectors were assembled according to either LUMINEU or LUCIFER standard schemes (see Table 3). The mechanical structure and the optical coupling to the crystal scintillators are designed to optimize the heat flow through the sensors and to maximize the light collection. The standard adopted by LUMINEU for the EDELWEISS-III set-up implies the use of a dedicated copper holder where the crystal scintillator is fixed by means of L- and S-shaped PTFE clamps [35-37]. The holder is completely covered internally by a reflector to improve the scintillation-light collection. For the prototype of the LUMINEU suspended tower, shown in Fig. 2, the holders were slightly modified to make the array structure able to pass through the holes in the copper plates of the EDELWEISS set-up. In case of the LUCIFER R&D standard, the crystal is fixed to a copper frame by S-shaped PTFE pieces and copper columns, as well as side-surrounded by a plastic reflective film (e.g. see in [8]). This frame is thermally anchored to the mixing chamber of the dilution refrigerator.

Table 3

The main construction elements of -containing heat detectors studied in the present work. Their IDs coincide with the scintillation crystal IDs defined above. Three types of reflectors were used: Radiant Mirror Film (RMF) VM2000/VM2002 and Enhanced Specular Reflector (ESR) film by 3M™and a thin silver layer (Ag) deposited on the holder. The masses of all used NTD sensors are . The enrZMO-t, enrZMO-b, and LMO-2 detectors were also equipped with a smeared source

Standard	Heat detector ID	Support		Reflector	NTD sensor type
		Copper	PTFE		No. 1	No. 2
LUMINEU	ZMO-t	Holder	L- and S-shaped	RMF	HR	HR
	ZMO-b				HR	LR
	LMO-1				LR	–
	LMO-3			Ag	LR	–
LUMINEU (tower)	ZMO-b	Holder	L- and S-shaped	ESR	LR	–
	enrLMO-t				LR	–
LUCIFER	enrZMO-t	Plate, columns	S-shaped	ESR	LR	LR
	enrZMO-b				LR	LR
	LMO-2				LR	LR
	enrLMO-b				HR	HR

Fig. 2

Photographs of a three-spring suspended tower (first column) and two LUMINEU scintillating bolometers (second column): the (top photo; ZMO-b, ) and the (bottom photo; enrLMO-t, ) bolometers together with two identical Ge light detectors (third column). A third detector of the tower (bottom in the left photo, not shown in details) is a 0.2 kg Ge bolometer

Thin bolometric light detectors (see Table 4) were coupled to the scintillating crystals to register the scintillation light. All of them are based on high purity Ge wafers and their typical size is 44–45 mm in diameter and 0.17–0.3 mm thickness, but two detectors have slightly lower area and tens thickness. Some light detectors were constructed according to the LUMINEU standard described in [58], with the additional deposition of a 70 nm SiO antireflecting coating on one surface of the Ge wafer to increase the light absorption [59]. Another type of light detectors used in the present study was developed by the LUCIFER group [60]. One bolometer was assembled according to CUPID-0 mounting standard [61]. In all these cases, the Ge wafer is held by PTFE clamps. The last type of used light detectors is the state-of-the-art optical bolometer developed at IAS (Orsay, France) [62]. The suspension of the Ge wafer is carried out by Nb–Ti wires in this case. All the light detectors were equipped with one NTD Ge thermistor.

Table 4

Information about Ge light detectors used in the present work. The detectors are grouped according to the mounting standard given in the first column

Standard	Light detector ID	Ge size (mm)	Coating	NTD sensor
				Type	Mass (mg)
LUMINEU	M1	⊘44×0.17	Yes	LR	20
	M3	⊘44×0.17	Yes	LR	9
	Lum11	⊘44×0.17	Yes	HR	5
	Lum12	⊘44×0.17	Yes	HR	5
LUCIFER	GeB	⊘45×0.30	No	LR	9
	GeT	⊘45×0.30	No	LR	9
CUPID-0	GeOld	⊘45×0.30	No	LR	9
IAS	B297	⊘40×0.045	No	LR	1
	B304	⊘25×0.030	No	LR	1

Underground cryogenic facilities

In the present investigations, we used two cryogenic set-ups: CUPID R&D and EDELWEISS-III located at Gran Sasso National Laboratories (LNGS, Italy) and Modane underground laboratory (LSM, France), respectively. The general description of these facilities is given in Table 5. Some features are related to the specific applications: the CUPID R&D is mainly oriented on the R&D of bolometers (including scintillating bolometers) for searches, with ROI at a few MeV, while the EDELWEISS-III set-up was conceived to perform direct dark-matter searches with the help of massive heat-ionization bolometers, with a ROI in the tens-of-keV range.

Table 5

A short description of the used underground cryogenic set-ups. The rock overburden is expressed in km of water equivalent (km w.e.). The base temperature indicates the minimal temperature of the cryostat. The sampling rate is given in kilo-samples per sec (kSPS )

	CUPID R&D [63–65]	EDELWEISS-III [35, 66, 67]
Location
Underground lab	LNGS (Italy)	LSM (France)
Rock overburden (km w.e.)	3.6	4.8
Cryostat
Dilution refrigerator	^3He/^4He	^3He/^4He
Type	Wet	Wet and dry
Geometry	Standard	Reversed
Experimental volume (L)	∼8	∼50
Outside mechanical decoupling	No	Yes
Inside mechanical decoupling	Yes	Yes (since 2016)
Base temperature (mK)	7	10
Shield (external)
Low activity lead (cm)	20	18
Roman lead (cm)	No	2
Polyethylene (cm)	8	55
Boron carbide (cm)	1	No
Anti-radon box	Yes	No
Muon veto	No	Yes
Shield (internal)
Roman lead (cm)	5.5	14
Polyethylene (cm)	No	10
Readout and DAQ
Electronics	Cold + Room-temp.	Cold
Dual readout channels	10 + 8	48
Bias	DC	AC
ADC digitization (bit)	18	16 or 14
Sampling rate (kSPS)	Up to 250/N_channels	Up to 1
Data taking mode	Trigger and/or stream	Trigger or stream
Calibration
Regular	^232Th	^133Ba
Exceptional	^40K, ^137Cs, AmBe	^232Th, ^40K, AmBe
In-situ source	Allowed	Prohibited
Pulser system	Yes	Yes (since 2015)

As one can see from Table 5, an efficient suppression of the cosmic-ray flux is provided by a deep underground location of both set-ups. The EDELWEISS-III is larger and can host up to 48 scintillating bolometers with a copper holder size of each. The reversed geometry of the EDELWEISS-III cryostat does not allow to decouple mechanically the detectors plate from the mixing chamber, as it was done by two-stage damping system inside the CUPID R&D set-up [63]. The external damping system (pneumatic dampers) of the EDELWEISS-III is adapted to the operation of tightly held massive EDELWEISS detectors, and not to scintillating bolometers. In particular, thin light detectors are very sensitive to the vibrations induced by the three thermal machines of the set-up. Therefore, an internal damping inside the EDELWEISS-III has been implemented through a mechanically-isolated suspended tower (see Fig. 2). Dilution refrigerators of both set-ups are able to reach a base temperature around 10 mK.

The EDELWEISS-III set-up is surrounded by a significantly massive passive shield against gamma and neutron background. The absence of an anti-radon system as that used in the CUPID R&D is somehow compensated by a deradonized (below ) air flow. The radon level is monitored continuously. An important advantage of the EDELWEISS-III set-up is a muon veto system with about 98% coverage (however, no clock synchronization with scintillating bolometers has been implemented yet).

All the EDELWEISS-III readout channels utilize a cold electronics stage, while only about half of those of CUPID R&D have this feature. The EDELWEISS-III readout system uses AC bolometer bias modulated at a frequency of up to 1 kHz, which is also kept for the demodulation procedure applied to the data sampled with a 100 kSPS rate (the modulated data can be also saved). Higher resolution without a significant enlargement of the data size is available for the CUPID R&D case, which envisages DC bolometer bias. In contrast to DC current, there are difficulties in the operation of high resistance NTDs with AC bias (e.g. unbalanced compensation of nonlinearities related to the differentiated triangular wave applied for NTD excitation – see details in [67, 68]).

An important difference between the set-ups is in the calibration procedure and the related policy. The CUPID R&D is well suited for a regular control of the detector’s energy scale in a wide energy range up to 2.6 MeV. On the contrary, a periodical calibration with the EDELWEISS-III set-up is available only with a source (’s with energies below 0.4 MeV), while the insertion of a source, as well as a few other available sources, requires lead/polyethylene shield opening, which is not supposed to be done frequently. Also, there is prohibition of the in-situ use of sources for the light detectors calibration, which is not the case for the CUPID R&D case. Finally, the control of the detector thermal response by a pulser system connected to the heaters is available for both set-ups.

Low-background measurements and data analysis

The list of the low-background bolometric experiments and the main technical details are given in Table 6. The bias currents of the order of a few nA were set to maximize the signal-to-noise ratio resulting in a working-point thermistor resistance of a few .

Table 6

General information about measurements with -containing bolometers operated at Modane and Gran Sasso underground laboratories. IDs of detectors used for the construction of double read-out hybrid bolometers correspond to the heat and light detectors ID defined above. denotes the base temperature of the cryostat

Set-up	Run ID	Detectors		Sampling (kSPS)		T_base (mK)	Data taking (h)
		Heat	Light	Heat	Light
EDELWEISS-III	Run308	ZMO-t	M1	1	1	18	5000
		ZMO-b	M3
	Run309	ZMO-b	M3	1	1	20	2966
	Run310	ZMO-b	Lum12	1	1	19–20	2090
		enrLMO-t	Lum11
CUPID R&D		LMO-1	B297	2	2	15	328
		LMO-2	GeB	1	1	19	201
		LMO-3	B304
		enrZMO-t	GeB	1	1	15	717
		enrZMO-b	GeT
		enrLMO-b	GeOld	2	4	12	487

An optimum filter technique [69, 70] was used to evaluate the pulse height and shape parameters. It relies on the knowledge of the signal template and noise power spectrum; both are extracted from the data by averaging about 2 MeV energy signals (40 individual pulses) and baseline waveforms (5000 samples), respectively. For those detectors that were equipped with two temperature sensors, the data of the thermistor with the best signal-to-noise ratio were analyzed. The light-detector signal amplitude is estimated at a fixed time delay with respect to the heat signals as described in [71]. Due to spontaneous temperature drifts, the amplitudes of the filtered signals from the crystal scintillator are corrected for the shift in thermal gain by using the heater pulses.3

The heat response of the scintillating bolometers is calibrated with quanta of (238.6, 338.3, 510.8, 583.2, 911.2, and 2614.5 keV), (1460.8 keV) and/or (356.0 keV) sources. The light detectors in the CUPID R&D set-up were calibrated with the X-ray doublet (5.9 and 6.5 keV) of the source.

Performance of -containing scintillating bolometers

Time profile of the pulses

The rising edge of the bolometric signal depends on the sensitivity of the sensor to athermal and/or thermal phonons created by a particle interaction and has a characteristic time ranging from microseconds (dominant athermal component) to milliseconds (dominant thermal component). Since the NTD sensors are sensitive mainly to thermal phonons, the rise time of the tested detectors given in Table 7 is within the expectation. The heat detectors have longer leading edge (tens of ms) than that of light detectors (few ms) due to the larger volume and therefore to the larger heat capacity of the absorber.

Table 7

Characteristics of -containing scintillating bolometers. The pulse-shape time constants are the rise () and decay () times defined as the time difference between the 10% and the 90% of the maximum amplitude on the leading edge and the time difference between the 90% and the 30% of the maximum amplitude on the trailing edge, respectively. The signal sensitivity is measured as the thermistor voltage change for a unitary energy deposition. The intrinsic energy resolution (FWHM baseline) is determined by noise fluctuations at the optimum filter output. The energy resolution (FWHM) of light detectors was measured with a X-ray source. The FWHM resolution of heat channels is obtained for quanta of , , and sources. and denote light yields for and , respectively. The quenching factor for particles and the discrimination power (above 2.5 MeV) are calculated according to the formulas given in the text

Scintillating bolometer	Scintillator	ZnMoO_4		Zn^100MoO_4		Li_2MoO_4			Li2100MoO_4
	Crystal ID	ZMO-b		enrZMO-t	enrZMO-b	LMO-1	LMO-2	LMO-3	enrLMO-t	enrLMO-b
	Size (⊘×h mm)	50×40		60×40	60×40	40×40	50×40	50×40	44×40	44×44
	Mass (g)	334		379	382	151	241	242	186	204
	Light detector ID	M3	Lum12	GeB	GeT	B297	GeB	B304	Lum11	GeOld
	Size (⊘×h mm)	44×0.17	44×0.17	45×0.30	45×0.30	40×0.045	45×0.30	25×0.030	44×0.17	45×0.30
Test	Underground lab	LSM	LSM	LNGS	LNGS	LNGS	LNGS	LNGS	LSM	LNGS
	T_base (mK)	18	19–20	15	15	15	19	19	19–20	12
Pulse-shape time constant (ms)	Light τ_R	5.2	4.6	2.4	2.9	2.3	4.0	2.7	3.5	5.1
	Light τ_D	23	24	12	14	2.6	8.5	6.4	13	13
	Heat τ_R	19	38	9.6	6.4	17	29	30	27	18
	Heat τ_D	200	204	37	18	67	339	414	169	88
Sensitivity (nV/keV)	Light detector	2200^a	2500^a	1047	1053	4030	850	15800	1900^a	2910
	Heat detector	48	26	73	39	166	11	23	32	89
Light FWHM (eV)	Baseline	∼140^a	∼60^a	∼490	∼230	∼42	∼420	∼18	∼70^a	∼140
	X-ray ^55Mn, 5.9 keV	–	–	787(3)	289(1)	334(4)	555(5)	504(4)	–	303(2)
Heat FWHM (keV)	Baseline	∼1.6	∼3.6	∼2.6	∼4.3	∼0.6	∼2.3	∼1.6	∼1.2	∼1.2
	γ ^133Ba, 356 keV	3.7(1)	5.1(1)	–	–	–	–	–	2.54(4)	–
	γ ^228Ac, 911 keV	4.4(7)	10(1)	5.6(7)	9(2)	2.0(3)	3.9(6)	3.1(6)	3.1(5)	3.1(2)
	γ ^40K, 1461 keV	–	7.9(2)	6.7(6)	14(1)	–	4.2(3)	4.4(3)	4.1(2)	–
	γ ^208Tl, 2615 keV	9(1)	12(1)	9.1(7)	22(2)	3.8(6)	6(1)	4.7(7)	6.3(6)	5.0(5)
	α ^210Po, 5407 keV	8.8(1)	9.0(2)	∼47	∼100	7(2)	9(1)	9(2)	5.4(3)	–
Response to γ(β) and α	LY_γ(β) (keV/MeV)	–	–	1.32(1)	1.20(1)	0.68(4)	0.99(1)	0.121(2)^b	–	0.775(4)
	LY_α (keV/MeV)	–	–	0.217(4)	0.148(2)	0.165(1)	0.203(4)	0.0236(3)	–	0.153(2)
	QF_α	0.15	0.17	0.17	0.13	0.23	0.20	0.17	0.22	0.19
	DP_α/γ(β)	12	21	7.8	11	16	8.7	11	18	12

Estimations are based on rough calibrations by scintillation light (see Sect. 4.2)

Low light yield is caused by non-optimal light collection conditions of the measurements (see Sect. 4.4)

The decaying edge time constant of the bolometric signal represents the thermal relaxation time, which is defined by the ratio of the heat capacity of the absorber to the thermal conductance to the heat bath. Therefore, it strongly depends on the material, the detector coupling to the heat bath and on the temperature. As one can see in Table 7, the variation of the decay time is even larger than that of the rise time, but again it has the typical values normally observed in light detectors (tens of ms) and massive bolometers (hundreds of ms). The improved coupling of the NTD sensor to the heat bath of B297 and B304 light detectors [62] leads to shorter decay time at the level of a few ms.

The only exception in the signal time constants of massive detectors is evident for both crystals, the largest of all tested samples, which exhibit signals faster by about a factor 2 than the other devices, in particular those tested in the same set-up and at similar temperatures. It is interesting also to note that the bottom crystal is twice faster than the top one, as it was observed also in the test of the detectors [42]. The fast response of the enriched bolometers has no clear explanation, but it is probably related to crystal quality. However, a fast detector response is crucial for a separation of the events pile-ups [24-26]. Thus, this feature of scintillators (e.g. below 10 ms) can lead to a better capability to discriminate random coincidences by heat pulse-shape analysis than that considered in Ref. [25].

Voltage sensitivity

In a bolometric detector, the response to a nuclear event is a temperature rise directly proportional to the deposited energy and inversely proportional to the detector heat capacity. A thermistor converts the temperature variations to a voltage output, digitized by the readout system. Therefore, a bolometric response is characterized by a voltage sensitivity per unit of the deposited energy.

A signal pulse height of the order of few tens to hundreds nV/keV is typical for NTD-instrumented massive bolometers. This figure corresponds to what is observed in all the tested crystals (see Table 7).

The reduced size of both the absorber and the sensor of light detectors (see Tables 3, 4) leads to lower heat capacities and therefore to higher sensitivites, which are in the range – for a good-performance detector. This is the case for all the tested light detectors4 (see in Table 7), except the ones with even smaller size (B297 and B304) and subsequently sensitivity enhanced by up to one order of magnitude.

Energy resolution

Most of the used light detectors have similarly good performance also in terms of energy resolution, in particular their baseline noise is –0.5 keV FWHM (see in Table 7). The only exceptions are detectors with enhanced sensitivity (B297 and B304) for which the baseline noise is below . However, they also exhibit a strong position-dependent response, therefore the energy resolution measured with an uncollimated source is near to that obtained with the other light detectors (FWHM 0.3–0.8 keV at 5.9 keV).

As it was mentioned above, better than 10 keV FWHM energy resolution at the ROI is one of the most crucial requirements for cryogenic double-beta decay detectors. This goal was successfully achieved with both natural and -enriched and based bolometers5 (see Table 7). Below we discuss the obtained results.

The detectors exhibit twice better energy resolution than the ones and the achieved values of 4–6 keV FWHM at 2615 keV are at the level of the best resolutions ever obtained with massive bolometers [9, 72, 73]. In particular, the energy resolution of bolometers is comparable to the performance of the cryogenic detectors of the CUORE-0 experiment (the effective mean FWHM at 2615 keV is 4.9 keV with a corresponding RMS of 2.9 keV [73]). This is mainly due to the fact that , as , demonstrates a low thermalization noise, i.e. a small deviation of the energy resolution from the baseline noise width. The results of the and detectors show possible improvement of the energy resolution by lowering of the temperature, as it is expected thanks to increased signal sensitivity. A dependence of the performance on the sample position in the boule, observed early with small [42] and now with large samples, is also evident. It could be related to the degradation of the crystal quality along the boule. Thanks to the higher crystal quality, no such effect is observed for crystals.

The energy spectra of a source measured by the -enriched bolometers (enrZMO-t and enrLMO-b) and the corresponding energy-dependance of the heat-channel resolution are illustrated in Fig. 3. The chosen data of and detectors represent the typical energy resolution for bolometers based on these materials in case of optimal experimental conditions (Table 7). Using the fitting parameters for the curves shown in Fig. 3 (right), the expected energy resolution of the enrZMO-t and enrLMO-b cryogenic detectors at of is and , respectively. Thus, the energy resolution of the detectors is acceptable but needs still an optimization, while bolometers already meet the resolution required for future generation bolometric experiments [7–9, 11].

Fig. 3

The energy spectra of the source measured by the (enrZMO-t; dashed histogram) and (enrLMO-b; solid histogram) bolometers over 64 h and 168 h, respectively, at LNGS (left figure). The energy bin is 5 keV. The 2615 keV peak of the quanta accumulated by the detectors is shown in the inset. The energy dependence of the energy resolution of the ZMO-t and enrLMO-b detectors (right figure). The fits to the data by a function FWHM ( and are free parameters) are shown by the dashed lines. The parameters of fits are and 0.0282(6) keV for the ZMO-t and and 0.0091(3) keV for the enrLMO-b. The dotted line indicates the value of (3034 keV)

Response to s and particle identification capability

Scintillation-assisted particle discrimination

By using coincidences between the heat and the light channels, one can plot a light-vs-heat scatter plot as the ones presented in Fig. 4. The heat channel of all data shown in Fig. 4 is calibrated by means of quanta of the calibration sources and it leads to 10% heat miscalibration for particles due to a so-called thermal quenching, common for scintillating bolometers (e.g. see results for different scintillators in [7, 34, 61, 74]). Therefore, in order to present the correct energy of the events, an additional calibration based on the peaks identification is needed.

Fig. 4

Scatter plots of light-versus-heat signals of the background data collected with scintillating bolometers based on (ZMO-b, top figures), (enrZMO-b, bottom left), and 186 g (enrLMO-t, bottom right) crystals over 2767, 1300, 593, and 1303 h, respectively. The detector was operated in the CUPID R&D cryostat, while the other data were accumulated in the EDELWEISS-III set-up (the light signals of the latter are in analog-to-digit units, ADU). The heat channels were calibrated with quanta. The () and events populations are distinguished in color by using the cuts on the heat energy and the light yield parameters (see the text). The particle identification capability of the detector affected by vibration noise (top left) was substantially improved in the suspended tower (top right). The features of the particle populations are discussed in the text

As it is seen in Fig. 4, the light-vs-heat scatter plot contains two separated populations: a band of ()’s and a distribution of events associated to decays. This is due to the fact that the amount of light emitted in an oxide scintillator by particles is quenched typically to 20% with respect to quanta ( particles) of the same energy (see, e.g., Ref. [22]). Therefore, the commonly used particle identification parameter for scintillating bolometers is the light yield (LY), that we will define as a ratio of the light-signal amplitude measured in keV to the heat-signal amplitude measured in MeV.

The data of all detectors with directly calibrated light channel (all measurements at LNGS) have been used to determine the and values of () and events selected in the heat-energy range 2.5–2.7 MeV and 2.5–7 MeV,6 respectively. In spite of the quite evident constancy of the in a wide energy range (as it is seen from the slop of s in Fig. 4), the event selection was applied above 2.5 MeV, because the same distributions have been used to calculate / discrimination power (see below) close to the ROI of . The LY values extracted from the present data are given in Table 7.

The light yields for events measured with both scintillating bolometers are in the range 1.2–1.3 keV/MeV, similar to the results of previous investigation of natural (see [23, 32] and references therein) and -enriched [42] detectors. Thanks to the progress in the development of high quality lithium molybdate scintillators – as documented in the present work and recently in Ref. [48] – the values for and scintillation bolometers, which lay in the range 0.7–1 keV/MeV, become comparable to the light yields of the detectors. The improvement of the LY with respect to the early investigations with detectors [33, 34] is of about a factor of 2. One bolometer (LMO-3) was viewed by a light detector with a significantly lower area implying a reduced light collection and consequently a rather small value of 0.12 keV/MeV.

The ratio of the LY parameters for s and ()s gives the quenching factor of the scintillation light signals for particles: . An absolute light detector calibration is not needed to calculate this parameter. As it is seen in Table 7, the results for ZnMoO and detectors are similar showing 20% quenching of the light emitted by particles with respect to the () induced scintillation.

The efficiency of discrimination between and populations can be characterized by the so-called discrimination power parameter defined as:where () denotes the average value (width) of the or () distribution. The value is estimated for () and events selected for the LY determination (see above).

As reported in Table 7, the achieved discrimination power for all the tested detectors is , which implies a high level of the /() separation: more than 99.9% rejection while preserving practically 100% signal selection efficiency. The separation efficiency is illustrated in Fig. 5 for the scintillating bolometer enrZMO-t with the lowest achieved due to the modest performance of the GeB light detector. It is to emphasize the 0.1 keV/MeV obtained with the LMO-3 detector, which would not allow effective particle identification by using a standard-performance light detector with 0.2–0.5 keV FWHM baseline noise.7 However, the performance of the B304 optical bolometer – which featured 0.02 keV FWHM baseline noise – was high enough to provide highly-efficient particle identification even with this detector ().

Fig. 5

The light yield distributions of particles from a smeared source and () events collected by a scintillating bolometer (enrZMO-t) over 593 h of background and 78 h of calibration measurements in the CUPID R&D set-up at LNGS (Italy). The distributions are fitted by Gaussian functions shown by solid lines. The corresponding discrimination power is . The intervals containing 99.9% of both event types and sigma interval of the band are also given

Peculiarities in particle identification

Figure 4 illustrates observed peculiarities of some detectors which could affect the particle identification capability. These peculiarities are originated either by a noise-affected detector performance or by a feature of the detector’s response to , which exhibits classes of events with more quenched or enhanced light signals. Below we will discuss briefly these observations and their impact on background in a decay experiment with .

High vibrational noise in a light detector affects the precision of the light-signal amplitude evaluation, especially for events with a low scintillation signal ( events and ()s below ). This was an issue of the measurements with detectors in Run308 and Run309, and this effect is apparent in Fig. 4 (top left). The problem can be solved by using a mechanically isolated system inside a cryostat (see Table 5 and e.g. Refs.  [63, 76, 77]). In particular, a stable and reliable light-channel performance of the scintillating bolometer in the suspended tower (Run310) is evident from Fig. 4 (top right).

The data of natural and -enriched bolometers contain some events that have more quenched light output and enhanced heat signals; e.g. see “dark and hot ” in Fig. 4. In the past, the same effect was observed in bolometric tests of small [31] and [78] crystals which also exhibit defects and macro inclusions. A major part of such events is distributed close to ( structures at around 6 MeV in electron-equivalent energy in Fig. 4), the main contamination of the investigated crystals. Only a short-range () interaction in the crystal bulk exhibits this anomaly, because it is not evident either for interactions at the crystal surface (energy-degraded events) or for () events which have longer mean path in the crystal than the bulk ’s. This phenomenon is probably related to the thermal quenching, as suggested by the pronounced anti-correlation between light and thermal signals in the response. The effect is more evident for the enriched crystals, which contain more inclusions than the natural ones: e.g. about 40% of events acquired by the enrZMO-b detector are attributed to the “dark and hot ”, while four times lower amount of such events is observed in the ZMO-b bolometer. Thereby, the origin of this anomaly in the response to interactions is probably related to the crystals imperfections. Taking into account that two electrons are expected in the of , we expect that the signal is unaffected by this anomaly. Furthermore, it does not affect the detector’s capability to identify and reject the -induced surface events, which constitute the most challenging background in a bolometric experiment without particle identification. Negative effects are only expected on the precision of the spectroscopy, which however is important not only to build a background model through radiopurity determination, but also for the off-line rejection of - delayed events from decays of - [7, 17, 32].

The light-vs-heat data of several scintillating bolometers contain also events with an enhanced light signal with respect to those of the prominent distribution. As it is seen in Fig. 4, these events belong to three families: BiPo events, surface events with escaped nuclear recoils hitting the light detector, and the so-called “bright ”.

The first family consists of unresolved coincidences in - decays (BiPo events in Fig. 4). Due to the slow bolometric response, the decays of () and Bi ( = 3272  keV) overlap with subsequent decays of (, ) and (, , respectively. Therefore, they are registered as a single event with a heat energy within 8–11 MeV range and a light signal higher than that of a pure event of the same energy. Since the BiPo events are distributed far away from 3 MeV, they have no impact on the ROI of .

In a case of an decay on a crystal surface, a nuclear recoil (or an particle) can escape from the scintillator and hit the light detector. Such events belong to the second family indicated in Fig. 4. Taking into account that only a few keV energy-degraded recoil can mimic a light signal of or bolometer, the heat energy release has to be close to the nominal -value additionally enhanced due to the thermal quenching. Therefore, independently on the surface activity of radionuclides from U/Th chains (4–9 MeV -values), they cannot populate the ROI of . Among other natural -active nuclides, a probable contaminant is ( [79]) due to the crystal growth in a platinum crucible. However, even in such case the expected heat signal is about 0.5 MeV away from the of , as well as the bulk contamination in the studied crystals is expected to be on the level of a few [35]. We can therefore conclude that also this class of events does not play a role in the search for the decay of .

The last family – consisting of “bright ” events in Fig. 4 – stem from the documented scintillation properties of the reflecting film. Specifically, an energy deposition in this film can take place for surface-originated decays, which can produce a heat and a light signal in the scintillating crystal but also a flash of scintillation light from the reflecting film, which adds up to that of the crystal scintillator. This results into an enhanced light signal. Consequently, the population of energy-degraded events can leak to the ROI of in the heat-light scatter plot, providing an unavoidable background. To check the scintillation response of the 3M film, we have performed a test using a photomultiplier and a source. The observed scintillation is at the level of 15–34% relatively to NE102A plastic scintillator (depending on the side of the film facing the photomultiplier). Therefore, such a feature of the reflector spoils the particle discrimination capability of the detector. In order to solve this issue, a reflecting material without scintillation properties has to be utilized or the reflecting film has to be omitted.8

Response to neutrons

The ZMO-b, LMO-1, and enrLMO-t detectors were also exposed to neutrons from an AmBe source. The results for Li-containing bolometers are illustrated in Figs. 6 and 7 (left). The () band exceeds the natural end-point because of the prompt de-excitation ’s following reaction. The cluster of events in the region is caused by the reaction (Q-value is 4784 keV [80, 81]). The has a natural abundance of 7.5% [82], and the large cross section for thermal neutrons ( barns [80, 81]) gives rise to the clear distribution at a heat energy of around 5 MeV. In the energy scale, the distribution is shifted by about with respect to the 4784 keV total kinetic energy released in the reaction. The energy resolution (FWHM) on the peak was measured as 7.7(3) and 5.9(2) keV for the LMO-1 and enrLMO-t detectors, respectively. This is an unprecedented result obtained with -containing detectors (e.g. compare with the results of Li-containing cryogenic detectors in Refs. [34, 83, 84] and references therein). A second structure at higher energy is attributed to the non-thermal neutrons, in particular to the resonant absorption of 240 keV neutrons. A linear fit to the less prominent lower band, ascribed to nuclear recoils induced by fast neutron scattering, gives a light yield of .

Fig. 6

The light-versus-heat data accumulated with the 186 g scintillating bolometer (enrLMO-t) in the EDELWEISS-III set-up (21–20 mK data) under neutron irradiation of an AmBe neutron source () over 33.5 h. Three populations ascribed to ()’s, events and nuclear recoils due to neutron scattering are well separated. (Inset) The thermal neutron capture peak, calibrated to the nominal energy of the reaction, together with a Gaussian fit. The energy resolution is

Fig. 7

Light-versus-heat scatter-plot obtained in a 20.5 h AmBe () calibration measurement with a 151 g scintillating bolometer (LMO-1; left figure). Rise and decay times as functions of the energy (right figures). The populations of and events used for the evaluation of the discrimination power are marked by black cycles and red triangles, respectively. The calculated discrimination power is 19, 5.4, and 8.1 by means of the light yield, the rise time, and the decay time parameters, respectively

Particle identification by heat signals

As it was shown before, particles exhibit a higher heat signal than ()s of the same energy. Even if a clear interpretation of this effect is lacking, this is probably related to the details of the phonon production mechanism in the particle interaction, which can lead to phonon populations with different features depending on the particle type. Therefore, one could expect some difference also in the shape of the heat signals between and () events and hence a pulse-shape discrimination capability of scintillating bolometers.9

Previous measurements with detectors demonstrated the possibility of pulse-shape discrimination by using only the heat channel [8, 31, 32]. However, the discrimination ability strongly depends on the experimental conditions and sometimes can fail [86]. No indication of this possibility has been claimed so far for bolometers.

A tiny difference between and ) heat pulses of the enrZMO-t detector (about 3% in the rising edge) allows us to perform an event-by-event particle identification using only the heat signals (e.g. was obtained for 2.5–3.5 MeV data). The data of the enrZMO-b bolometer, more affected by noise, show partial pulse-shape discrimination. It is worth noting that these results were obtained in spite of a low sampling rate (1 kSPS) and in one of the worst noise conditions among all the tested detectors.

-based bolometers also demonstrate the possibility of the pulse-shape discrimination by a heat-signal shape analysis. Unfortunately, the data of most detectors were acquired with a low sampling rate (1 kSPS) and/or do not contain a large statistics of and radiation in the same energy range, essential condition to investigate precisely this remarkable feature. However, significant results have been obtained by the analysis of the neutron calibration data (2 kSPS) of the LMO-1 detector. An example of the tiny difference in the time constants of ) and heat pulses (less than , i.e. a bin for the 2 kSPS sampling) is reported in Fig. 7 (right). By exploiting the rise and decay time parameters, we evaluated a between s in the 2.5–7 MeV and -triton events in the 5–7 MeV range as 5.4 and 8.1, respectively. These results could probably be improved by using other pulse-shape parameters, as it was demonstrated with detectors [8, 31]. However, due to a few per mille difference of the thermal signals induced by )s and , the pulse-shape discrimination of scintillating bolometers is expected to be less efficient in comparison to the light-assisted particle identification which exploits an about 80% difference in response (an exception for has been reported in [8]). This is also the case for the LMO-1 detector, for which the double read-out allows to reach about twice better discrimination power. However, the requirement of 99.9% rejection of -induced background (with a acceptance larger than 90%) is achieved even for 3, therefore pulse-shape discrimination with the heat signals only could allow to simplify the detector structure and to avoid doubling the read-out channels in a CUPID-like experiment.

Backgrounds and radiopurity of -containing scintillating bolometers

Alpha background

The spectrum measured by the ZMO-b detector in Run308 can be found in [35, 37], therefore the illustration of other spectra of the bolometers is omitted. The background spectra of events accumulated by the natural and all the enriched detectors are shown in Figs. 8 and 9, respectively. The anomaly (“dark and hot ”) in the response to in the bolometers was corrected by using the results of the fit to the events distribution in the LY-vs-heat data. The calibration data (168 h) of the enrLMO-b detector were combined with the background data to increase the statistics.

Fig. 8

The background energy spectra measured with the LMO-1 (over 237.5 h), LMO-2 (135 h), and LMO-3 (135 h) scintillating bolometers in the CUPID R&D set-up. The energy bin is 10 keV. The events in red are selected by the LY parameter (the events of the smeared source for the LMO-2 detector are not shown below 3.25 MeV). An internal potassium contamination of the LMO-1 crystal generates the continuum up to and the de-excitation peak at 1464 keV. The line visible in the LMO-1 data can be ascribed to the thorium contamination of the set-up. The spectrum of in the data of the LMO-2 detector is due to the presence of the smeared source. The peaks of (common for all the crystals) and with daughters (in LMO-3) are caused by the contamination of the crystals

Fig. 9

The energy spectra of events detected by the (left) and the (right) scintillating bolometers. The energy bin is 20 keV and 10 keV, respectively. The data of enrZMO-t and enrZMO-b (both over 593 h), and enrLMO-b (487 h) detectors were collected in the low-background measurements in the CUPID R&D cryostat at LNGS. The enrLMO-t (1303 h) bolometer has been operated in the EDELWEISS-III set-up at LSM

All the crystals exhibit a contamination by , however we cannot distinguish precisely a surface pollution from a bulk one. Furthermore, most likely the observed is due to contamination of the crystals, as this is the case for the scintillator (ZMO-b) [35]. The LMO-3 crystal, produced from the compound strongly polluted by (see Table 2), is contaminated by too. There is also a hint of a contamination of the other natural and crystals (ZMO-t, ZMO-b, LMO-1, and LMO-2), but the low statistics does not allow to estimate activity in the crystals. In addition to and , both crystals demonstrate a weak contamination by and .

The spectra were analyzed to estimate the activity of radionuclides from the U/Th chains and . Determination of activity in the detectors operated with the smeared sources (enrZMO-t, enrZMO-b, and LMO-2) is difficult. We assumed that the energy resolution of the peaks searched for is the same as the resolution of the peak present in the spectra of all detectors. The area of the peaks was determined within energy interval, where is a standard deviation of the Po peak. If no peak observed, the Feldman–Cousins approach [87] was applied to determine upper limits at 90% C.L. A summary of the radioactive contamination of the natural and -enriched and crystals is given in Table 8.

Table 8

Radioactive contamination of and crystal scintillators. The errors of the activities are estimated at 68% C.L., the upper limits are given at 90% C.L. The contamination of the LMO-3 sample is due to high activity of this radionuclide in the powder used for the crystal growth (see text)

Scintillator	ZnMoO_4		Zn^100MoO_4		Li_2MoO_4			Li2100MoO_4
Mo sublimation	Single		Double		Single	Single	Single	Double
Mo recrystallization	Double		Double		Double	Double	Double	Double
Boule crystallization	Double		Single		Single	Double	Single	Triple
Crystal ID	ZMO-t	ZMO-b	enrZMO-t	enrZMO-b	LMO-1	LMO-2	LMO-3	enrLMO-t	enrLMO-b
Position in boule	Top	Bottom	Top	Bottom	–	–	–	Top	Bottom
Crystal mass (g)	336	334	379	382	151	241	242	186	204
Radiputity test at	LSM	LSM	LNGS	LNGS	LNGS	LNGS	LNGS	LSM	LNGS
Time of measurements (h)	1540	1300	593	593	237	135	135	1303	487
Activity(μBq/kg)
^232Th	≤7.3	≤1.4	≤8.0	≤9.0	≤18	≤21	≤18	≤2.7	≤11
^228Th	≤26	≤4.6	≤8.0	≤21	≤18	≤21	≤18	≤8.4	≤6.2
^238U	≤13	≤2.6	10 (4)	39 (7)	≤18	≤37	≤48	≤4.9	≤11
^234U	≤20	≤3.0	11 (6)	43 (10)	≤18	≤21	≤46	≤6.7	≤11
^230Th	≤28	≤1.4	≤17	≤24	≤18	≤21	≤18	≤2.7	≤11
^226Ra	≤26	≤6.2	14(3)	23(4)	≤44	≤37	130 (19)	≤6.7	≤11
^210Po	575(18)	1320(30)	809(32)	2390(50)	139(33)	195(41)	76(25)	230(20)	60(10)
^235U	≤19	≤2.6	≤13	≤19	≤18	≤21	≤18	≤4.9	≤6.2
^231Pa	–	≤1.4	≤8.0	≤36	≤18	≤21	≤18	≤2.7	≤6.2
^227Th	–	≤2.6	≤8.5	≤12	≤18	≤21	≤18	≤4.9	≤6.2
^40K	–	–	–	–	62,000(2000)	≤12,000	≤3200	≤3500	≤3500
^190Pt	≤4.4	2.6(13)	–	–	≤18	–	≤18	≤2.7	≤11

The measured activity of in the crystals is 0.1 to 2 mBq/kg (see Table 8). If the contamination of samples is originated by , one can expect a growth of the activity up to the . The limits on the activity of other radionuclides from the U/Th families have been set on the level of 0.001–0.05 mBq/kg (a few exceptions of contamination by and/or will be discussed below). A hint for a content on the level is evident only for the ZMO-b sample, while for other crystals it is below 4–. It should be stressed that the sensitivity of the most measurements in the CUPID R&D set-up at LNGS is limited by the low exposure, while the constrains on radioactive contamination of the ZMO-t crystal are affected by the vibrational noise induced poor energy resolution of the bolometer (e.g. FHWM at 5407 keV of ).

The efficient segregation of thorium and radium in the growing process is evident from the comparison of the radioactive contamination of the sample LMO-3 (Table 8) and the powder (Alfa Aesar in Table 2) used for the crystal growth: the latter exhibits a clear pollution by (12 mBq/kg) and (705 mBq/kg), while no indication of () and a significantly reduced activity of (0.13 mBq/kg) were observed in the crystal. In addition, there is a clear sign of segregation of and its daughters along the crystal boule, because their concentration in the samples produced from the bottom part of the boule is around 2–4 times larger than that in the samples cut from the top part. Similar segregation has been reported, e.g., for CsI(Tl) [88], [89] and [90-92]. The results of the present work and Refs. [48, 53] show that the mechanism of segregation in the crystal growth process is less clear and further study would be useful to clarify this item. In general, it is expected [89, 92] that crystals produced by a double crystallization should be less contaminated. This is indeed observed for a boule (, and content at the level of few tens of ) in comparison to the recrystallized and scintillators (only limits below ten ). In summary, it is evident that the radiopurity level of the -enriched and crystals satisfies the demands of a next-generation bolometric experiment [7-9].

Surface radioactive contamination

As one can see in Fig. 4, the counting rate of energy-degraded s (that are expected due the surface contamination of the detector) in Run310 is lower than in previous runs. In particular, the rate in the energy range 2.7–3.9 MeV, excluding the region of , was reduced from 0.7(1) and 0.52(6) counts/year/kg/keV in Run308 (for ZMO-t and ZMO-b, respectively) to 0.20(6) and 0.09(5) counts/year/kg/keV in Run310 (for ZMO-b and enrLMO-t, respectively). These results are comparable to the 0.110(1) counts/year/kg/keV rate measured in CUORICINO [20]), but a factor 5 worse than the purity achieved in CUORE-0 (0.016(1) counts/year/kg/keV [21]). However, it is worth noting that no special efforts were dedicated to surface cleaning in the and detectors, while a significantly reduced amount of copper structure, a special surface treatment [93] and a dedicated mounting system [94] were adopted in CUORE-0.

Neutron background

The data acquired with the detectors were used to estimate the thermal neutron flux inside the EDELWEISS-III and CUPID R&D set-ups by exploiting the signature of neutron captures by . The data of the enrLMO-t detector do not contain any evidence of such events, while one event is found in the enrLMO-b data. The expected background in the region (the same used for the radiopurity analysis) is 0.054 (0.24) counts for enrLMO-t (enrLMO-b). According to Ref. [87], the number of events which can be excluded at 90% C.L. is 2.39 (4.11) counts for the enrLMO-t (enrLMO-b). Assuming 100% detection efficiency for such large-volume -containing detectors (e.g. see in Ref. [83]) and taking into account the total live time (1303/487 h) and the surface area (85.7/) of the enrLMO-t / enrLMO-b crystals, we estimate the following upper limits on the thermal neutron flux inside the EDELWEISS-III and CUPID R&D set-ups: and at 90% C.L., respectively. The constraint for the EDELWEISS-III is comparable to the limit of reported in [95] for the thermal neutron flux inside the lead/polyethylene shielding of the EDELWEISS-II set-up. The shield of both configurations of the set-up was the same except for 150 kg of polyethylene recently installed outside the thermal screens and inside the cryostat close to the detector volume. The limit for CUPID R&D is by order of magnitude improved to the one, which can be extracted in the same way from the data of previous measurements with a scintillating bolometer in this set-up [34]. It is worth noting that the deduced results are affected by an uncertainty which is difficult to estimate without the Monte Carlo simulations of the neutron propagation in the low temperature environment. It concerns a possible competition between the detectors and the neighbor materials in the capture of cold neutrons, further thermalized thermal neutrons as a result of interactions with a cold moderator of the set-up (e.g. 1 K polyethylene shield of the EDELWEIS-III).

Gamma(beta) background

background below 2615 keV

The background spectra of ) events measured by the ZMO-b detector in the Runs 308–310 are shown in Fig. 10. The data acquired at different positions of the detector inside the cryostat are superimposed. Few of the peaks present in the spectra are caused by the contamination of the set-up [96] and the detector components by K, Th and U. The natural isotopic abundance of molybdenum contains the isotope of at the level of 9.7% [82], therefore the decay of this nucleus gives a dominant background above 1.5 MeV even for the non-enriched detector (see Fig. 10).

Fig. 10

The normalized energy spectra of ) events accumulated in low-background measurements with the scintillating bolometer in the EDELWEISS-III set-up. A Monte-Carlo-simulated energy spectrum of the decay of with half-life (measured in the present work, see Sect. 5.4.3) is shown (upper panel). The energy bin is 10 keV. The same data in the 2.5–3.5 MeV energy interval (lower panel)

The ) background accumulated by three detectors in the CUPID R&D set-up is shown in Fig. 8. The region below 1.5 MeV of the LMO-1 detector is dominated by due to potassium contamination of the crystal. The main decay mode (branching ratio BR = 89.3% [97]) is a decay with [79]. The 1460.8 keV de-excitation -quanta following the electron capture in (; K-shell electron binding energy is 3.2 keV) is also clearly visible with a total energy of 1464 keV. The peak in the LMO-1 data is due to the thorium contamination of the set-up. The ) spectra of the LMO-2 and LMO-3 detectors contain only the peak caused by the potassium contamination of the set-up. The background of the LMO-2 bolometer is dominated by the spectrum of originated from the smeared source. The activity in the LMO-1 and the limits for the LMO-2 and LMO-3 crystals are given in Table 8. The comparison of the content in the LMO-1 and LMO-2 crystals demonstrates a segregation of potassium in the crystal growth process by at least a factor of 5.

Figure 11 shows the ) background of the -enriched detectors dominated above 1 MeV by the decay of with an activity of . Some difference in the background counting rate for several peaks measured by the bolometers (Fig. 11, left) indicates a position-dependent background inside the CUPID R&D set-up. In addition, one can see in Fig. 11 (right) the excess of events below 0.8 MeV for the enrLMO-b data, which indicates a higher external background in the CUPID R&D set-up in comparison to the EDELWEISS-III set-up.

Fig. 11

The energy spectra of ) events measured by the (left) and (right) scintillating bolometers. The energy bin is 10 keV. The data of both detectors and one (enrLMO-b) detector were accumulated in the CUPID R&D cryostat (593 and 319 h of data taking, respectively), while the enrLMO-t bolometer was measured in the EDELWEISS-III set-up (over 1303 h). The origin of the most intensive peaks is marked

background above 2615 keV

The background spectra of the detectors (except natural samples) contain events above the 2615 keV peak of , see Figs. 10 and 11). An event-by-event analysis excludes that they are due to random coincidences. Also they cannot be explained by decay of from the internal thorium contamination of the crystals, because the activity in the scintillators is low enough and no evidence of decays was found. The present surface purity of the detectors (see in Sect. 5.2) does not play a role on the surface-induced ) background, whose contribution is expected to be about two orders of magnitude lower that that of the surface [9]. These events can be originated by the muon-induced background, because no dedicated muon counter is available for the CUPID R&D set-up, while the scintillating bolometers operated in the EDELWEISS-III did not have a synchronization with the available muon veto. However, at least for detector ZMO-b the background above the 2615 keV peak cannot be completely ascribed to muons because of a clear run-dependent difference in the counting rate of events in the 2.65–3.5 MeV energy range: 0.14(3), 0.08(3), and 0.02(2) counts/day in Runs 308, 309 and 310, respectively (see Fig. 10). The ZMO-b crystal was kept at sea level before Runs 308 and 310 over about 60 days. Therefore, cosmogenic activation, relevant for the tested crystals due to a rather short cooling period underground (typically, less than one month), cannot be the origin of the observed decrease of rate in time. A crucial difference in the bolometer design of Run310 in comparison to early measurements is related to the absence of - connectors with CuBe press-fit contacts that were previously placed on the external lateral surface of the detector holder. According to [96], a considerable part of the background of the EDELWEISS-III set-up originates from radioactive contamination of the press-fit contacts (10(2) Bq/kg of ; the total mass of the press-fit is 40 mg per connector). The Monte Carlo simulations of the connector-induced background of the ZMO-b bolometer show that 0.57% of all decays of populate the 2.65–3.5 MeV energy region, corresponding to a rate of 0.07 counts/day, comparable to the ones measured in Runs 308 and 309. Therefore, we can conclude that for the detectors tested in the EDELWEISS-III set-up the main source of events above 2615 keV is the detector’s connectors. Moreover, as it is seen in Fig. 11 (right), the ) background rate inside the CUPID R&D set-up is even higher than that in the EDELWEISS-III, that can be explained by the radioactive contamination of the set-up. Therefore, special attention should be focused on selection of radiopure materials, in particular nearby the detectors, to realize a background-free decay experiment.

Double-beta decay of

To extract the decay half-life, the energy spectrum of the () events accumulated by the enrLMO-t detector in the EDELWEISS-III set-up was fitted by a simplified background model (Fig. 12). Taking into account a high crystal radiopurity, only two components of the internal background – the decay of and the bulk decay – are expected to give a significant contribution to the measured spectrum. The response function of the detector has been simulated with the help of the GEANT4-based code [98] and the DECAY0 event generator [99, 100]. A model of the residual background (assuming it was caused by external quanta from radioactive contamination of the materials surrounding the crystal) was built from an exponential function and the distributions of and (built from the calibration data). The best fit () obtained in the energy interval 160–2700 keV gives () decays of the of and () decays of the internal . Taking into account the mass of the crystal , the live time (, and the pulse-shape discrimination efficiency %, the bulk activity of in the enrLMO-t crystal is estimated to be ( (or ). The enriched crystal (% of ) contains nuclei of and, therefore, the decay half-life of is (statistical uncertainty only). Systematic uncertainties are related to the Monte Carlo simulations (5% – corresponds to the uncertainty of the GEANT4 modeling of electromagnetic interactions [101]) and to the fit by the background model (0.4%). To estimate the latter value, we made a fit with the described model in different energy intervals: the left side was varied in the 160–300 keV range with a 10 keV step, while the right side was varied within the interval 2650–2750 keV with a 50 keV step. The fits gave a number of decays in the range of 8821–8888. The combination of all the contributions results 5.4% as a total systematic uncertainty. So, the half-life of the decay of is measured toThe obtained value is in a good agreement with the most accurate results achieved by NEMO-3 experiment, [102] and [6.93 ± 0.04(stat.)] 10 [103],10 and the bolometric measurements with crystals [7.15 ± 0.37(stat.) ± 0.66(syst.)] [86], as well as with the average value [5].

Fig. 12

The () background spectrum accumulated over 1303 h with the 186 g -based detector (enrLMO-t) in the EDELWEISS-III set-up together with the fit by a simplified background model built from the distribution of (), internal (), and external quanta represented by exponential background (ext ), external and . The signal-to-background ratio above 1.5 MeV is 8:1

Because of the large mass and the relatively long measurement ( 600 h), the largest exposure was accumulated with the two detectors enrZMO-t and enrZMO-b operated in the CUPID R&D set-up at the Gran Sasso laboratory, containing nuclei. The usage of the smeared source, which also emits electrons with an end-point (see above the case of the LMO-2 detector), prevents us from getting a more precise half-life value of than that obtained from the analysis of the enrLMO-t background. However, these data were used to search for neutrinoless double-beta decay of and no counts were observed in the region of interest around 3034 keV. Considering an efficiency of in a 10 keV energy window, we set a limit on decay of Mo of  years at 90% C.L. Of course, this result is by far inferior to that achieved by NEMO-3 with 6.914 kg of Mo over the live time of 4.96 years ( years at 90% C.L. [103]), but – given the low sensitive mass and the short duration of the test – it shows the high potential of scintillating bolometers approach.

Down selection of -based scintillating bolometers technology

According to the results of the present work and some early related investigations, the state-of-the-art of the and scintillating bolometer technology is summarized below:In conclusion, the advanced crystal production process and better detector performance are crucial advantages of with respect to . For these reasons, was selected for the realization of a 10-kg-scale experiment (CUPID-0/Mo) aiming at demonstrating the viability of the LUMINEU scintillating bolometer technology for CUPID. Mass production of twenty enriched crystals with a size of has been recently completed for the first phase of this experiment, to be performed in the EDELWEISS set-up at LSM (France). The start of CUPID-0/Mo phase-I data taking is planned in early 2018 and the full-scale operation is expected by the end of the year. A second phase will follow aiming at a full use of the available 10 kg of .

In spite of about 10% higher concentration of molybdenum (55 vs. 44% weight), the unit volume of contains 5% less Mo because of the lower density (3.04 vs. ; see the properties of the materials in [41] and [23] respectively).

Naturally occurring Zn and Li do not contain radioactive isotopes. The only radioactive isotope in Mo natural composition is itself and its comparatively “fast” decay rate ( activity in the enriched crystal) requires fast detector’s response and pulse shape discrimination to avoid populating the decay ROI of by pile-uped events.

A rather low melting point (705 vs. ) and the absence of phase transitions are compatible with comparatively easier crystallization with respect to , and lower losses of the enriched material during the crystal growth process (0.1 vs. 0.6%).

Highly purified Mo-enriched molybdenum oxide [23] is usable for both materials. No special purification is needed for use with commercially available high purity zinc oxide and lithium carbonate. However, there is an issue with high contamination of Li-containing powder due to chemical affinity of lithium and potassium. Therefore, pre-screening measurements and purification are required to reduce potassium contamination in the crystal scintillators.

Double crystallization is an efficient approach to produce high optical quality radiopure and crystal scintillators.

The established technology of crystal growth allows the use of most of the material for the production of scintillation elements. In a case of , the crystalline material quality along the boules is not stable enough to reach the same high level of the ready-to-use scintillation elements production.

The hygroscopicity of is weak enough not to require a strict handling for the production of scintillation elements, mounting and operation of the detectors. The necessity of further improvement of the crystal surface purity is not presently evident but it would require the development of special mechanical/chemical treatment. is not hygroscopic and therefore an acid etching could be applied to improve the surface purity if it is needed.

The time response of and bolometers equipped with an NTD Ge thermistor (order of one to few tens ms) is comparable with an efficient suppression of the background caused by pile-ups of the decay of Mo. A further improvement (below is expected with faster temperature sensors, e.g. Metallic Magnetic Calorimeters [30], or with light detectors with enhanced signal-to-noise ratio (e.g. exploiting the signal amplification by Neganov–Luke effect [26]).

The energy resolution of bolometers satisfies the CUPID requirement. The resolution of does not meet this requirement by a factor of 2. Moreover, an addition degradation by a factor of 2 is observed for the -enriched crystals produced from the bottom part of the boule with presently available quality.

Typical light yield of is about 30% lower than that of . However, light-assisted alpha rejection at satisfactory high level ( and more) is achieved by detectors based on both materials.

The imperfections of crystals affect the bolometric response to bulk events, a fraction of which is characterized by more quenched light and enhanced thermal signals. The observed anomaly does not spoil the rejection capability, but affects the quality of the spectroscopy. This is not an issue of bolometers thanks to a significantly higher crystal’s quality.

The ability to perform heat-pulse-shape discrimination is a feature of both and detectors which allows a substantial simplification of the detector structure. The reproducibility of alpha particles rejection at the level of about has to be demonstrated but it is not mandatory for the scintillating bolometer technique.

High thermal-neutron cross section of () leads to reaction, which can be exploited to suppress neutron-induced background. The lack of a similar feature in ZnMoO would not suppress such background due to () reactions on Zn, Mo, and O isotopes which produce quanta with energies up to 7 MeV [104].

A possible cosmogenic activation of is expected to be much less significant than that of because no cosmogenically activated isotopes, with the decay energy high enough to contribute to the ROI, can be produced from lithium natural isotopes (in contrast to zinc isotopes). Therefore, a cosmogenic-originated background of would be only associated to molybdenum.

Very low contamination of both materials by U/Th completely satisfies the radiopurity demands even in the case of a single crystallization. The second crystallization further improves the crystals radiopurity thanks to the observed segregation of radioactive impurities.

Conclusions

A technology suitable for mass production of massive (), high-optical-quality zinc and lithium molybdate crystal scintillators from highly purified molybdenum enriched in has been established. The required performance and radiopurity of scintillating bolometers based on large-volume (50–) and crystals (including -enriched) have been demonstrated in low-background measurements at the Modane and Gran Sasso underground laboratories.

The detectors show an excellent energy resolution (in particular, 4–6 keV FWHM of detectors at the quanta of ), which is among the best resolution ever achieved with massive bolometers. The exploited heat-light dual read-out provides an efficient particle discrimination between () and events, which is compatible with more than 99.9% rejection while preserving approximately 100% selection efficiency of a signal. Furthermore, we demonstrated the possibility of pulse-shape discrimination by using the heat channel only, which is an important step towards detector simplification for the CUPID experiment. The scintillating bolometers are also found to be excellent neutron low-counting detectors. Their operation in the Modane and Gran Sasso cryogenic set-ups has allowed us to set very stringent limits, on the level of , on the thermal neutron flux in the EDELWEISS-III and CUPID R&D facilities.

The radioactive contamination of the developed -enriched crystal scintillators is very low. The activity of () and is below (down to a few of in case of single boule crystallization). The total bulk activity of U/Th is below a few mBq/kg. The activity of in the samples is less than 4 mBq/kg. The () background of the enriched detectors is dominated by the decay with a .

By utilizing the data accumulated over about 50 days with a 0.2 kg detector, the half-life of relative to the decay to the ground state of is measured with up-to-date highest accuracy: . The sensitivity to decay half-life on the order of has been reached with -enriched detectors operated over less than one month. It is far to be competitive to the most stringent limits in the field (e.g. limit deduced by NEMO-3 for ), but this sensitivity was achieved by using 3–4 orders of magnitude less accumulated statistics in comparison with the leading experiments. These results definitely demonstrate a potential of scintillating bolometers to perform high sensitivity searches.

Taking into account the reproducible technology to grow large-mass, high-optical quality crystals and their high bolometric performance together with low radioactive contamination, -based scintillating bolometers have been chosen for the realization of a cryogenic experiment with of enriched (CUPID-0/Mo) to prove the viability of this approach for CUPID. The first batch of twenty 0.2-kg crystal scintillators has been produced to carry out a first phase of the experiment in the EDELWEISS-III set-up at Modane (France).