Abstract

We report the development of a CaMoO4 crystal low temperature detector for the AMoRE neutrinoless double beta decay search experiment. The prototype detector cell was composed of a 216 g CaMoO4 crystal and a metallic magnetic calorimeter. An overground measurement demonstrated FWHM resolution of 6–11 keV for full absorption gamma peaks. Pulse shape discrimination was clearly demonstrated in the phonon signals, and 7.6  of discrimination power was found for the and separation. The phonon signals showed rise-times of about 1 ms. It is expected that the relatively fast rise-time will increase the rejection efficiency of two-neutrino double beta decay pile-up events which can be one of the major background sources in searches.

1. Introduction

Recent neutrino oscillation experiments have been unveiling the properties of neutrinos [1, 2]. Their experimental evidences strongly suggest that neutrinos are massive and encounter flavor mixing of mass eigenvalues. The mixing angles and the differences between the square masses have been estimated. However, those observations do not provide a direct measurement of the absolute mass and do not answer the question of whether neutrino is its own antiparticle (Majorana-type) or not (Dirac-type).

Search for neutrinoless double beta decay () is a key experiment to reveal unanswered nature of neutrinos [37]. The double beta decay () that accompanies the simultaneous emission of two electrons and two antineutrinos is a rare process that is an allowed transition in the standard model. Another type of double beta decay that does not emit any neutrinos, , can occur if neutrino is massive Majorana particles (i.e., it is its own antiparticle). In the process, the full available energy of the decay is carried by the two electrons and the recoiled daughter which has a very small amount of energy compared with that of the electrons. Therefore, while the electron sum energy spectrum in the is continuous up to the available energy release (), in the , the spectrum should have a sharp peak at [2].

The observation of would clearly demonstrate that neutrino is not Dirac-type but rather is Majorana-type particle. In that case, physical processes that do not conserve lepton number would be allowed. Moreover, the absolute mass scale, so-called “effective Majorana neutrino mass”, can be estimated using the half life ,where is the electron mass, is the kinematic phase-space factor calculable with reasonable precision, and is the model dependent nuclear matrix element. Here, the Majorana mass is defined aswhere the ’s are the mass eigenstates of the neutrino and ’s are the elements of the mixing matrix between the flavor states and mass eigenstates.

Experimentally, the measurement limit of half life is often used as the sensitivity to probe the rare event [7]. In a measurement with nonnegligible backgrounds, the sensitivity becomeswhere is the concentration of isotope in the detector, is the detection efficiency, is the detector mass, is the measurement time, is the background rate per unit mass and energy, and is the energy resolution of the detector, in other words, the region of interest (ROI) of the energy window at the value. However, in a case of a zero-background experiment that observes no event in ROI during the measurement time, the sensitivity becomes proportional to the detector mass and the measurement time,

To increase the detection sensitivity, it is essential to have a detector with high concentration of the isotope of interest, detection efficiency, energy resolution, and efficient background rejection capability as well as to minimize backgrounds from internal and external sources in the region of interest. High energy resolution and detection efficiency experiment can be realized with crystal detectors containing the isotope of interest. The detector performance of recently developed low temperature detectors (LTDs) that operate at sub-Kelvin temperatures can perfectly meet the requirements by utilizing state-of-the-art detector technologies with extreme energy sensitivity, such as neutron transmutation doped (NTD) Ge thermistors, superconducting transition edge sensors (TESs), or metallic magnetic calorimeters (MMCs) [8].

The AMoRE (Advanced Mo-based Rare process Experiment) project is an experiment to search for of 100Mo [9]. AMoRE uses CaMoO4 crystals as the absorber and MMCs as the sensor [10, 11]. CaMoO4 is a scintillating crystal that has the highest light output at room and low temperatures among Mo-containing crystals (molybdates) [12, 13].

Choosing of 100Mo as the source of is advantageous. The nucleus has a high value of 3034.40(17) keV [14] that is above the intensive 2615 keV gamma quanta from 208Tl decay (232Th family). The natural abundance of 100Mo is 9.8% [15], which is comparatively high. Furthermore, enriched 100Mo can be produced by centrifugation method in amount of tens of kilograms per year with a reasonable price. Also, the theoretically estimated half life of 100Mo is relatively shorter than that of other candidates [16, 17]. However, of 48Ca with keV (despite rather low concentration of the isotope ) can be an irremovable background source in the ROI of 100Mo. The AMoRE collaboration has successfully grown 40Ca100MoO4 crystals using 100Mo enriched and 48Ca depleted materials. Three of the doubly-enriched crystals with masses in the range of 0.2–0.4 kg were tested in a low background 4 veto system to determine their internal backgrounds [18].

