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

We compute the nuclear matrix elements (NMEs) corresponding to the neutrinoless double beta (0v beta beta) decays of nuclei which attract current experimental interest. We concentrate on ground-state-to-ground-state decay transitions mediated by light (1-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-exchangeGmatrix. Both the 1-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 (g(pp)) is decomposed into isoscalar (g(pp)(T=0)) and isovector (g(pp)(T= 1)) pp) parts. A detailed analysis of the 1-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.