  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 of Ge76, Mo100, Cd116, Te128 and Xe136 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 behaviour, 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.
