A flexible scintillation light apparatus for rare event searches

Compelling experimental evidences of neutrino oscillations and their implication that neutrinos are massive particles have given neutrinoless double beta decay a central role in astroparticle physics. In fact, the discovery of this elusive decay would be a major breakthrough, unveiling that neutrino and antineutrino are the same particle and that the lepton number is not conserved. It would also impact our efforts to establish the absolute neutrino mass scale and, ultimately, understand elementary particle interaction unification. All current experimental programs to search for neutrinoless double beta decay are facing with the technical and financial challenge of increasing the experimental mass while maintaining incredibly low levels of spurious background. The new concept described in this paper could be the answer which combines all the features of an ideal experiment: energy resolution, low cost mass scalability, isotope choice flexibility and many powerful handles to make the background negligible. The proposed technology is based on the use of arrays of silicon detectors cooled to 120 K to optimize the collection of the scintillation light emitted by ultra-pure crystals. It is shown that with a 54 kg array of natural CaMoO4 scintillation detectors of this type it is possible to yield a competitive sensitivity on the half-life of the neutrinoless double beta decay of 100Mo as high as ~10E24 years in only one year of data taking. The same array made of 40CaMoO4 scintillation detectors (to get rid of the continuous background coming from the two neutrino double beta decay of 48Ca) will instead be capable of achieving the remarkable sensitivity of ~10E25 years on the half-life of 100Mo neutrinoless double beta decay in only one year of measurement.


Introduction
In this paper we propose a novel approach to neutrinoless double beta decay searches which uses high performance solid state detectors to read the scintillation light emitted by large mass crystals with high energy resolution.
Neutrinoless double beta decay is an extremely rare phenomenon hypothesized long time ago but never observed. Its discovery would be a major breakthrough in astroparticle physics, since it would demonstrate the lepton number non-conservation and unveil the Majorana character of the neutrino, i.e. prove that the neutrino is equal to its own anti-particle. At the same time it would allow to assess the absolute neutrino mass scale with high sensitivity and would help point us towards the proper extension of the Standard Model of Particle Physics. This justifies the large number of experimental programs devoted to the search of this decay with different techniques. On the other hand, all these programs are facing with the technical and financial challenge of increasing the experimental mass while maintaining the contribution of the spurious background at incredibly low levels.
The technique proposed in this paper is the combination in a single device of all the demanding features needed by next generation experiments: high energy resolution, low cost mass scalability, flexibility in the choice of the double beta decaying isotope, and many powerful handles to keep the background negligible. Its novelty is the enhancement and optimization of the collection of the scintillation light emitted by ultra-pure crystals through the use of arrays of high performance silicon photodetectors cooled to 120 K. This would provide scintillation detectors with 1% level energy resolution and a µs time resolution. We show that with a small array (∼ 54 kg) of such detectors it is possible to yield in only one year of measurement a competitive physics result in the search for the neutrinoless double beta decay of 100 Mo using natural CaMoO 4 scintillating crystals (i.e. without enriching in 100 Mo or depleting in 48 Ca isotopes). This result would pave the way to a new class of ton scale experiments to search not only for neutrinoless double beta decay with potentially zero background, but also for dark matter and neutrino charged current reactions.

