A lab scale experiment for keV sterile neutrino search

We developed a simple small-scale experiment to measure the beta decay spectrum of $^{3}$H. The aim of this research is to investigate the presence of sterile neutrinos in the keV region. Tritium nuclei were embedded in a 1$\times$1$\times$1 cm$^3$ LiF crystal from the $^6$Li(n,$\alpha$)$^3$H reaction. The energy of the beta electrons absorbed in the LiF crystal was measured with a magnetic microcalorimeter at 40 mK. We report a new method of sample preparation, experiments, and analysis of $^3$H beta measurements. The spectrum of a 10-hour measurement agrees well with the expected spectrum of $^3$H beta decay. The analysis results indicate that this method can be used to search for keV-scale sterile neutrinos.


Introduction
Low-temperature detectors (LTDs) have become one of the key detector technologies in astroparticle physics applications [1,2,3]. For neutrino physics in particular, the fundamental properties of neutrinos have been investigated in several outstanding LTD experiments of neutrinoless double beta decay [4,5,6,7]. Moreover, the direct detection of the neutrino mass has been intensively carried out with promising sensitivities from an accurate end-point measurement of the 163 Ho electron capture spectrum using low temperature microcalorimeters [8,9]. On the other hand, the microcalorimeters together with the source manipulation technologies established decay energy spectroscopy that has been adopted in radionuclide analysis and the study of beta decay spectra [10,11,12,13].
Sterile neutrinos are hypothetical right-handed neutrinos that undergo no fundamental interactions other than gravity in the standard model [14]. In recent years, the existence of sterile neutrinos was implied by neutrino oscillation experimental results unexpected from the three-flavor model of ordinary neutrinos [15,16]. On the other hand, if sterile neutrinos exist with their mass scale in the keV range and a very weak interaction strength with ordinary matter, they are viable dark matter candidates [17,18]. In the search for sterile neutrinos in the mass region of several hundred keV, the most stringent experimental limit was obtained by a low-temperature experiment based on superconducting tunneling junction devices [19]. However, for the mass range of 1-100 keV, only conventional detection methods have been applied with the measurement of beta spectra of 3 H and 63 Ni [20,21,22]. Their experimental bounds have not yet been reached at the level of the low-temperature experiment in this mass region.
This present experiment is a feasibility study for detecting the presence of sterile neutrinos using a sensitive LTD technology. We employed a LiF crystal with a magnetic microcalorimeter (MMC) [23] readout to measure the beta decay events of 3 H ions that are embedded inside the crystal. We report a new method of source preparation, experimental procedure, and analysis to study sterile neutrinos from the 3 H beta spectrum.

