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
arXiv: 2202.11970
DOI: 10.1103/PhysRevD.106.014502

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

Submitted: 24 February 2022

Subjects:

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>.