The present experimental work aims to test the low temperature detection concept with a CaMoO4 crystal and an MMC that is suitable for a high resolution experiment to search for of 100Mo. A 216 g natural CaMoO4 crystal with an MMC phonon sensor was employed in this experiment, which was performed in an overground measurement facility. The energy resolution and linearity of the detector setup, particle, and randomly coinciding events discrimination by pulse shape analysis for background rejection in a experiment are discussed in this report.

2. Experimental Details

The detector setup was structured in a cylindrical shape with copper support details as shown in Figure 1. A CaMoO4 crystal, 4 cm in diameter and 4 cm height, was mounted inside the copper structure using metal springs. The mass of the crystal was 216 g. It was grown with natural Ca and Mo elements at the Bogoroditsk plant in Russia. A patterned gold film was evaporated on one side of the crystal to serve as a phonon collector. An MMC device, the primary sensor for detecting the phonon signals absorbed in the gold film, was placed on a semicircular copper plate over the crystal. The thermal connection between the gold film and the MMC was made using annealed gold wires. Details regarding the measurement principle and the detector structure of the MMC device were presented in previous reports [19, 20].

When a particle hits a dielectric material, most of the energy deposited into the absorber is converted into the form of phonons. High energy phonons with frequencies that are close to the Debye frequency are generated initially. However, they quickly decay to lower frequency phonons via anharmonic processes. When their energy becomes 20–50 K, they can travel ballistically in the crystal [21]. The major down-conversion processes of these athermal phonons are isotope scattering, inelastic scatterings by impurities and lattice dislocations, and inelastic scatterings at crystal surfaces [22]. These excess phonons eventually change the equilibrium thermal phonon distribution, thereby causing temperature increase.

In the detector setup with the CaMoO4 crystal and the gold phonon collector film, the ballistic athermal phonons can hit the crystal and gold interface, transmit into the gold film, and transfer their energy to the electrons in the film [23]. The electron temperature of the gold film increases quickly via electron-electron scatterings. This temperature change is measured by the MMC sensor that is thermally connected with the gold wires. The size of the gold film and number of gold wires were chosen based on a thermal model study that considered the efficient athermal heat flow process [11]. Consequently the gold film had a diameter of 2 cm, a thickness of 200 nm, and an additional gold pattern of 200 nm thickness on top of the gold film to increase the lateral thermal conductivity of the gold film.

The detector assembly was installed in a dilution refrigerator in an overground laboratory at KRISS (Korea Research Institute of Standards and Science). The refrigerator was surrounded by a 10 cm thick lead shield (except the top surface) to reduce environmental gamma-ray background. The detector with an MMC operates well in the temperature range of 10–50 mK. The signal size increases at lower temperatures since the MMC sensitivity enhanced and the heat capacities decreased. However, the signals have slower rise and decay times at lower temperatures as thermal conductances become poorer. Larger signal size improves the energy threshold and baseline energy resolution of the detector. The energy resolution of the detector measured for particle absorption events can be worse than the baseline resolution because of any uncorrelated mechanism that affects the signal size and shape. Examples of such mechanisms include temperature fluctuations due to instrumental instability or frequent event rates, position dependence of signal shapes, or scintillation processes that are associates with phonon generations in an inhomogeneous way. Therefore, larger signal size does not guarantee a better energy resolution at certain temperatures. At the present experimental condition including the background rate from cosmic muons and external gamma-rays, 40 mK was selected as the main measurement temperature. At this temperature, about 1 ms rise-time was obtained for the 2.6 MeV gamma line without degrading the energy resolution. A typical signal of 2.6 MeV gamma-ray full absorption events is shown in Figure 2. The rise-time of the DC coupled signal is 1.1 ms, which is somewhat slower than that of earlier measurements for which shorter gold wires were used [11].

3. Pulse Shape Analysis

