Uncertainties in nuclear transition matrix elements for neutrinoless ββ decay within the projected-Hartree-Fock-Bogoliubov model

The nuclear transition matrix elements M-(0 nu) for the neutrinoless beta beta decay of Zr-94,Zr-96, Mo-98,Mo-100, Ru-104, Pd-110, Te-128,Te-130, and Nd-150 isotopes in the case of 0(+) -> 0(+) transition are calculated using the projected-Hartree-Fock-Bogoliubov wave functions, which are eigenvectors of four different parametrizations of a Hamiltonian with pairing plus multipolar effective two-body interaction. Employing two (three) different parametrizations of Jastrow-type short-range correlations, a set of eight (twelve) different nuclear transition matrix elements M-(0 nu) is built for each decay, whose averages in conjunction with their standard deviations provide an estimate of the model's uncertainties.