Pion: QCD, Structure & Precision Tests

• Pion is an isotriplet bosonic hadron—a bound state of light quark-antiquark pairs—that acts as a pseudo-Goldstone boson from chiral symmetry breaking.
• It underpins precise QCD tests by linking nonperturbative methods, decay constant measurements, and form factor determinations to Standard Model validations.
• Studies of the pion integrate lattice QCD, Bethe–Salpeter methods, and neutrino scattering models to enhance our understanding of hadron structure and medium modifications.

The pion is an isotriplet of bosonic hadrons, understood in QCD as bound states of quarks and anti-quarks of the two lightest flavours and, simultaneously, as pseudo-Goldstone bosons of spontaneously broken approximate axial-vector symmetries. They are the lightest particles of the strong-interaction spectrum, with experimental masses $m_{\pi^\pm}=139.57061\pm0.00024$ MeV and $m_{\pi^0}=134.9768\pm0.0005$ MeV, and the two-flavour normalization of the decay constant gives $f_\pi=92.09\pm0.25$ MeV (Ananthanarayan et al., 2018). Because the pion is both a Nambu-Goldstone boson and a quark-antiquark bound state, it occupies a distinctive position in the study of confinement, dynamical chiral symmetry breaking, hadron structure, precision Standard Model tests, and nuclear and neutrino interactions (Horn et al., 2016).

1. QCD origin and chiral dynamics

In the limit $m_u,m_d\to0$, the QCD Lagrangian acquires a global $SU(2)_L\times SU(2)_R$ symmetry, and the vacuum spontaneously breaks it to $SU(2)_V$, yielding three massless Goldstone bosons. For small quark masses these become light, with

$$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$

The leading-order two-flavour chiral Lagrangian is

$$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$

with $U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$ and $M_q=\mathrm{diag}(m_u,m_d)$ (Ananthanarayan et al., 2018).

A continuum-QCD description emphasizes dynamical chiral symmetry breaking (DCSB). Nonperturbative dressing of quarks generates a momentum-dependent mass function

$m_{\pi^0}=134.9768\pm0.0005$0

even when the current-quark mass $m_{\pi^0}=134.9768\pm0.0005$1. The axial-vector Ward–Green–Takahashi identity then guarantees, in the chiral limit, the appearance of a massless $m_{\pi^0}=134.9768\pm0.0005$2 bound state whose Bethe–Salpeter amplitude satisfies the quark-level Goldberger–Treiman relation

$m_{\pi^0}=134.9768\pm0.0005$3

For nonzero light-quark masses one obtains the exact Gell-Mann–Oakes–Renner relation,

$m_{\pi^0}=134.9768\pm0.0005$4

with

$m_{\pi^0}=134.9768\pm0.0005$5

Using $m_{\pi^0}=134.9768\pm0.0005$6 MeV and typical light-quark masses $m_{\pi^0}=134.9768\pm0.0005$7 MeV yields $m_{\pi^0}=134.9768\pm0.0005$8, which explains the pion’s exceptional lightness (Horn et al., 2016).

Within this framework, the pion is not merely a light meson but a direct manifestation of the mechanism that generates nearly all hadronic mass. This dual role explains why pion observables are routinely used to test chiral effective theory, lattice QCD, and nonperturbative continuum methods (Horn et al., 2016).

2. Electromagnetic, transition, and partonic structure

The charged-pion elastic form factor is defined by

$m_{\pi^0}=134.9768\pm0.0005$9

with $f_\pi=92.09\pm0.25$0 and

$f_\pi=92.09\pm0.25$1

The neutral-pion transition form factor $f_\pi=92.09\pm0.25$2 appears in

$f_\pi=92.09\pm0.25$3

with the Abelian anomaly imposing

$f_\pi=92.09\pm0.25$4

At asymptotically large spacelike $f_\pi=92.09\pm0.25$5, perturbative factorization gives

$f_\pi=92.09\pm0.25$6

At the hadronic scale $f_\pi=92.09\pm0.25$7 GeV, however, one finds a broad, concave leading-twist pion PDA,

$f_\pi=92.09\pm0.25$8

rather than the conformal-limit form $f_\pi=92.09\pm0.25$9 (Horn et al., 2016).

A Minkowski-space Bethe–Salpeter treatment makes the same point in a different language. In the model of de Melo et al., the dressed-quark propagator is