Neutrinoless double beta decay
In recent years, the discovery of neutrino oscillations unambiguously proved that neutrinos are massive particles and thus provided the first confirmation of the existence of new physics beyond the Standard Model [1]. However oscillation experiments cannot shed light on the fundamental open issues concerning the absolute mass values or the quantum nature (Dirac or Majorana fermions) of neutrinos. Neutrinoless double beta decay is a lepton number violating process with the unique feature of being the only feasible means to provide insight into both the above questions. In fact its existence would imply that neutrinos are massive Majorana fermions and could put important constraints on the absolute mass scale [2,3].
Double beta decay (ββ) is a very rare nuclear process in which a nucleus (A, Z) decays into its (A, Z+2) isobar. As the expected decay rates are extremely low, the choice of the parent nuclei is limited to those which are more bounded than the intermediate ones (to avoid the background caused by the sequence of two single beta decays): this condition is met for a number of even-even nuclei. In the Standard Model of Particle Physics ββ is allowed with the contemporary emission of 2 electrons and 2 antineutrinos (ββ2ν), and it has indeed been observed experimentally in a dozen of isotopes with half-lives of the order 10 18 − 10 21 y. Non-standard decay channels appear whenever the Majorana character of the neutrino is allowed. In this case the lepton number is not conserved and neutrinoless decay modes are possible. Neutrinoless double beta decay (ββ0ν) can proceed via different mechanisms, the simplest one being the virtual exchange of a light Majorana neutrino between the two nucleons. In this case the decay rate is proportional to the square of the so-called effective Majorana mass | m ββ |: where T 0ν 1/2 is the decay half-life, G 0ν is the two-body phase-space integral, M 0ν is the ββ0ν Nuclear Matrix Element (NME), and m e is the electron mass. The product F 0ν N = G 0ν |M 0ν | 2 includes all the nuclear details of the decay and it is usually referred to as nuclear factor of merit. While G 0ν can be calculated with reasonable accuracy, the NME value is strongly dependent on the nuclear model used for its evaluation. Significant improvements towards compatibility among various theoretical models have been achieved in recent years, even if discrepancies of about a factor 2-3 still exist. Fig. 1 summarizes the current situation [4,5,6,7,8,9,10,11,12,13].
The effective Majorana mass is a coherent sum over mass eigenstates with complex Majorana phases and it can be parametrized as a function of neutrino oscillation parameters. The expected allowed ranges for | m ββ | as a function of the lightest neutrino mass are depicted in Fig. 2, where the solid lines are obtained by taking into account the errors of the oscillation parameters (at 3σ level).
Experimentally, ββ0ν searches rely on the measurement of the two electron signal. In the so-called homogeneous or calorimetric experiments the active volume of the detector contains the double beta decaying isotope (source = detector). Since the energy of the recoiling nucleus is negligible, the sum kinetic energy of the two electrons is equal to the Q-value of the transition (Q ββ ). This  The left panel shows the expected ranges for | m ββ | as a function of the lightest neutrino mass. The two sets of solid lines correspond to the Normal (lighter line with horizontal hatch) and Inverted (darker line with diagonal hatch) Hierarchies, which merge for | m ββ |≥ 0.1 eV into the degenerate mass pattern. The width of the hatched areas is due to the unknown Majorana phases and thus irreducible. The central and right panels show the dependencies of | m ββ | as a function of the summed neutrino mass and of the kinematical mass | m β | (reproduced from [1]). monochromatic signal, is the main signature used by experiments. On the other hand, given the rarity of the process, the gathering of the ββ0ν counts is hindered by the background events in the energy region under investigation. The sensitivity of a given ββ0ν experiment is usually expressed in term of the detector factor of merit F 0ν D , defined as the process half-life corresponding to the maximum signal that could be hidden by a background fluctuation n B at a given statistical Confidence Level (C.L.). At 1σ level this is given by: with n B = √ BT M ∆, where B is the background level per unit mass, energy, and time, M is the detector mass, T is the measuring time, ∆ is the FWHM energy resolution, N ββ is the number of ββ decaying nuclei under observation, x is the stoichiometric multiplicity of the element containing the ββ candidate, η is the ββ candidate isotopic abundance, N A is the Avogadro number, A is the compound molecular mass, and is the detection efficiency. Finally, f (∆) is the "analysis efficiency", i.e. the fraction of signal events that fall in an energy window equal to the FWHM ∆ centered on Q ββ (Region Of Interest or ROI). The experimental sensitivity F 0ν D is then translated into an effective Majorana mass sensitivity as | m ββ | −1 = F 0ν D · F 0ν N showing both the slow dependence (only with the fourth root) of | m ββ | on the crucial experimental parameters and how the isotope choice plays a role in the game from the nuclear model point of view.
Apart from one controversial claim of evidence (in tension with several recent experimental results [14,15,16]), no ββ0ν signal has been detected so far. Experimental half-life lower limits exist for many isotopes. Current best sensitivities have been obtained for 76 Ge, 82 Se, 100 Mo, 130 Te, and 136 Xe, and are at the level of 10 24 − 10 25 y, which translate to | m ββ | values between 0.2 and 0.7 eV, depending on the ββ0ν candidate and on the corresponding NME (see [17] for a recent review).
To improve these results several experimental programs are currently on the way with a variety of challenging techniques and with the common goal of approaching the Inverted Hierarchy (IH) region of the mass spectrum (see Fig. 2) in the next years. They require not only very strong efforts towards the reduction of the radioactive background in the region of interest, below the intrinsic limit set by the ββ2ν continuous spectrum, but also good energy resolutions (approximately better than 2% at Q ββ [17]). In fact, an approximate expression for the ββ0ν signal S to ββ2ν background B ratio in the ROI can be written as [2]: where R = ∆/Q ββ is the fractional FWHM resolution. It is then clear that a bad energy resolution can spoil all background reduction efforts, because of the broadening of the upper tail of the ββ2ν continuous spectrum. Table 1 reports the predicted count rates due to ββ2ν in a ROI of 50 keV (therefore assuming a resolution of 1.7% at 3 MeV) as well as the number of expected ββ0ν signal events for two different values of | m ββ | (50 and 10 meV, 3σ boundaries of the IH region): the signal S to background B ratio is safely bigger than one for all candidate isotopes. Finally, energy resolution plays a unique role also when discussing the discovery potential of any given experimental approach, since it is a crucial parameter to single-out the ββ0ν peak over an almost flat background. When the background rate is so low that the expected number of spurious events in the region of interest is close to zero (ZB -Zero Background condition), the detector factor of merit becomes: where the background is now described by n L , a constant term which represents the maximum number of counts compatible, at a given confidence level, with no counts observed [1]. In this case the sensitivity to the ββ0ν signal goes linearly not only with the isotopic fraction and the detection efficiency, but also with the detector mass and the measuring time, a particularly appealing regime for next generation experiments. From the above discussion the strategy for future experimental programs aiming at exploring the IH region is clear: • minimization of the environmental and cosmogenic radioactive background (experiments placed underground with heavy shielding, careful selection of radiopure materials to build the experimental set-up, detector technology with active event identification capability and/or event position reconstruction); • large detector mass, on the tonne scale or larger, and mass scalability; • a well performing detector (good energy resolution and time stability), with flexibility in the choice of the ββ0ν isotope; • a favourable ββ0ν candidate isotope (high Q ββ value, high natural isotopic abundance, short half-life, large ratio of ββ0ν/ββ2ν in the ROI) • a reliable and easy to operate detector technology requiring little service on long run times; Any of the proposed new projects [17] is a compromise among these often conflicting requests (see also Table 2). We show that large arrays of inorganic scintillators could address all the above items in a single approach.

