The CUPID-Mo experiment for neutrinoless double-beta decay: performance and prospects

CUPID-Mo is a bolometric experiment to search for neutrinoless double-beta decay (0νββ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0\nu \beta \beta $$\end{document}) of 100Mo\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{100}\hbox {Mo}$$\end{document}. In this article, we detail the CUPID-Mo detector concept, assembly and installation in the Modane underground laboratory, providing results from the first datasets. The CUPID-Mo detector consists of an array of 20 100Mo\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{100}\hbox {Mo}$$\end{document}-enriched 0.2 kg Li2MoO4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Li}_2\hbox {MoO}_4$$\end{document} crystals operated as scintillating bolometers at ∼20mK\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 20\hbox { mK}$$\end{document}. The Li2MoO4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Li}_2\hbox {MoO}_4$$\end{document} crystals are complemented by 20 thin Ge optical bolometers to reject α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} events by the simultaneous detection of heat and scintillation light. We observe a good detector uniformity and an excellent energy resolution of 5.3 keV (6.5 keV) FWHM at 2615 keV, in calibration (physics) data. Light collection ensures the rejection of α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} particles at a level much higher than 99.9% – with equally high acceptance for γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}/β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} events – in the region of interest for 100Mo\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{100}\hbox {Mo}$$\end{document}0νββ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0\nu \beta \beta $$\end{document}. We present limits on the crystals’ radiopurity: ≤3μBq/kg\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\le 3~\mu \hbox {Bq/kg}$$\end{document} of 226Ra\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{226}\hbox {Ra}$$\end{document} and ≤2μBq/kg\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\le 2~\mu \hbox {Bq/kg}$$\end{document} of 232Th\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{232}\hbox {Th}$$\end{document}. We discuss the science reach of CUPID-Mo, which can set the most stringent half-life limit on the 100Mo\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{100}\hbox {Mo}$$\end{document}0νββ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0\nu \beta \beta $$\end{document} decay in half-a-year’s livetime. The achieved results show that CUPID-Mo is a successful demonstrator of the technology developed by the LUMINEU project and subsequently selected for the CUPID experiment, a proposed follow-up of CUORE, the currently running first tonne-scale bolometric 0νββ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0\nu \beta \beta $$\end{document} experiment.

uniformity and an excellent energy resolution of 5.3 keV (6.5 keV) FWHM at 2615 keV, in calibration (physics) data. Light collection ensures the rejection of α particles at a level much higher than 99.9% -with equally high acceptance for γ /β events -in the region of interest for 100 Mo 0νββ. We present limits on the crystals' radiopurity: ≤ 3 μBq/kg of 226 Ra and ≤ 2 μBq/kg of 232 Th. We discuss the science reach of CUPID-Mo, which can set the most stringent halflife limit on the 100 Mo 0νββ decay in half-a-year's livetime. The achieved results show that CUPID-Mo is a successful demonstrator of the technology developed by the LUMINEU project and subsequently selected for the CUPID experiment, a proposed follow-up of CUORE, the currently running first tonne-scale bolometric 0νββ experiment.

