While the discovery of non-zero neutrino masses is one of the most important accomplishments by physicists in the past century, it is still unknown how and in what form these masses arise. Lepton number-violating neutrinoless double beta decay is a natural consequence of Majorana neutrinos and many BSM theories, and several experimental efforts are involved in the search for these processes. Understanding how neutrinoless double beta decay would manifest in nuclear environments is key for understanding the underlying mechanism behind any observed signals. In this talk I will present an overview of a set of one- and two-body matrix elements relevant for experimental searches for neutrinoless double beta decay, along with preliminary lattice QCD results.