The benefit of scintillation detectors
In realm of ββ0ν calorimetric experiments, many techniques are available. So far Ge diode detectors and bolometers have provided excellent energy resolutions [18,19,20,21], at the level of few parts per thousand, and both offer promising solutions for background reduction. On the other hand, both techniques suffer from limitations in the achievable mass increase (mainly for cost reasons in the case of Ge diodes, and also technological issues in the case of bolometers). On the contrary, large scale organic scintillation detectors provide mass scalability in low-background environment without apparent limitations, but suffer from a seemingly irreducible restriction in the attainable energy resolution [15,22]. Finally, Xe (liquid or gas) based TPC [14,23] are bound to a trade off between source mass and energy resolution, although they may implement powerful, but complex, background abatement concepts. The possibility to implement in a single approach all the critical features -high energy resolution, background abatement techniques and mass scalability - Table 1. Properties of several ββ0ν candidates: Q-value, isotopic abundance, ββ2ν half-life. R2ν is the rate due to ββ2ν in a FWHM of 50 keV centered around Q ββ and is calculated using the Primakoff and Rosen approximation [24]. R0ν is the range of ββ0ν count rates (from the different NME values of Fig. 1) expected for two values of | m ββ | (the effect of the analysis efficiency f (∆) has been included). Both R2ν and R0ν are expressed in counts per year per ton of ββ0ν candidate isotope.
is the true challenge for next generation projects. Indeed, large arrays of inorganic scintillators may be the solution to that. Inorganic scintillating crystals can be grown with high level of intrinsic radiopurity and good scintillation properties (required for high energy resolution). Moreover, there exist many scintillators containing ββ0ν candidate nuclei, thus allowing a high flexibility in the choice of the scintillation detector (i.e. a wide range of ββ0ν candidates) without the limitation imposed by other experimental techniques. Furthermore, thanks to Quenching Factors 1 (QF) for α particles as small as 0.2 -which practically remove all α background from the energy region around Q ββand intrinsic α/β discrimination ability (Pulse Shape Discrimination or PSD) -useful to tag α − β coincidences -scintillators offer powerful tools for background minimization. Finally, large arrays of scintillating crystals can be assembled with quite simple technological skills and in rather low-cost experimental setups. Despite all those features, this approach has not exploited its full potential so far, the main reason being the relatively poor energy resolution achievable with scintillating crystals optically coupled to photomultiplier tubes (see next Section). We discuss in this paper an innovative method to read the scintillation light emitted by large mass scintillators containing the ββ0ν candidate nuclei without degrading 1 The quenching factor is here defined as the ratio of the light yield for an alpha particle or a nuclear recoil to that of an electron or a gamma. the energy resolution. In fact, the limit imposed to the scintillator resolution by the carrier statistics alone is not so crucial. Assuming to have a crystal that emits 10,000 photons per MeV -a value definitely reasonable for many scintillators [25] -the relative FWHM resolution R stat at an energy of 3 MeV (close to the Q ββ of many isotopes) in the hypothesis of being able to collect all the photons can be estimated as R stat (FWHM) = 2.355/ √ 3 · 10000 = 1.36%, which corresponds to a FWHM energy resolution ∆ of 41 keV at 3 MeV. Even if this value is not comparable to the resolution of Ge diodes or bolometers (which is of the order of 4-6 keV at 3 MeV), it is not a limiting factor in the case of a low background ββ0ν experiment (see also Table 1). In the following we will demonstrate that it is possible to maintain this resolution performance even with the complete readout chain. The great advantage of this new method is the availability of a homogeneous detector with the strongest potentialities towards background reduction and the possibility to indefinitely scale the experimental mass to boost the sensitivity towards unexplored regions of the neutrino mass spectrum. This detection concept is compared to the other techniques in Table 2.