Introduction
Two-neutrino double-beta decay (2νββ) is one of the rarest processes in nature. Initially proposed by Maria Goeppert-Meyer in [1], it has since been observed for 11 nuclei with typical half-lives ranging from 10 18 to 10 24 year [2,3].
Numerous extensions of the Standard Model predict that double-beta decay could occur without neutrino emission (e.g., see [4][5][6][7][8]). This hypothetical transition, called neutrinoless double-beta decay (0νββ), is a lepton-number violating process. Its signature is a peak in the electron sum-energy spectrum at the Q-value of the transition (Q ββ ). Its observation could help explain the cosmological baryon number excess [9], and would prove that neutrinos are Majorana fermions (i.e., their own antiparticles) [5,10].
The current leading 0νββ decay experiments have a sensitivity on the 0νββ half-life of 10 25 -10 26 year [11][12][13][14][15]. At present, there is no confirmed observational evidence for 0νββ decay, which implies that next-generation experiments have to further increase their discovery potential by at least one order of magnitude.
One of the most promising technologies for 0νββ decay searches are cryogenic calorimeters, historically also referred to as bolometers [16]. These detectors are sensitive to the minute temperature rise induced by energy deposited in a crystal cooled to cryogenic temperatures (∼ 10 mK). Key benefits of bolometers are some of the best energy resolutions in the field ( E(FWHM)/E ∼ 0.2%), high detection efficiency and the possibility to grow radiopure crystals with a large degree of freedom in the choice of the material. Dual readout devices i.e. scintillating bolometers further allow for particle identification and thus yield the prospect of studying multiple 0νββ decay candidate isotopes in the background free regime [17,18].
The CUORE (Cryogenic Underground Observatory for Rare Events) experiment [14,19], currently collecting data at Laboratori Nazionali del Gran Sasso (LNGS, Italy), demon-strates the feasibility of a tonne-scale detector based on this technology. The success of this experiment is the starting point of CUPID (CUORE Upgrade with Particle IDentification), which aims to increase the mass of the 0νββ decay isotope via isotopic enrichment while decreasing the background in the region of interest.
According to the CUORE background model, the dominant background in the 0νββ decay region originates from α particles emitted by radioactive contamination of the crystals or nearby materials [20]. CUPID aims to identify and suppress this background using scintillating crystals coupled to light detectors [21]. A further background suppression can be attained by choosing 0νββ decay emitters with a Q ββ well above the 2.6 MeV line of 208 Tl, which is typically the end-point of natural γ radioactivity. Bolometers containing isotopes such as 100 Mo (Q ββ = 3034.40 ± 0.17 keV [22]), 82 Se (Q ββ = 2997.9 ± 0.3 keV [23]) or 116 Cd (Q ββ = 2813.50 ± 0.13 keV [24] satisfy this condition. The dual readout concept, where both the heat and light signals are recorded, has been implemented in two mediumscale CUPID demonstrators: CUPID-0, taking data at LNGS since 2017, and CUPID-Mo, which started the physics datataking at the beginning of 2019 in the Laboratoire Souterrain de Modane (LSM, France). Following the CUPID strategy, both experiments make use of enriched crystals (24 Zn 82 Se crystals for CUPID-0 [25] and 20 Li 2 100 MoO 4 crystals for CUPID-Mo) to search for 0νββ decay.
With an exposure of ∼ 10 kg × year, CUPID-0 proved that dual readout could suppress the dominant α background to a negligible level, obtaining the lowest background level for a bolometric experiment to date [26]. Nevertheless, the radiopurity and energy resolution (20.05±0.34 keV FWHM at Q ββ ) of the CUPID-0 crystals [26,27] do not meet the requirements of CUPID and demand further R&D activity if Zn 82 Se were to be used.
Conversely, the Li 2 100 MoO 4 crystals chosen by CUPID-Mo have demonstrated excellent radiopurity and energy resolution in the tests performed within the LUMINEU (Luminescent Underground Molybdenum Investigation for NEUtrino mass and nature) experiment [28,29].
The primary goal of CUPID-Mo is to demonstrate on a larger scale the reproducibility of detector performance in terms of the high energy resolution and efficient α rejection power combined with high crystal radiopurity. Given the high number of 100 Mo emitters contained in enriched crystals and the favorable 0νββ transition probability for 100 Mo, CUPID-Mo also enables a competitive 0νββ decay search.
The present work describes the CUPID-Mo experimental setup, currently operating in the EDELWEISS-III [30,31] cryostat at LSM. The detector was constructed in the clean rooms of the Laboratoire de l'Accélélarateur Linéaire (LAL) and the Centre de Sciences Nucléaires et de Sciences de la Matiére (CSNSM, Orsay, France) in the fall of 2017 and then moved to LSM and installed in the cryostat in January 2018. The detector was successfully operated through the summer of 2018 (Commissioning I). The fall of 2018 was devoted to cryostat maintenance, after a severe cryogenic failure, and detector upgrades. After optimization of the cryogenic system and detectors over the winter of 2019 (Commissioning II), the experiment has been collecting data in a stable configuration since the end of March 2019 (Physics run). In this paper we present the CUPID-Mo detector concept and construction (Sect. 2), the operation and initial performance of the first Physics run dataset (Sect. 3), and the prospects of the experiment in 0νββ decay search (Sect. 4).

Experimental setup
CUPID-Mo consists of an array of 20 scintillating bolometer modules arranged in five towers, each with four modules, as shown in Fig. 1. Each module contains one Li 2 100 MoO 4 crystal and one germanium wafer assembled inside a single-piece copper housing, instrumented with Neutron Transmutation Doped (NTD) Ge thermistors. All the materials used for the towers' construction were carefully selected, and additionally cleaned as needed to minimize radioactive contamination. The detector construction, transportation, and assembly into the underground cryogenic facility were performed in a clean environment. The key ingredients of the detector, its assembly, and the cryogenic apparatus are detailed below.   [28,29]. An additional sixteen new Li 2 100 MoO 4 crystals were fabricated with the identical procedure as that employed by LUMINEU [28,32,33]. All crystals have a cylindrical shape with ∼ 44 mm diameter and ∼ 45 mm height, and a mass of ∼ 0.2 kg. The crystals were produced at the Nikolaev Institute of Inorganic Chemistry (NIIC, Novosibirsk, Russia) as follows: • purification of the ∼ 97% enriched molybdenum [32], previously used in the NEMO-3 experiment [34]; • selection of lithium carbonate with low U/Th and 40 K content [28] and purified 100 Mo oxide [33]; • crystal growth via a double crystallization process using the low-thermal-gradient Czochralski technique [28,33]; • slicing of the scintillation elements, and treatment of their surfaces with radio-pure SiO powder.
The total mass of the 20 Li 2 100 MoO 4 crystals used in CUPID-Mo is 4.158 kg, corresponding to 2.264 kg of 100 Mo.

Ge slabs
The high-purity Ge wafers (Umicore Electro-Optical Material, Geel, Belgium), used as absorbers for the scintillation light, have a diameter of 44.5 mm and a 175 µm thickness. A ∼ 70 nm SiO coating was evaporated on both sides of the Ge wafer to make them opaque, thus increasing the light collection by ∼ 35% [35]. A small part of the wafer surface was left uncoated to ease the gluing of a temperature sensor.

Sensors
CUPID-Mo employs NTD Ge thermistors [36] as thermal sensors. These thermistors were provided by the Lawrence Berkeley National Laboratory (LBNL, Berkeley, USA) and come from a single production batch. The NTDs used for the Li 2 MoO 4 (LMO) bolometers are 3.0 × 3.0 × 1.0 mm 3 in dimension, and have a temperature-dependent resistance given by R = R 0 · e (T 0 /T ) 0.5 where the average values for the parameters are T 0 = 3.8 K and R 0 = 1.5 . Given the lower heat capacity of the Ge absorbers for the bolometric light detector (LD), we opted to better match and reduce the heat capacity of their sensors by dicing the NTDs into multiple pieces.
In Commissioning I, we produced three sensors with 3.0 × 0.8 × 0.4 mm 3 dimensions from the slicing of a single NTD in two directions. The LDs with these sensors showed an unexpectedly high noise with a strong 1/ f [Hz] component reaching frequencies up to hundreds of Hz. For Commission-ing II and beyond, we replaced all but two sensors with new ones with 3.0 × 0.8 × 1.0 mm 3 dimensions, avoiding the horizontal cut of the original NTDs.
In addition to the thermistor, each Li 2 100 MoO 4 crystal is instrumented with a silicon-based resistive chip [37] operated as a heater. This heater allows us to periodically inject a constant power and generate a pulse of constant energy. The resulting reference pulses can be used in the offline analysis to monitor and correct for a change of the signal gain due to temperature drifts of the bolometer [38].

Sensor coupling
The NTDs were glued on the Li 2 100 MoO 4 crystals using the dedicated tool shown in Fig. 2, similar to the one used by CUPID-0 [25]. The glue is a two-component epoxy resin (Araldite® Rapid) well-tested for cryogenic applications and demonstrated to have acceptable radiopurity [20]. The gluing tool features a part for holding the NTD, and it can be moved along the vertical axis and fixed at any level. The performance of the bolometer is strongly dependent upon the quality of gluing, and we obtain optimal results when separate glue spots connect the NTD to the crystals. This helps in compensating the different thermal contractions of the involved materials. To maintain separate glue spots during and after the epoxy curing, the NTD is kept 50 µm from the Li 2 100 MoO 4 crystal, which is positioned on the top surface of the gluing tool. NTDs of five crystals (LMO-1-4,15), used in the LUMINEU experiment and/or the CUPID-Mo single tower test, were glued with six spots, while nine spots were applied for the remaining crystals. The heaters were glued with a single glue spot using a 50 µm Mylar mask to provide a gap between the crystal surface and the chip.
The gluing of NTDs to the Ge wafer was also performed with the two-component epoxy resin described above. However, instead of the six or nine spot matrix, we applied a uniform veil of glue. This choice was motivated by the small size of the sensor and the less pronounced effect of thermal contraction expected for the Ge-glue-Ge interface. We used the manipulator of an ultrasonic bonding machine to provide a controlled force to attach the NTD to the Ge wafer surface and provide better reproducibility.

Detector structure
The CUPID-Mo single module and tower structure were designed by the Service de Physique de l'Etat Condensé (SPEC) at CEA (Commissariat à l'Énergie Atomique et aux énergies alternatives, Gif-sur-Yvette, France) according to the following requirements: -the single module structure should be compact but permit the housing of four scintillating crystals in a single tower, taking into account the restricted space in the experimental set-up; -the towers should be suspended by dedicated springs to mitigate the vibrational noise of the set-up [28]; -the design should allow a simple installation inside the cryostat (see Fig. 1 and Sect. 2.8).
The mechanical workshop of LAL (Orsay, France) fabricated the detector support structure. Each detector module (see Fig. 3) is a single-piece holder, made of highly radiopure NOSV™ copper from Aurubis (Hamburg, Germany). It contains both a Li 2 100 MoO 4 scintillation element and a Ge wafer. The bolometers are kept in place by small Polytetrafluoroethylene (PTFE) holders which decouple them from the thermal bath. The Li 2 100 MoO 4 crystals are supported with three PTFE elements on the top and bottom, while the LDs are clamped with three PTFE pieces. In Commissioning I, we did not install any reflecting foil around the crystals, because previous measurements demonstrated efficient particle identification performance despite a factor of 2 lower light collection efficiency [18]. Commissioning I was characterized by sub-optimal LD performance (see Sect. 2.3), hence we decided to surround the crystals' lateral side with reflecting foil (3M Vikuiti™) in addition to the replacement of the LDs' NTDs.