A two dimensional scatter plot of the pulse heights and mean-times of signals obtained in a 95 h background measurement is shown in Figure 3. The pulse height is the difference between the maximum value and baseline level of a signal. The maximum value is found using a quadratic polynomial fit to the region of the signal near the pulse maximum. The baseline level is the average voltage value in the time region before the signal rises. The mean-time parameter is defined aswhere is the measured voltage value at time subtracting the baseline level, is the time when it reaches 10% of the pulse height, and and indicate the time length of the signal toward left and right directions from , respectively, to calculate the mean-time. is set to reach the time at the baseline level, while is a free parameter that was selected to achieve the most efficient particle discrimination. Here, the value was set to not include the negative part of an AC coupled signal (see Figure 2), which was used for the energy spectrum because it was recorded with finer digitizer resolution (i.e., bigger gain is used) than the DC coupled signal.

The pulse-shape discrimination (PSD) between and particles can be realized with the mean-time as a pulse shape parameter due to the difference in the rise and decay times of the two types of events. A similar tendency of PSD was reported for other low temperature scintillating detectors [24]. The separation into two groups of the events in the energy region between 4 and 5 MeV of alpha-equivalent energy can be readily observed from the distribution of mean-time values shown in Figure 4. The energy of induced events was determined as described in the following section. Although the two peaks have noticeable right-hand tails toward higher mean-time values, normal Gaussian functions were used to fit the distributions. We interpret the right-hand tails as a result of signal pile-up. A parameter of discrimination power (DP) is defined as where are the mean values and are the standard deviations of the Gaussian distributions for and events. DP was found to be 7.6.

The averaged pulse shapes for the two groups of events are compared in Figure 5. Templates of events were obtained by averaging out the 232Th -decay events pulse profiles with energy release in the crystal 4082 keV (due to the contamination of the crystal by thorium), whereas the events caused by cosmic muons with the same pulse height were selected for the template of induced events. The normalized pulses of alpha and beta templates are aligned at the time of their maximum values, as shown in Figure 5. Both the rise and decay times of the and signals are clearly different. The signals induced by events have faster rise and faster decay than those of the events.

According to scintillation measurements of a CaMoO4 crystal at 7–300 K [25, 26], the scintillation decay time of the CaMoO4 crystal reaches hundreds of at 7 K. The crystal also shows different light outputs for and events [27]. A slowly decaying scintillation mechanism would cause slow generation of phonon in the CaMoO4 crystal. and particle events may have different fractions of slow component for phonon generation. This difference in slow phonon generation at mK temperatures may induce different pulse shapes for and particle events.

4. Energy Spectrum

Because athermal phonon absorption in the phonon collector significantly contributes to the signal size [11], the signals have some degrees of position dependence for their pulse height and shape. In this detector, the signals with faster rise-times show bigger pulse heights for the same energy events. This effect can be observed in Figure 3. The distribution of mean-time values and pulse heights for alpha and gamma-ray full absorption lines has anticorrelated slopes. This negative slope more dominantly appears in the alpha signals, which is likely because the full-energy peaks originate from multiple Compton scatterings in the crystal that smear the position dependence on the event location. The pulse height is not an optimal parameter to obtain a high resolution. In the optimal filtering method [28] that is often adopted in high resolution microcalorimeters, the signals are assumed to have one shape but with different amplitudes. Thus, the optimal filtering method is not applicable to provide a high resolution spectrum for the present signals.

In the present analysis, a new parameter LA (left area) was used as an amplitude parameter to reduce the position dependence effect of the large crystal detector. It is defined aswhere the variables are as defined in (5). LA is a partial integration for the leading part of a signal. Because the integration range of the LA parameter is set by a shape-dependent parameter, the mean-time, this parameter is less influenced by the pulse shape. For instance, the correlation coefficient between the pulse height and mean-time was −0.62 for 4082 keV signals, but it was 0.08 between the LA and mean-time.

In a calibration run, a thoriated tungsten rod was used as an external gamma-ray source. The source was placed in the gap between the cryostat and the lead shield.

When a linear energy calibration was applied to gamma-ray peaks, deviations from the linear calibration of less than 0.4% were found for low-energy peaks. A quadratic function with no constant term was used for the calibration of electron-equivalent energy for 511, 583, 911, and 2615 keV peaks for the spectrum shown in Figure 6. The corresponding energy resolutions of the peaks are listed in Table 1.

