---
title: Pseudoscalar Dominance in QCD
url: https://www.emergentmind.com/topics/pseudoscalar-dominance
type: topic
---

# Pseudoscalar Dominance in QCD

Pseudoscalar dominance denotes a family of saturation hypotheses for amplitudes, form factors, and susceptibilities carrying pseudoscalar quantum numbers. In its most developed contemporary form, it is a statement about the nucleon pseudoscalar density: the pion pole supplies the dominant contribution, while the lowest excited pseudoscalar states provide the minimal corrections required by analyticity, chiral symmetry, and pQCD short-distance behavior. In related settings, the same phrase can denote pion-pole saturation of the induced pseudoscalar response, low-temperature saturation of pseudoscalar susceptibilities by the pion, or more loosely the direct dynamical involvement of the pseudoscalar sector in effective interactions [2306.05267].

## 1. Formal meaning and basic relations

In the nucleon axial sector, the standard decomposition of the isovector axial current introduces the axial and induced pseudoscalar form factors,
$$
\langle N(p',s') \vert A_\mu\vert N(p,s) \rangle
=
\bar{u}_N(p',s')
\bigg[
\gamma_\mu G_A(Q^2)-\frac{Q_\mu}{2m_N}G_P(Q^2)
\bigg]\gamma_5 u_N(p,s),
$$
while the pseudoscalar density defines
$$
\langle N(p',s') \vert P \vert N(p,s) \rangle
=
G_5(Q^2)\,\bar{u}_N(p',s')\gamma_5 u_N(p,s).
$$
These quantities are related by the nucleon-level PCAC identity
$$
G_A(Q^2)-\frac{Q^2}{4m_N^2}G_P(Q^2)=\frac{m_q}{m_N}G_5(Q^2),
$$
so any claim of pseudoscalar dominance is a claim about how the pseudoscalar channel saturates this relation and the associated pole structure [2502.07583].

Two related but distinct notions are standard. “Pion-pole dominance” (PPD) is the approximation
$$
G_P(Q^2)\simeq \frac{4m_N^2G_A(Q^2)}{Q^2+m_\pi^2},
$$
where the longitudinal axial response is controlled by a single pion pole. A broader pseudoscalar-dominance statement uses
$$
G_5(Q^2)=\frac{m_\pi^2}{Q^2+m_\pi^2}\frac{f_\pi}{m_q}G_{\pi NN}(Q^2),
$$
so that the pion-nucleon form factor \(G_{\pi NN}\), the residue of the pseudoscalar channel at the pion pole, and the Goldberger-Treiman relation are treated within one common framework. In this usage, the pion pole is dominant but not necessarily exhaustive, because regular non-pole pieces and excited pseudoscalar states can generate controlled deviations [2103.02054].

## 2. Dispersive nucleon formulation and Extended PCAC

