Calculation Of Nuclear Matrix Elements Research Articles

First Study of the PIKACHU Project: Development and Evaluation of High-Purity Gd3Ga3Al2O12:Ce Crystals for 160Gd Double Beta Decay Search

Abstract Uncovering neutrinoless double beta decay (0ν2β) is crucial for confirming neutrinos’ Majorana characteristics. The decay rate of 0νββ is theoretically uncertain, influenced by nuclear matrix elements that vary across nuclides. To reduce this uncertainty, precise measurement of the half-life of neutrino-emitting double beta decay (2ν2β) in different nuclides is essential. We have launched the PIKACHU (Pure Inorganic scintillator experiment in KAmioka for CHallenging Underground sciences) project to fabricate high-purity Ce-doped Gd3Ga3Al2O12 (GAGG) single crystals and use them to study the double beta decay of 160Gd. Predictions from two theoretical models on nuclear matrix element calculations for 2ν2β in 160Gd show a significant discrepancy in estimated half-lives, differing by approximately an order of magnitude. If the lower half-life estimation holds true, detecting 2ν2β in 160Gd could be achievable with a sensitivity enhancement slightly more than an order of magnitude compared to prior investigations using Ce-doped Gd2SiO5 (GSO) crystal. We have successfully developed GAGG crystals with purity levels surpassing previous standards through refined purification and selection of raw materials. Our experiments with these crystals indicate the feasibility of reaching sensitivities exceeding those of earlier studies. This paper discusses the ongoing development and scintillator performance evaluation of High-purity GAGG crystals, along with the anticipated future prospects of the PIKACHU experiment.

Impact of nuclear matrix element calculations for current and future neutrinoless double beta decay searches

Nuclear matrix elements (NME) are a crucial input for the interpretation of neutrinoless double beta decay data. We consider a representative set of recent NME calculations from different methods and investigate the impact on the present bound on the effective Majorana mass mββ by performing a combined analysis of the available data as well as on the sensitivity reach of future projects. A crucial role is played by the recently discovered short-range contribution to the NME, induced by light Majorana neutrino masses. Depending on the NME model and the relative sign of the long- and short-range contributions, the current 3σ bound can change between mββ< 40 meV and 600 meV. The sign-uncertainty may either boost the sensitivity of next-generation experiments beyond the region for mββ predicted for inverted mass ordering or prevent even advanced setups to reach this region. Furthermore, we study the possibility to distinguish between different NME calculations by assuming a positive signal and by combining measurements from different isotopes. Such a discrimination will be impossible if the relative sign of the long- and short-range contribution remains unknown, but can become feasible if mββ ≳ 40 meV and if the relative sign is known to be positive. Sensitivities will be dominated by the advanced 76Ge and 136Xe setups assumed here, but NME model-discrimination improves if data from a third isotope is added, e.g., from 130Te or 100Mo.

Probing the mechanism of neutrinoless double-beta decay in multiple isotopes

A large experimental program is being mounted to search for neutrinoless double-beta decay over the next decade. Multiple experiments using different target isotopes are being prepared to explore the whole parameter space allowed for inverted-ordered light neutrinos, and have the potential to make discoveries in several other scenarios, including normal-ordered light neutrinos and other exotic mechanisms. We investigate to what extent long-range and exotic short-range contributions may be distinguished by combining measurements of the decay half-life across isotopes in the framework of a global Bayesian analysis. We demonstrate how measurements in two isotopes will constrain the parameter space up to a two-fold degeneracy, and how a further measurement in a third isotope removes such a degeneracy. We also discuss the impact of uncertainties and correlations in nuclear matrix element calculations. Our work motivates an experimental program measuring neutrinoless double-beta decay in more than one isotope, as this would break parameter degeneracies and advance our understanding of particle physics beyond the Standard Model.

Search for the Majorana Nature of Neutrinos in the Inverted Mass Ordering Region with KamLAND-Zen.

The KamLAND-Zen experiment has provided stringent constraints on the neutrinoless double-beta (0νββ) decay half-life in ^{136}Xe using a xenon-loaded liquid scintillator. We report an improved search using an upgraded detector with almost double the amount of xenon and an ultralow radioactivity container, corresponding to an exposure of 970kgyr of ^{136}Xe. These new data provide valuable insight into backgrounds, especially from cosmic muon spallation of xenon, and have required the use of novel background rejection techniques. We obtain a lower limit for the 0νββ decay half-life of T_{1/2}^{0ν}>2.3×10^{26} yr at 90%C.L., corresponding to upper limits on the effective Majorana neutrino mass of 36-156meV using commonly adopted nuclear matrix element calculations.

Neutrinoless double- β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

