Comparison between variational Monte Carlo and shell model calculations of neutrinoless double beta decay matrix elements in light nuclei

Benchmark comparisons between many-body methods are performed to assess the ingredients necessary for an accurate calculation of neutrinoless double beta decay matrix elements. Shell model and variational Monte Carlo (VMC) calculations are carried out for A=10 and 12 nuclei. Different variational wavefunctions are used to evaluate the uncertainties in the ab initio calculations, finding fairly small differences between the VMC double beta decay matrix elements. For shell model calculations, the role of model space truncation, radial wavefunction choices, and short-range correlation are investigated and all found to be important. Based on the detailed comparisons between the VMC and shell model approaches, we conclude that accurate descriptions of neutrinoless double beta decay matrix elements require a proper treatment of both long-range and short-range correlations.