We approach the calculation of the nuclear matrix element of the neutrinoless double-β decay process, considering the light-neutrino-exchange channel, by way of the realistic shell model. To this end, we start from a realistic nucleon-nucleon potential and then derive the effective shellmodel Hamiltonian and 0νββ-decay operator within the many-body perturbation theory. We focus on investigating the perturbative properties of the effective shell-model operator of such a decay process, aiming to establish the degree of reliability of our predictions. The contributions of the socalled short-range correlations and of the correction of Pauli-principle violations to the effective shellmodel operator, the latter introduced in many-valence nucleon systems, are also taken into account. The subjects of our study are a few candidates to the 0νββ-decay detection, in a mass interval ranging from A = 48 up to A = 136, whose spin-and spin-isospin-dependent decay properties we have studied in previous works. Our results evidence that the effect of the renormalization of the 0νββ-beta decay operator on the values of the nuclear matrix elements is less relevant than what we have obtained in previous studies of the effective single-body GT transitions operating also in the two neutrinos double-beta decay PACS numbers: 21.60.Cs, 21.30.Fe, 27.60.+j, 23.40-s 1/2 −1/2 . Moreover, combining the calculated nuclear structure factor with neutrino mixing parameters [10] and limits on m ν from current experiments, one may extract an estimation of the half-life an experiment should measure in order to be sensitive to a particular value of the neutrino effective mass [1]. All the above considerations evidence that reliable calculations of M 0ν are of paramount importance, and, currently, various nuclear structure models arXiv:2001.00890v2 [nucl-th] 3 Apr 2020