---
title: Pion–Nucleon Sigma Term
url: https://www.emergentmind.com/topics/pion-nucleon-sigma-term
type: topic
---

# Pion–Nucleon Sigma Term

The pion–nucleon sigma term, \(\sigma_{\pi N}\), is a low-energy QCD observable defined by the scalar light-quark matrix element in the nucleon and, equivalently, by the Feynman–Hellmann derivative of the nucleon mass with respect to the average light-quark mass. It quantifies the contribution of explicit chiral-symmetry breaking to the nucleon mass, coincides with the scalar form factor at \(t=0\), and sits at the intersection of \(\pi N\) scattering, pionic atoms, chiral effective theory, lattice QCD, and scalar nucleon couplings relevant for dark-matter phenomenology [1301.3067][1602.07688].

## 1. Operator definition and physical interpretation

In the isospin-symmetric limit, the sigma term is conventionally written as
\[
\sigma_{\pi N} = \hat m\,\langle N|\bar u u + \bar d d|N\rangle,
\qquad
\hat m = \frac{m_u+m_d}{2},
\]
or, in equivalent notation,
\[
\sigma_{\pi N} = m_q\,\langle N|\bar u u + \bar d d|N\rangle.
\]
By the Feynman–Hellmann theorem,
\[
\sigma_{\pi N} = \hat m\,\frac{\partial m_N}{\partial \hat m}
              = M_\pi^2\,\frac{\partial m_N}{\partial M_\pi^2},
\]
so \(\sigma_{\pi N}\) measures the response of the nucleon mass to a variation of the light-quark masses [1301.3067][1710.07164].

This definition admits two closely related interpretations. First, \(\sigma_{\pi N}\) is the nucleon’s isoscalar scalar form factor at zero momentum transfer, \(\sigma(t=0)\). Second, it is the part of the nucleon mass generated by explicit, rather than spontaneous, chiral-symmetry breaking. The quantity therefore fixes how strongly Higgs-generated light-quark masses feed into the nucleon mass, while the remainder is dominated by gluonic dynamics and dynamical chiral symmetry breaking [1602.07688][2508.11435].

The sigma term is also tied to the scalar flavor structure of the nucleon. In phenomenological discussions it is often combined with the strange sigma term, \(\sigma_s = m_s\langle N|\bar s s|N\rangle\), or with the strangeness ratio
\[
y = \frac{2\langle N|\bar s s|N\rangle}{\langle N|\bar u u+\bar d d|N\rangle},
\]
because the size of \(\sigma_{\pi N}\) constrains how much strange scalar density can be accommodated in the nucleon [1301.3067][1911.00846].

## 2. Relation to \(\pi N\) scattering and the Cheng–Dashen theorem

Phenomenological determinations of \(\sigma_{\pi N}\) proceed through the isoscalar \(\pi N\) amplitude. The key low-energy theorem is the Cheng–Dashen relation, which connects the Born-subtracted isoscalar amplitude at the unphysical Cheng–Dashen point, \((\nu,t)=(0,2M_\pi^2)\), to the scalar form factor. In the dispersive formulation used in modern analyses,
\[
\sigma_{\pi N}
=
\Sigma_d + \Delta_D - \Delta_\sigma - \Delta_R,
\qquad
\Sigma_d = F_\pi^2\left(d_{00}^+ + 2M_\pi^2 d_{01}^+\right),
\]
where \(d_{00}^+\) and \(d_{01}^+\) are subthreshold coefficients of the Born-subtracted isoscalar amplitude [1506.04142].

The practical difficulty is that the Cheng–Dashen point lies outside the physical \(\pi N\) region. Roy–Steiner equations solve this problem by combining analyticity, unitarity, and crossing symmetry into a controlled extrapolation from physical data to the subthreshold region. In that framework, \(\sigma_{\pi N}\) becomes tightly correlated with the \(S\)-wave \(\pi N\) scattering lengths. One representative Roy–Steiner relation is
\[
\sigma_{\pi N} = (59.1\pm 3.1)\,\text{MeV}
+ \sum_{I_s} c_{I_s}\big(a^{I_s}-\bar a^{I_s}\big),
\]
with
\[
c_{1/2} = 0.242\,\text{MeV}\times 10^3,\qquad
c_{3/2} = 0.874\,\text{MeV}\times 10^3,
\]
and reference scattering lengths
\[
\bar a^{1/2}=(169.8\pm 2.0)\times 10^{-3}M_\pi^{-1},\qquad
\bar a^{3/2}=(-86.3\pm 1.8)\times 10^{-3}M_\pi^{-1}.
\]
Using pionic-atom input, this program yields
\[
\sigma_{\pi N}=(59.1\pm 3.5)\,\text{MeV}
\]
[1602.07688][1506.04142].

