The nuclear matrix elements that govern the rate of neutrinoless double beta decay must be accurately calculated if experiments are to reach their full potential. Theorists have been working on the problem for a long time but have recently stepped up their efforts as ton-scale experiments have begun to look feasible. Here we review past and recent work on the matrix elements in a wide variety of nuclear models and discuss work that will be done in the near future. Ab initio nuclear-structure theory, which is developing rapidly, holds out hope of more accurate matrix elements with quantifiable error bars.
Accurate nuclear matrix elements (NMEs) for neutrinoless double beta decays of candidate nuclei are important for the design and interpretation of future experiments. Significant progress has been made in the modeling of these NMEs from first principles. The NME for 48Ca shows a good agreement among three different ab initio calculations starting from the same nuclear interaction constructed within the chiral EFT and the same decay operator. These studies open the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as 76Ge, 130Te, and 136Xe. The ultimate goal is the computation of NMEs in many-body calculations with controllable approximations, using nuclear interactions and weak transition operators derived consistently from chiral EFT. We are expecting more progress towards this goal in the near future.
The nuclear matrix elements of neutrinoless double-$beta$ decay for nuclei $^{76}$Ge, $^{82}$Se, $^{100}$Mo, $^{130}$Te, and $^{150}$Nd are studied within the triaxial projected shell model, which incorporates simultaneously the triaxial deformation and quasiparticle configuration mixing. The low-lying spectra and the $B(E2:0^+rightarrow2^+)$ values are reproduced well. The effects of the quasiparticles configuration mixing, the triaxial deformation, and the closure approximation on the nuclear matrix elements are studied in detail. For nuclei $^{76}$Ge, $^{82}$Se, $^{100}$Mo, $^{130}$Te, and $^{150}$Nd, the nuclear matrix elements are respectively reduced by the quasiparticle configuration mixing by 6%, 7%, 2%, 3%, and 4%, and enhanced by the odd-odd intermediate states by 7%, 4%, 11%, 20%, and 14%. Varying the triaxial deformation $gamma$ from $0^circ$ to $60^circ$ for the mother and daughter nuclei, the nuclear matrix elements change by 41%, 17%, 68%, 14%, and 511% respectively for $^{76}$Ge, $^{82}$Se, $^{100}$Mo, $^{130}$Te, and $^{150}$Nd, which indicates the importance of treating the triaxial deformation consistently in calculating the nuclear matrix elements.
We calculate basis-space converged neutrinoless $beta beta$ decay nuclear matrix elements for the lightest candidates: 48Ca, 76Ge and 82Se. Starting from initial two- and three-nucleon forces, we apply the ab initio in-medium similarity renormalization group to construct valence-space Hamiltonians and consistently transformed $beta beta$-decay operators. We find that the tensor component is non-negligible in 76Ge and 82Se, and resulting nuclear matrix elements are overall 25-45% smaller than those obtained from the phenomenological shell model. While a final matrix element with uncertainties still requires substantial developments, this work nevertheless opens a path toward a true first-principles calculation of neutrinoless $beta beta$ decay in all nuclei relevant for ongoing large-scale searches.
A new generation of neutrinoless double beta decay experiments with improved sensitivity is currently under design and construction. They will probe inverted hierarchy region of the neutrino mass pattern. There is also a revived interest to the resonant neutrinoless double-electron capture, which has also a potential to probe lepton number conservation and to investigate the neutrino nature and mass scale. The primary concern are the nuclear matrix elements. Clearly, the accuracy of the determination of the effective Majorana neutrino mass from the measured 0 ubetabeta-decay half-life is mainly determined by our knowledge of the nuclear matrix elements. We review recent progress achieved in the calculation of 0 ubetabeta and 0 u ECEC nuclear matrix elements within the quasiparticle random phase approximation. A considered self-consistent approach allow to derive the pairing, residual interactions and the two-nucleon short-range correlations from the same modern realistic nucleon-nucleon potentials. The effect of nuclear deformation is taken into account. A possibility to evaluate 0 ubetabeta-decay matrix elements phenomenologically is discussed.