Neutrino Masses and OscillationsView this Special Issue
Current Status and Future Perspectives of the LUCIFER Experiment
In the field of fundamental particle physics, the neutrino has become more and more important in the last few years, since the discovery of its mass. In particular, the ultimate nature of the neutrino (if it is a Dirac or a Majorana particle) plays a crucial role not only in neutrino physics, but also in the overall framework of fundamental particle interactions and in cosmology. The only way to disentangle its ultimate nature is to search for the neutrinoless double beta decay. The idea of LUCIFER is to combine the bolometric technique proposed for the CUORE experiment with the bolometric light detection technique used in cryogenic dark matter experiments. The bolometric technique allows an extremely good energy resolution while its combination with the scintillation detection offers an ultimate tool for background rejection. The goal of LUCIFER is not only to build a background-free small-scale experiment but also to directly prove the potentiality of this technique. Preliminary tests on several detectors containing different interesting DBD emitters have clearly demonstrated the excellent background rejection capabilities that arise from the simultaneous, independent, double readout of heat and scintillation light.
1. State of the Art
The double beta transition, which is the rarest nuclear weak process, takes place between two even-even isobars, when the decay to the intermediate nucleus is energetically forbidden due to the pairing interaction, which opens a gap between the even-even and the odd-odd mass parabolas in a given isobaric chain. The two-neutrino decay conserves the lepton number and was originally proposed by Goeppert-Mayer in 1935 . It is a second-order weak process—that explains its low rate—and it has been observed for a dozen of nuclei, with lifetimes in the range of 1018–1022 y .
Besides the two-neutrino decay, a much more intriguing process, the so-called neutrinoless double beta decay 0-DBD [3–6], was proposed by Furry  shortly after the Majorana theory of the neutrino . In this case, the simultaneous transformation of two neutrons into two protons is accompanied by the emission of two electrons and nothing else. The main feature of 0-DBD is just the violation of the lepton number. In the modern (standard model) perspective, this is as important as the violation of the baryon number. In full generality, we can imagine this process as a mechanism capable of creating electrons in a nuclear transition.
It is remarkable that 0-DBD is not necessarily due to the exchange of Majorana neutrinos (mass mechanism) as a leading contribution, even though its detection would demonstrate that neutrinos are self-conjugate particles . In spite of the wide variety of mechanisms possibly inducing the 0-DBD, after the discovery of neutrino flavor oscillations (which prove that neutrinos are massive particles), the mass mechanism occupies a special place. It relates neatly the 0-DBD to important parameters of neutrino physics, fixes clear experimental targets, and provides a clue to compare on equal footing experiments which present considerable differences from the methodological and technological points of view. In fact, the lifetime of the 0-DBD decay is related to the so-called effective Majorana neutrino mass , a crucial parameter that contains the three neutrino masses , , and , the elements of the first row of the neutrino mixing matrix, and the unknown CP-violating Majorana phases, which make the cancellation of terms possible; could be smaller than any of the three neutrino masses.
Majorana mass and weak isospin selection rules make it possible to find a natural explanation to the smallness of neutrino mass. The pattern of neutrino masses and mixing admits an elegant solution, the so-called see-saw mechanism. At the same time, the only viable explanation of the matter-antimatter asymmetry available today, the justification for our very existence, is based on the leptogenesis, which again requires a Majorana nature for neutrinos . This well motivates why 0-DBD plays a central role in particle physics and cosmology.
Thanks to the information we have from oscillations and assuming the standard three active neutrino scenarios, it is useful to express the effective Majorana neutrino mass in terms of three unknown quantities: the mass scale, represented by the mass of the lightest neutrino , and the two Majorana phases. It is then common to distinguish three mass patterns: normal hierarchy, where , inverted hierarchy where, , and the quasi-degenerate spectrum, where the differences between the masses are small with respect to their absolute values. We ignore the neutrino mass ordering at the moment, and the 0-DBD has the potential to provide this essential information.
In the standard interpretation of 0-DBD in terms of mass mechanism, three challenges are in front of the experimentalists who study this process. The first one consists in scrutinizing the much debated 76Ge claim ; recent experimental results [12, 13] have practically accomplished this task. The second one consists in approaching and then covering the inverted hierarchy region of the neutrino mass pattern. The third and ultimate goal is to explore the direct hierarchy region.
While scrutinizing the 76Ge claim can be done in principle with only 10 kg of isotope, we need typically 1 ton of isotope mass in order to explore the inverted hierarchy region, just to accumulate a few signal counts. The direct hierarchy region seems for the moment out of the reach of the present technologies, since one would need sources of the order of 1 Mmol (typically 100 tons).
Just having a large source is not enough; in order to appreciate such tiny signal rates, the experimentalists are obliged to operate in conditions of almost zero background. Acceptable background rates are of the order of 1–10 counts/y/ton if the goal is just to approach or touch the inverted hierarchy region, whereas one needs at least one order of magnitude lower values to explore it fully, around or even less than 1 count/y/ton.
The experimental strategy pursued to investigate the 0-DBD consists of the development of a proper nuclear detector, with the purpose of revealing the two emitted electrons in real time and collecting their sum energy spectrum as minimal information. Indeed, the shape of the two electron sum energy spectrum enables us to distinguish between the two discussed decay modes. In case of -DBD, this spectrum is a continuum between 0 and the -value with a maximum around . For 0-DBD, the spectrum is just a peak at the transition energy, enlarged only by the finite energy resolution of the detector.
The ideal desirable features of this nuclear detector are: (i) high energy resolution, as a peak must be identified over an almost flat background; (ii) large source, in order to monitor many candidate nuclides; and (iii) low background, which requires underground detector operation to shield from cosmic rays, very radiopure materials, well-designed passive and/or active shielding against local environmental radioactivity and possibly event tracking andtopology capability. This last point is useful not only to reject the background but also to provide additional kinematical information on the emitted electrons . Typically, the listed features cannot be met simultaneously in a single detection method.
The searches for 0-DBD can be further classified into two main categories: the so-called calorimetric technique, in which the source is embedded in the detector itself and which provides extremely high efficiency, and the external-source approach, in which source and detector are two separate systems, allowing excellent event reconstruction.
Which are the best isotopes to search for 0-DBD? Experimental practice shows that the following three factors weigh the most in the design of an experiment: (i) the -value; (ii) the isotopic abundance together with the ease of enrichment; and, last but not least, (iii) the compatibility with an appropriate detection technique. The -value is probably the most important criterion. It influences both the phase space and the background. On the basis of -value selection, at the moment there are only 9 experimentally relevant isotopes. The transition energies of all these isotopes are larger than 2.4 MeV, with the important exception of 76Ge (-value = 2.039 MeV). There are two important energy markers in terms of background which have to be compared with the -value: (i) the 2615 keV line represents the endpoint of the bulk of the natural gamma radioactivity; (ii) the 3270 keV line represents the -value of the 214Bi -decay, which, among the 222Rn daughters, is the one releasing the highest-energy ’s and ’s. The 9 candidates are divided according to these two markers into three groups of three isotopes each.
The first group (76Ge, 130Te, and 136Xe) has to cope with some gamma background and with the radon-induced one; however, very sensitive experiments can be performed with these nuclides since they are particularly suitable to be studied with a calorimetric approach. Germanium semiconductor diodes are excellent devices for the isotope 76Ge [11, 13]; TeO2 bolometers are at the forefront due to their high content of 130Te , while the last of these three isotopes can be easily embedded in gaseous or liquid TPC [16, 17] or in large volumes of a liquid scintillator .
The second group (82Se, 100Mo, and 116Cd) is out of the reach of the bulk of the gamma environmental background; as shown in the next section, they belong to the realm of the scintillating bolometers, which are almost perfectly suited to investigate these isotopes. LUCIFER’s experimental results described below show that at least two of them, namely, 82Se and 100Mo, can be studied with high sensitivity.
The candidates of the third group (48Ca, 96Zr, and 150Nd) are in the best position to realize a background-free experiment. For a sort of conspiracy of nature, however, these three golden-plated elements cannot be enriched at low cost and high throughput.
We are now (September 2013) at a turning point in the experimental search for 0-DBD decay. The 76Ge claim is strongly disfavored after the results provided by EXO-200 , KamLAND-Zen , and especially GERDA (phase I) , which investigates the same isotope as the claim and is therefore free from the systematics induced by the calculation of the nuclear matrix elements. The follow-up of these searches and others which are in an advanced construction phase (CUORE  and SNO+ ) promises to go well below 0.2 eV. Therefore, the first of the three aforementioned challenges is very close to be achieved, and it will be in the near future.
However, there is no univocal strategy to deal with the second challenge, namely, to explore deeply the inverted hierarchy band of the neutrino mass pattern. Searches in preparation can only approach the onset of this region at eV.
The place of LUCIFER in this context is very clear. This experiment does not foresee at the moment a sensitivity comparable to those of the most advanced experiments. In fact, LUCIFER is essentially an R&D activity, aiming at building a technology demonstrator at the 10–15 kg scale. However, it will face and hopefully solve all the issues related to a final search. Therefore, the LUCIFER results will have a major impact on the investigation for 0-DBD. All the experimental indications collected so far and reported in detail in this work show not only that the LUCIFER technology is viable, but also that the experiments which will adopt it are among the best positioned to achieve a background of the order of 1 count/y/ton, combined with the high energy resolution of the bolometric technique, the sensitive mass scalability guaranteed by the detector modular structure, and the high efficiency property of the calorimetric approach. They would have therefore all the numbers to achieve the second challenge in the search for 0-DBD decay.
2. Bolometers and Scintillating Bolometers
Cryogenic particle detectors (see, e.g., ) have been used for many years for the study of rare events such as 0-DBD and dark matter (DM) searches.
A cryogenic bolometer is simply made of a suitable absorber and a thermometer coupled with it. When a particle releases energy within the absorber, a tiny temperature rise can be measured with an appropriate thermometer, provided that the heat capacity of the absorber is reasonably small. This can be fulfilled even with a multi-kg absorber  provided that it is made of a dielectric crystal kept at cryogenic temperatures of the order of tens of mK.
Bolometers offer a wide choice of absorber materials, and at the same time, they are able to achieve an energy resolution competitive with that of Ge diodes (viz. of the order of 5 keV FWHM at 3 MeV). The freedom in the choice of the absorber provides the unique opportunity of selecting the DBD isotope without the limitations usually induced by the experimental technique (e.g., with semiconductor detectors).
The CUORE experiment represents the most advanced stage in the use of bolometers for 0-DBD search. CUORE will consist of an array of 988 crystals for a total mass of ≃1 ton of TeO2. Its data taking is expected to begin in 2015.
Unlike other solid-state devices, however, bolometers are not ionization detectors but phonon detectors. As a consequence, they are almost equally sensitive to any kind of particle, despite the way energy is released. In other words electrons, particles, and nuclear recoils—depositing the same amount of energy in the detector—produce a pulse with the same amplitude and shape. The Cuoricino  data demonstrated that the background in the region of interest is dominated by radioactive contaminations on the surfaces facing the detectors. particles produced by these contaminants can lose a fraction of their energy in the host material and the rest in the detector, thus producing a flat background from the energy of the decay (several MeV) down to the 0-DBD region .
The background induced by surface -emitters represents, in fact, the limiting background in the region of interest for this technique .
Scintillating bolometers, developed in recent years for DM searches [24, 25] and proposed also for 0-DBD surveys , allow overcoming this disadvantage by providing the possibility to distinguish interactions (background only) from interactions (background and signal). A scintillating bolometer is obtained by coupling a scintillating crystal, operating as a cryogenic bolometer (as described above), to a proper light detector (see Section 5).
When a particle traverses the scintillating crystal, a large fraction of the deposited energy is converted into heat (thus inducing a temperature rise), while the remaining small fraction is spent to produce scintillation light.
The interesting feature of scintillating bolometers is that the ratio between the two signals (light/heat) depends on the particle mass and charge. Particles like ’s and ’s have the same light emission (conventionally referred to as the light yield (LY), i.e., the fraction of particle energy emitted in photons) which is typically different from the light emission of particles or neutrons. Consequently, the simultaneous readout of the heat and light signals allows particle discrimination.
If the scintillating crystal contains a DBD candidate, the 0-DBD signal (i.e., the energy deposition produced by the two electrons emitted after the decay) can be distinguished from an signal, and only ’s and ’s can give a sizable contribution to the background that limits the experimental sensitivity. The feasibility of this technique is today widely proven. Scintillating bolometers containing Ca, Mo, Cd, and Se have been successfully tested, coupled to a thin Ge wafer operating as bolometer for the light readout [27–30]. Presently, those that look most promising for a 0-DBD large scale experiment are ZnSe and .
3. The LUCIFER Layout
The LUCIFER setup will consist of an array of individual single module detectors, arranged in a tower-like structure, as schematized in Figure 1. The single module (see Section 6) will consist of a cylindrical crystal of 45 mm diameter and 55 mm height, equipped with a Ge-crystal light detector (see Section 5). The crystal is fastened by means of four S-shaped Teflon pieces fixed to the two cylindrical Cu frames. The frames are held together through two Cu columns. The crystal is surrounded (without being in thermal contact) with a reflecting plastic foil (3M VM2002) in order to increase the light collection. Individual single modules allow us to test in advance each single detector before assembling the entire tower experiment.
The tower will be installed in the same dilution cryostat that hosted the Cuoricino experiment. The cryogenic layout and the readout electronics (see Section 4) are based on the experience gained in the Cuoricino experiment and implemented in the CUORE project.
The key point of the experiment will be the crystal growth made with enriched isotope, as discussed in Section 7. With the enriched material being extremely expensive, the losses during the various stages of the production have to be minimized. This R&D is still ongoing and could slightly change some parameter of the layout (number of crystals and their final dimensions).
The evaluation of the background and of the sensitivity foreseen for this layout, as explained in Section 9, takes advantage of the Cuoricino experiment that was able to characterize and simulate the various dangerous contaminants present in the cryogenic facility. As explained in Section 1, the main source of background arises from 226Ra that will produce 214Bi whose rare, but high, energies s can travel enough to reach the tower. With respect to 228Th, producing the 208Tl whose decay represents the real “nightmare” for all the DBD experiments, the situation is rather different. Even if the -value of the 208Tl -decay is the highest among the natural radioactive chains, in order to produce background above 2615 keV, the contamination must be present in the proximity of the detectors, so that summing effects can take place. For the LUCIFER project, all the materials (Cu holders, and thermal shields of the cryostat) within the first 5 cm from the detectors will be replaced with new ones that will be produced with already selected and tested Cu, whose 228Tl activity is found to be negligible.
4. The Electronic System
The thermometers used in LUCIFER are thermistors (see Section 8). These devices show a resistivity that changes exponentially with the temperature being below 100 mK. The readout scheme is, in principle, very simple; the thermistor is polarized by a simple bias circuit, consisting of a voltage generator closed on two load resistors; when a particle releases energy in the absorber, the temperature rises and the resistance of the thermistors changes, thus producing a measurable voltage drop. Given that the study of rare events such as 0-DBD is performed over a long period of time, the overall stability of the readout system plays a crucial role.
To fully exploit this stringent requirement and to match the DAQ dynamic range, a dedicated system has been developed, as schematized in Figure 2; the bolometer, Bol, is coupled with the thermistor, . It is inside the fridge, at a temperature between 10 and 20 mK.
Signals from bolometers (depending on their mass) show a signal bandwidth from a few tens of Hz up to a few hundreds of Hz. In case of relatively fast detector response (e.g., light detectors; see next section), a unity gain buffer and the load resistors are located inside the fridge [31–34], at a larger temperature this time, 120 K. This is the box Option in Figure 2, which is close to the detectors (a few tens of cm in terms of signal wire), to minimize parasitic capacitance. The box Option is not present for the readout of the slow (massive) detectors, and, in this case, the thermistor is connected directly to the box Con. Block, at room temperature, through a small thermal conductance electrical link , and, in this case also, the load resistors are at room temperature (Bias and Load Res.). All the electrical connections of the detectors consist of twisted pairs of wires, so that all the possible sources of common mode disturbances and crosstalk between nearby channels are minimized [36, 37].
The Con. Block of Figure 2 allows us to connect the preamplifier, PR, to the detector, to the ground, or to a test pulse. This block, together with Bias and Load Res., that can switch between 2 different sets of load resistances and set the level and polarity of the detector bias voltage, is exploited for DC characterization of the detectors to optimize their working points, detector by detector . The detector load resistances may have an impact on the thermal and low-frequency parallel noise. The minimization of their low-frequency noise, otherwise proportional to the square of the applied voltage, has been extensively studied  and addressed. Other noise sources originate from the PR. Also in this case, the input transistors, a pair of JFETs, have been selected with a semicustom approach . Moreover, in order to decrease the thermal drift of the preamplifiers, a special circuit was developed [41, 42] that ensures a maximum drift of the order of a fraction of μV/°C.
Remote programmability  is available for many parts of the electronic system, marked with a star in Figure 2. DC coupling is necessary not only for what concerns signal analysis but also for detector monitoring.
The second-stage amplifier (SS) buffers the outputs of PR with a programmable gain. The outputs of SS are connected to the DAQ system after the antialiasing filter BF that consists of a 6-pole-roll-off Thomson or Bessel filter. The antialising filter has a few programmable frequency bandwidths . The blocks before the BF are all located on the top of the cryostat, just a few cm away from it. The BF is instead located close to the DAQ, in a remote location, so as to reject every possible disturbance from the long connecting links.
The power supply of LUCIFER, as well as parts of the electronics system previously described, is an upgraded version of . The DC input to the power supply is derived from a 2-stage AC/DC system based on an AC (main line)/48 V DC and a 48 V/±11 V/±6 V/6 V.
5. Light Detectors
The first light/heat measurement was performed with a thermal bolometer and a silicon photodiode in 1992  but was no longer pursued due to the difficulties of running a “standard photodevice” at cryogenic temperatures. The use of a bolometer as a light detector (LD) was first developed in 1997  and further optimized for DM searches [24, 25] after a few years. The first developments for 0-DBD searches began in 2006 .
The LUCIFER light detectors consist of thin germanium slabs, operating as bolometers, facing the main scintillating crystal. The choice of such material is driven by several considerations. Being opaque semiconductors, they are sensitive over an extremely wide range of photon wavelengths and, moreover, satisfy the very stringent radiopurity requirements of rare event searches. Their overall quantum efficiency can be as good as the one of photodiodes. To further increase the light collection, a dark layer  is deposited on the surface of the Ge crystal that faces the bolometer. Each LUCIFER light detector consists of a disk-shaped pure Ge crystal ( = 44 mm, thickness = 180 μm) grown using Czochralski technique and purchased from Umicore. The temperature sensor is a neutron transmutation doped (NTD) Ge thermistor (see Section 8).
The scintillation photons, as well as other particles, interact in the LD and deposit their energy within its volume. The energy deposition is eventually converted into phonons, causing a temperature rise in the system and, in particular, in the NTD thermistor.
To allow proper calibration of the signal, a 55Fe source, producing two X-rays at 5.9 and 6.5 keV, is faced to the light detector. The choice of such low-energy calibration lines is due to the fact that light signals produced in scintillating bolometers have typical energies of the order of 10 keV.
As shown in , LUCIFER light detector performances have a clear dependence on the applied bias current, both in terms of energy resolution and signal time development (see Figure 3). As will be shown in Section 6.2, a very efficient particle discrimination can be obtained using the time development of the scintillation signal. Therefore, very often, the working point of the LD is chosen as a compromise between good energy resolution and fast detector response, obtained by tuning the polarization current.
6. Bolometric Investigation of the Crystals
6.1. Experimental Setup
As explained in Section 2, bolometers offer a wide choice in terms of the absorber material. An extensive R&D activity was carried out over these last years in order to understand the features of the various crystals, and several prototypes of different masses and production lines were investigated. All the analyzed crystals were operated as bolometers in the deep underground facilities of Laboratori Nazionali del Gran Sasso, Italy.
The detectors were installed in a 3He/4He dilution refrigerator (operating at about 10 mK) properly shielded from external radiations (’s/’s and neutrons) by means of lead and polyethylene shields, inside a Faraday cage continuously flushed with nitrogen. The experimental setup consists of a copper structure for the housing of the crystals and light detectors. In order to avoid the noise induced by the vibration of the cryogenic facility, the structure is mechanically decoupled  from the cryostat by using a two stage damping system . The thermal stability of the holder (whose variations are within a few tens of μK) is ensured by a completely custom proportional, integrative, and derivative (PID) feedback loop . In order to further increase the stability of each single detector, a resistor of kΩ, realized with a heavily doped meander on a 3.5 mm3 silicon chip, is attached to each bolometer and acts as a heater to stabilize the gain of the detector  by injecting “monochromatic” voltage pulses [54, 55] across the heater.
The heat and light signals produced by interacting particles into the absorbers are transformed into voltage pulses by NTD thermistors; then the signals are amplified and fed into an 18 bit NI-6284 PXI ADC unit. Software triggers ensure that every thermistor pulse is recorded with a kHz sampling rate over a s duration (depending on the crystal). Moreover, when a trigger is generated by the main crystal, the corresponding waveform from the LD is recorded, irrespectively of its trigger. The amplitude and the shape of the voltage pulse are determined by the offline analysis that makes use of the optimum filter technique [56, 57]. The signal amplitudes are computed as the maximum of the filtered pulse. The amplitude of the light signal is estimated from the value of the filtered waveform at a fixed time delay with respect to the signal of the bolometer, as described in detail in .
The energy calibration of the crystals is performed using removable sources placed outside the cryostat. The heat channel is calibrated, attributing to each identified -peak the nominal energy of the line, as if all the energy is converted into heat. Consequently, this calibration does not provide an absolute evaluation of the heat deposited in the crystal. It has to be remarked that in scintillating bolometers, the energy scale of / and particles is different. This is simply induced by the different scintillation yield; in standard scintillators, for example, / particles scintillate more than ’s with the same impinging energy. This energy, therefore, “escapes from the crystal” and, as a consequence, implies a different energy scale for a different type of particles . For this reason, the energy scale is, usually, labelled as keVee (electron equivalent), meaning that the energy scale is calibrated on / particles. The energy difference between the two scales is normally at the few % level, but for ZnSe it can reach larger values  that cannot be explained in terms of escaping light.
Because of its high content of Se (56%), as well as its good bolometric and scintillating properties, ZnSe has always been an interesting candidate for the search of the 0-DBD of 82Se ( keV).
The largest ZnSe bolometer ever realized, a 431 g crystal grown at ISMA, Ukraine, was recently characterized in terms of energy resolution, internal contaminations, and particle identification capabilities.
The FWHM energy resolution was found to be keV at 1461 keV and keV at 2615 keV. However, starting from the consideration that light and heat are two correlated estimates of the particle energy, we developed an analysis algorithm that allowed us to considerably improve the energy resolution by removing this correlation. The application of the algorithm, whose details are described in , resulted in a FWHM energy resolution of keV at 1461 keV and keV at 2615 keV.
Many runs were performed in order to assess the scintillation properties of ZnSe at cryogenic temperatures. In Figure 4(a), the detected scintillation light is reported as a function of the heat released in the crystals. In order to study the discrimination power between and events, a uranium source covered with a thin mylar foil was placed close to the ZnSe surface. The mylar foil degrades the energy of the emitted ’s, producing a continuum at low energy (red dots in the plot). Using a 228Th -source, (outside the cryostat) and the “internal” source, we were able to investigate the discrimination potential in the DBD energy region.
Looking at Figure 4(a), we observed that the events produced by the -source (blue) lie in a different region with respect to particles with the same energy. The LY of events resulted in keV/MeV, irrespective of the energy deposit. The LY of particles was studied separately for surface and internal contaminations, resulting in keV/MeV and keV/MeV. Several tests showed that the difference in the LY of bulk/surface events cannot be ascribed to self-absorption in the crystal or to energy dependence of the effects. Therefore, the origin of this behavior must reside in an effectively larger light production on the crystal surface or in nonuniformities in the light collection.
The larger LY of events with respect to in ZnSe crystals is not yet understood and, in principle, could affect the discrimination capability. A poor efficient light collection for particles, indeed, can generate event leaking in the band, increasing the background in the region of interest (see Figure 4(a)).
In order to overcome this problem, we developed an analysis algorithm which is sensitive to the difference in the pulse shape of ’s and ’s (see details in ). Each pulse was fitted using a model that takes into account the development of the phonon signals, the thermistor and electronics response, and the scintillation process, which cannot be considered instantaneous. Applying this algorithm to the heat and light pulses of ZnSe, we discovered that a poor particle discrimination can be obtained using the heat channel alone. The shape of the light pulses, on the contrary, is very sensitive to the type of interacting particles (see Figure 4(b)).
Thanks to the pulse shape discrimination, it is possible to identify and reject the interaction with very high efficiency, thus strongly reducing the background in the energy region of interest.
Finally, we performed a 524-hour background run to investigate the internal contaminations of the ZnSe crystal.
The study of the spectrum highlighted a contamination of μBq/kg in 232Th and of μBq/kg in 238U, both in equilibrium with their daughters.
The energy spectrum of this measurement, obtained using cuts on the LY and pulse shape, is reported in Figure 5. Besides the 40K and 208Tl peaks, due to the natural contaminations of the environment, only contaminations in 75Se ( = 119.8 d, -value = 863.6 keV) and 65Zn ( = 244 d, -value = 1359.9 keV) were found. The presence of these isotopes, due to the activation of 74Se and 64Zn, respectively, does not affect the background in the DBD region because of their short half-lives and low -values.
Above the 2615 keV -line, a single event can be observed. This event was in time coincidence with high energy ’s detected by other bolometers in the same setup, and, therefore, it was likely produced by a muon interaction.
In order to further improve the background suppression, a feasibility study of a muon veto for LUCIFER is ongoing.
A possible alternative to ZnSe is , which is an interesting candidate for the search of 0-DBD of 100Mo ( keV). This compound, synthesized in the last years, gained a lot of attention because of the excellent bolometric performances and discrimination capabilities.
In this paper, we summarize the results obtained by measuring a 330 g crystal , grown by the Nikolaev Institute of Inorganic Chemistry (NIIC, Novosibirsk, Russia), starting with high purity ZnO (produced by Umicore) and , synthesized by NIIC.
The FWHM energy resolution achieved with this bolometer ranged from keV at 583 keV to keV at 2615 keV. Similar results were achieved by other large mass crystals, demonstrating the reproducibility of the bolometric performances.
The internal contaminations of the crystals were investigated in a 524-hour background run, showing that can reach excellent radiopurity levels even without dedicating too much effort to the radiopurity controls. Only internal contaminations in 226Ra ( μBq/kg) and 210Pb ( μBq/kg) were found.
For the isotopes of 232Th, 228Th, 238U, 234U, and 230Th, upper limits of <8 Bq/kg, <6 μBq/kg, <6 μBq/kg, <11 μBq/kg, and <6 Bq/kg, respectively, were set.
The internal contaminations measured in other prototypes (even if with low exposure) were found to be compatible with the previous results [61, 62], showing that the radiopurity levels achievable with crystals match the requirements for a low background 0-DBD detector.
Finally, the scintillation properties of the bolometer were investigated in several calibration runs performed with a 228Th -source, a smeared source (similar to the ones described in the previous section), and an Am-Be neutron source to provide high energy ’s.
In Figure 6(a), the relative LY is reported as a function of the heat released in the crystal, showing that an excellent particle discrimination can be achieved. In this run, we measured keV/MeV, which is much larger with respect to keV/MeV (evaluated on the internal line of 210Po). The QF for particles obtained with this crystal is , in agreement with the values measured in smaller samples .
An interesting feature of can be observed in Figure 6(b), where the shape of the heat pulses is reported as a function of the energy deposited in the crystal. This plot shows that the background rejection provided by the pulse shape discrimination on the heat channel is very similar to the one given by the LY (Figure 6(a)). This feature is induced by the “long” decay constant of the scintillation signal  and has several advantages in view of a large mass experiment. Besides the reduction of costs and complexity of the setup, discarding the LDs will also reduce the mass of inert materials next to the detectors, thus possibly reducing the background.
7. Isotopic Enrichment and Crystal Growth
The first choice material for LUCIFER is ZnSe crystals grown from 82Se enriched raw material. ZnSe crystals are well known for their extended infrared (IR) transmission and for their scintillating properties and are currently produced for the use in IR optics, optoelectronic devices, and security control systems . ZnSe crystals’ growth, however, is quite difficult, and the constraints imposed on their use in DBD experiments add further complications to the obtainment of relatively large crystal samples. Preliminary cryogenic tests show that only ZnSe crystals grown from the melt can have the crystal perfection needed for good bolometric performances. The growth of ZnSe crystals from the melt is hampered by several factors.(i)The 1525°C melting point is relatively high and the high total vapor pressure of ZnSe (2 Bar at 1525°C) causes containment problems due to the high volatility of components (Zn and Se) resulting in vaporization and dissociation according to the following reaction: where the different vapor pressures of zinc and selenium lead to deviation from the stoichiometric compositions of the melt during melt growth. The growth from off-stoichiometric melt in nonreactive high pressure enclosures is applied in order to circumvent this problem.(ii)Melt grown crystals are prone to constitutional supercooling, resulting in local deviations from stoichiometry and formation of inclusions of the excess component. This effect is mainly due to the very sharp homogeneity range in the ZnSe phase diagram . (iii)The low thermal conductivity of ZnSe (15 W/m·°C) leads to a very difficult control of the growth interface and therefore to the formation of stresses in the crystal. This is further worsened by a relatively low stacking fault energy (10 erg/cm2), leading to the easy formation of stacking faults (change in the stacking sequence over few interplanar distances) and twinning (stacking faults over several atomic spacings). (iv)The solid-solid phase transition from wurtzite (hexagonal symmetry) to sphalerite (cubic symmetry) occurring in the proximity of 1400°C [65, 66] is another cause of extended defects (high dislocation density) in melt grown ZnSe crystals. The application of low cooling rates in the proximity of the phase transition reduces the formation of defects. The postgrowth annealing further contributes to the recovery of the crystal perfection.
The challenging obtainment of ZnSe crystal for DBD application is further hampered by the necessity to make a dedicated synthesis of the ZnSe compound starting from elemental Zn and Se (82Se enriched). The use of 82Se enriched material is required by the relatively low isotopic abundance (8.73%) in natural selenium. The enrichment, starting from natural selenium, is ongoing at URENCO, Stable Isotope Group in Almelo, The Netherlands. The starting material is natural produced by an external supplier. The technique is based on centrifuges, in a cascade recently updated and fully separated from 235U production in order to avoid any U and/or Th contamination of the final product. In addition, the centrifuges were flushed before use with fluoride, thus ensuring very good cleaning of the entire cascade. The enriched is further chemically processed to obtain elemental selenium . Care and control are applied at every step of the process from the procurement of natural to the delivery of enriched Se in order to guarantee the chemical and radiochemical purity of the final product . The chemical conversion of enriched is the most challenging step in the process of enriched Se production. The selected chemical reaction and corresponding equipment had to be compliant with very strict rules concerning radiocontamination prevention. In addition, the process has to be conducted in such a way to reduce the contamination by impurities defined as critical for the scintillation performance of ZnSe crystals (Fe, Cr, V, Ni, As, Cu, Mo, Si, and S). The chosen technology required the construction of a dedicated area at the URENCO site in Almelo where the dedicated to 82Se conversion rig will be operated.
Presently, the overall chemical purity turns out to be better than 99.8% on trace metal base; in particular, the concentrations of 238U and 232Th fall below 10−10 g/g and the critical impurities have concentrations below the accepted limits for a good scintillation performance of ZnSe crystals . The distribution of different phases of the very complex process of ZnSe crystal production is currently fixed, and the related production contracts are in progress or under discussion. The diagram of the production of enriched crystals is given in Figure 7.
The crystals discussed in this work were grown from the melt, in graphite crucibles using the Bridgman technique in a vertical furnace under inert gas (Argon) pressured up to 15 Bar in a temperature gradient with the maximum temperature reaching 1575°C in the melt zone. The growth rate (variable along crucible length) was 2–5 mm/hour. After growth, crystals were annealed in the same furnace and then cut and polished to the standard (LUCIFER) dimensions: 45 mm diameter and 55 mm length for a total volume of 87.4 cm3 and a mass of 460.7 g in the case of natural ZnSe. The expected mass of ZnSe crystals grown from enriched material (enriched selenium with 95% 82Se content) with standard dimensions is 469.98 g corresponding to 248.1 g content of 82Se. Table 1 gives the mass balance of ZnSe crystals grown from natural and enriched raw material.
The amount of enriched Se foreseen for the experiment is 15 kg. For an estimated yield of ZnSe synthesis of 97% (3% irrecoverable loss of enriched Se, i.e., 14.55 kg effective enriched Se embodied in 26.2 kg of enriched ZnSe powder) and an estimated yield of crystal growth and processing of 65% (35% irrecoverable loss of enriched ZnSe), it comes out that the total expected enriched crystal mass is 16.9 kg (8.9 kg content of 82Se) embodied in a total number of 36 crystals with standard dimensions. The relatively low yield of crystal production comes from the irrecoverable loss in evaporated material during the growth process and the losses during the mechanical processing (cutting, shaping, and polishing) of the as grown crystal to the final crystal sample.
8. Thermistor Production
Neutron transmutation doped (NTD) thermistors are semiconductor devices doped close to the metal to insulator transition. Melt-doped Ge crystals cannot achieve the necessary uniformity due to a variety of dopant segregation effects. The only technique available for producing a uniform doping is to produce the dopant by bombarding the Ge crystal with thermal neutrons. The most important aspect of this process is that 70Ge transmutes into 71Ga, an acceptor, and 74Ge transmutes into 75As, a donor, the primary active dopants in NTD Ge. The optimization of the detector performance is strictly related to the characteristics of the NTDs.
The thermal signal from ZnSe crystals and light detectors will be read out using NTDs that must be optimized for such specific application. The variation of resistivity of heavily doped thermistors as function of temperature, , is well represented by the formula where is a parameter that normally shows a value around 1/2 and and are directly correlated with the doping profile of the NTD. In particular, the more critical parameter is , due to its extremely strong dependence on the total neutron fluence.
Starting from our previous experience , there are precise indications for the characteristics of sensors that will give us the best performance in our experimental temperature range. In fact, we have chosen a of the order of 3.8 K that, together with the fixed dimensions of the thermistors (3 × 3 × 1 mm3), gives a working resistance of few MΩ at 15 mK.
In order to produce thermistors with the adequate features, we proceeded in two independent ways.
8.1. New Production
A completely new production has begun in 2012, starting from nine HPGE Ge wafers. The shape of the single wafer is close to a circular shape with mm and 3.2 mm thickness (see Figure 8).
All the neutron irradiations (ten independent runs) were performed at MIT Nuclear Reactor Laboratory between February 2012 and March 2013. Considering that a change in the total dose of neutrons of only 2.5% will change the value of by 1 unity (our acceptance window is between 3.3 and 4.3 K) and that the whole process of production (irradiation + cooling time) takes 2 years, we asked for four different doses: 3 wafers at −7% of the nominal dose, 2 wafers at −3% of the nominal does, 3 wafers at the nominal dose, and 1 wafer at the +2% of the nominal dose. From the single wafer, more than 500 thermistors can be obtained. The first cryogenic characterizations of the irradiated wafers will take place in October 2013.
8.2. NTD Fine Tuning
We have a previous “old” NTD wafer that is slightly underdoped, with K. This value, which is considered slightly too high for LUCIFER, can be corrected by reirradiating the sample with a total fluence of ~0.9·1017 n/cm2. This is not a large value but the new irradiation phase needs a very precise integral flux control. Due to the fact that MIT NRL cannot provide this fine tuning, we found an alternative nuclear reactor. To fulfill the previously mentioned requests, we characterized a selected neutron irradiation channel in the TRIGA Mark II nuclear reactor of the LENA (Laboratorio Energia Nucleare Applicata) Laboratory of Pavia University. After a long characterization of the neutron spectrum of the various irradiation channels  (in terms of flux and energy spectrum), the irradiation of the “old” NTD wafer was started in July 2013. Due to the cooling time, we will make the first test of those thermistors by the beginning of 2014.
9. Background and Sensitivity
The sensitivity of a given 0-DBD experiment is defined as the half-life corresponding to the signal that could be emulated by a background fluctuation of a chosen significance level, expressed in numbers of the Gaussian standard deviations () where is the stoichiometric coefficient of the DBD candidate, is the DBD candidate isotopic abundance, is Avogadro’s number, is the molecular weight of the active mass, is the background rate per unit mass and energy, is the detector mass, is the live time, is the FWHM energy resolution, is the detection efficiency, and is the fraction of signal events that fall in an energy window around the -value.
The sensitivity from (3) is then translated into an effective Majorana mass sensitivity where is the electron mass and is the nuclear factor of merit, defined as . The quantities and represent, respectively, the two-body phase-space factor and the 0-DBD nuclear matrix element (NME). While can be calculated with reasonable accuracy, the NME value is strongly dependent on the nuclear model used for its evaluation.
Equation (3) holds if the number of background counts is large enough so that its distribution can be considered to be Gaussian. In the limit of zero background counts in the region of interest (see e.g., ), the above formula reduces to
To calculate the sensitivity of LUCIFER to 0-DBD, we need to define the assumptions on the experimental parameters and quantities of (3), as well as the setup of the experiment. These assumptions are summarized in Table 2. The value of for is 0.76. On the basis of what was discussed in Section 3, we assume a tower-like detector of 36 crystals (ZnSe or ) arranged in a modular configuration of 9 floors with 4 crystals each, held together by a copper structure. Light detectors, consisting of Ge ultrapure wafers of about 180 μm thickness and 44 mm diameter, will be faced to each of the scintillating crystals, which are wrapped in a reflecting sheet to improve the light collection efficiency.
The detection efficiency () of ZnSe and crystals was estimated by means of GEANT4 Monte Carlo simulations and represents the fraction of 0-DBD decays, where the total energy of the two emitted electrons is contained within the crystal. The energy resolution (FWHM) is assumed to be 10 keV for ZnSe crystals and 5 keV for crystals.
The LUCIFER tower will be inserted into a dilution refrigerator made of few nested cylindrical copper thermal shields. A 10 cm thick lead disk, placed just above the tower, provides supplementary shielding against the radioactivity of the various components of the refrigerator located above the detectors: the dilution unit, the pumping lines, and the cabling system. Outside the refrigerator, a 20 cm (minimum) thick lead shield and a 10 cm thick borated polyethylene shield will be used to absorb ’s and neutrons from the environmental background of the LNGS experimental hall.
Generally, the background components for an underground DBD experiment can be divided into the so-called near sources (radioactive contaminations of crystals and copper mounting structures), far sources (radioactive contaminations of refrigerator and its external shields), and environmental sources (muons, neutrons, and gamma rays fluxes at the experimental site).
GEANT4 simulations of the various background components for a detector setup similar to what was described here have shown  that, exploring the background discrimination potential of the scintillating bolometers discussed in this work and given the radiopurity of the selected crystals and the detector’s materials presently available, a background index of 10−4 counts/kg/keV/y in the ROI is within reach.
However, as pointed out in Section 3, the detector will be mounted in the Cuoricino cryostat. This cryostat, unfortunately, was constructed in the late 1980s, and at that time, the materials’ selection requirements were not so stringent as nowadays. This implies that, even with changing the most internal Cu thermal and vacuum shields of the cryostat, the ultimate background will be dominated by the liquid helium Dewar in which the dilution unit is placed. The unavoidable background induced by 214Bi ’s generated in the Dewar is evaluated to be ≤1.5 × 10−3 counts/kg/keV/y.
For crystals, the dominant background contribution comes from an additional source with respect to the ones mentioned above: the accidental pileup of -DBD decays , given the slow time response of bolometric detectors and the relatively “fast” -DBD decay of 100Mo.
It has been shown, however, that a faster time response of light detectors can improve the pileup discrimination potential and reduce the background induced by 100Mo -DBD to well below 10−4 counts/kg/keV/y .
From the numbers in Table 2 and assuming a background index in the ROI of 10−3 counts/kg/keV/y, we can estimate the number of background counts in the ROI to be ~1 (~2) for a live time of 5 (10) yr. This means we cannot use the Gaussian approximation for the background fluctuations, but indeed we need to use the zero-background approximation.
In Table 3, we report the LUCIFER half-life background-fluctuation sensitivity at 90% C.L. in 5 or 10 yr, calculated using (5). From the value of , we extracted the corresponding limit on the effective Majorana neutrino mass that could be set by LUCIFER, considering the spread on most recent NME calculations [72–78] and the phase space factor from .
Moreover, LUCIFER can be considered a demonstrator of the scintillating bolometer technology with a significant mass and a full test of all the critical experimental issues of this technique and could lead the way towards a new generation of 0-DBD experiments with scintillating bolometers able to scrutinize the inverted hierarchy region of the Majorana neutrino mass spectrum.
In this paper, we have exposed a set of relevant results which show how the LUCIFER technology is now mature not only for the fabrication of a LUCIFER demonstrator at the kg scale (immediate objective of the LUCIFER program) but also for a much larger next-generation experiment.
The activity described here has allowed us to set the bases for the procurement of a considerable amount of 82Se, including for the first time a European company in the extremely exclusive group of potential producers of enriched isotopes for 0-DBD search. We have studied in detail the production chain which allows us to embed the enriched isotope into the final ZnSe crystals, basic elements of the 0-DBD detectors. The technical performances of the ZnSe scintillating bolometers, in terms of energy resolution and alpha-particle rejection factor, are fully compatible with an extremely sensitive 0-DBD experiment. The present indications on the background are also very promising and show that a background index of at least 10−3 counts/kg/keV/y is reachable.
We have presented also results on an alternative isotope and compound, that is, 100Mo embedded in crystals, which competes with the 82Se way in view of a next-generation experiment, offering better energy resolution thanks to superior crystal features. The radiopurity of the crystals studied so far is also quite promising.
Coming to more technical aspects, we have fully determined the structure of the light detectors and fabricated successful prototypes of these devices. We have also started the production of a considerable amount of neutron transmutation doped thermistors (the temperature sensors used both in the heat and in the light channels), largely sufficient to equip all the scintillating crystals and the light detectors foreseen in the LUCIFER demonstrator.
A detailed study of the background and of the sensitivity, assuming detector configuration and performance based on the features—reported here—exhibited by working prototypes, shows that the LUCIFER demonstrator can approach the inverted hierarchy region both in the ZnSe and in the versions. In the former case, the required amount of isotope is available. A well-founded extrapolation to the sensitivities of larger setups, limited essentially by the cost of the enrichment, indicates clearly that arrays of ZnSe and/or scintillating bolometers are viable candidates for a next-generation 0-DBD experiment, capable of fully exploring the inverted hierarchy region of the neutrino mass pattern.
This work was made in the frame of the LUCIFER experiment, funded by the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007–2013)/ERC Grant Agreement no. 247115. Thanks are due to E. Tatananni, A. Rotilio, A. Corsi, B. Romualdi, and M. Guetti for continuous and constructive help in the overall setup construction.
J. D. Vergados, H. Ejiri, and F. Simkovic, “Theory of neutrinoless double-beta decay,” Reports on Progress in Physics, vol. 75, Article ID 106301, 2012.View at: Google Scholar
M. Fukugita and T. Yanagida, “Barygenesis without grand unification,” Physics Letters B, vol. 174, no. 1, pp. 45–47, 1986.View at: Google Scholar
A. Gando, Y. Gando, H. Hanakago et al., “Limit on neutrinoless ββ decay of 136Xe from the first phase of KamLAND-Zen and comparison with the positive claim in 76Ge,” Physical Review Letters, vol. 110, no. 6, Article ID 062502, 2013.View at: Google Scholar
M. Agostini, M. Allardt, E. Andreotti et al., “Results on neutrinoless double beta decay of 76Ge from GERDA Phase I,” Physical Review Letters, vol. 111, Article ID 122503, 2013.View at: Google Scholar
M. Auger, D. J. Auty, P. S. Barbeau et al., “Search for neutrinoless double-beta decay in 136Xe with EXO-200,” Physical Review Letters, vol. 109, no. 3, Article ID 032505, 2012.View at: Google Scholar
V. Álvarez, F. I. G. M. Borges, S. Cárcel et al., “Initial results of NEXT-DEMO, a large-scale prototype of the NEXT-100 experiment,” Journal of Instrumentation, vol. 8, Article ID P04002, 2013.View at: Google Scholar
A. Alessandrello, C. Brofferio, D. V. Camin et al., “Measuring thermistor resistance with very low D.C. power dissipation,” Cryogenics, vol. 37, no. 1, pp. 27–31, 1997.View at: Google Scholar
G. Pessina, “Low-noise, low drift, high precision linear bipolar (±10 V) voltage supply/reference for cryogenic front-end apparatus,” Review of Scientific Instruments, vol. 70, no. 8, pp. 3473–3478, 1999.View at: Google Scholar
A. Alessandrello, V. Bashkirov, C. Brofferio et al., “Development of a thermal scintillating detector for double beta decay of 48Ca,” Nuclear Physics B, vol. 28, no. 1, pp. 233–235, 1992.View at: Google Scholar
C. Bobin, I. Berkes, J. P. Hadjout, N. Coron, J. Leblanc, and P. De Marcillac, “Alpha/gamma discrimination with a CaF2(Eu) target bolometer optically coupled to a composite infrared bolometer,” Nuclear Instruments and Methods in Physics Research A, vol. 386, no. 2-3, pp. 453–457, 1997.View at: Publisher Site | Google Scholar
J. W. Beeman, A. Gentils, A. Giuliani et al., “Effect of SiO2 coating in bolometric Ge light detectors for rare event searches,” Nuclear Instruments and Methods in Physics Research A, vol. 709, pp. 22–28, 2013.View at: Google Scholar
C. Arnaboldi, A. Giachero, C. Gotti, and G. Pessina, “A very high performance stabilization system for large mass bolometer experiments,” Nuclear Instruments and Methods in Physics Research A, vol. 652, pp. 306–309, 2011.View at: Google Scholar
V. Radeka and N. Karlovac, “Least-square-error amplitude measurement of pulse signals in presence of noise,” Nuclear Instruments and Methods, vol. 52, no. 1, pp. 86–92, 1967.View at: Google Scholar
J. W. Beeman, F. Bellini, S. Capelli et al., “ZnMoO4: a promising bolometer for neutrinoless double beta decay searches,” Astroparticle Physics, vol. 35, no. 12, pp. 813–820, 2012.View at: Google Scholar
L. Gironi, “Pulse shape analysis with scintillating bolometers,” The Journal of Low Temperature Physics, vol. 167, pp. 504–509, 2012.View at: Google Scholar
H. Okada, T. Kawanaka, and S. Ohmoto, “Study on the ZnSe phase diagram by differential thermal analysis,” Journal of Crystal Growth, vol. 165, no. 1-2, pp. 31–36, 1996.View at: Google Scholar
M. P. Kulakov, V. D. Kulakovskii, I. B. Savchenko, and A. V. Fadeev, “Phase transition in zinc selenide crystals,” Soviet Physics—Solid State, vol. 18, pp. 526–527, 1976.View at: Google Scholar
P. Benetti, I. Dafinei, F. Ferroni et al., “Enriched isotopes for double beta decays experiments: Se82 and Lucifer,” in Proceedings of the 12th International Workshop on Separation Phenomena in Liquids and Gases (SPLG '12), pp. 434–837, Paris, France, June 2012.View at: Google Scholar
P. Benetti, I. Dafinei, F. Ferroni et al., “Synthesis of 82Se enriched ZnSe for double beta decay experiments,” in Proceedings of the 12th International Workshop on Separation Phenomena in Liquids and Gases (SPLG '12), pp. 438–843, Paris, France, June 2012.View at: Google Scholar
A. Borio di Tigliole, A. Cammi, M. Clemenza, V. Memoli, L. Pattavina, and E. Previtali, “Benchmark evaluation of reactor critical parameters and neutron fluxes distributions at zero power for the TRIGA Mark II reactor of the University of Pavia using the Monte Carlo code MCNP,” Progress in Nuclear Energy, vol. 52, no. 5, pp. 494–502, 2010.View at: Publisher Site | Google Scholar
D. M. Chernyak, F. A. Danevich, A. Giuliani, E. Olivieri, M. Tenconi, and V. I. Tretyak, “Random coincidence of 2ν2β decay events as a background source in bolometric 02β decay experiments,” The European Physical Journal C, vol. 72, p. 1989, 2012.View at: Google Scholar
J. Barea, J. Kotilia, and F. Iachello, “Nuclear matrix elements for double-B decay,” Physical Review C, vol. 87, Article ID 014315, 2013.View at: Google Scholar
P. K. Rath, R. Chandra, K. Chaturvedi, P. K. Raina, and J. G. Hirsch, “Uncertainties in nuclear transition matrix elements for neutrinoless ββ decay within the projected-Hartree-Fock-Bogoliubov model,” Physical Review C, vol. 82, no. 6, Article ID 064310, 2010.View at: Publisher Site | Google Scholar
T. R. Rodriguez and G. Martinez-Pinedo, “Energy density functional study of nuclear matrix elements for neutrinoless ββ decay,” Physical Review Letters, vol. 105, Article ID 252503, 2010.View at: Google Scholar