Detector assembly
We performed all activities related to the detectors' assembly in a cleanroom environment. All the used detector compo-

Fig. 3
All components used to assemble a CUPID-Mo single detector module: a copper holder, a Li 2 100 MoO 4 crystal with glued NTD and heater, a Ge LD with NTD, the copper screws, the PTFE spacers and fixing elements, the Kapton foil with golden pads. Note that the reflecting film is not shown here

Wiring
A dedicated wiring scheme was designed and implemented for the CUPID-Mo experiment as the existing EDELWEISS-III readout could not accommodate the additional 20 dualreadout modules required for CUPID-Mo.
We bonded gold wires from the NTDs to flat Kapton pads with gold contacts to provide the electrical readout connection as well as the weak thermal link to the heat bath. In fact, the pads are glued on the copper detector holder, which is in excellent thermal contact with the coldest part of the dilution refrigerator. Silk-covered constantan twisted wires were sol-dered on the other side of the pads and run up each tower to a larger Kapton pad with gold contacts glued at the top of the tower.
On this pad the constantan wires and copper wires (connection to the cold electronics) were soldered. This connection provides a link to Si-JFET (junction gate field-effect transistor) based pre-amplifiers at 100 K through the copper plate inside the EDELWEISS cryostat (see Sect. 2.8).

