$I=rac{3}{2}$ $πK$ $s$-wave scattering length from lattice QCD

The $I=\frac{3}{2}$ $πK$ $s$-wave scattering phase shift is computed by lattice quantum chromodynamics with $N_f=3$ flavors of Asqtad-improved staggered fermions. The energy-eigenvalues of $πK$ systems at one center of mass frame and six moving frames using moving wall source technique are used to get phase shifts by Lüscher’s formula and its extensions. The calculations are good enough to acquire effective range expansion parameters: scattering length $a$, effective range $r$, and shape parameter $P$, which are in good agreement with our explicit analytical predictions in three-flavor chiral perturbation theory at next-to-leading order. All results are fairly consistent with experimental measurements, phenomenological studies, and lattice estimations. Numerical computations are implemented at a fine ($a\approx0.082$ fm, $L^3 T = 40^3 96$) lattice ensemble with physical quark masses.

💡 Research Summary

In this work the authors present a state‑of‑the‑art lattice QCD determination of the I = 3/2 πK s‑wave scattering parameters directly at the physical point. The calculation uses the MILC fine ensemble with three flavors of Asqtad‑improved staggered fermions (a ≈ 0.082 fm, spatial volume 40³, temporal extent 96) where the up/down and strange quark masses are tuned to their physical values (m_π ≈ 140 MeV, m_K ≈ 494 MeV).

To extract the scattering information the authors compute πK two‑point correlators in the center‑of‑mass frame and in six moving frames with total momenta P = (0,0,1), (0,1,1), (1,1,1), (0,0,2), (0,0,3) and (0,0,4) (in units of 2π/L). The moving‑wall source technique is employed: a wall source with a plane‑wave phase factor is placed on every time slice, which dramatically increases statistics (signal‑to‑noise improves as 1/√(N_slice L³)). The correlators are built from the two quark‑line diagrams that survive in the I = 3/2 channel (direct and crossed), and the staggered taste factor N_f is accounted for.

Energy levels are extracted from exponential fits to the correlators, and the corresponding relative momentum k is obtained from the relativistic dispersion relation. The Lüscher formalism and its extensions to moving frames (including the Lorentz boost factor γ) are used to relate k and the finite‑volume phase shift δ(k) via
 k cot δ(k) = 2 γ L √π Z_{d00}(1; q²).
The generalized zeta functions Z_{d00} are evaluated numerically.

Having obtained δ(k) at fourteen distinct k‑values (two levels in the CM frame and one level in each moving frame, with an extra excited level in two frames), the authors fit the data to the effective‑range expansion (ERE)
 k cot δ = 1/a + ½ r k² + P k⁴ + …,
thereby extracting simultaneously the scattering length a, the effective range r, and the shape parameter P.

On the theoretical side the paper derives explicit NLO SU(3) chiral perturbation theory (χPT) expressions for a, r and P. While the scattering length has been known analytically, the effective‑range and shape‑parameter formulas are new. Using the low‑energy constants L_i (i = 1…8) from the BE14 fit, the χPT predictions at the physical point are
 m_π a = −0.0595(62), m_π r = 16.92(4.01), P = −8.76(1.91).

The lattice fit yields
 m_π a = −0.0588(28), m_π r = 21.54(6.90), P = −6.98(3.80).
These numbers are in excellent agreement with the χPT NLO estimates, with experimental determinations based on Roy‑Steiner analyses, and with previous lattice studies (e.g., PACS‑CS, RBC‑UKQCD). The relatively large value of r compared with earlier calculations is consistent with the expectation that the effective range grows as one approaches the physical point.

The authors discuss systematic uncertainties: finite‑volume effects are suppressed by the large spatial extent (L ≈ 3.3 fm), discretization errors are expected to be small at the fine lattice spacing, and the use of multiple moving frames reduces the model dependence associated with truncating the ERE. The moving‑wall source also mitigates excited‑state contamination, as demonstrated by the stability of the extracted energies under variations of fit ranges.

In summary, this study demonstrates that a single lattice ensemble at the physical point, combined with moving‑frame Lüscher analysis and the moving‑wall source, can deliver not only the πK scattering length but also the effective range and shape parameter with competitive precision. The work validates SU(3) χPT at NLO in the strange sector and provides a benchmark for future investigations of higher partial waves, coupled‑channel dynamics, and the inclusion of isospin‑breaking effects.