The measurement technique
The basic idea of the proposed new technique is to optically couple suitable scintillating crystals to Silicon Drift Detectors (SDD) at a working temperature around 120 K. Indeed, many scintillating crystals have been used at low temperatures showing improved light emission [26] with respect to room temperature. On the other hand, in the recent years there have been many advances in the light detection using solid state devices [27]. The innovative idea is to combine these features to provide the winning technology for future experiments. Table 3 lists a selection of scintillating crystals containing interesting ββ0ν candidate nuclei with transition energy above the 208 Tl γ peak (the most energetic and intense γ line from natural radioactivity) -crucial issue to fulfill the ZB condition. All of them have fairly well known scintillation properties at low temperatures [26,25] and offer Light Yields (LYs) well above the 10000 ph/MeV value considered in the previous section. Therefore, all of them are potentially good candidates for high resolution detectors. It is worth stressing here that the development of inorganic scintillators is a field of continuous innovation that keeps on offering new crystals with competitive characteristics for many applications -see e.g. [25].
The key point is to avoid a degradation of the resolution (i.e. a loss of carriers) when the crystal is coupled to a proper device to detect the scintillation signal. For example, the use of commercial phototubes introduces strong limitations because of low quantum efficiency Qof the order of 40% in the best cases -and rather narrow spectral response, which implies that many of the crystals of Table 3 do not emit light within the wavelength intervals required to match the phototube characteristics, thus causing a further deterioration of the energy resolution. Moreover it is very difficult to obtain commercial photubes with a low intrinsic radioactivity of the glass envelope. The use of SDDs as photodetectors represents therefore an appealing possibility, since these devices are characterized by a Q larger than 80% in a wide range of wavelenghts [27]. In addition, silicon is one of the cleanest materials from the radioactivity point of view. To evaluate the energy resolution potentially achievable by coupling SDDs to scintillating crystals, we must consider the additional contribution of the noise of the measurement chain: where ENC e is the equivalent noise charge of the electronic chain and N ph is total number of photons emitted by the scintillating crystal after an energy deposition. A competitive detector therefore requires to reduce as much as possible the term R noise . In practice, one should keep this term at least within a tenth of the statistical one: ENC e = Q · N ph /10 = √ 0.8 · 30000/10 = 16 e − rms, where a LY of 10,000 ph/MeV at an energy of 3 MeV has been considered. This value of ENC e is quite demanding but is within reach of presently available technology, as we will discuss in the following Section.
Another possible approach that can be adopted to readout the signal emitted by the scintillating crystals is the use of a different solid state device: the Silicon Photomultiplier (SiPM). The optimization of this photodetector had a strong development in the last few years. The noise induced by such type of device operated at low temperatures is negligible; on the other hand the main limitations, at the moment, are given by the lower quantum efficiency with respect to the SDDs (due to the larger dead layers that are actually needed to realize a "large" surface detector) and by the strong non linearity that can be observed when measuring a large number of photons (due to the intrinsic avalanche signal production in a single cell).
Since the optimization of SiPM is expected to continue in the future, it is not excluded that also these devices will became suitable for a rare event search with the same performances that we point out in this paper for the SDDs.