Experiment
We developed a small low-temperature thermal calorimeter with a LiF crystal with an MMC readout. A 1×1×1 cm 3 LiF crystal was employed as a hosting material for a 3 H source and a target absorber for the energetic electrons emitted from 3 H beta decay. 3 H can be embedded in LiF from the 6 Li(n,α) 3 H reaction. Because the natural abundance of 6 Li is 7.6% and 6 Li has a large cross-section for thermal neutrons, a considerable amount of 3 H can be created when a LiF crystal is exposed to a neutron flux. From the reaction, the 3 H ions are embedded in a LiF crystal and remain in the crystal until they decay to 3 He with a half-life of 12.32 years via a β -decay process [24,25]: Electron emission results in a continuous energy spectrum with an end-point at the Q value of 18.59 keV. We placed a pure LiF crystal in a neutron flux in laboratory source storage at the Korea Research Institute of Standards and Science (KRISS). The storage hosts three neutron sources (two AmBe and one 252 Cf) surrounded by polyethylene bricks. The neutron emissions are expected to be (1-2)×10 7 n/s from the three sources. The storage has a small sample loading volume approximately 17 cm away from the neutron sources, which is a sufficient distance for the neutrons to be thermalized in the polyethylene bricks. A 10-day exposure resulted in a rate of 22 3 H beta decays per second. Moreover, a Monte Carlo simulation taking neutron thermalization and neutron capture in the LiF crystal into account for 3 H production showed a similar activity we measured with the present experiment. The mean distance for thermal neutrons to travel in a LiF crystal with natural 6 Li abundance is 2.3 mm. This indicates that 3 H implantation was relatively uniform over the LiF crystal with 1×1×1 cm 3 dimensions.
The detector composed of the LiF crystal and an MMC sensor was configured as shown in Fig. 1. Similar detector setups were developed to survey various molybdate crystals for the AMoRE project [26]. The crystal positioned in the copper frame was clamped with two Teflon pieces. At the four bottom corners, thin Teflon sheets (not shown in Fig. 1) were placed between the crystal and the copper frame to minimize heat losses from the crystal to the copper sample holder. A 5 mm×5 mm×300 nm gold film was e-beam evaporated on the surface of the crystal after the neutron activation procedure. The gold film is an intermediate thermal component that serves as a phonon collector to make a strong thermal connection to an MMC sensor with gold bonding wires [27,28]. The MMC sensor used in the present experiment was fabricated in the Au:Er (1000 ppm) batch [29]. A 55 Fe source was located at the side of the crystal as an internal calibration source. The detector setup was attached to an adiabatic demagnetization refrigerator (ADR). The ADR has a 5-cm-thick lead shield around the cryostat to reduce environmental gamma background. We used another 241 Am gamma-ray source located between the cryostat and the lead shield. This source was used only for calibration purposes and was removed when taking the 3 H beta spectrum.

