# Lattice QCD Calculation of \<math display\="inline"\>\<msup\>\<mi\>π\</mi\>\<mn\>0\</mn\>\</msup\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>e\</mi\>\<mo\>\+\</mo\>\</msup\>\<msup\>\<mi\>e\</mi\>\<mo\>\-\</mo\>\</msup\>\</math\> Decay

**INSPIRE:** [2133490](https://inspirehep.net/literature/2133490)
**arXiv:** [2208\.03834](https://arxiv.org/abs/2208.03834)
**DOI:** [10\.1103/PhysRevLett\.130\.191901](https://doi.org/10.1103/PhysRevLett.130.191901)

**Authors:** Christ, Norman, Feng, Xu, Jin, Luchang, Tu, Cheng, Zhao, Yidi

**Submitted:** 7 August 2022

**Subjects:**
- hep\-lat
- Lattice

**Journal reference:** Phys\.Rev\.Lett\. 130 \(19\) 191901 \(2023\)

**Citations:** 27

## Abstract

We extend the application of lattice QCD to the two\-photon\-mediated, order \<math display\="inline"\>\<msup\>\<mi\>α\</mi\>\<mn\>2\</mn\>\</msup\>\</math\> rare decay \<math display\="inline"\>\<msup\>\<mi\>π\</mi\>\<mn\>0\</mn\>\</msup\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>e\</mi\>\<mo\>\+\</mo\>\</msup\>\<msup\>\<mi\>e\</mi\>\<mo\>\-\</mo\>\</msup\>\</math\>\. By combining Minkowski\- and Euclidean\-space methods we are able to calculate the complex amplitude describing this decay directly from the underlying theories \(QCD and QED\) which predict this decay\. The leading connected and disconnected diagrams are considered; a continuum limit is evaluated and the systematic errors are estimated\. We find \<math display\="inline"\>\<mrow\>\<mi\>Re\</mi\>\<mi mathvariant\="script"\>A\</mi\>\<mo\>\=\</mo\>\<mn\>18\.60\</mn\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>1\.19\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>1\.05\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>eV\</mi\>\</mrow\>\</math\>, \<math display\="inline"\>\<mi\>Im\</mi\>\<mi mathvariant\="script"\>A\</mi\>\<mo\>\=\</mo\>\<mn\>32\.59\</mn\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>1\.50\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>1\.65\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>eV\</mi\>\</math\>, a more accurate value for the ratio \<math display\="inline"\>\<mrow\>\<mo stretchy\="false"\>\(\</mo\>\<mi\>Re\</mi\>\<mi mathvariant\="script"\>A\</mi\>\<mo\>/\</mo\>\<mi\>Im\</mi\>\<mi mathvariant\="script"\>A\</mi\>\<mo stretchy\="false"\>\)\</mo\>\<mo\>\=\</mo\>\<mn\>0\.571\</mn\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>10\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>4\</mn\>\<mo stretchy\="false"\>\)\</mo\>\</mrow\>\</math\>, and a result for the partial width \<math display\="inline"\>\<mrow\>\<mi mathvariant\="normal"\>Γ\</mi\>\<mo stretchy\="false"\>\(\</mo\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mn\>0\</mn\>\</mrow\>\</msup\>\<mo stretchy\="false"\>→\</mo\>\<mi\>γ\</mi\>\<mi\>γ\</mi\>\<mo stretchy\="false"\>\)\</mo\>\<mo\>\=\</mo\>\<mn\>6\.60\</mn\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>0\.61\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>0\.67\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>eV\</mi\>\</mrow\>\</math\>\. Here the first errors are statistical and the second systematic\. This calculation is the first step in determining the more challenging, two\-photon\-mediated decay amplitude that contributes to the rare decay \<math display\="inline"\>\<mi\>K\</mi\>\<mo stretchy\="false"\>→\</mo\>\<msup\>\<mi\>μ\</mi\>\<mo\>\+\</mo\>\</msup\>\<msup\>\<mi\>μ\</mi\>\<mo\>\-\</mo\>\</msup\>\</math\>\.