Nuclear matrix elements for 0νββ decays with light or heavy Majorana-neutrino exchange

We compute the nuclear matrix elements (NMEs) corresponding to the neutrinoless double beta (0νββ) decays of nuclei which attract current experimental interest. We concentrate on ground-state-to-ground-state decay transitions mediated by light (l-NMEs) or heavy (h-NMEs) Majorana neutrinos. The computations are done in realistic single-particle model spaces using the proton-neutron quasiparticle random-phase approximation (pnQRPA) with two-nucleon interactions based on the Bonn one-boson-exchange G matrix. Both the l-NMEs and the h-NMEs include the appropriate short-range correlations, nucleon form factors, and higher-order nucleonic weak currents. In addition, both types of NMEs are corrected for the isospin symmetry by the recently proposed method in which the particle-particle proton-neutron interaction parameter (gpp) is decomposed into isoscalar (gT =0 pp ) and isovector (gT =1 pp ) parts. A detailed analysis of the l-NMEs and the h-NMEs is performed to benchmark our computer code and to compare with other recent calculations which produce h-NMEs that are in tension with each other.peerReviewe