The same dispersive machinery can be anchored directly to low-energy \(\pi N\) scattering rather than pionic atoms. A Roy–Steiner-based fit to low-energy \(\pi N\) cross sections extracted scattering lengths consistent with the pionic-atom determination and obtained
\[
\sigma_{\pi N}=58(5)\,\text{MeV},
\]
thereby reinforcing the “large-\(\sigma\)” phenomenology from an independent data set [1706.01465].

## 3. Pionic atoms, in-medium amplitudes, and nuclear observables

Pionic atoms provide a distinct route to \(\sigma_{\pi N}\) because the strong interaction shifts and broadens Coulombic \(\pi^-\) atomic levels in a way that depends on the \(\pi N\) threshold amplitudes. In the Ericson–Ericson optical potential, the crucial quantity is the in-medium isovector \(s\)-wave amplitude \(b_1(\rho)\). Using the Tomozawa–Weinberg relation and a finite-density Gell-Mann–Oakes–Renner relation, Friedman and Gal wrote
\[
b_1(\rho)
=
b_1^{\rm free}
\left(
1-\frac{\sigma_{\pi N}}{m_\pi^2 f_\pi^2}\rho
\right)^{-1},
\]
so that a decrease of \(f_\pi^2(\rho)\) and of the quark condensate with density enhances \(|b_1(\rho)|\) [1901.03130].

A global fit to 98 pionic-atom level shifts and widths across the periodic table, including deeply bound states in Sn isotopes and \(^{205}\)Pb, gave
\[
\sigma_{\pi N}=57\pm 7~\text{MeV}.
\]
That extraction was reported to be robust against variations of neutron-density parameters, \(p\)-wave modeling, and other optical-potential terms. Higher-order corrections to the leading-density relation were found to involve partial cancellations, implying only a few percent overall systematic uncertainty [1901.03130].

The pionic-atom method probes in-medium chiral dynamics rather than only free \(\pi N\) scattering. In that setting, \(\sigma_{\pi N}\) governs the leading density dependence of
\[
\frac{\langle\bar q q\rangle_\rho}{\langle\bar q q\rangle}
\simeq
1-\frac{\sigma_{\pi N}}{m_\pi^2 f_\pi^2}\rho,
\]
and therefore the partial restoration of chiral symmetry in nuclei. Because pionic atoms mainly sample densities around \(0.5\)–\(0.6\,\rho_0\), the linear-density approximation is comparatively stable in practice [1901.03130].

Deeply bound pionic atoms sharpen the same idea at the level of selected observables. A later analysis of Sn isotopes found that the \(1s\)–\(2p\) binding-energy gap and the \(1s\) width are particularly sensitive to \(\sigma_{\pi N}\). For realistic present-day experimental errors, the gap \(B_\pi(1s)-B_\pi(2p)\) was estimated to constrain \(\sigma_{\pi N}\) at the level of roughly \(2.5\) MeV, provided neutron densities and optical-potential parameters are sufficiently controlled. The same study also emphasized a strong correlation between \(\sigma_{\pi N}\), the real part of \(B_0\), and neutron-density profiles, so the method is intrinsically coupled to nuclear-structure systematics [2204.09211].

## 4. Effective-field-theory extractions and lattice-QCD strategies

