Neutrinoless double beta decay searches are currently among the major foci of experimental physics. The observation of such a decay will have important implications in our understanding of the intrinsic nature of neutrinos and shed light on the limitations of the Standard Model. The rate of this process depends on both the unknown neutrino effective mass and the nuclear matrix element associated with the given neutrinoless double-beta decay transition. The latter can only be provided by theoretical calculations, hence the need of accurate theoretical predictions of the nuclear matrix element for the success of the experimental programs. This need drives the theoretical nuclear physics community to provide the most reliable calculations of the nuclear matrix elements. Among the various computational models adopted to solve the many-body nuclear problem, the shell model is widely considered as the basic framework of the microscopic description of the nucleus. Here, we review the most recent and advanced shell-model calculations of the nuclear matrix elements considering the light-neutrino-exchange channel for nuclei of experimental interest. We report the sensitivity of the theoretical calculations with respect to variations in the model spaces and the shell-model nuclear Hamiltonians.
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.
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 (NMEs) for two-neutrino double-beta decay ($2 ubetabeta$) are studied in the framework of the relativistic nuclear energy density functional. The properties of nuclei involved in the decay are obtained using relativistic Hartree-Bogoliubov model and relevant nuclear transitions are described using the relativistic proton-neutron quasiparticle random phase approximation (pn-ReQRPA). Three effective interactions have been employed, including density-dependent meson-exchange and point coupling interactions, supplemented with nuclear pairing correlations. The $2 ubetabeta$ matrix elements and half-lives are calculated for several nuclides experimentally known to undergo this kind of decay: $^{48}$Ca, $^{76}$Ge, $^{82}$Se, $^{96}$Zr, $^{100}$Mo, $^{116}$Cd, and $^{128}$Te. The model dependence of the NMEs and their sensitivity on the isoscalar pairing strength $V_0$ is investigated, and the optimized value of this parameter is determined. The results of the present study represent an important benchmark for the future applications of the relativistic framework in studies of neutrinoless double-beta decay.