We devise a framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of experimentally relevant nuclei, including $^{76}$Ge, $^{130}$Te, and $^{136}$Xe. In light nuclei, we validate our nuclear matrix element calculations by comparing against accurate variational Monte Carlo results. Due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, in heavier systems, we obtain long-range nuclear matrix elements that are about 30% smaller than previous shell-model results. We also evaluate the recently recognized short-range nuclear matrix element estimating its coupling by the charge-independence-breaking term of the Argonne $v_{18}$ potential used in the Monte Carlo calculations. Our results indicate an enhancement of the total nuclear matrix element by around 30%.

Discovering neutrinoless double-beta decay in the era of precision neutrino cosmology

We evaluate the discovery probability of a combined analysis of proposed neutrinoless double-beta decay experiments in a scenario with normal ordered neutrino masses. The discovery probability strongly depends on the value of the lightest neutrino mass, ranging from zero in case of vanishing masses and up to 80--90% for values just below the current constraints. We study the discovery probability in different scenarios, focusing on the exciting prospect in which cosmological surveys will measure the sum of neutrino masses. Uncertainties in nuclear matrix element calculations partially compensate each other when data from different isotopes are available. Although a discovery is not granted, the theoretical motivations for these searches and the presence of scenarios with high-discovery probability strongly motivates the proposed international, multi-isotope experimental program.

Majorana neutrino mass constraints in the landscape of nuclear matrix elements

We discuss up-to-date constraints on the Majorana neutrino mass $m_{\beta\beta}$ from neutrinoless double beta decay ($0\nu\beta\beta$) searches in experiments using different isotopes: KamLAND-Zen and EXO ($^{136}$Xe), GERDA and MAJORANA ($^{76}$Ge) and CUORE ($^{130}$Te). Best fits and upper bounds on $m_{\beta\beta}$ are explored in the general landscape of nuclear matrix elements (NME), as well as for specific NME values obtained in representative nuclear models. By approximating the likelihood of $0\nu\beta\beta$ signals through quadratic forms, the analysis of separate and combined isotope data becomes exceedingly simple, and allows to clarify various aspects of multi-isotope data combinations. In particular, we analyze the relative impact of different data in setting upper bounds on $m_{\beta\beta}$, as well as the conditions leading to nonzero $m_{\beta\beta}$ at best fit, for variable values of the NMEs. Detailed results on $m_{\beta\beta}$ from various combinations of data are reported in graphical and numerical form. Implications for future $0\nu\beta\beta$ data analyses and NME calculations are briefly discussed.

High-precision Q-value measurement and nuclear matrix element calculations for the double-beta decay of ^{98}Mo

The ^{98}Mo double-beta decay Q-value has been measured, and the corresponding nuclear matrix elements of neutrinoless double-beta (0nu beta beta ) decay and the standard two-neutrino double-beta (2nu beta beta ) decay have been provided by nuclear theory. The double-beta decay Q-value has been determined as Q_{beta beta }=113.668(68) keV using the JYFLTRAP Penning trap mass spectrometer. It is in agreement with the literature value, Q_{beta beta }=109(6) keV, but almost 90 times more precise. Based on the measured Q-value, precise phase-space factors for 2nu beta beta decay and 0nu beta beta decay, needed in the half-life predictions, have been calculated. Furthermore, the involved nuclear matrix elements have been computed in the proton–neutron quasiparticle random-phase approximation (pnQRPA) and the microscopic interacting boson model (IBM-2) frameworks. Finally, predictions for the 2nu beta beta decay are given, suggesting a much longer half-life than for the currently observed cases.

Prospects for searching new Double Beta Decay physics with KamLAND-Zen 800

Two-neutrino double beta decay (2υββ) is a rare second-order weak radioactive decay process. It has been observed in neutrinoless double beta decay (0υββ) search experiments. Precise observation of 2υββ is essential to reduce the theoretical uncertainty in the calculation of nuclear matrix elements required to obtain the effective Majorana mass from a lifetime of 0υββ. Also, 2υββ itself is interesting because new physics could be hiding in the energy spectrum, such as for example, Majoron emission mode 0υββ and 2υββ with neutrino self-interaction. KamLAND-Zen 800 is an experiment to search for 0υββ of 136Xe with a large ultra-pure liquid scintillator detector, KamLAND. KamLAND-Zen 800 has been observing 745 kg of xenon gas at 91% enriched in 136Xe since 2019, providing a high statistics 2υββ decay sample. We describe the potential for new physics searches in 2υββ with KamLAND-Zen 800.

Interacting Shell Model Calculations for Neutrinoless Double Beta Decay of 82Se With Left-Right Weak Boson Exchange