$m_u,m_d\to0$0

with

$m_u,m_d\to0$1

and the full Bethe–Salpeter amplitude admits a Nakanishi integral representation. In that construction, a $m_u,m_d\to0$2 variation of the quark self-energy parameters changes $m_u,m_d\to0$3 from $m_u,m_d\to0$4 to $m_u,m_d\to0$5 GeV, $m_u,m_d\to0$6 from $m_u,m_d\to0$7 to $m_u,m_d\to0$8 fm, and $m_u,m_d\to0$9 from $SU(2)_L\times SU(2)_R$0 to $SU(2)_L\times SU(2)_R$1 MeV; the resulting band encloses essentially all of the Jefferson Lab data up to $SU(2)_L\times SU(2)_R$2 GeV$SU(2)_L\times SU(2)_R$3 (Melo et al., 2019).

Lattice QCD has progressively moved from low moments to $SU(2)_L\times SU(2)_R$4-dependent observables. Earlier calculations with $SU(2)_L\times SU(2)_R$5 $SU(2)_L\times SU(2)_R$6-improved Wilson fermions and pion masses down to $SU(2)_L\times SU(2)_R$7 MeV found

$SU(2)_L\times SU(2)_R$8

with first determinations of $SU(2)_L\times SU(2)_R$9 and $SU(2)_V$0 at nonzero momentum transfer (Bali et al., 2013). A later physical-mass LaMET calculation, using clover valence quarks on $SU(2)_V$1 HISQ configurations with $SU(2)_V$2 fm, $SU(2)_V$3 fm, and $SU(2)_V$4 GeV, extracted the zero-skewness valence-quark GPD $SU(2)_V$5 directly. At $SU(2)_V$6 the distribution peaks near $SU(2)_V$7–$SU(2)_V$8; from its first moment one obtains $SU(2)_V$9, with the numerical example $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$0 versus experiment $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$1. The same study predicted higher generalized form factors,

$$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$2

and found that the transverse profile is broad at $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$3, with $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$4 fm, but narrows sharply by $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$5, where $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$6 fm (Lin, 2023).

The neutral-pion transition form factor remains a notable point of interpretation. One approach matches a model-independent low-energy monopole form to a QCD expression at $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$7, assumes a flat pion distribution amplitude $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$8, evolves it with ERBL evolution, and adds a small $$m_\pi^2 = 2 B_0 \hat m + O(m_q^2), \qquad \hat m=(m_u+m_d)/2.$$9 twist-3 term. In that framework, the scenario $$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$0 GeV, $$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$1 GeV, $$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$2 GeV$$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$3 gives an excellent fit over the published BaBar range up to $$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$4 GeV$$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$5 (Noguera et al., 2010). By contrast, continuum-DSE analyses describe $$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$6 as approaching $$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$7 from below with logarithmic damping and regard the BaBar excess above $$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$8 for $$L_2 = \frac{f_\pi^2}{4}\,\mathrm{Tr}[\partial_\mu U\,\partial^\mu U^\dagger] +\frac{f_\pi^2 B_0}{2}\,\mathrm{Tr}[M_q(U+U^\dagger)],$$9 GeV$U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$0 as inconsistent with Belle and with that theoretical framework (Horn et al., 2016).

3. Precision Standard Model tests and rare decays

The neutral-pion decay $U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$1 is fixed at low energy by the chiral anomaly: $U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$2 Using $U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$3 MeV, $U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$4, and $U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$5 MeV gives

$U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$6

including a $U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$7 estimate of higher corrections and isospin-violating $U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$8–$U(x)=\exp(i\,\tau\!\cdot\!\pi(x)/f_\pi)$9 mixing. The PrimEx result at JLab is

$M_q=\mathrm{diag}(m_u,m_d)$0

which agrees at the $M_q=\mathrm{diag}(m_u,m_d)$1 level (Ananthanarayan et al., 2018).

Charged-pion leptonic decays provide a separate precision arena. The tree-level ratio

$M_q=\mathrm{diag}(m_u,m_d)$2

is helicity suppressed, so radiative corrections are comparatively prominent: $M_q=\mathrm{diag}(m_u,m_d)$3 with $M_q=\mathrm{diag}(m_u,m_d)$4. The Standard Model prediction is

$M_q=\mathrm{diag}(m_u,m_d)$5

