We report a lattice QCD calculation of the pion scalar form factor and its
associated radii, with fully controlled systematic uncertainties. The
computation utilizes 17 gauge ensembles with
improved sea quarks, spanning a lattice spacing range of
to
the challenging quark-disconnected contributions enables an unprecedented
momentum resolution, particularly on large and fine ensembles near the
physical quark mass.
To ensure reliable extraction of ground-state matrix elements at both zero and
non-zero momentum transfer, we use a broad range of source-sink separations,
from 1 fm to 3.25 fm. For the first
time, the scalar radii are obtained using a
the
approximation at low momentum transfer.
The physical extrapolation is guided by the next-to-leading order expressions
of the scalar radii from chiral perturbation theory. The pertinent low-energy
constants are fitted from lattice data, leading to the first lattice
determination of
extraction, form-factor parametrization, and the physical extrapolation are
rigorously quantified using model averaging based on the Akaike Information
Criterion (AIC).