Low background cryogenic facility
The CUPID-Mo detector array is installed (see Fig. 7) in the EDELWEISS-III cryogenic set-up [28,30], located in LSM. This site is among the deepest underground laboratories in the world; the 1700 m (4800 m water equivalent) rock overburden, provided by the Frejus mountain, reduces the cosmic muon flux to 5 muons/m 2 /day [39].
The EDELWEISS cryostat is a custom dilution refrigerator with a reversed geometry [30], developed by Institut Néel (Grenoble, France). During a cryogenic run, this set-up requires periodic refilling of the liquid helium (LHe) bath every 10 days. The consumption of LHe is minimized by the use of a cold vapor reliquefaction system based on three Gifford-MacMahon cryocoolers. The cryocoolers are responsible for most of the vibrational noise in the set-up, thus necessitating the use of the suspension to achieve highperformance operation of the scintillating bolometers [28]. The passive shielding of the set-up against environmental radiation consists of lead (20 cm thickness) and polyethylene (55 cm thickness). The inner part of the lead shield is made The position and types of calibration sources are also indicated of 2 cm thick low 210 Pb radioactivity (< 0.12 Bq/kg) lead recovered from sunken Roman-era galleys (hereafter called "Roman lead"). An additional internal Roman lead (14 cm) and polyethylene (10 cm) shield at the 1K-plate is used to protect the detectors from radioactivity from the cryostat components. A muon veto system is surrounding the whole cryostat providing 98% geometrical coverage. The muon veto is constructed from 46 individual plastic scintillator modules with a total surface of 100 m 2 and provides a detection efficiency of 97.7% for muons passing through a central sphere with 1 m radius [39]. With a two module trigger rate of O(1 mHz) an efficient rejection of muon-induced events with negligible dead-time is expected for the 0νββ search program. The setup is located inside a class 10000 cleanroom with a depleted radon air supply (∼ 30 mBq/m 3 of 222 Rn).
The experimental volume of EDELWEISS-III contains four floors (detector plates) with twelve slots each (see Fig. 8). The CUPID-Mo towers were inserted through the slots T3, T4, T10, T11, and T12 (see Table 1) and mechanically decoupled from the EDELWEISS-III detector plate with three metal springs for each tower. The remaining experimental space is partially occupied (tower slots T2, T5, T7, and T8) by 11 Ge bolometers for the EDELWEISS dark matter search program [40] and a cadmium tungstate based scintillating bolometer for CUPID R&D [41].
The construction of the CUPID-Mo detector has implied the introduction of two new elements in proximity of the Li 2 MoO 4 crystals: solder for the connection of the constantan wires to the Kapton pads and springs for the suspension of the towers. In prospects, these elements can contribute to the background budget of the experiment. High purity Ge counting revealed that the solder is not Pb-free (the activity We added two mixed U/Th sources made of thorite mineral to the EDELWEISS-III automatic source deployment system [30] (see Fig. 8). These sources complement the already available γ -calibration sources of 133 Ba (∼ 1 kBq) and 60 Co (∼ 100 kBq) for periodic calibration of the CUPID-Mo detectors. The activities of the sources are ∼ 50 Bq of 232 Th, ∼ 100 Bq of 238 U, and few Bq of 235 U. The 133 Ba source emits γ s with energies up to 0.4 MeV and was only used during the commissioning stage. The high-activity 60 Co γ source is used to eliminate space charges in the dual readout heat-ionisation Ge bolometers of EDELWEISS-III [30] and also to calibrate the CUPID-Mo Ge LDs via source-induced X-ray fluorescence [42,43] (see Sect. 3.3). The 60 Co source is used mainly during and just after each LHe refill (every 10 days) while a regular ∼ 2-days-long Th/U calibration is scheduled for each period between subsequent LHe refills.
The detector readout in the EDELWEISS-III setup is based on AC-biased cold electronics [30], which restricts the use of high-resistivity thermistors to at most a few M resistance at a given bias current (working point) [28]. Custom made room-temperature electronics modules called bolometer boxes (BBs) are mounted just outside of the cryostat to ensure short cables to limit noise pick-up. These BBs contain the electronics for the cold Si-JFET pre-amplifiers' biasing, Digital to Analog Converters (DACs) for the detectors' biasing, post-amplification, anti-aliasing filter, and ADCs to record the CUPID-Mo NTDs [30]. All LMOs and five LDs are operated with BBs containing 16-bit ADCs, while the signal digitization for the remaining LDs is done with 14bit ADCs. The pulser system, used to inject a constant Joule power through the heaters, is based on a 4-channel pulse generator with a typical injection periodicity of a few minutes. The data acquisition system [28,30] can record both online triggered and stream data; the triggered data is used only for monitoring purposes.

CUPID-Mo detector operation
Of the 20 LMO and LD pairs, only a single LD was lost due to a hardware issue, resulting in 39 out of 40 active channels. Additionally, 18 out of 20 heaters are available to inject pulses. The optimal working point of the LMO detectors was chosen to maximize the signal amplitude. LDs, instrumented with smaller, more resistive, sensors operate in an over-biased regime to obtain an NTD resistance of ∼ 1 M mitigating the impact of AC biasing (see details in [28]). The modulation frequency of 500 Hz was chosen to reduce the pick-up of cryocooler-induced high-frequency noise 1 . The nominal base temperature of the empty EDELWEISS cryostat is 11.5 mK. In the present, densely populated, cryogenic setup, an additional heat load increases this base temperature to ∼ 20 mK and we could stably operate at 20.7 mK with a few μW of regulation power. This temperature is considerably higher than the operating temperature in the LUMINEU predecessor [28] and it is expected to have an adverse effect on the detector performance.
Nevertheless, the following analysis of a ∼ 2 week period with 11.1 days of physics data, 2.2 days of mixed Th/U source calibration, and 1.6 days of 60 Co irradiation provides a robust confirmation of the bolometric performance achieved within LUMINEU [18,28]. The data were acquired between March 24th 2019 and April 6th 2019 and correspond to a physics exposure of 0.1 kg × yr Li 2 100 MoO 4 . This early data is comparable with the prior exposure presented in [18], and emphasizes the reproducibility of Li 2 MoO 4 detectors using a total of 20 detectors.

Data processing
Two independent analysis frameworks, both exploiting the optimum filter technique [44], are used for the data processing: one, called DIANA [45,46], is adapted from the CUORE [14] and CUPID-0 [46] experiments and the other was developed at CSNSM [47] and used for the analysis of the LUMINEU data [28]. The CSNSM code has been developed (using the MATLAB MULTI Integrated Development Environment) specifically for the analysis of scintillating bolometer data. It is more nimble and readily adapted to different experimental setups. In contrast, DIANA is a much broader framework, including analysis packages for larger detector arrays (in particular allowing for the analysis of coincident events). It is object-oriented C++ code with a PostgreSQL [48] database interface to track detector and electronics settings. The use of DIANA allows for comparison between different CUPID project demonstrator experiments with effectively the same analysis tools, and DIANA is expected to be used as the primary package for CUPID-Mo in the future. Therefore, all results presented below are based on the use of DIANA, while it is noted that very similar results were obtained with the CSNSM code, providing a cross-check of the DIANA processing.

Performance of bolometric Ge light detectors
Characteristic pulse shape parameters such as the rise-and decay-times, defined as 10-90% of the rising edge and 90% to 30% of the trailing edge of the LD pulse shape have been investigated (see Table 2). We estimate typical (median) riseand decay-times of 4.2 ms and 9.2 ms respectively from an averaged pulse, triggered and aligned on events recorded in an associated LMO crystal.
Averaging of pulses was necessary since Li 2 MoO 4 has a moderate Relative Light Yield (RLY) which does not exceed 1 keV/MeV relative to the heat signal (see [28,29] and Sect. 3.4), and estimates from individual light pulses are subject to bias from noise fluctuations. We note that in particular for the rise-time both the 500 Hz sampling and the alignment of the average pulse become limiting factors for a more precise estimate. At a previous surface test at CSNSM with a similar temperature, working point, and a 10 kHz sampling rate, a factor of 3 faster rise-time (0.96 ms) was observed in LD 4.
To estimate the performance of the Ge LDs, we perform an in situ calibration. We employ the X-ray fluorescence of Mo or Cu that is generated when the crystals and setup are exposed to a higher intensity γ source [42,43]. For Mo we expect characteristic peaks from the K α1 (17.48 keV, intensity I = 100%), K α2 (17.37 keV, I = 52%), and K β1 (19.61 keV, I = 15%) lines [49]. The Cu X-rays can give additional peaks from K α1 (8.05 keV, I = 100%), K α2 (8.03 keV, I = 51%), and K β1 (8.91 keV, I = 17%). Figure 9 shows a typical X-ray spectrum obtained during the 60 Co source irradiation. The prominent features are a sum K α peak from Cu and both a sum K α and a distinct K β peak from Mo. The intensity of the Cu X-rays is much lower than those associated with Mo, as the Cu is only facing the LDs on the side. Also, the statistics in the Cu K α peak are very low for detectors far from the 60 Co source, and we chose to omit this peak from the LD calibration.
With a stable operating temperature of 20.7 mK and a strong NTD polarization for the Ge LDs, negligible nonlinearity is expected. We use a Gaussian fit to the most intense Energy (keV) 5 1 0 1 5 2 0 2 5 3 0 Entries / (0.2 keV) 10 2 10 X-ray α K Mo X-ray β K Mo X-ray α K Cu Fig. 9 Energy spectrum of CUPID-Mo light detector LD2 after a 33 h 60 Co irradiation in the EDELWEISS-III set-up. The Mo K α X-rays are used for the LD calibration peak, the Mo K α X-rays, and perform a first-order polynomial calibration with zero intercepts. The 1.4 g Ge LDs are instrumented with small-size NTDs that achieve a typical sensitivity of 1.1 μV/keV with an RMS of ∼ 40% (see Table 2). Uncertainties in the individual sensitivity estimates are dominated by the gain in the analog chain with typical uncertainties of order 10% for several of the operational amplifiers in the amplification chain.
The LDs sensitivity is limited by a comparatively high regulation temperature of the detector plate and the strong NTD polarization. 2 We estimate the baseline resolution for all detectors from a set of forced random trigger events injected every 101 s. We exclude one detector instrumented with a different NTD sensor (used in Commissioning I), see Table 2 and runs with atypical noise performance, resulting in 183/209 (LD-bolometer, run) pairs. The median of these estimates yields a typical baseline resolution of 148 eV FWHM in agreement with the channel based estimate in Table 2. We see good reproducibility with individual channel estimates ranging from 66 up to 368 eV.
A resulting scatter plot of the correlation between the sensitivity and the achieved baseline resolution is shown in Fig.  10. We note that the spread in detector performance is only slightly higher than for the NTDs on the Li 2 MoO 4 crystals (see Sect. 3.4). We want to emphasize the uniformity that lends itself to applications in larger cryogenic detector arrays.
For reference, we list the performance characteristics on an individual LD basis in Table 2. The reported performance in terms of the baseline resolution exceeds the requirements to achieve a better than 99.9% rejection of α events at 99.9% acceptance of γ /βs as is discussed in detail in Sect. 3. Several improvements can be pursued for the full-size CUPID experiment. DC-biased electronics, higher sampling rate, and the  implementation of additional analysis and de-noising techniques can improve the quoted performance. Futhermore, lower noise NTD Ge sensors and a lower operational temperature resulting in a higher detector sensitivity can also yield a significantly better LD performance, as demonstrated with a 20 eV FHWM baseline resolution in [51].

Performance of Li 2 100 MoO 4 bolometers
The time constants of LMO bolometers are much longer than those of LDs. We obtain median values of 24 ms for the risetime and 299 ms for the decay-time with a significant spread of 208 ms in the decay-times and a smaller spread of 8 ms in the rise-time (see Table 3). These values are consistent with previously reported values [28], and in the typical range for macroscopic cryogenic bolometers operated in the tens of mK range. We calibrate with a mixed Th/U source, with a most prominent peak at 2615 keV ( 208 Tl) and negligible gamma continuum, see Fig. 13. This is the closest observable γ -line, ∼ 415 keV lower than the Q ββ -value of 100 Mo.
The calibration data were acquired over a short period (2.2 days), resulting in limited statistics of the detected γ peaks. We neglect nonlinearities in the detector response and fit using zero and the 2615 keV 208 Tl line. Additionally, we use this peak to correct for changes in thermal gain due to slow temperature drifts in the experimental setup. The resulting correction is a linear scaling factor obtained from the optimum filter (OF) amplitude versus baseline dependence in the calibration data and is applied to both calibration and background data.
The detector sensitivity at the 20.7 mK operation temperature has a median value of 17 nV/keV with an RMS of about  Table 3). For unknown reasons, the detector LMO 2 shows very low sensitivity in comparison to the results of the CUPID-Mo Commissioning I (6 nV/keV at 20.5 mK) and LUMINEU (47 nV/keV at 17 mK [42]). As in the case of LDs, larger sensitivity of LMO bolometers is expected at colder temperatures (e.g., compare results given in [28]). The same method utilized for the investigation of LDs' baseline resolution (see Sect. 3.3) is also applied for Li 2 100 MoO 4 bolometers. We obtained characteristic (median) values of 1.96 keV FWHM for the baseline resolution with the spread of the distribution given in Fig. 11 and individual detector based resolutions presented in Table 3. The baseline noise versus sensitivity data is also illustrated in Fig. 12.
For further analysis, we utilize a preliminary set of analysis cuts. First, periods of atypical noise and temperature spikes of the cryostat are rejected, removing ∼11% of the data from the commissioning period. A large part of the loss of livetime is caused by a suboptimal setting of the cryo- All the major peaks have been labeled. The inset shows a fit of the 208 Tl γ peak at 2614.5 keV stat suspension, and improved stability has been observed in more recent data. 3 We exclude pile-up events with another trigger in a (−1, +2) s window which incurs an inefficiency of ∼ 4% in physics data, dominated by the frequent injection of our heater pulses during the commissioning period. We require a baseline slope consistent with the typical behavior of the channel, and require both the rise-time as well as the optimum filter peak position to be within 5 median-absolute deviations (MAD) of the mean range as defined by the overall distribution of these values. We further select γ /β events by requiring events to have a RLY (see Sect. 3.5) within 4σ of the mean amplitude incident in a LD associated with a LMO bolometer.
The resulting calibration data are presented as a summed spectrum in Fig. 13. The 2615 keV 208 Tl resolution is 5.3 keV FWHM estimated with an unbinned extended maximum likelihood (UELM) fit shown in the inset. The fit model includes a Gaussian function and two components, a smeared step function for multi-Compton events and a locally flat background. We note a potential bias on the resolution since we perform the thermal gain stabilization on this gamma peak and are in a low statistics limit. A toy Monte-Carlo (MC) with a typical value of 20 counts per detector resulted in an estimated bias (underestimate of the 208 Tl peak width) of 0.3 keV.
In addition to the good energy resolution, we highlight the linearity and uniformity of the data. The maximum residual between observed peak position and expected peak position in the summed calibration spectrum was 3 keV for the 1120 keV line from 214 Bi. Similarly, we observe an excess width for all γ peaks of at most 5 keV due to not yet accounted for individual detector non-linearities. Table 3 Performance of 100 Mo-enriched Li 2 MoO 4 bolometers of the CUPID-Mo experiment operated at 20.7 mK in the EDELWEISS setup at LSM (France). This table contains the following information: the crystal size and mass, the NTD resistance at the working point (R W ork ), the rise-time (τ R ), the decay-time (τ D ), the voltage sensitivity (A Signal ), the baseline noise resolution (FWHM Noise ), the scintilla-tion light yield (RLY) measured by top LD (RLY T op ) and bottom LD (RLY Bottom ), and the light yield quenching for alpha particles (Q F α ). The omission of a measured parameter due to lack of statistics or insufficient performance/non-operational light detector is indicated by "-". The median value for RLY Bottom is given for scintillators coupled to two LDs ( a ) and for single LD ( b ); see text enrLMO Size