Covariant baryon chiral perturbation theory (BChPT) in the EOMS scheme provides a direct link between \(\pi N\) amplitudes and the sigma term through the low-energy constant \(c_1\). At \(\mathcal{O}(p^3)\), the nucleon mass behaves as
\[
M_N = M_N^{(0)} -4c_1 M_\pi^2 + \dots,
\]
so that
\[
\sigma_{\pi N}\approx -4c_1 M_\pi^2
\]
at leading order, with a calculable loop correction at the same chiral order. Fits of the covariant EOMS \(\pi N\) amplitude with explicit \(\Delta(1232)\) to modern partial-wave analyses gave
\[
\sigma_{\pi N}=59(7)\,\text{MeV},
\]
while the older Karlsruhe analysis yielded a smaller value near \(43\) MeV. The higher value was associated with modern meson-factory and pionic-atom input [1301.3067][1111.4934].

Lattice QCD accesses \(\sigma_{\pi N}\) in two main ways. The direct method computes the scalar three-point function and extracts
\[
g_S^{u+d}=\langle N|\bar u u+\bar d d|N\rangle,
\qquad
\sigma_{\pi N}=m_{ud}\,g_S^{u+d},
\]
whereas the Feynman–Hellmann method fits the nucleon mass as a function of \(M_\pi\) and differentiates with respect to \(M_\pi^2\) [2203.13862].

These strategies have not always agreed numerically. A covariant SU(2) EOMS analysis of ETMC nucleon masses obtained
\[
\sigma_{\pi N}=50.2(1.2)(2.0)\,\text{MeV},
\]
and found that including a virtual \(\Delta(1232)\) does not change the result qualitatively [1710.07164]. By contrast, a heavy-baryon SU(3) \(\mathcal{O}(p^3)\) analysis that fit \(\pi N\) phase shifts together with octet-baryon masses reported
\[
\sigma_{\pi N}=(34.57\pm 11.85)\,\text{MeV},
\]
together with a very small nucleon strangeness content, \(y_N=0.05\pm0.05\) [1911.00846].

A more recent two-loop EOMS treatment of the nucleon mass, applied to \(N_f=2+1\) lattice data, reported
\[
\sigma_{\pi N}=55.9(2.5)\,\text{MeV}.
\]
That analysis identified intermediate \(\pi\pi\) rescattering effects, which enter only at two-loop order, as the mechanism by which the long-standing lattice–dispersive tension can be naturally resolved [2508.11435].

## 5. Tension, systematics, and disputed determinations

For much of the last decade, the central controversy concerned the gap between phenomenological extractions near \(59\) MeV and several direct lattice calculations near the physical point. A Roy–Steiner review collected representative lattice results such as \(38(3)(3)\) MeV from BMW, \(44.4(3.2)(4.5)\) MeV from \(\chi\)QCD, \(37.22(2.57)^{+0.99}_{-0.63}\) MeV from ETMC, and \(35(6)\) MeV from RQCD, and characterized the mismatch with \((59.1\pm 3.5)\) MeV as a tension of roughly \(3\sigma\) or more [1602.07688].

A major proposed explanation is excited-state contamination in lattice nucleon correlators. A ChPT-guided reanalysis argued that direct lattice calculations of \(\sigma_{\pi N}\) are contaminated by multihadron \(N\pi\) and \(N\pi\pi\) states, whose contributions are large and negative at realistic source–sink separations. In that picture, standard fits return values near \(40\) MeV, but fits that incorporate the multihadron spectrum move the result toward \(\sim 60\) MeV [2203.13862]. A subsequent study summarized the same conclusion more explicitly: \(N\pi\) and \(N\pi\pi\) states each contribute about \(10\) MeV to \(\sigma_{\pi N}\), and including the \(\Delta\) as an explicit degree of freedom does not alter that conclusion [2301.07885].

A second systematic issue is the convention used for the isospin limit. A recent analysis pointed out that phenomenology conventionally defines the isospin limit with the charged pion mass, while lattice QCD usually adopts the neutral pion mass. In ChPT this mismatch induces
\[
\Delta \sigma_{\pi N}=3.1(5)\,\text{MeV},
\]
which should be included when comparing lattice and phenomenological results. The same work updated the Roy–Steiner-plus-pionic-atom benchmark, using the latest pionic-hydrogen width, to
\[
\sigma_{\pi N}=59.0(3.5)\,\text{MeV}
\]
[2305.07045].