Results and Discussions
The beta spectra were measured three times at 40 mK. Each measurement lasted 10 hours. The duration was limited by the ADR hold time. The 10-hour data stream was stored in a continuous data acquisition system with a 1 MS/s sampling rate. We developed an offline trigger using a Butterworth filter with a pass band of [30 kHz 300 kHz] [30]. The highfrequency combination was set to the trigger filter to minimize the influence of the lowfrequency noise on the trigger signals. Fig. 2   set at 4.3 times the standard deviation for the filtered baselines with no signal pulses. We found the trigger efficiency in a pseudo experiment that simulates realistic signal pulses by adding noise signals to a noiseless template scaled to the energy. The template pulse was found from the average signal for the 6 keV events from the 55 Fe source and 3 H beta decay. The trigger logic showed 50% efficiency at 0.95 keV, 95% efficiency at 1.45 keV and nearly 100% efficiency at 2 keV. The pulse amplitude was determined as a scale factor of the time-domain least-squares method in the comparison of the pulse template and the signals filtered with another combination of pass bands [310 Hz 1.09 kHz]. The pass band frequencies used to determine the pulse amplitude were obtained to maximize the signal-to-noise ratio (the pulse height vs. the rms value of the baseline). We also applied a pulse-shape cut using the rms of the difference between the signals and the scaled template.
The energy spectrum of the electron emission from 3 H beta decay is well understood when no other decay branch is associated, as in the sterile neutrino decay mode. The spectral shape dn/dE e of the electrons can be expressed as [31] dn where F is the relativistic Fermi function, Z = 2 is the atomic number of the daughter, Θ is the unit step function, Q is the total associated energy of the decay, K e is the kinetic energy of the electron, and p i , E i , and m i are the momentum, energy, and mass of the electron and the neutrino respectively. In thermal calorimetric detection, the detector measures the total decay energy from the 3 H beta decay except that taken by the neutrino. The recoil energy of the daughter nucleus has its maximum kinetic energy of 1.72eV [31] and can be considered negligible. The measured energy was set to the kinetic energy of the electrons in the analysis. Eq.
(2) was fitted to the measured histogram of the pulse amplitudes in the energy region between 2 keV and 40 keV. We fitted the total spectrum by including the events from 3 H beta decay, 55 Fe source, unresolved pileups, and constant backgrounds in the energy region, as shown in Fig. 4. The energy calibration and the resolution were considered as fitting parameters to be determined in the spectrum fit. Fig. 3 shows the fit results of energy calibration and resolution.
The fit result of the energy calibration showed a linear behavior, although a general quadratic relation was subject to the fitting parameters. The reference points are the results of the individual analysis near the 6 keV and 60 keV regions. They lie in the 95% confidence interval of the calibration curve. Similarly, the fit result of the energy resolution, assumed to be linear, to the energy agrees well with the individual analysis. To add 60 keV reference data, the fit analysis included the measured spectrum near 60 keV acquired with the external 241 Am source. We found a noticeable difference in the low-energy spectra with and without the 241 Am source due to excess Compton scattering events.
The spectral shape of the unresolved pileups was obtained by considering the probability of pulse shape discrimination for any two events in the spectrum occurring in various time intervals. Approximately 260 µs timing resolution was found with 99% rejection probability for the two events with the mean energy of the spectrum. The total fit result includes a constant level of background events. It influences events with energies larger than 25 keV and has little effect on the 3 H beta spectrum.
The resulting spectrum from the fitting analysis matches well with the measured data in the overall energy range between 2 keV and 40 keV, as shown in Fig. 4. The χ 2 goodness of fit, χ 2 /N, is equal to 1.047 where N is the degrees of freedom, 423 in our case. Moreover, the p-value of χ 2 is 0.24 for N. This suggests that the model fits the observed data at an acceptable level of error.
The spectral shape of the 3 H decay model used in the fitting analysis was set to follow Eq. (2), ignoring the possible decay branch associated with the emission of a massive sterile neutrino. This was mainly because the number of 10-hour events was less than that of the experiment that provides the present experimental bound for sterile neutrinos in the 1-20 keV region [20,21].
Considering the presence of the beta decay branch with a sterile neutrino, the spectral shape depends on the mixing angle θ between the massive neutrino and the ordinary threeflavor neutrinos: dλ dE e = sin 2 θ dn dE e m νs + cos 2 θ dn dE e m ν where m ν s is the mass of the sterile neutrino. Here, the mass of ordinary neutrinos can be set to zero in a keV sterile neutrino search. The mass of the electron neutrino is found to be less than 1 eV in the recent experiment of a 3 H end-point measurement [32]. The resulting spectral shape of Eq. (3) is subject to characteristic distortion at the end point (i.e., Q − m ν s ) of the beta spectrum of the sterile branch. The presence of sterile neutrinos can be explored with the 10-hour dataset used in this analysis using Eq. (3). A limit bound is expected for the mixing parameters of sin 2 θ and m ν s for the null result. We plan to expand this method for the other two sets of data. With all the datasets, the analysis is expected to reach a level close to the present up-to-date limit in the 1-16 keV region.

Conclusion
The energy spectrum of the electron emission in 3 H beta decay can be accurately measured with a LiF detector with an MMC readout. Because the detector is sensitive to almost the full energy range of the decay, the detection method allows for the search of sterile neutrinos in a similar mass range. Moreover, it was a successful demonstration of 3 H implantation in a LiF crystal, having a reasonable decay rate of 3 H ions. A larger number of crystals can be prepared with the control of the activity.
As a next step of the experiment, we plan to increase the channel number and the event rate in a crystal. Furthermore, future experiments should be carried out in a dilution refrigerator. Data acquisition does not have to be interrupted as it was in the ADR from the magnet cycles but can be settled for long-term measurement. Measurements performed at a lower temperature increase the signal amplitudes. This will increase the sensitivity to reject the pileup signals as well as the energy resolution. As a rule of thumb, a measurement for 12 channel·month with 40 Bq/channel reaches the mixing amplitude sin 2 θ of 10 −4 -10 −3 in the keV mass region.