Neutrinoless double beta decay from lattice QCD: The short-distance π- → π plus e- e- amplitude

This work presents a determination of potential short-distance contributions to the unphysical 71-- -> 71-+e-e- decay through lattice QCD calculations. The hadronic contributions to the transition amplitude are described by the pion matrix elements of five Standard Model effective field theory operators, which are computed on five ensembles of domain-wall fermions with Nf 1/4 2 + 1 quark flavors with a range of heavier-than-physical values of the light quark masses. The matrix elements are extrapolated to the continuum, physical light-quark mass, and infinite volume limit using a functional form derived in chiral effective field theory (,yEFT). This extrapolation also yields the relevant low-energy constants of ,yEFT, which are a necessary input for ,yEFT calculations of neutrinoless double beta decay of nuclei.