and the theory error on $M_q=\mathrm{diag}(m_u,m_d)$6 is $M_q=\mathrm{diag}(m_u,m_d)$7 ($M_q=\mathrm{diag}(m_u,m_d)$8) (Collaboration et al., 2022).

The PIONEER program is designed to exploit that sensitivity. Its design includes a $M_q=\mathrm{diag}(m_u,m_d)$9 sr, $m_{\pi^0}=134.9768\pm0.0005$00 radiation length liquid-xenon calorimeter, a segmented low gain avalanche detector stopping target, and a positron tracker, using intense pion beams at the PSI ring cyclotron. For the $m_{\pi^0}=134.9768\pm0.0005$01 measurement, the statistical goal is $m_{\pi^0}=134.9768\pm0.0005$02 $m_{\pi^0}=134.9768\pm0.0005$03 events, corresponding to $m_{\pi^0}=134.9768\pm0.0005$04 ‰ per three years at $m_{\pi^0}=134.9768\pm0.0005$05 uptime; the total systematic uncertainty is projected at $m_{\pi^0}=134.9768\pm0.0005$06 ‰, giving $m_{\pi^0}=134.9768\pm0.0005$07, which matches the Standard Model error (Collaboration et al., 2022).

Pion beta decay,

$m_{\pi^0}=134.9768\pm0.0005$08

is theoretically pristine because its rate can be written as

$m_{\pi^0}=134.9768\pm0.0005$09

with $m_{\pi^0}=134.9768\pm0.0005$10 and $m_{\pi^0}=134.9768\pm0.0005$11. Current data imply

$m_{\pi^0}=134.9768\pm0.0005$12

while a $m_{\pi^0}=134.9768\pm0.0005$13 improvement in branching-ratio precision yields $m_{\pi^0}=134.9768\pm0.0005$14, and a $m_{\pi^0}=134.9768\pm0.0005$15 improvement yields $m_{\pi^0}=134.9768\pm0.0005$16 (Collaboration et al., 2022).

The same experimental program also targets non-Standard-Model channels. PIONEER projects sensitivity at the level of $m_{\pi^0}=134.9768\pm0.0005$17 for heavy-neutrino searches in the electron sector, $m_{\pi^0}=134.9768\pm0.0005$18 in the muon sector, $m_{\pi^0}=134.9768\pm0.0005$19 for $m_{\pi^0}=134.9768\pm0.0005$20 channels involving axions, ALPs, or light dark photons, and $m_{\pi^0}=134.9768\pm0.0005$21 for lepton-flavor-violating modes (Collaboration et al., 2022).

4. Scattering, proton structure, and the hadronic environment

Pion-pion scattering is one of the most tightly constrained low-energy hadronic processes. The partial-wave amplitudes

$m_{\pi^0}=134.9768\pm0.0005$22

lead to the S-wave scattering lengths through

$m_{\pi^0}=134.9768\pm0.0005$23

Roy-equation analyses yield

$m_{\pi^0}=134.9768\pm0.0005$24

and the corresponding phase shifts agree with $m_{\pi^0}=134.9768\pm0.0005$25 and $m_{\pi^0}=134.9768\pm0.0005$26 data at the percent level (Ananthanarayan et al., 2018).

An extended linear sigma model provides a complementary dynamical interpretation of these numbers. In the two-flavour eLSM supplemented by a light scalar four-quark field $m_{\pi^0}=134.9768\pm0.0005$27 associated predominantly with $m_{\pi^0}=134.9768\pm0.0005$28 and by a scalar glueball $m_{\pi^0}=134.9768\pm0.0005$29, the tree-level $m_{\pi^0}=134.9768\pm0.0005$30 amplitude receives a four-pion contact term, $m_{\pi^0}=134.9768\pm0.0005$31-channel exchange of $m_{\pi^0}=134.9768\pm0.0005$32, $m_{\pi^0}=134.9768\pm0.0005$33, and $m_{\pi^0}=134.9768\pm0.0005$34, and $m_{\pi^0}=134.9768\pm0.0005$35-channel $m_{\pi^0}=134.9768\pm0.0005$36 exchange. A global fit to scalar masses, decay widths, and scattering lengths gives

$m_{\pi^0}=134.9768\pm0.0005$37

