# Long\-distance nuclear matrix elements for neutrinoless double\-beta decay from lattice QCD

**INSPIRE:** [2758519](https://inspirehep.net/literature/2758519)
**arXiv:** [2402\.09362](https://arxiv.org/abs/2402.09362)
**DOI:** [10\.1103/PhysRevD\.109\.114514](https://doi.org/10.1103/PhysRevD.109.114514)

**Authors:** Davoudi, Zohreh, Detmold, William, Fu, Zhenghao, Grebe, Anthony V\., Jay, William, Murphy, David, Oare, Patrick, Shanahan, Phiala E\., Wagman, Michael L\.

**Submitted:** 14 February 2024

**Subjects:**
- hep\-lat
- nucl\-th
- Lattice
- Theory\-Nucl

**Journal reference:** Phys\.Rev\.D 109 \(11\) 114514 \(2024\)

**Citations:** 16

## Abstract

Neutrinoless double\-beta \(\<math display\="inline"\>\<mn\>0\</mn\>\<mi\>ν\</mi\>\<mi\>β\</mi\>\<mi\>β\</mi\>\</math\>\) decay is a heretofore unobserved process which, if observed, would imply that neutrinos are Majorana particles\. Interpretations of the stringent experimental constraints on \<math display\="inline"\>\<mn\>0\</mn\>\<mi\>ν\</mi\>\<mi\>β\</mi\>\<mi\>β\</mi\>\</math\>\-decay half\-lives require calculations of nuclear matrix elements\. This work presents the first lattice quantum chromodynamics \(LQCD\) calculation of the matrix element for \<math display\="inline"\>\<mn\>0\</mn\>\<mi\>ν\</mi\>\<mi\>β\</mi\>\<mi\>β\</mi\>\</math\> decay in a multinucleon system, specifically the \<math display\="inline"\>\<mi\>n\</mi\>\<mi\>n\</mi\>\<mo stretchy\="false"\>→\</mo\>\<mi\>p\</mi\>\<mi\>p\</mi\>\<mi\>e\</mi\>\<mi\>e\</mi\>\</math\> transition, mediated by a light left\-handed Majorana neutrino propagating over nuclear\-scale distances\. This calculation is performed with quark masses corresponding to a pion mass of \<math display\="inline"\>\<msub\>\<mi\>m\</mi\>\<mi\>π\</mi\>\</msub\>\<mo\>\=\</mo\>\<mn\>806\</mn\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>MeV\</mi\>\</math\> at a single lattice spacing and volume\. The statistically cleaner \<math display\="inline"\>\<msup\>\<mi mathvariant\="normal"\>Σ\</mi\>\<mo\>\-\</mo\>\</msup\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi mathvariant\="normal"\>Σ\</mi\>\<mo\>\+\</mo\>\</msup\>\<mi\>e\</mi\>\<mi\>e\</mi\>\</math\> transition is also computed in order to investigate various systematic uncertainties\. The prospects for matching the results of LQCD calculations onto a nuclear effective field theory to determine a leading\-order low\-energy constant relevant for \<math display\="inline"\>\<mn\>0\</mn\>\<mi\>ν\</mi\>\<mi\>β\</mi\>\<mi\>β\</mi\>\</math\> decay with a light Majorana neutrino are investigated\. This work, therefore, sets the stage for future calculations at physical values of the quark masses that, combined with effective field theory and nuclear many\-body studies, will provide controlled theoretical inputs to experimental searches of \<math display\="inline"\>\<mn\>0\</mn\>\<mi\>ν\</mi\>\<mi\>β\</mi\>\<mi\>β\</mi\>\</math\> decay\.