Silicon Drift Detectors as photodetectors
The keystone for the success of the proposed technique is the achievable noise performance of the SDD detectors. Considering a solid state device with capacity C d coupled to a charge preamplifier with input capacitance C i and series noise e 2 w , the equivalent noise charge can be written as [28]: where C p is the parasitic capacitance between the detector and the amplifying circuit, I l is the leakage current both at the input of the gate of the FET preamplifier and of the solid state detector, q is the electron charge (1.6×10 −19 C), and τ is the shaping time of the acquired signal. k 1 , k 2 and k 3 are parameters related to the type of shaping form used in the electronic chain (Gaussian, triangular, trapezoidal, etc.) -at this level all of them can be assumed of unit value, also in view of a possible application of an optimum filter procedure [28] as it is actually done in many experiments. The three terms in the ENC expression are the so called series noise, 1/f noise, and parallel noise, respectively. As the above formula shows, the achievement of low levels of ENC requires a careful design and selection of all components: the detector, the amplifier front-end, and the coupling between the two. In particular: • The contribution of the 1/f noise can be kept at negligible levels compared to the other terms through a proper choice of the input device for the preamplifier (for example, a JFET) [29]. • SDDs are characterized by a very low capacitance C d of the electrode collecting the signal charge [30], typically of the order of 0.5-1 pF/cm 2 ; moreover, the frontend electronics can be integrated directly on the same wafer of the photodetector [31], thus minimizing C p . • The leakage current is one of the most critical parameters. Its value is related to the size of the solid state detector and is approximatly proportional to the area of the detector. Typical values at room temperature are of the order of 1 nA/cm 2 for many silicon devices.
Since a desirable order of magnitude for the detector surface -to collect the light from a scintillating crystal devoted to ββ0ν search -is of tens of square centimeters, this would result in a leakage current of several nA that would spoil the ENC. On the other hand, measurements on SDDs [27] show a marked decrease in the leakage current with decreasing temperature. In principle this decrease should have an exponential trend: therefore the leakage current could become negligible by reducing the operating temperature, although carrier freeze-out sets a lower limit around 77 K [32]. An extrapolation below −25 • C of the experimentally observed behaviour of I l in SDDs shows that a working temperature around 120 K should be enough to have a leakage current as low as 10 −14 A per cm 2 , value observed for many semiconductor devices at this temperature. Moreover, at low temperatures also the leakage current of the input JFET of the preamplifier becomes negligible: therefore it is advantageous to place this component at the cold stage of the experimental setup too. Another advantage of placing both the detector and the preamplifier in close vicinity is that the parasitic capacitance C p can be reduced to small values (C p ∼ 0.5 pF) therefore diminishing its contribution to the first two terms of ENC. • The series noise of the preamplifier is given by: e 2 w = 2k B T α/g m where k B is the Boltzmann constant, T is the operating temperature, g m is the transconductance of the input JFET, and α=0.7 for an ideal JFET. The value of the transconductance increases with the JFET area, therefore by increasing the size of the device it is possible to achieve a sensible reduction of e 2 w . On the other hand this would cause an increase of both the leakage current and of the input capacitance C i . Also in this case a low operating temperature [28] helps the optimization of the device geometry keeping a small area to reduce I l and C i while controlling e 2 w with a proper choice of the working temperature, thus achieving a reduction of the total ENC. Typical values of e 2 w at temperatures of about 120-150 K are of the order of 10 −18 V 2 /Hz while for the input capacitance we can assume C i ∼ 1 pF.
In conclusion, an SDD of 1 cm 2 operated at a temperature of about 120-150 K and coupled to a JFET closely placed are a viable solution for reading out the light emitted by a scintillating crystal with high resolution. To evaluate the expected total noise ENC e associated to such an electronic chain we can choose a shaping time of about 50 µs, an acceptable value for the typical decay times of inorganic scintillators, especially at low temperatures [26]. With this selection of τ , together with the values of the capacitances, e 2 w and I l just inferred, the contributions of the series and of the parallel noises to the total ENC e are comparable: This level of noise is for one solid state photodetector of 1 cm 2 of total area. On the other hand, for applications which need to cover large areas with light sensitive detectors -like the one proposed in this paper -bigger size SDDs would be necessary. A possible solutionwhich may have additional advantages for background reduction, as discussed later -is to use several independent small area SDDs to cover completely the surface of interest. The use of SDDs for segmented readout has already given interesting results for γ and X-ray detection [33,34,35]. A segmented optical readout has the additional benefit of keeping small drift times within each SDD, i.e. the time needed by the charge generated inside the SDD by an optical photon to reach the collecting anode. This parameter influences mainly the time response of the SDD. On the other hand, the number of SDDs (N SDD ) used to readout the emitted scintillation light has to be included in the evaluation of the overall energy resolution: assuming ∼ 40 SDDs per detector (see later), the additional contribution to the resolution amounts to √ N SDD ×ENC e = √ 40×3.1 ∼ 19.6 e − rms, very close the initial goal of 16 e − rms, and looks indeed extremely promising.

Scintillating crystal choice
Among the different scintillating crystals listed in Table 3, CaMoO 4 , CdWO 4 , and CaF 2 look very promising to build large scintillation detector arrays devoted to rare event searches. Their main physical properties are listed in Table 4. In order to give relevance to the discussion, in the following we will refer to CaMoO 4 crystals, even if it is clear that similar considerations are valid also for CdWO 4 , CaF 2 , and in general for any scintillating crystal with adequate physical properties.

CaMoO 4
CaMoO 4 is an intrinsic scintillator with a high molybdenum content (48% in mass), thus extremely attractive for the search of the ββ0ν of 100 Mo, which in turn is one of the most interesting isotopes because of its high transition energy and relatively high natural isotopic abundance (Table 1). Moreover the relatively low Z of CaMoO 4 reduces its efficiency for γ detection, with the largest γ contribution to counts in the ROI coming from the 3143 keV 214 Bi  [41]. Another important feature of CaMoO 4 is a measured QF of 0.2 [42], which prevents the dangerous alpha background from contributing to the ROI. Moreover, CaMoO 4 offers a pulse-shape discrimination (PSD) capability between γ/β events and α particles, which is very useful in the understanding and reduction of background spectra (see later): in [42] the demonstrated discrimination capability was at the level of ∼ 90%, but it was recently measured to be higher than 99.9% [41]. A drawback in the use of natural CaMoO 4 as detector for the search of the ββ0ν of 100 Mo is the presence of the isotope 48 Ca, which is also a ββ0ν candidate with a high transition energy (see Table 1). Despite its very low isotopic abundance, the expected count rate due to the continuous spectrum of the ββ2ν of 48 Ca in the ROI (∼ 50 keV) of 100 Mo is at the level of ∼ 10 −2 counts/(keV·kg·y). This must be taken into account when conceiving a ββ0ν experiment aiming at exploring the IH region of the neutrino mass spectrum. In this respect, recent developments in the production of 40 Ca 100 MoO 4 crystals, isotopically depleted in 48 Ca and enriched in 100 Mo, with scintillation properties very similar to those of the natural crystals [43] already offer a viable solution to this problem.