Performance of light-vs-heat dual readout
We estimate the RLY from events in the 2-3 MeV region, close to the Q-value for 0νββ of 100 Mo. We create a distribution of light/heat energies and fit a Gaussian to this distribution to obtain the RLY μ γ for γ /β events.   Fig. 15 Light yield versus heat signal scatter-plot using 11 days of physics data from LMO 1. The presented detector has the highest 210 Po contamination to illustrate best the distributions of α and γ events and the scintillation light quenching for α events for downward going light. The obtained results are consistent with previous observations [28,29]. The summed light collected from two adjacent LDs is the closest estimate we have for ideal light collection. It is as high as 1.44 keV/MeV with a median value of 1.35 keV/MeV. The uncertainty for individual RLY estimates has been quantified from the spread in RLY estimates of three distinct 60 Co plus 208 Tl datasets. We observed a ∼ 4% spread around the mean (RMS), with a maximum deviation of 16% for a single detector.
For this analysis, we opt to use the LD in the same detector module just below the crystal by default. In cases where the lower LD is unavailable or performs significantly worse (LMO 1, 3, 6 and 7) we switched to associating the upper LD to this crystal (see also Tables 2 and 3).
Taking into account the measured RLY (Table 3) and the LD performance (Table 2), all detectors achieve better than 99.9% discrimination of α events (see Sect. 3.6) with a typical example of the discrimination power given in Fig. 15. The preliminary γ /β selection by RLY (blue) defined before eliminates a significant population of α events with ∼20% of the RLY of γ /β events and a few remaining events at higher light yield than expected (red). This particular crystal is characterized by the highest contamination level of 210 Po with ∼ 0.5 mBq/kg and hence best exemplifies the alpha discrimination power achieved for a scintillating bolometer with typical performance values of 0.67 keV/MeV RLY and 0.18 keV FWHM N oise of a coupled LD. We observe that the 210 Po α events misreconstructed at ∼7% higher energy at 5.8 MeV instead of 5.4 MeV. This shift is much larger than nonlinearities in the γ region would suggest, but we note that a similar difference in the detector response for α particles has been observed previously with lithium molybdate based detectors [28,29]. Events at higher light yield than γ /β events can be observed due to noise spikes and misreconstructed amplitude estimates in the LD, as well as due to close β contaminations with a coincident γ depositing energy in the Li 2 100 MoO 4 crystal. We estimate a scintillation light quenching of α-particles with respect to γ /β particles of (19.7 ± 1.0)% across the detectors (see Table 3). These results are also within expectations for this scintillation material [18,28,29].
3.6 Extrapolated α discrimination of Li 2 100 MoO 4 scintillating bolometers We systematically evaluate the α discrimination level following Refs. [28,29,52] and report the discrimination of α versus γ /β events in terms of the discrimination power (DP) at the Q-value for 0νββ in 100 Mo The parameters in the definition of the DP are the mean RLYs μ α , μ γ for α and γ /β events respectively, and resolutions σ α , σ γ . We obtain detector based values μ α = Q F α ·μ γ from the measured μ γ and approximate the very uniform light quenching of α events with Q F α = 0.2 (see Table 3). The expected LD resolutions σ α and σ γ at the endpoint of the 100 Mo decay are extrapolated by adding the baseline resolution and a statistical photon noise component with an average photon energy of 2.07 eV [52] in quadrature. The resulting median discrimination power is 15.0, with the worstperforming detector having a discrimination power of 6.3. Hence all detectors are expected to achieve better than 99.9% α rejection with more than 99.9% γ /β acceptance.
We note that this model calculation does not take into account additional sources of uncertainty such as variation associated with the position of the incident particle interaction and subsequent light propagation. However, the validity of the model is supported by the excellent agreement between the predicted and achieved discrimination in neutron calibration data in previous measurements [29]. The computed discrimination level exceeds the requirements for CUPID, and we plan to study adverse effects due to non-Gaussian tails with larger statistics in the future. If multiple alpha peaks emerge in individual detectors we will also be able to study the α energy scale and the energy dependence of the α discrimination from data. It should be noted that we are only using a single of the two LDs, typically the one at the bottom of each detector module (see Sect. 3.4). An optimized selection of the better performing LD, or a combined light estimate using both LDs will further improve the quoted discrimination. In addition it is expected that information from the combination of the LDs could potentially be relevant to break degenerecies if non-gaussian tails related to contamination at the NTDs or the LDs were encountered. We apply an additional anticoincidence cut with a time coincidence window of 100 ms between Li 2 100 MoO 4 detectors, a so called multiplicity one (M1) cut to reject multi-Compton and muon shower events and obtain the background spectrum shown in Fig. 16. The γ /β spectrum of Li 2 100 MoO 4 bolometers above ∼1 MeV is dominated by the 2ν2β decay of 100 Mo with an activity of 10 mBq/kg [28]. In 11.1 days of background data, we observe no event compatible with the RLY of γ /β events above 3034 keV, the Q-value for double-beta decay in 100 Mo. The estimate for the resolution at 2615 keV, (6.5 ± 1.4) keV, is compatible with the prediction from the calibration data, albeit this is subject to considerable uncertainty due to the limited statistics.
By relaxing the pulse shape cuts, and removing the γ /β RLY cuts and the M1-cut, we can investigate the α region with a set of basic cuts that are designed to have better than 99% acceptance. The only clear contaminants seen in the spectrum (see Fig. 17) are a bulk and surface peak by 210 Po ( 210 Pb), as observed in the LUMINEU studies [28,29]. For a few of the nuclei in the Th-( 232 Th, 228 Th, 224 Ra, 212 Bi) and the U-chain ( 238 U, 234 U+ 226 Ra, 230 Th, 222 Rn, 218 Po), that have a clean decay signature we start to see first hints of a contamination. However, the number of events in a ±30 keV window around each nominal decay energy (4082 keV, 5520 keV, 5789 keV, 6207 keV and 4270 keV, 4770 keV, (4871 keV and 4858 keV), 5590 keV, 6115 keV) is compatible with zero at the 2σ level for all of the decay signatures. We thus place a conservative upper limit at a level of 2 μBq/kg (Th-series) and 3 μBq/kg (U-series) (90% CL) on the activity in the U-/Th-chains using the largest observed event count for any of the decays in the U-/Th-chain.
We look for possible backgrounds from surface contaminants in the 3-4 MeV region. Excluding a potential 190 Pt alpha bulk contribution in a ±30 keV window around 3269 keV we observe 14 events, which is equivalent to a background of (0.14 ± 0.04) counts/(keV×kg×yr) in degraded alpha events before the rejection by RLY.

