Origin of a sensitive dependence of calculated ββ-decay amplitudes on the particle-particle residual interaction

In the present work the sensitivity of calculated ββ-decay amplitudes to a realistic residual interaction is analyzed in the framework of the approach of O. A. Rumyantsev and M. H. Urin, Phys. Lett. BPYLBAJ0370-269310.1016/S0370-2693(98)01291-X 443, 51 (1998). and V. A. Rodin, M. H. Urin, and A. Faessler, Nucl. Phys. ANUPABL0375-947410.1016/j.nuclphysa.2004.10.010 747, 297 (2005). Both the Gamow-Teller (GT) and Fermi (F) matrix elements M^2ν for two-neutrino ββ decay (2νββ decay), along with the monopole transition contributions to the total matrix elements M^0ν of neutrinoless ββ decay (0νββ decay), are calculated within the quasiparticle random-phase approximation (QRPA). In the aforementioned approach decompositions of M^2ν and M^0ν can be obtained in terms of the corresponding energy-weighted sum rules S. It is shown that in most of the cases almost the whole dependence of M^2ν and M^0ν on the particle-particle (p-p) renormalization parameter g_pp is accounted for by the g_pp dependence of the corresponding sum rules S. General expressions relating S to a realistic residual particle-particle interaction are derived, which show a pronounced sensitivity of S to the singlet-channel interaction in the case of F transitions and to the triplet-channel interaction in the case of GT transitions. Thus, the sensitivity of M^2ν and M^0ν to the SU(4)-symmetry-breaking part of the p-p residual interaction is dictated by the generic structure of the ββ-decay amplitudes. Therefore, a choice of this part in a particular calculation needs a special caution. Finally, a better isospin-consistent way of renormalization of a realistic residual p-p interaction to use in QRPA calculations is suggested.