If one sets $m_{\pi^0}=134.9768\pm0.0005$38, thereby removing the four-quark field, then $m_{\pi^0}=134.9768\pm0.0005$39 drops to $m_{\pi^0}=134.9768\pm0.0005$40, which is too small. In that model the inclusion of $m_{\pi^0}=134.9768\pm0.0005$41 improves the description of both $m_{\pi^0}=134.9768\pm0.0005$42 and $m_{\pi^0}=134.9768\pm0.0005$43 scattering lengths, whereas the heavy glueball has only a marginal effect on low-energy $m_{\pi^0}=134.9768\pm0.0005$44 scattering (Lakaschus et al., 2018).

The pion also appears as an explicit degree of freedom in proton structure. In a light-front description, the proton state is expanded as

$m_{\pi^0}=134.9768\pm0.0005$45

and the pion-cloud contributions enter deep-inelastic scattering and Drell–Yan cross sections through convolutions of splitting functions with pion and baryon PDFs. For the integrated isovector asymmetry,

$m_{\pi^0}=134.9768\pm0.0005$46

the theory gives $m_{\pi^0}=134.9768\pm0.0005$47, compared with the E866 result $m_{\pi^0}=134.9768\pm0.0005$48. In the SeaQuest kinematics, the theory band reproduces the measured $m_{\pi^0}=134.9768\pm0.0005$49 ratio: at $m_{\pi^0}=134.9768\pm0.0005$50, $m_{\pi^0}=134.9768\pm0.0005$51 versus $m_{\pi^0}=134.9768\pm0.0005$52–$m_{\pi^0}=134.9768\pm0.0005$53; at $m_{\pi^0}=134.9768\pm0.0005$54, $m_{\pi^0}=134.9768\pm0.0005$55 versus $m_{\pi^0}=134.9768\pm0.0005$56–$m_{\pi^0}=134.9768\pm0.0005$57; and at $m_{\pi^0}=134.9768\pm0.0005$58, $m_{\pi^0}=134.9768\pm0.0005$59 versus $m_{\pi^0}=134.9768\pm0.0005$60–$m_{\pi^0}=134.9768\pm0.0005$61, with overall $m_{\pi^0}=134.9768\pm0.0005$62 across the measured range (Alberg et al., 2021).

These results place the pion in two distinct but related roles: as the key low-energy carrier of long-range hadronic dynamics, including the long-range part of the $m_{\pi^0}=134.9768\pm0.0005$63 force, and as a nonperturbative component of the proton sea that affects deep-inelastic and Drell–Yan observables (Alberg et al., 2021).

5. Nuclear matter, hot dense media, and neutrino interactions

In symmetric nuclear matter, the quark-meson coupling model predicts a substantial modification of pion structure. At the hadron level the nucleon effective mass is written as

$m_{\pi^0}=134.9768\pm0.0005$64

while at the quark level

$m_{\pi^0}=134.9768\pm0.0005$65

Using a light-front Bethe–Salpeter wave function for the pion, the model finds that $m_{\pi^0}=134.9768\pm0.0005$66 falls from $m_{\pi^0}=134.9768\pm0.0005$67 MeV in vacuum to about $m_{\pi^0}=134.9768\pm0.0005$68 MeV at $m_{\pi^0}=134.9768\pm0.0005$69 fm$m_{\pi^0}=134.9768\pm0.0005$70. Over the same density range, $m_{\pi^0}=134.9768\pm0.0005$71 decreases from $m_{\pi^0}=134.9768\pm0.0005$72 MeV to $m_{\pi^0}=134.9768\pm0.0005$73 MeV, the charge radius $m_{\pi^0}=134.9768\pm0.0005$74 increases from $m_{\pi^0}=134.9768\pm0.0005$75 fm to $m_{\pi^0}=134.9768\pm0.0005$76 fm, and the valence probability $m_{\pi^0}=134.9768\pm0.0005$77 increases from $m_{\pi^0}=134.9768\pm0.0005$78 to $m_{\pi^0}=134.9768\pm0.0005$79 (Melo et al., 2014).

In hot dense matter, leptonic pion decay can be blocked kinematically by chemical potentials. The in-medium decay widths obey

$m_{\pi^0}=134.9768\pm0.0005$80

with $m_{\pi^0}=134.9768\pm0.0005$81. At zero temperature, a charged pion is stable whenever $m_{\pi^0}=134.9768\pm0.0005$82; in particular, $m_{\pi^0}=134.9768\pm0.0005$83 is stable if $m_{\pi^0}=134.9768\pm0.0005$84, and $m_{\pi^0}=134.9768\pm0.0005$85 is stable if $m_{\pi^0}=134.9768\pm0.0005$86. Under $m_{\pi^0}=134.9768\pm0.0005$87-equilibrium, the corresponding metastability thresholds are

