First direct lattice calculation of the chiral perturbation theory low-energy constant ℓ7

INSPIRE: 1892216
arXiv: 2107.11895
DOI: 10.1103/PhysRevD.104.074513

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

Submitted: 25 July 2021

Subjects:

Journal reference: Phys.Rev.D 104 (7) 074513 (2021)

Citations: 10

Abstract

We evaluate by means of lattice QCD calculations the low-energy constant <math display="inline"><msub><mo>ℓ</mo><mn>7</mn></msub></math> which parametrizes strong isospin effects at next-to-leading order (NLO) in SU(2) chiral perturbation theory. Among all low-energy constants at NLO, <math display="inline"><msub><mo>ℓ</mo><mn>7</mn></msub></math> is the one known less precisely, and its uncertainty is currently larger than 50%. Our strategy is based on the RM123 approach in which the lattice path-integral is expanded in powers of the isospin breaking parameter <math display="inline"><mi mathvariant="normal">Δ</mi><mi>m</mi><mo>=</mo><mo stretchy="false">(</mo><msub><mi>m</mi><mi>d</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo stretchy="false">)</mo><mo>/</mo><mn>2</mn></math>. In order to evaluate the relevant lattice correlators we make use of the recently proposed rotated twisted-mass (RTM) scheme. Within the RM123 approach, it is possible to cleanly extract the value of <math display="inline"><msub><mo>ℓ</mo><mn>7</mn></msub></math> from either the pion mass splitting <math display="inline"><msub><mi>M</mi><msup><mi>π</mi><mo>+</mo></msup></msub><mo>-</mo><msub><mi>M</mi><msup><mi>π</mi><mn>0</mn></msup></msub></math> induced by strong isospin breaking at order <math display="inline"><mi mathvariant="script">O</mi><mrow><mo stretchy="false">(</mo><mo stretchy="false">(</mo><mi mathvariant="normal">Δ</mi><mi>m</mi><msup><mo stretchy="false">)</mo><mn>2</mn></msup><mo stretchy="false">)</mo></mrow></math> (mass method), or from the coupling of the neutral pion <math display="inline"><msup><mi>π</mi><mn>0</mn></msup></math> to the isoscalar operator <math display="inline"><mrow><mo stretchy="false">(</mo><mover accent="true"><mi>u</mi><mo stretchy="false">¯</mo></mover><msub><mi>γ</mi><mn>5</mn></msub><mi>u</mi><mo>+</mo><mover accent="true"><mi>d</mi><mo stretchy="false">¯</mo></mover><msub><mi>γ</mi><mn>5</mn></msub><mi>d</mi><mo stretchy="false">)</mo></mrow><mo>/</mo><msqrt><mn>2</mn></msqrt></math> at order <math display="inline"><mi mathvariant="script">O</mi><mo stretchy="false">(</mo><mi mathvariant="normal">Δ</mi><mi>m</mi><mo stretchy="false">)</mo></math> (matrix element method). In this pilot study we limit the analysis to a single ensemble generated by the Extended Twisted Mass Collaboration (ETMC) 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, which corresponds to a lattice spacing <math display="inline"><mi>a</mi><mo>≃</mo><mn>0.095</mn><mtext> </mtext><mtext> </mtext><mi>fm</mi></math> and to a pion mass <math display="inline"><msub><mi>M</mi><mi>π</mi></msub><mo>≃</mo><mn>260</mn><mtext> </mtext><mtext> </mtext><mi>MeV</mi></math>. We find that the matrix element method outperforms the mass method in terms of resulting statistical accuracy. Our determination, <math display="inline"><msub><mo>ℓ</mo><mn>7</mn></msub><mo>=</mo><mn>2.5</mn><mo stretchy="false">(</mo><mn>1.4</mn><mo stretchy="false">)</mo><mo>×</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>3</mn></mrow></msup></math>, is in agreement and improves previous calculations.