Outlook
Based on these first physics data we are confident that the Li 2 100 MoO 4 cryogenic detectors possess a high degree of reproducibility and are well suited to scale to a much larger CUORE sized detector array in CUPID [50]. We expect the current CUPID-Mo experiment to be able to set significant limits on 0νββ in 100 Mo. Consequently, we evaluate the CUPID-Mo sensitivity using the Bayesian method for limit setting for a counting experiment in a ±2σ region of interest around Q ββ . At present, several steps of the data analysis procedure have not been fully optimized, leaving room for improvement. We expect to achieve a 5 keV energy resolution (FWHM) at Q ββ with a dedicated optimization of the energy reconstruction algorithms. With an average containment efficiency for 0νββ decay events of 75%, and assuming a ∼ 90% analysis efficiency, as demonstrated in CUORE [14], CUPID-0 [26,27] and LUMINEU [28] for the combined trigger efficiency, multiplicity, pulse shape analysis (PSA) and RLY cut efficiencies we obtain the exclusion sensitivity curves reported in Fig. 18. If we demonstrate a background index of 10 −2 counts/(keV×kg×yr) with increased statistics, CUPID-Mo reaches a sensitivity superior to the most recent limit on the 100 Mo half-life set by NEMO-3 [34] in just 6 months of accumulated livetime. Figure 18 also reports the sensitivity for the more optimistic scenario where the background level is 10 −3 counts/(keV×kg×yr); in this case, the experiment is practically background free for a total of 1 year of livetime reaching a final sensitivity of T 0νββ 1/2 = 2.43 × 10 24 year. The exclusion sensitivity has a very minor dependence on the detector energy resolution  Fig. 18 Bayesian exclusion sensitivity at 90% CI for 5 keV resolution (FWHM) at Q ββ and different background levels and decreases by ∼ 10% for a factor two worse resolution and a background index of 10 −2 counts/(keV×kg×yr).