In the present work, the λ mechanism (left-right weak boson exchange) and the light neutrino-exchange mechanism of neutrinoless double beta decay is studied. In particular, much attention is paid to the calculation of nuclear matrix elements for one of the neutrinoless double beta decaying isotopes 82Se. The interacting shell model framework is used to calculate the nuclear matrix element. The widely used closure approximation is adopted. The higher-order effect of the pseudoscalar term of nucleon current is also included in some of the nuclear matrix elements that result in larger Gamow-Teller matrix elements for the λ mechanism. Bounds on Majorana neutrino mass and lepton number violating parameters are also derived using the calculated nuclear matrix elements.

Ab initio benchmarks of neutrinoless double- β decay in light nuclei with a chiral Hamiltonian

We report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-beta ($0\nu\beta\beta$) decays in light nuclei with mass number ranging from $A=6$ to $A=22$. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two-plus-three-nucleon interaction from chiral effective field theory is employed, and both isospin-conserving ($\Delta T=0$) and isospin-changing ($\Delta T=2$) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of $\Delta T=0$ transitions are in good agreement among all calculations, at the level of 10%. For $\Delta T=2$, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps towards improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross-checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the $0\nu\beta\beta$ decays in heavier candidate nuclei, where quasi-exact benchmarks are not available.

Analysis of light neutrino exchange and short-range mechanisms in0νββdecay

Neutrinoless double beta decay ($0\nu\beta\beta$) is a crucial test for lepton number violation. Observation of this process would have fundamental implications for neutrino physics, theories beyond the Standard Model and cosmology. Focussing on so called short-range operators of $0\nu\beta\beta$ and their potential interplay with the standard light Majorana neutrino exchange, we present the first complete calculation of the relevant nuclear matrix elements, performed within the interacting boson model (IBM-2). Furthermore, we calculate the relevant phase space factors using exact Dirac electron wavefunctions, taking into account the finite nuclear size and screening by the electron cloud. The obtained numerical results are presented together with up-to-date limits on the standard mass mechanism and effective $0\nu\beta\beta$ short-range operators in the IBM-2 framework. Finally, we interpret the limits in the particle physics scenarios incorporating heavy sterile neutrinos, Left-Right symmetry and R-parity violating supersymmetry.

Comparative analysis of muon-capture and 0νββ -decay matrix elements

Average matrix elements of ordinary muon capture (OMC) to the intermediate nuclei of neutrinoless double beta ($0\ensuremath{\nu}\ensuremath{\beta}\ensuremath{\beta}$) decays of current experimental interest are computed and compared with the corresponding energy and multipole decompositions of $0\ensuremath{\nu}\ensuremath{\beta}\ensuremath{\beta}$-decay nuclear matrix elements (NMEs). The present OMC computations are performed using the Morita-Fujii formalism by extending the original formalism beyond the leading order. The $0\ensuremath{\nu}\ensuremath{\beta}\ensuremath{\beta}$ NMEs include the appropriate short-range correlations, nuclear form factors, and higher-order nucleonic weak currents. The nuclear wave functions are obtained in extended no-core single-particle model spaces using the spherical version of the proton-neutron quasiparticle random-phase approximation with two-nucleon interactions based on the Bonn one-boson-exchange $G$ matrix. Both the OMC and $0\ensuremath{\nu}\ensuremath{\beta}\ensuremath{\beta}$ processes involve 100-MeV-range momentum exchanges and hence similarities could be expected for both processes in the feeding of the $0\ensuremath{\nu}\ensuremath{\beta}\ensuremath{\beta}$ intermediate states. These similarities may help improve the accuracy of the $0\ensuremath{\nu}\ensuremath{\beta}\ensuremath{\beta}$ NME calculations by using the data from the currently planned OMC experiments.

Lattice QCD and neutrino-nucleus scattering

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD in neutrino-oscillation physics, which inevitably entails nucleon and nuclear structure. In addition to discussing pertinent lattice-QCD calculations of nucleon and nuclear matrix elements, the interplay with models of nuclei is discussed. This program of lattice- QCD calculations is relevant to current and upcoming neutrino experiments, becoming increasingly important on the timescale of LBNF/DUNE and HyperK.

Study of the Effect of Newly Calculated Phase Space Factor on β-Decay Half-Lives