The 2023 dispersive formulation starts from the nucleon matrix element of the pseudoscalar quark density,
$$
\langle N(p')|\bar q\,\{\tau,\hat m\}\, i\gamma_5\, q|N(p)\rangle
=
\bar u(p')\,\{\tau,\hat m\}\, i\gamma_5\, u(p)\, F_P(t),
\qquad t=q^2,
$$
and emphasizes the analytic structure of \(F_P(t)\): it is real in the spacelike region, has a pion pole at \(t=m_\pi^2\), branch cuts beginning at \(3\pi,5\pi,\dots\), and must fall at large Euclidean momentum in a way compatible with pQCD. The corresponding dispersion relation makes the pion pole explicit,
$$
\hat F_P(t)=\frac{2 f_\pi m_\pi g_{\pi NN}}{m_\pi^2-t}
+\frac{1}{\pi}\int_{(3m_\pi)^2}^{\infty}ds\,\frac{\mathrm{Im}\,\hat F_P(s)}{s-t},
$$
with the residue defining \(g_{\pi NN}\). Neglect of the continuum yields the Goldberger-Treiman relation \(M_N g_A=f_\pi g_{\pi NN}\), while the observed nonzero discrepancy motivates a structured contribution beyond the pion alone [2306.05267].

The 2025 dispersive analysis recasts the same content in terms of
$$
m_q F_P(s)=\frac{M_\pi^2}{M_\pi^2-s}\,D(s),
$$
with
$$
D(0)=M_N g_A,\qquad D(M_\pi^2)=F_\pi g_{\pi NN},
$$
so that
$$
\Delta_{\rm GT}=1-\frac{M_N g_A}{F_\pi g_{\pi NN}}
=1-\frac{D(0)}{D(M_\pi^2)}.
$$
This is presented as a modern realization of extended PCAC: the divergence of the axial current contains not only the pion field but also higher pseudoscalar states,
$$
\partial^\mu A_\mu \sim F_\pi M_\pi^2 \phi_\pi + F_{\pi'} M_{\pi'}^2 \phi_{\pi'}+\cdots,
$$
which in the large-\(N_c\) limit becomes an infinite tower of narrow poles. The practical claim is not unrestricted pole proliferation, but minimal pseudoscalar saturation subject to analyticity, PCAC, and pQCD asymptotics [2507.22726].

## 3. Spectral constraints and the Goldberger-Treiman discrepancy

A central result of the modern dispersive program is that pure pion dominance is insufficient once the Goldberger-Treiman discrepancy and QCD asymptotics are imposed simultaneously. The 2025 analysis derives superconvergent sum rules from the deep-Euclidean behavior
$$
m_q F_P(-Q^2)\sim \frac{m_q\Lambda_N^4}{Q^4}\alpha_s^2(Q^2),
$$
and concludes that the spectral density must have at least one zero. Low-energy ChPT fixes the threshold behavior but contributes only
$$
\Delta_{\rm GT}^{\rm ChPT}\sim 0.02\%-0.3\%,
$$
far below the observed discrepancy. The dominant correction is instead assigned to the first excited isovector pseudoscalar, \(\pi(1300)\), with a negative residue \(Z_{\pi'}<0\), while a Regge-like tail provides the intermediate/high-energy continuation. In this framework,
$$
g_{\pi NN}=13.14(^{+7}_{-4}),\qquad
\Delta_{\rm GT}=1.26(^{+51}_{-34})\%,
$$
and the approximate relations
$$
g_{\pi NN}\simeq \frac{g_A F_\pi}{M_N}\frac{1}{1-M_\pi^2/M_{\pi'}^2},
\qquad
\Delta_{\rm GT}\simeq \frac{M_\pi^2}{M_{\pi'}^2}
$$
make explicit that the correction is controlled primarily by the first excited pseudoscalar state [2507.22726].

The 2023 large-\(N_c\)-inspired minimal ansatz implements the same idea with the pion and the two lightest excited isovector pseudoscalars,
$$
F_P(t)=2M_N g_A
\frac{1}{\left(1-\frac{t}{m_\pi^2}\right)
\left(1-\frac{t}{m_{\pi'}^2}\right)
\left(1-\frac{t}{m_{\pi''}^2}\right)}.
$$
From the pion-pole residue it obtains
$$
g_{\pi NN}=13.21^{+0.11}_{-0.06},
$$
to be compared with the Granada-2013 scattering benchmark \(g_{\pi^+np}=13.25(5)\), and equivalently
$$
\Delta_{\rm GT}=1.8^{+0.9}_{-0.4}\%.
$$
The result is interpreted as evidence that the strong pion-nucleon-nucleon vertex is almost flat: the pion pole dominates the form factor, while higher pseudoscalar resonances generate the small but necessary correction [2306.05267].

## 4. Lattice QCD and continuum tests

Direct physical-point lattice calculations have produced both tension and support for pseudoscalar dominance, depending on how cutoff effects and excited states are controlled. A 2020 study found that \(G_P(Q^2)\) and \(G_5(Q^2)\) show the expected strong low-\(Q^2\) enhancement, with \(G_P(0)\sim 237\) and \(G_5(0)\sim 313\) on the main ensemble, but the diagnostic ratios
$$
r_{\rm PCAC}=
\frac{\frac{m_q}{m_N}G_5(Q^2)+\frac{Q^2}{4m_N^2}G_P(Q^2)}{G_A(Q^2)},
\qquad
r_{\rm PPD}=
\frac{G_P(Q^2)}{\frac{4m_N^2}{m_\pi^2+Q^2}G_A(Q^2)}
$$
show sizeable deviations from unity at small \(Q^2\), even after an improved excited-state treatment. The interpretation given there is that low-energy \(\pi N\) states enter the three-point functions in a manner not fully captured by simpler two-state analyses. At the same time, indirect use of \(G_A\) together with PCAC/PPD gives phenomenologically reasonable values,
\(g_P^\ast=8.69(14)\), \(g_{\pi NN}=13.48(27)(2)\), and \(\Delta_{GT}=0.0276(38)(17)\) [2011.13342].

A later continuum-limit calculation with three \(N_f=2+1+1\) twisted-mass ensembles at the physical point identifies the dominant obstruction as a lattice artifact. Using lattice spacings \(a=0.07957,\ 0.06821,\ 0.05692\ \mathrm{fm}\), multiple source-sink separations from about \(0.5\) fm to \(1.5\) fm, two-state and three-state fits, and AIC model averaging, it finds significant \(\mathcal O(a^2)\) effects in both \(r_{\rm PCAC}\) and \(r_{\rm PPD}\) at finite lattice spacing, but restoration of both relations in the continuum limit:
\(r_{\rm PCAC}(Q^2)\to 1\) and \(r_{\rm PPD}(Q^2)\to 1\). The quoted continuum results are
$$
g_A = 1.245(28)(14),\qquad
\langle r_A^2\rangle = 0.339(48)(06)\ \mathrm{fm}^2,
$$
$$
g_P^* = 8.99(39)(49),\qquad
g_{\pi NN} = 13.25(67)(69),
$$
which align the lattice data with the standard pseudoscalar-dominance picture once discretization effects are removed [2502.07583].

A symmetry-preserving continuum quark+diquark calculation reaches a similar conclusion from a different direction. There, exact satisfaction of PCAC requires consistent construction of the external axial and pseudoscalar currents, including the axial and pseudoscalar seagull couplings. The directly computed \(G_P(Q^2)\) is then extremely close to the PPD estimate,
with the ratio \(R_{\rm PPD}=\frac{4m_N^2G_A}{(Q^2+m_\pi^2)G_P}\) staying near unity and differing by only about \(1\%\) near \(Q^2\simeq 0\). The same framework gives
\(g_A=1.25(3)\), \(g_p^\ast=8.80(23)\), \(g_{\pi NN}/m_N=14.02(33)/{\rm GeV}\), and \(\Delta_{\rm GT}=0.030(1)\), supporting the view that pion-pole dominance is an excellent approximation when current identities are implemented exactly [2103.02054].

## 5. Other QCD realizations of pseudoscalar-channel saturation

In finite-temperature QCD, pseudoscalar dominance appears in the low-temperature saturation of the pseudoscalar susceptibility by the pion pole. One-loop \(N_f=2\) ChPT gives
$$
\chi_P^{\rm ChPT}(T)= -\frac{\langle\bar q q\rangle_T}{m_q}+\mathcal O(F^{-2}),
$$
so the pseudoscalar susceptibility tracks the chiral order parameter throughout the low-\(T\) regime. This behavior is then linked to scalar-pseudoscalar partner degeneration near the chiral transition, with \(\chi_S(T)\simeq \chi_P(T)\) around restoration and the scalar channel described by saturation with the thermal \(\sigma/f_0(500)\) resonance [1304.3356].

At the quark level, a Landau-gauge Schwinger-Dyson analysis finds the dressed pseudoscalar density to be large for light quarks and nonconvergent in the chiral limit. The quoted renormalized zero-momentum values are
\(P(0)_{\mu=2\,\mathrm{GeV}}=177\) for \(m_q=2.2\) MeV,
\(86.3\) for \(m_q=4.8\) MeV, and
\(7.76\) for \(m_q=95\) MeV.
This behavior is explicitly interpreted as pion-pole dominance: the pseudoscalar vertex is enhanced by the ladder-generated Nambu-Goldstone mode and inherits the \(1/m_\pi^2\) singularity expected from current algebra and PCAC [1401.2852].

In the \(\Delta(1232)\) sector, lattice QCD yields an analogous but richer structure because spin-\(\tfrac32\) matrix elements involve four axial and two pseudoscalar form factors. Two Goldberger-Treiman-type relations are derived,
$$
f_\pi G_{\pi\Delta\Delta}(0)=m_\Delta g_1(0),\qquad
f_\pi H_{\pi\Delta\Delta}(0)=m_\Delta h_1(0),
$$
and the form factors \(g_3,h_3,\tilde g,\tilde h\) are tested against pion-pole-dominance expectations. The reported outcome is qualitative and often quantitative support for PCAC and pion-pole dominance in the \(\Delta\) sector, although subdominant form factors remain statistically challenging [1304.4614].

The excited pseudoscalar states that enter nucleon-dominance analyses also admit hadronic dynamical realizations. A unitarized scalar-pseudoscalar scattering study generates resonances associated with \(K(1460)\), \(\pi(1300)\), \(\pi(1800)\), \(\eta(1475)\), and \(X(1835)\) from rescattering involving the dynamically generated \(f_0(980)\) and \(a_0(980)\). This suggests a concrete hadronic mechanism for the intermediate pseudoscalar spectrum that appears as \(\pi',\pi''\), or specifically \(\pi(1300)\), in dispersive saturation schemes [1011.5843].

## 6. Boundaries of the concept and relation to other dominance models

The term is not interchangeable with every departure from vector meson dominance. In an NJL-type analysis with gauge-covariant diagonalization of axial-vector–pseudoscalar mixing, the correct treatment generates direct photon–pion–quark vertices and therefore departures from naive VMD. However, the authors explicitly reject a simple “pseudoscalar-dominance” replacement of VMD: some amplitudes acquire essential non-VMD contributions, but the pion electromagnetic form factor remains unchanged at the order considered, and anomalous decays such as \(f_1(1285)\to\gamma\pi^+\pi^-\) and \(a_1(1260)\to\gamma\pi^+\pi^-\) retain a surface-term ambiguity that must be fixed experimentally [1807.07742].

Electromagnetic form factors of pseudoscalar mesons also illustrate the limited domain of hadronic pole dominance. A full-lattice study of a light pseudoscalar meson finds that low-\(Q^2\) behavior is well approximated by a vector-meson pole, but the form factor peels away from that pole for \(Q^2\gtrsim 1~\mathrm{GeV}^2\), and \(Q^2F(Q^2)\) is nearly flat between \(Q^2\approx 2\) and \(6~\mathrm{GeV}^2\), with slope \(-0.014(8)\) at \(Q^2=6~\mathrm{GeV}^2\), indicating crossover toward partonic scaling. By contrast, matched simulations with \(N_c=3,4,5\) show that the spacelike pseudoscalar meson vector form factor is nearly independent of \(N_c\) in shape and is consistent with VMD across those color numbers. These results delimit dominance ideas: they work well in specific low-energy domains, but not as universal descriptions across all momentum scales [1701.04250] [2412.14143].

The same caution appears in proposed experimental extractions of the nucleon pseudoscalar form factor. A study of charged-current quasielastic neutrino scattering argues that polarization observables could provide the first extraction of \(F_P\) from neutrino data without imposing PCAC or pion-pole dominance, especially with low-energy \(\nu_\mu,\bar\nu_\mu\) beams or with \(\nu_\tau,\bar\nu_\tau\). In a broader conceptual sense, reviews of emergent hadron mass treat pseudoscalar mesons as the most revealing laboratory for QCD mass generation, with the identity
$$
f_P^0\,\mathbb{E}_P^0(q;P=0)=\frac{M_0(q^2)}{Z_0(q^2)}
$$
expressing the connection between one-body quark dressing and the two-body Goldstone bound state. This usage elevates the pseudoscalar sector because of its symmetry content, not because a single universal dominance ansatz applies to every observable [2008.03527] [2403.00629].

Source: https://www.emergentmind.com/topics/pseudoscalar-dominance