Uncertainties in nuclear transition matrix elements for neutrinoless   $\beta \beta $ decay within the PHFB model
  The nuclear transition matrix elements $M^{(0\nu)}$ for the neutrinoless double beta decay of $^{94,96}$Zr, $^{98,100}$Mo, $^{104}$Ru, $^{110}$Pd, $^{128,130}$Te and $^{150}$Nd isotopes in the case of $0^{+}\rightarrow 0^{+}$ transition are calculated using the PHFB wave functions, which are eigenvectors of four different parameterizations of a Hamiltonian with pairing plus multipolar effective two-body interaction. \QCOM{35}{In addition, the consideration of} Employing two (three) different parameterizations of Jastrow-type short range correlations, \QCOM{19}{provides us with} 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 uncertainties. 
P. K. Rath, R. Chandra, K. Chaturvedi, P. K. Raina and J. G. Hirsch
nucl-th