We present results for β-decay half-lives based on a new recipe for calculation of phase space factors recently introduced. Our study includes fp-shell and heavier nuclei of experimental and astrophysical interests. The investigation of the kinematics of some β-decay half-lives is presented, and new phase space factor values are compared with those obtained with previous theoretical approximations. Accurate calculation of nuclear matrix elements is a prerequisite for reliable computation of β-decay half-lives and is not the subject of this paper. This paper explores if improvements in calculating the β-decay half-lives can be obtained when using a given set of nuclear matrix elements and employing the new values of the phase space factors. Although the largest uncertainty in half-lives computations come from the nuclear matrix elements, introduction of the new values of the phase space factors may improve the comparison with experiment. The new half-lives are systematically larger than previous calculations and may have interesting consequences for calculation of stellar rates.

Single Particle Levels and ββ-Decay Matrix Elements in The Interacting Boson Model

Recently a new method to calculate the occupancies of single particle levels in atomic nuclei was developed in the context of the microscopic interacting boson model, in which neutron and proton degrees of freedom are treated explicitly (IBM-2). The energies of the single particle levels constitute a very important input for the calculation of the occupancies in this method, and further they play important role in the calculation of double beta decay nuclear matrix elements. Here we discuss how the 0νββ, 0νhββ, and 2νββ-decay nuclear matrix elements (NMEs) are affected when the energies of single particle levels are changed.

First Result on the Neutrinoless Double-β Decay of ^{82}Se with CUPID-0.

We report the result of the search for neutrinoless double beta decay of ^{82}Se obtained with CUPID-0, the first large array of scintillating Zn^{82}Se cryogenic calorimeters implementing particle identification. We observe no signal in a 1.83kg yr ^{82}Se exposure, and we set the most stringent lower limit on the 0νββ ^{82}Se half-life T_{1/2}^{0ν}>2.4×10^{24} yr (90% credible interval), which corresponds to an effective Majorana neutrino mass m_{ββ}<(376-770) meV depending on the nuclear matrix element calculations. The heat-light readout provides a powerful tool for the rejection of α particles and allows us to suppress the background in the region of interest down to (3.6_{-1.4}^{+1.9})×10^{-3} counts/(keV kg yr), an unprecedented level for this technique.

Search for double β-decays of 124Xe with XENON100 & XENON1T

The rare nuclear process of two-neutrino double electron capture, where two electrons are simultaneously captured from the atomic shell, has not yet been observed for 124Xe. A detection of this decay would provide a new reference for nuclear matrix element calculations. Moreover, if a neutrinoless mode were discovered, it would prove a Majorana nature of neutrinos and would shed light on the effective neutrino mass. The XENON dark matter project, with its dual-phase xenon time projection chambers XENON100 and XENON1T, is well suited for this rare event searches with signatures in the keV-region. The search with the XENON100 detector, containing 29 g of 124Xe, is explained as well as the outlook of its successor XENON1T, which contains 2 kg of the isotope in its active volume.

New limit for the half-life of double beta decay of 94Zr to the first excited state of 94Mo

Neutrinoless Double Beta Decay is a phenomenon of fundamental interest in particle physics. The decay rates of double beta decay transitions to the excited states can provide input for Nuclear Transition Matrix Element calculations for the relevant two neutrino double beta decay process. It can be useful as supplementary information for the calculation of Nuclear Transition Matrix Element for the neutrinoless double beta decay process. In the present work, double beta decay of $^{94}$Zr to the $2^{+}_{1}$ excited state of $^{94}$Mo at 871.1 keV is studied using a low background $\sim$ 230 cm$^3$ HPGe detector. No evidence of this decay was found with a 232 g.y exposure of natural Zirconium. The lower half-life limit obtained for the double beta decay of $\rm^{94}Zr$ to the $2^{+}_{1}$ excited state of $\rm^{94}Mo$ is $T_{1/2} (0\nu + 2\nu)> 3.4 \times 10^{19}$ y at 90% C.L., an improvement by a factor of $\sim$ 4 over the existing experimental limit at 90\% C.L. The sensitivity is estimated to be $T_{1/2} (0\nu + 2\nu) > 2.0\times10^{19}$ y at 90% C.L. using the Feldman-Cousins method.

Search for two-neutrino double electron capture of Xe124 with XENON100

Two-neutrino double electron capture is a rare nuclear decay where two electrons are simultaneously captured from the atomic shell. For $^{124}$Xe this process has not yet been observed and its detection would provide a new reference for nuclear matrix element calculations. We have conducted a search for two-neutrino double electron capture from the K-shell of $^{124}$Xe using 7636 kg$\cdot$d of data from the XENON100 dark matter detector. Using a Bayesian analysis we observed no significant excess above background, leading to a lower 90 % credibility limit on the half-life $T_{1/2}>6.5\times10^{20}$ yr. We also evaluated the sensitivity of the XENON1T experiment, which is currently being commissioned, and find a sensitivity of $T_{1/2}>6.1\times10^{22}$ yr after an exposure of 2 t$\cdot$yr.