Systematic analysis of the uncertainty in the $0νββ$ decay nuclear matrix elements

The nuclear matrix elements $M^{0\nu}$ of the neutrinoless double beta decay ($0\nu\beta\beta$) of most nuclei with known $2\nu\beta\beta$-decay rates are systematically evaluated using the Quasiparticle Random Phase Approximation (QRPA) and Renormalized QRPA (RQRPA). The experimental $2\nu\beta\beta$-decay rate is used to adjust the most relevant parameter, the strength of the particle-particle interaction. New results confirm that with such procedure the $M^{0\nu}$ values become essentially independent on the size of the single-particle basis. Furthermore, the matrix elements are shown to be also rather stable with respect to the possible quenching of the axial vector strength parameterized by reducing the coupling constant $g_A$, as well as to the uncertainties of parameters describing the short range nucleon correlations. Theoretical arguments in favor of the adopted way of determining the interaction parameters are presented.