Within the squark-neutrino mechanism of RpRp-violating SUSY, sets of 12 nuclear transition matrix elements (NTMEs) are calculated for the neutrinoless double-ββ decay (0+→0+transition)(0+→0+transition) of 94,9694,96Zr, 100100Mo, 110110Pd, 128,130128,130Te and 150150Nd isotopes. Specifically, four sets of HFB wave functions generated with four different parametrizations of the pairing plus multipolar two-body interactions, dipole form factor and three different parametrizations of the Jastrow short-range correlations are employed in the calculation of NTMEs with two possible prescriptions for the hadronization, namely the two-nucleon mode and the pionic mode. Without (with) Miller–Spencer parametrization of short-range correlation, uncertainties in average NTMEs ¯M(̃qν)2N¯¯¯¯¯¯M(˜qν)2N (QBM), ¯M(̃qν)2N (NRQM), ¯M(̃qν)2N (FF3) and ¯M(̃qν)π turn out be 11–18% (29–37%), 11–16% (27–31%), 5–12% (13–17%) and 3–13% (9–15%), respectively.