Figure 7 shows the linearities of the electron and alpha signals. The LA values of the electron and alpha peaks divided by the linear calibration of the gamma-ray peaks are plotted in the upper figure. The LA/energy ratios for the alpha peaks are about 6% larger than those for the gamma-ray peaks. The quadratic fit functions are shown as dotted and dashed lines for gamma and alpha peaks, respectively. The residuals in the lower figure indicate deviations from the quadratic functions for the two groups. With the quadratic calibration, a very small deviation is expected near 3 MeV for electron measurements. It implies that this method can supply an accurate energy calibration at the value of the of 100Mo.

Alpha events can be separated from events using the mean-time parameter. The energy spectrum of alpha events is shown in Figure 8. An energy calibration for this spectrum was performed with a quadratic function for the alpha peaks as discussed above. These backgrounds were bulk events of alpha decays in the crystal, mainly from radionuclides of U and Th decay chains, most of them are identified as shown in Figure 8. Because this crystal was developed to investigate the scintillation properties of CaMoO4, its internal background was not exceptionally low. The AMoRE collaboration has developed 40Ca100MoO4 growing technology with low radioimpurities. Internal alpha activities of about 80 Bq/kg of 226Ra and 70 Bq/kg of 228Th were found in a 196 g 40Ca100MoO4 [18].

Not only alpha signals can be rejected efficiently, but also decays of 212Bi, 214Bi, and 208Tl can be tagged and eliminated from the data by using information about associated alpha-emitting nuclides. Random pile-ups of events (first of all from the ) could be a substantial source of background of LTD to search for due to the poor time resolution [29]. Relatively fast response of the MMC among the LTDs is a certain advantage to discriminate the background.

5. Conclusion

LTDs made of crystal scintillators containing isotopes of interest have distinct advantages in searches. Such LTDs make it possible to provide a high detection efficiency to the . Taking the advantage of high resolution sensor technologies, these dielectric detectors in the heat (phonon) measurement can have comparable energy resolution to those of HPGe detectors. The comparison in heat/light measurement channels makes unambiguous separation of alpha events from electron events. As discussed in this report, pulse shape discrimination is also possible using only the phonon measurement. This PSD capability of the phonon sensor increases discrimination power for alpha background signals or can simplify the detector cell design by using only one phonon sensor without a photon sensor that is commonly used for particle discrimination [30] and will reduce the number of measurement channels.

In comparison with signals from PMTs or conventional semiconductor detectors, phonon signals from a crystal detector in the LTD concept are typically slow. Even though the energy resolution of these detectors does not suffer from the slow signal in the low activity environment of an underground experiment, random coincidence of events leads to an unavoidable background because of the slow rise-time [29]. Two consecutive electron events that occur in a time interval that is much shorter than the signal rise-time can be regarded as a single event. Such randomly coincident events, particularly of , are an unavoidable source of backgrounds in the taking into account that the event rate of 100Mo is expected to be about 10 mBq in 1 kg 40Ca100MoO4.

In the present experiment, the phonon signals had rise-times of 1.1 ms, which is much faster than the rise-times of LTDs with NTD Ge thermistors. The fast rise-time provides efficient rejection possibility for randomly coincident events [31]. Moreover, a photon detector composed of a 2 inch Ge wafer and an MMC sensor showed a temperature independent rise-time of about 0.2 ms with reasonable energy resolution [32]. Simultaneous measurements with the photon detector will increase the discrimination power not just for alpha events but also for randomly coincident events.

For the AMoRE project, 40Ca100MoO4 crystals will be used as the detector material together with MMCs in phonon and photon measurement setups. We aim to reach a zero background with improved energy resolution of a few keV. The first stage experiment is expected to be constructed with a 10 kg prototype detector by 2016. We plan to perform a large scale experiment with 200 kg 40Ca100MoO4 crystals in the next 5-6 years. The sensitivity of the experiment to the effective Majorana neutrino mass is estimated to be in the range of 20–50 meV, which corresponds to the inverted scheme of the neutrino mass.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was funded by Grant no. IBS-R016-G1 and partly supported by the National Research Foundation of Korea Grant funded by the Korean Government (NRF-2011-220-C00006 and NRF-2013K2A5A3000039). F. A. Danevich was supported in part by the Space Research Program of the National Academy of Sciences of Ukraine.