Influence of pairing on the nuclear matrix elements of the neutrinoless ββ decays

We study in this Letter the neutrinoless double beta decay nuclear matrix elements (NME's) in the framework of the interacting shell model. We analyze them in terms of the total angular momentum of the decaying neutron pair and as a function of the seniority truncations in the nuclear wave functions. This point of view turns out to be very adequate to gauge the accuracy of the NME's predicted by different nuclear models. In addition, it gives back the protagonist role in this process to the pairing interaction, the one which is responsible for the very existence of double beta decay emitters. We show that low seniority approximations, comparable to those implicit in the quasiparticle RPA in a spherical basis, tend to overestimate the NME's in several decays.