Expected detector performance
A prototype detector single module to prove this new measuring technique can be conceived as made of a cylindri-cal scintillating crystal with the two flat circular surfaces optically coupled to two arrays of SDD photodetectors. Scintillating crystal indicative dimensions could be 5 cm of diameter and 6 cm of height. We imagine each crystal to be laterally covered with a reflective sheet or paint to improve the light collection, and coupled to the SDD array by means of a suitable optical grease. Given the above crystal dimensions, about 40 exagonal SDDs of ∼1 cm 2 of area are needed to completely cover the two circular surfaces of the cylinder 2 .
To evaluate the expected energy resolution of a single detector module, besides the contributions of the statistical term and of the electronics noise, it is necessary to take into account also the capability to collect all the generated photons on the entrance window of the solid state photodetector. For this purpose, a dedicated Monte Carlo simulation based on the GEANT4 toolkit [44] was developed with the main goal of understanding the dependence of the photon collection efficiency on the type of coupling to the scintillating crystal. In the simulation, the scintillator is described as a cylindrical crystal of CaMoO 4 with the above dimensions and the characteristics listed in Table 4. The lateral crystal surfaces can be chosen to be polished or covered with reflecting materials to vary the probability that a photon is refracted outside the crystal and then lost. On the two circular surfaces of the cylinder it is possible to put a coupling layer of 0.1 mm of air (refraction index 1.003) or optical grease (refraction index 1.6) and, additionaly, a SDD detector with a refraction index of 1.45 and the Q shown in Fig. 3. The preliminary simulation results are reported in Table 5. They show that the fraction of collected light is strongly enhanced by the presence of a grease interface between the crystal and the silicon device and by the covering of the lateral crystal surfaces with a diffused reflector. The fraction of photons generated inside the CaMoO 4 crystal which reach the SDD photodetector seems very close to 60% (including the SDD Q ), roughly confirming data available in the litera- ture [45]. These results can be further optimized by using a more precise design of the photodetectors, of the optical coupling to the scintillator, and of the crystal lateral surface treatment. Finally, the FWHM resolution attainable at an energy of 3 MeV with the conceived single detector module operated at 120 K can be evaluated as: where α ph ∼ 70% is the light collection efficiency just discussed, N ph ∼ 75,000 is the total number of photons generated inside a CaMoO 4 crystal operated at 120 K after an energy deposition of 3 MeV, Q ∼ 80% is the SDD quantum efficiency, and √ N SDD ×ENC e ∼ 20 e − rms is the equivalent noise charge value calculated above. In the general expression of the total FWHM resolution it has been introduced also the contribution of R int which is the intrinsic resolution of the scintillator and is related to several effects, like inhomogeneities due to the local variation of the light output in the scintillating crystal, variations of the reflectivity of the diffuse reflector surrounding the scintillator, as well as non-proportionalities of the scintillator response [46,47]. At the energies of interest for ββ0ν searches (∼ 3 MeV) the contribution of R int can be completely neglected, but this assumption must be verified at energies lower than 1 MeV. Therefore the predicted FWHM resolution at 3 MeV of a CaMoO 4 detector coupled to 40 SDDs which read the scintillation light is of about 37 keV: this is a very good value indeed, that allows the minimization of the 100 Mo ββ2ν contribution to the ROI as well as a real discrimination capability of the expected ββ0ν sharp peak over the background events.
The attainable time resolution of such a detector module is directly connected to the scintillation decay time constant and the total photon emission (Hyman theory [48]) through: t 2 = 2τ r τ d N (t)/N phe where N (t) is the total number of photoelectrons needed to overcome the trigger threshold, N phe is the total number of photoelectrons collected by the photodetector, τ d is the scintillator decay time constant, and τ r is the photodetector rise time.
In the case of SDDs, τ r is the time needed to drift the electrons towards the anode (∼ 1µs). Therefore, using the the data reported in Table 4 for CaMoO 4 at a temperature of 120 K, and assuming a threshold six time larger than the evaluated rms noise of 20 e − , a time resolution of ∼ 1µs can be expected (for 1 MeV energy deposition), completely dominated by the SDD drift time.
Finally, an additional potentiality of this type of detectors is the possibility of reconstructing the event topology by exploiting the segmentation of the light readout. In principle the photodetector matrix placed on the scintillator circular surfaces should be able to distinguish a single-site event (like the energy deposition by electrons, as in the case of ββ0ν) from a multiple-site event (like an energetic γ interaction) of similar energy. Preliminary Monte Carlo simulations -which do not include the optical photon propagation yet -showed that among the γ interactions of initial energy equal to 3 MeV depositing the whole energy in a single CaMoO 4 crystal, only 30% would be labelled as single-site event if the threshold distance for distinguishing a multiple-site interaction is set equal to 10 mm; this fraction becomes as low as 14% if the minimum detectable inter-site distance is 5 mm. These results are referred to interactions of single γs of 3 MeV: the abatement fraction of the multiple-site events may increase when the 3 MeV energy is deposited by coincident γs (e.g. coming from the 232 Th and 238 U chains or from showers caused by muons interacting in the external shielding). These are very preliminary remarks: if confirmed they would open new background discrimination opportunities for this type of detectors.