$m_{\pi^0}=134.9768\pm0.0005$88

In this regime, the antineutrino emissivity behaves as $m_{\pi^0}=134.9768\pm0.0005$89 at high $m_{\pi^0}=134.9768\pm0.0005$90, and the additional neutrino emission can accelerate the early cooling of protoneutron-star cores for temperatures above $m_{\pi^0}=134.9768\pm0.0005$91 MeV (Loewe et al., 2011).

Pion production and propagation are equally important in accelerator neutrino physics. In the few-GeV region relevant for DUNE, the total charged-current $m_{\pi^0}=134.9768\pm0.0005$92-nucleon cross section is decomposed as

$m_{\pi^0}=134.9768\pm0.0005$93

Once produced, pions undergo absorption, charge exchange, and elastic or inelastic scattering in the residual nucleus. For $m_{\pi^0}=134.9768\pm0.0005$94-Ar CC interactions in one million simulated events, GENIE gives primary multiplicities $m_{\pi^0}=134.9768\pm0.0005$95 $m_{\pi^0}=134.9768\pm0.0005$96, $m_{\pi^0}=134.9768\pm0.0005$97 $m_{\pi^0}=134.9768\pm0.0005$98, $m_{\pi^0}=134.9768\pm0.0005$99 $f_\pi=92.09\pm0.25$00, which become final-state $f_\pi=92.09\pm0.25$01 $f_\pi=92.09\pm0.25$02, $f_\pi=92.09\pm0.25$03 $f_\pi=92.09\pm0.25$04, $f_\pi=92.09\pm0.25$05 $f_\pi=92.09\pm0.25$06 after FSI. NuWro gives primary $f_\pi=92.09\pm0.25$07 $f_\pi=92.09\pm0.25$08, $f_\pi=92.09\pm0.25$09 $f_\pi=92.09\pm0.25$10, $f_\pi=92.09\pm0.25$11 $f_\pi=92.09\pm0.25$12, becoming final $f_\pi=92.09\pm0.25$13 $f_\pi=92.09\pm0.25$14, $f_\pi=92.09\pm0.25$15 $f_\pi=92.09\pm0.25$16, $f_\pi=92.09\pm0.25$17 $f_\pi=92.09\pm0.25$18. For single-pion primaries, an initial $f_\pi=92.09\pm0.25$19 in GENIE emerges as $f_\pi=92.09\pm0.25$20 in $f_\pi=92.09\pm0.25$21 of cases, is absorbed in $f_\pi=92.09\pm0.25$22, converts to $f_\pi=92.09\pm0.25$23 in $f_\pi=92.09\pm0.25$24, and to $f_\pi=92.09\pm0.25$25 in $f_\pi=92.09\pm0.25$26; in NuWro the corresponding fractions are $f_\pi=92.09\pm0.25$27, $f_\pi=92.09\pm0.25$28, $f_\pi=92.09\pm0.25$29, and $f_\pi=92.09\pm0.25$30. For an initial $f_\pi=92.09\pm0.25$31, GENIE gives $f_\pi=92.09\pm0.25$32 survival as $f_\pi=92.09\pm0.25$33, $f_\pi=92.09\pm0.25$34 absorption, and $f_\pi=92.09\pm0.25$35 conversion to $f_\pi=92.09\pm0.25$36, whereas NuWro gives $f_\pi=92.09\pm0.25$37, $f_\pi=92.09\pm0.25$38, and $f_\pi=92.09\pm0.25$39 (Devi et al., 2022).

These FSI have direct oscillation consequences. “Pionless” events arising from RES/DIS plus absorption mimic true QE topologies, biasing the reconstructed neutrino energy downward by the undetected pion’s energy, typically $f_\pi=92.09\pm0.25$40–$f_\pi=92.09\pm0.25$41 MeV. Charge exchange, such as $f_\pi=92.09\pm0.25$42, injects neutral pions into final states relevant to $f_\pi=92.09\pm0.25$43–$f_\pi=92.09\pm0.25$44 separation. The systematic impact on the energy scale is at the $f_\pi=92.09\pm0.25$45–$f_\pi=92.09\pm0.25$46 level, FSI-rate uncertainties are $f_\pi=92.09\pm0.25$47–$f_\pi=92.09\pm0.25$48, and DUNE systematics typically assign $f_\pi=92.09\pm0.25$49 MeV energy-reconstruction biases from FSI (Devi et al., 2022).