Conclusion
The first physics data of the CUPID-Mo experiment validates and extends the previously reported bolometric performance for cryogenic Li 2 MoO 4 crystals [28,29] on a much larger array of 20 detectors. We find that crystal growth and detector assembly can be well controlled to obtain excellent uniformity in performance and radiopurity. In particular, the summed energy resolution was 5.3 keV (6.5 keV) FWHM at 2615 keV in calibration (physics) data of 19 out of 20 detectors. The measured light yield for γ /β events (0.6-0.9 keV/MeV), the quenching of the scintillation light for α particles (20%) with respect to γ /βs and the achieved baseline resolution of bolometric Ge light detectors (146 eV FWHM) are compatible with full α to γ /β separation (median discrimination power value of 15). The Li 2 100 MoO 4 crystals also exhibit a high level of radiopurity, particularly ≤ 3 μBq/kg of 226 Ra and ≤ 2 μBq/kg of 232 Th. The results indicate the prospect to surpass the sensitivity of NEMO-3 with ∼6 months of physics data in the current demonstrator. The technology is scalable, and the first results presented in this article strengthen the choice of Li 2 100 MoO 4 as the baseline option for application in the CUPID next-generation cryogenic 0νββ experiment. Additional data from the current and future demonstrators is essential to develop a detailed background model, investigate and optimize the performance in the region of interest for 0νββ, and to further strengthen the projections for CUPID.