A prototype array of high performance scintillation detectors for ββ0ν searches
In this Section we investigate the expected performance of an array of CaMoO 4 detectors with SDDs' light readout for the search of the ββ0ν of 100 Mo. We suppose the full detector to be composed by 108 CaMoO 4 cylindrical crystals, 5 cm of diameter and 6 cm of height, arranged in four layers of 27 single modules for a total mass of 54 kg. We will consider two options: i) an array of 108 natural CaMoO 4 crystals, and ii) an array of 108 40 Ca nat MoO 4 crystals, isotopically depleted in 48 Ca. Each scintillator is laterally covered with a reflective sheet or paint to improve the light collection. The two circular surfaces of every cylindric crystal are optically coupled to an array of ∼ 20 SDD light sensors each. The 27 single detectors of any layer are lying side by side and mechanically supported by a thin copper grid with teflon gaskets. The tentative design is shown in Fig. 4. The full structure can be enclosed in a 5 cm thick cylindrical copper cryostat of ∼ 45 cm of diameter and ∼ 50 cm of height. The Cu vessel is vacuum tight and the operating temperature of 120 K can be provided by an external cryocooler thermally coupled to the detector array. This very simple experimental set-up already reflects many of the shrewdnesses required by a careful background control strategy. In fact, the only materials in close contact with the detectors are high-purity copper, teflon, and silicon, besides the wiring. The whole set-up can be surrounded by at least 20 cm of lead to shield the experiment from the environmental radioactivity. We imagine this experimental set-up to be located at the Laboratori Nazionali del Gran Sasso, Italy, where an average rock overburden of ∼ 3600 m.w.e. shields the experiments from cosmic rays.
In order to reach a competitive sensitivity in a ββ0ν search, as already discussed, it is extremely important to remove any spurious source of counts from the ROI. This is a tricky task, which absorbs most efforts of any experimental programs. Main sources of background are: i) radioactive contaminations of the detector and of the experimental set-up; ii) environmental muons, neutrons and gammas; iii) cosmogenic activation of the detector and set-up materials. The main strategies for background control include a careful selection of the construction materials, special treatments for reducing surface contaminations of the elements close to the detectors, adequate shielding against environmental sources, and limited exposure to cosmic rays during construction and transportation. The above conceived experimental set-up already addresses many of these requests. A GEANT4-based Monte Carlo simulation of the whole experimental set-up was performed to investigate the effect of the possible background sources on the ROI count rate. In the code, all detector elements can act as a radioactive source with the bulk and surface contaminations independently simulated; for the surfaces, different contamination depths can be chosen. The following possible sources have been considered in the simulation: , respectively -are of even greater concern, since the "quenched" α energy from the Po isotopes sums to the energy of the β emitted by the parent Bi (because of the very short Po half-life), thus contributing with a continuum spectrum to the ROI. The assumed bulk contamination limits for copper are available in [49]. For the CaMoO 4 crystals bulk, we have assumed the contamination levels reported in [42]. This is a very conservative assumption since there are promising purification methods of the raw materials needed for the crystal growing which allow to reduce the intrinsic radioactivity of the scintillators, as already demonstrated in [41]. • 232 Th and 238 U surface contaminations of the CaMoO 4 crystals and of the Cu (holder + tank). Surface contaminations very often exceed the background coming from the bulk of the materials, even in the most radiopure ones: this is supposedly caused by the mechanical and/or chemical treatments of the materials which may contaminate the surfaces or leave them in a highly reactive chemical state. The results of the simulation are presented in Table 6: they have been obtained with the detectors operated in anticoincidence 3 and by assuming a QF=0.2 for the α particles in CaMoO 4 . The time resolution of the detectors has been fixed at 1 µs, the FWHM resolution at (conservatively) 50 keV and the threshold for anticoincidence rejection at 50 keV.
The limits reported for the crystals in the last column of Table 6 are obtained with the aid of the α/β pulse shape discrimination in CaMoO 4 or by exploiting the technique of delayed coincidences within 232 Th and 238 U radioactive chains. In particular: 1. 208 Tl β decay (τ 1/2 = 3.05 m) induced background can be strongly reduced by tagging the α particle emitted by the parent nucleus 212 Bi (α decay branching ratio is 36%): the dead time resulting by vetoing of the detector for ten half-lives of 208 Tl after every 212 Bi α decay is negligible (< 0.5%). 2. The background induced by 214 Bi β decay (branching ratio = 99.98%) can be strongly reduced by tagging the α particle emitted by its daughter nucleus 214 Po: the dead time resulting by the removal from the recorded spectrum, of the data segment acquired during the ten 214 Po half-lives previous to each α decay is completely neglectable. 3. The α+β pile-up events of the 214 Bi → 214 Po cascades, given the good time resolution of the detectors, can be resolved in 99% of the cases. 4. Since the 212 Po half-life is much shorter than the time resolution of the detectors, we have assumed to be able to reject 99% of the α + β pile-up events of the 212 Bi → 212 Po cascades thanks to the α/β pulse shape discrimination.
As Table 6 shows, with the highly performing detectors proposed in this paper the overall background expected in a ROI of 50 keV centered around the Q ββ of 100 Mo is at the extremely encouraging level of 10 −4 counts/(keV·kg·y). This implies that in the hypothesis of an array of 54 kg of natural CaMoO 4 detectors, the limiting background is the unavoidable ββ2ν continuous spectrum of 48 Ca (about 10 −2 counts/(keV·kg·y)). Even in this case, it is thus possible to reach a sensitivity of ∼ 10 24 y on the half-life of the ββ0ν of 100 Mo in only one year of data taking at the Gran Sasso Underground Laboratory -to be compared with the presently available limit on 100 Mo (1.1 × 10 24 y [50]). In the hypothesis of an array of 54 kg of 40 Ca nat MoO 4 crystals (isotopically depleted in 48 Ca) operated in the same conditions, the estimated sensitivity on the half-life of 100 Mo ββ0ν after one year of data taking is ∼ 10 25 y.
This shows the tremendous potential of the proposed technique.

