Abstract
A large Hilbert space is used for the calculation of the nuclear matrix elements governing the light neutrino mass mediated mode of neutrinoless double beta decay (Ovββ-decay) of76Ge,100Mo,116Cd,128Te, and136Xe within the proton-neutron quasiparticle random phase approximation (pn-QRPA) and the renormalized QRPA with proton-neutron pairing (full-RQRPA) methods. We have found that the nuclear matrix elements obtained with the standard pn-QRPA for several nuclear transitions are extremely sensitive to the renormalization of the particle-particle component of the residual interaction of the nuclear hamiltonian. Therefore the standard pn-QRPA does not guarantee the necessary accuracy to allow us to extract a reliable limit on the effective neutrino mass. This behavior already known from the calculation of the two-neutrino double beta decay matrix elements, manifests itself in the neutrinoless double-beta decay but only if a large model space is used. The full-RQRPA, which takes into account proton-neutron pairing and considers the Pauli principle in an approximate way, offers a stable solution in the physically acceptable region of the particle-particle strength. In this way more accurate values on the effective neutrino mass have been deduced from the experimental lower limits of the half-lifes of neutrinoless double beta decay.