# Lattice calculation of short\-range contributions to neutrinoless double\-beta decay \<math display\="inline"\>\<msup\>\<mi\>π\</mi\>\<mo\>\-\</mo\>\</msup\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>π\</mi\>\<mo\>\+\</mo\>\</msup\>\<mi\>e\</mi\>\<mi\>e\</mi\>\</math\> at physical pion mass

**INSPIRE:** [2957151](https://inspirehep.net/literature/2957151)
**arXiv:** [2508\.01900](https://arxiv.org/abs/2508.01900)
**DOI:** [10\.1103/pf3g\-kmcr](https://doi.org/10.1103/pf3g-kmcr)

**Authors:** Boyle, Peter, Erben, Felix, Feng, Xu, Flynn, Jonathan M\., Garron, Nicolas, Izubuchi, Taku, Jin, Luchang, Mukherjee, Rajnandini, Tsang, J\. Tobias, Tuo, Xin\-Yu

**Submitted:** 3 August 2025

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

**Journal reference:** Phys\.Rev\.D 113 \(5\) 054504 \(2026\)

**Citations:** 1

## Abstract

Neutrinoless double\-beta \(\<math display\="inline"\>\<mn\>0\</mn\>\<mi\>ν\</mi\>\<mi\>β\</mi\>\<mi\>β\</mi\>\</math\>\) decays provide an excellent probe for determining whether neutrinos are Dirac or Majorana fermions\. The short\-range matrix elements associated with the \<math display\="inline"\>\<msup\>\<mi\>π\</mi\>\<mo\>\-\</mo\>\</msup\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>π\</mi\>\<mo\>\+\</mo\>\</msup\>\<mi\>e\</mi\>\<mi\>e\</mi\>\</math\> process contribute at leading order in the \<math display\="inline"\>\<mn\>0\</mn\>\<mi\>ν\</mi\>\<mi\>β\</mi\>\<mi\>β\</mi\>\</math\> decay channel \<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\> through pion exchange between nucleons\. However, current lattice calculations show notable discrepancies in predicting these short\-range contributions\. To address this issue, we perform a lattice QCD calculation of the \<math display\="inline"\>\<msup\>\<mi\>π\</mi\>\<mo\>\-\</mo\>\</msup\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>π\</mi\>\<mo\>\+\</mo\>\</msup\>\<mi\>e\</mi\>\<mi\>e\</mi\>\</math\> matrix elements using domain wall fermion ensembles at the physical pion mass generated by the RBC/UKQCD Collaboration\. To mitigate contamination from around\-the\-world effects, we develop a new method to reconstruct and subtract them directly from lattice data\. We then perform nonperturbative renormalization in the regularization\-independent symmetric momentum\-subtraction scheme \(RI/SMOM\), using the \<math display\="inline"\>\<mo stretchy\="false"\>\(\</mo\>\<msub\>\<mi\>γ\</mi\>\<mi\>μ\</mi\>\</msub\>\<mo\>,\</mo\>\<msub\>\<mi\>γ\</mi\>\<mi\>μ\</mi\>\</msub\>\<mo stretchy\="false"\>\)\</mo\>\</math\> and \<math display\="inline"\>\<mo stretchy\="false"\>\(\</mo\>\<menclose notation\="updiagonalstrike" other\="updiag1"\>\<mi\>q\</mi\>\</menclose\>\<mo\>,\</mo\>\<menclose notation\="updiagonalstrike" other\="updiag1"\>\<mi\>q\</mi\>\</menclose\>\<mo stretchy\="false"\>\)\</mo\>\</math\> projectors\. Compared with previous studies, this work reduces the uncertainties in the matrix elements and provides an independent cross\-check that helps to reconcile the discrepancies among previous lattice calculations\.