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
arXiv: 2508.01900
DOI: 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:

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.