Other physics measurements with CaMoO 4
A ββ0ν experiment capable of reaching a high enough sensitivity to deeply explore the IH region of the neutrino mass spectrum will require a mass of the candidate isotope of the order of (or higher than) 1 ton. With an array of CaMoO 4 detectors of this size, interesting topics other than ββ0ν studies may be at reach. Dark matter may be investigated through the study of the modulation signal. The achievable energy threshold can be evaluated by considering the global light emission of the CaMoO 4 crystals operated at low temperatures coupled to the low electronic noise obtainable with the proposed SDD approach. As reported in the previous paragraphs, we estimate a global ENC around 20 e − rms in normal operating conditions, which means an estimated charge threshold around 100 e − considering a Gaussian distribution of the noise fluctuations. Taking into account the photon collection efficiencies reported in Section 5, it is possible to estimate an energy threshold of about 6-7 keV. This threshold value together with the large detector mass and the low background requested by the ββ0ν experiment, make the proposed CaMoO 4 detector suitable for a competitive Dark Matter experiment.
Neutrino oscillations and neutrino-nucleus interactions can also be studied through the scattering off Mo nuclei. Molybdenum isotopes, in fact, are a good target for neutrino interactions. The first suggestion of using a charged current reaction on 100 Mo to detect solar neutrinos dates back to the year 2000 [51] with a proposal of using a 100 Mo based detector (MOON experiment) both for ββ0ν and real-time solar neutrino studies. Recently charged current neutrino reactions on different Mo isotopes have been discussed [52] in view of their use for the detection of supernova neutrinos or for the study of neutrino-nucleus interactions, as suggested in [53]. As an example of the potentialities, with the cross sections calculated in [52] for neutrinos of ∼ 1 MeV scattering off 100 Mo nuclei, an interesting interaction rate of ∼ 0.4 d −1 is expected already on the 54 kg natural CaMoO 4 crystals array (discussed in Section 6) using the nominal rate for a 51 Cr ν e source of 370 PBq -like the one proposed in [54] -at 1 m distance. Obviously this rate will scale linearly with the number of 100 Mo nuclei, so a proportionally stronger signal is expected by increasing the detector mass and/or by isotopically enriching the CaMoO 4 crystals.

Conclusions
In this paper we have discussed a new powerful technology for ββ0ν searches with the potentialities of achieving zero background in the region of interest. Adequate radiopure scintillating crystals coupled to low noise SDD photodetectors result in a detector module with all the demanding characteristics required by future ββ0ν projects: high Table 6. Monte Carlo simulated contribution of the various background sources to the energy region centered on 100 Mo Q ββ for the full detector array (108 CaMoO4 detectors). The limits in the last column are obtained by applying off-line rejection analyses, as explained in the text.