# Lattice calculation of the pion mass difference \<math display\="inline"\>\<mrow\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mo\>\+\</mo\>\</mrow\>\</msup\>\</mrow\>\</msub\>\<mo\>\-\</mo\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mn\>0\</mn\>\</mrow\>\</msup\>\</mrow\>\</msub\>\</mrow\>\</math\> at order \<math display\="inline"\>\<mrow\>\<mi mathvariant\="script"\>O\</mi\>\<mo stretchy\="false"\>\(\</mo\>\<msub\>\<mrow\>\<mi\>α\</mi\>\</mrow\>\<mrow\>\<mi\>em\</mi\>\</mrow\>\</msub\>\<mo stretchy\="false"\>\)\</mo\>\</mrow\>\</math\>

**INSPIRE:** [2037634](https://inspirehep.net/literature/2037634)
**arXiv:** [2202\.11970](https://arxiv.org/abs/2202.11970)
**DOI:** [10\.1103/PhysRevD\.106\.014502](https://doi.org/10.1103/PhysRevD.106.014502)

**Authors:** Frezzotti, R\., Gagliardi, G\., Lubicz, V\., Martinelli, G\., Sanfilippo, F\., Simula, S\.

**Submitted:** 24 February 2022

**Subjects:**
- hep\-lat
- hep\-ph
- Lattice
- Phenomenology\-HEP

**Journal reference:** Phys\.Rev\.D 106 \(1\) 014502 \(2022\)

**Citations:** 12

## Abstract

We present a lattice calculation of the charged/neutral pion mass difference \<math display\="inline"\>\<mrow\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mo\>\+\</mo\>\</mrow\>\</msup\>\</mrow\>\</msub\>\<mo\>\-\</mo\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mn\>0\</mn\>\</mrow\>\</msup\>\</mrow\>\</msub\>\</mrow\>\</math\> at order \<math display\="inline"\>\<mrow\>\<mi mathvariant\="script"\>O\</mi\>\<mo stretchy\="false"\>\(\</mo\>\<msub\>\<mrow\>\<mi\>α\</mi\>\</mrow\>\<mrow\>\<mi\>em\</mi\>\</mrow\>\</msub\>\<mo stretchy\="false"\>\)\</mo\>\</mrow\>\</math\> using the gauge configurations produced by the Extended Twisted Mass Collaboration with \<math display\="inline"\>\<msub\>\<mi\>N\</mi\>\<mi\>f\</mi\>\</msub\>\<mo\>\=\</mo\>\<mn\>2\</mn\>\<mo\>\+\</mo\>\<mn\>1\</mn\>\<mo\>\+\</mo\>\<mn\>1\</mn\>\</math\> dynamical quark flavors at three values of the lattice spacing \(\<math display\="inline"\>\<mrow\>\<mi\>a\</mi\>\<mo\>≃\</mo\>\<mn\>0\.062\</mn\>\<mo\>,\</mo\>\<mn\>0\.082\</mn\>\<mo\>,\</mo\>\<mn\>0\.089\</mn\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>fm\</mi\>\</mrow\>\</math\>\) and pion masses in the range \<math display\="inline"\>\<mrow\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\</msub\>\<mo\>≃\</mo\>\<mn\>250\</mn\>\<mi\>–\</mi\>\<mn\>500\</mn\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>MeV\</mi\>\</mrow\>\</math\>\. We employ the RM123 method and expand the path integral around the isospin symmetric point at leading order in the electromagnetic coupling \<math display\="inline"\>\<mrow\>\<msub\>\<mrow\>\<mi\>α\</mi\>\</mrow\>\<mrow\>\<mi\>em\</mi\>\</mrow\>\</msub\>\</mrow\>\</math\>\. Making use of the recently proposed rotated twisted\-mass scheme, we evaluate the full \<math display\="inline"\>\<mrow\>\<mi mathvariant\="script"\>O\</mi\>\<mo stretchy\="false"\>\(\</mo\>\<msub\>\<mrow\>\<mi\>α\</mi\>\</mrow\>\<mrow\>\<mi\>em\</mi\>\</mrow\>\</msub\>\<mo stretchy\="false"\>\)\</mo\>\</mrow\>\</math\> contribution, with the inclusion of the disconnected diagram\. At the physical point, after performing the continuum and infinite volume extrapolation, we obtain the value \<math display\="inline"\>\<mrow\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mo\>\+\</mo\>\</mrow\>\</msup\>\</mrow\>\</msub\>\<mo\>\-\</mo\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mn\>0\</mn\>\</mrow\>\</msup\>\</mrow\>\</msub\>\<mo\>\=\</mo\>\<mn\>4\.622\</mn\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>95\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>MeV\</mi\>\</mrow\>\</math\> which is in good agreement with the experimental result \<math display\="inline"\>\<mrow\>\<mo stretchy\="false"\>\[\</mo\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mo\>\+\</mo\>\</mrow\>\</msup\>\</mrow\>\</msub\>\<mo\>\-\</mo\>\<msub\>\<mrow\>\<mi\>M\</mi\>\</mrow\>\<mrow\>\<msup\>\<mrow\>\<mi\>π\</mi\>\</mrow\>\<mrow\>\<mn\>0\</mn\>\</mrow\>\</msup\>\</mrow\>\</msub\>\<msup\>\<mrow\>\<mo stretchy\="false"\>\]\</mo\>\</mrow\>\<mrow\>\<mi\>exp\</mi\>\</mrow\>\</msup\>\<mo\>\=\</mo\>\<mn\>4\.5936\</mn\>\<mo stretchy\="false"\>\(\</mo\>\<mn\>5\</mn\>\<mo stretchy\="false"\>\)\</mo\>\<mtext\> \</mtext\>\<mtext\> \</mtext\>\<mi\>MeV\</mi\>\</mrow\>\</math\>\.