A separate controversy concerns the low-energy CHAOS data set. One analysis of CHAOS \(\pi^\pm p\) data extracted
\[
\sigma_{\pi N}=(44\pm 12)\,\text{MeV},
\]
placing the sigma term near the lower edge of the historical range [1211.1148]. A later comment argued that the angular distribution of the CHAOS \(\pi^+p\) differential cross sections is incompatible in shape with the rest of the modern low-energy \(\pi^+p\) data and that this problem must be resolved before any extrapolation into the unphysical Cheng–Dashen region can be trusted [1309.3469].

The range of representative determinations illustrates both the methodological diversity and the historical spread:

| Approach | Result | Source |
|---|---:|---|
| Roy–Steiner + pionic atoms | \(59.1\pm 3.5\) MeV | [1506.04142] |
| Low-energy \(\pi N\) scattering + Roy–Steiner | \(58(5)\) MeV | [1706.01465] |
| Covariant EOMS BChPT fit to \(\pi N\) PWAs | \(59(7)\) MeV | [1301.3067] |
| Global pionic-atom fit | \(57\pm 7\) MeV | [1901.03130] |
| Covariant SU(2) EOMS fit to lattice \(m_N\) | \(50.2(1.2)(2.0)\) MeV | [1710.07164] |
| HB SU(3) fit to \(\pi N\) phase shifts and octet masses | \(34.57\pm 11.85\) MeV | [1911.00846] |
| Two-loop EOMS extrapolation of \(N_f=2+1\) lattice data | \(55.9(2.5)\) MeV | [2508.11435] |

## 6. Phenomenological consequences and current outlook

Because \(\sigma_{\pi N}\) controls the scalar light-quark content of the nucleon, it propagates directly into scalar nucleon couplings. In the Roy–Steiner framework, a sigma term near \(59\) MeV leads to
\[
f_u^p = (20.8 \pm 1.5)\times 10^{-3},\quad
f_d^p = (41.1 \pm 2.8)\times 10^{-3},
\]
\[
f_u^n = (18.9 \pm 1.4)\times 10^{-3},\quad
f_d^n = (45.1 \pm 2.7)\times 10^{-3},
\]
with immediate implications for Higgs-mediated and other scalar interactions in dark-matter direct detection [1506.04142]. More generally, both phenomenological and lattice-oriented discussions identify \(\sigma_{\pi N}\) as a dominant hadronic input to spin-independent WIMP–nucleon scattering and to scalar contributions in \(\mu\to e\) conversion and EDM analyses [1602.07688][2203.13862].

In nuclear matter, \(\sigma_{\pi N}\) fixes the leading density dependence of the quark condensate and thus the strength of partial chiral restoration. Pionic-atom analyses near \(\sigma_{\pi N}\approx 57\) MeV imply a noticeable decrease of \(\langle\bar q q\rangle_\rho\) and \(f_\pi^2(\rho)\) already at densities around \(0.5\,\rho_0\), which feeds directly into the in-medium renormalization of the isovector \(\pi N\) amplitude \(b_1(\rho)\) [1901.03130]. Continuum DSE studies that infer \(\sigma_{\pi N}\) from the density dependence of the chiral condensate likewise find values near \(62\) MeV, reinforcing the link between the sigma term and in-medium chiral dynamics [1910.08298].

The present status is more convergent than the older “small-\(\sigma\)” versus “large-\(\sigma\)” dichotomy suggests. Phenomenological analyses based on Roy–Steiner equations, low-energy \(\pi N\) scattering, pionic atoms, and covariant BChPT continue to cluster near \(58\)–\(60\) MeV [1706.01465][2305.07045]. Several recent developments on the lattice side—explicit treatment of multihadron excited states, correction of isospin-limit conventions, and two-loop chiral extrapolations with \(\pi\pi\) rescattering—move the preferred range upward and reduce the discrepancy [2203.13862][2305.07045][2508.11435].

The sigma term therefore remains a benchmark quantity for the mutual consistency of hadronic phenomenology, effective field theory, and lattice QCD. Its importance lies not only in the number itself, but in the fact that \(\sigma_{\pi N}\) encodes, in a single scalar matrix element, how chiral symmetry breaking propagates from the QCD Lagrangian into nucleon mass, \(\pi N\) amplitudes, nuclear observables, and scalar couplings beyond the Standard Model.

Source: https://www.emergentmind.com/topics/pion-nucleon-sigma-term