Using the renormalized quasiparticle random phase approximation (RQRPA), we calculate the light neutrino mass mediated mode of neutrinoless double beta decay (0νββ-decay) of 76 Ge, 100 M o, 128 T e and 130 T e. Our results indicate that the simple quasiboson approximation is not good enough to study the 0νββ-decay, because its solutions collapse for physical values of g pp . We find that extension of the Hilbert space and inclusion of the Pauli Principle in the QRPA with proton-neutron pairing, allows us to extend our calculations beyond the point of collapse, for physical values of the nuclear force strength.