A recent generator development, Achilles, integrates the electroweak single-pion production vertex from a Dynamical Coupled-Channels model with realistic hole spectral functions and a semi-classical intranuclear cascade. In that framework the collision probability is

$f_\pi=92.09\pm0.25$50

and pion absorption can be treated either through an Oset–Salcedo optical potential or by explicit resonance propagation. Validation studies report agreement with $f_\pi=92.09\pm0.25$51-C reaction data to better than $f_\pi=92.09\pm0.25$52 for $f_\pi=92.09\pm0.25$53 MeV, reproduction of $f_\pi=92.09\pm0.25$54-absorption cross sections on C and Ar to within $f_\pi=92.09\pm0.25$55–$f_\pi=92.09\pm0.25$56 mb in the $f_\pi=92.09\pm0.25$57-peak region, JLab $f_\pi=92.09\pm0.25$58 cross sections on $f_\pi=92.09\pm0.25$59C and $f_\pi=92.09\pm0.25$60Ar described to better than $f_\pi=92.09\pm0.25$61, and neutrino-nucleus observables typically reproduced at the $f_\pi=92.09\pm0.25$62–$f_\pi=92.09\pm0.25$63 level, with some regions underpredicted by $f_\pi=92.09\pm0.25$64–$f_\pi=92.09\pm0.25$65 where MEC or DIS dominate (Isaacson et al., 26 Aug 2025).

6. Outstanding tensions and current directions

Several active directions follow directly from current pion research. On the structure side, experiments planned at JLab 12 aim to measure $f_\pi=92.09\pm0.25$66 up to $f_\pi=92.09\pm0.25$67–$f_\pi=92.09\pm0.25$68 GeV$f_\pi=92.09\pm0.25$69, while $f_\pi=92.09\pm0.25$70 measurements up to $f_\pi=92.09\pm0.25$71 GeV$f_\pi=92.09\pm0.25$72 are expected to sharpen tests of the onset of hard scattering, factorization theorems, and the flavour dependence of DCSB. The same program is intended to clarify the transition between nonperturbative and perturbative regimes and to settle the remaining disagreement between the BaBar and Belle trends for the neutral-pion transition form factor (Horn et al., 2016).

On the lattice side, the first physical-mass $f_\pi=92.09\pm0.25$73-dependent pion GPD already provides transverse-plane tomography, but its stated next steps are simulations at multiple lattice spacings, larger $f_\pi=92.09\pm0.25$74, and nonzero skewness $f_\pi=92.09\pm0.25$75. These are directed at reducing discretization and power-correction uncertainties and at extending the framework to kaon and nucleon GPDs (Lin, 2023).

In weak-interaction and precision-frontier studies, PIONEER is explicitly motivated by the gap between the Standard Model prediction for $f_\pi=92.09\pm0.25$76 and current experimental precision, by the possibility of extracting $f_\pi=92.09\pm0.25$77 from pion beta decay with competitive or superior accuracy, and by rare-decay searches for sterile neutrinos, axions, dark photons, and lepton-flavor violation (Collaboration et al., 2022).

In nuclear and neutrino physics, the current bottleneck is the control of pion final-state interactions. The DUNE-motivated GENIE versus NuWro comparison identifies strong sensitivity to absorption and charge exchange, and its summary states that future improvements hinge on better $f_\pi=92.09\pm0.25$78–Ar scattering data and unified generator tuning (Devi et al., 2022). The Achilles program points in the same direction by combining a unitary DCC description of the electroweak vertex with a cascade that uses the same meson-baryon amplitudes during propagation, while indicating that incorporation of two-body currents and DIS remains necessary for a complete neutrino-oscillation-grade generator (Isaacson et al., 26 Aug 2025).

Across these domains, the pion remains a uniquely constraining object: light because of chiral symmetry and its breaking pattern, structurally rich because of nonperturbative QCD, and experimentally central because its decays, form factors, scattering amplitudes, and in-medium propagation connect hadron structure, Standard Model precision tests, and neutrino- and nuclear-physics systematics.