In 1990 Guagnelli, Marinari and Parisi proposed extracting two-hadron scattering lengths from configuration-by-configuration fluctuations of hadron masses in finite-volume lattice QCD, converting the fluctuation-induced energy shifts with Lüscher's relation, in quenched simulations with staggered and Wilson quarks. The paper itself flagged four limitations: the quenched approximation omits quark-exchange/annihilation ('eye') diagrams and sea-quark effects; the pion-nucleon channel was too noisy for its scattering length to yield a controlled determination of the nucleon sigma term; Wilson-fermion results suffer discretization artefacts from fluctuations of the critical hopping parameter; and the isospin decomposition remained within a four-different-flavour 'toy' world, with chiral extrapolation to the physical point deferred to future unquenched calculations. The meson-meson half of the program is closed: full four-point functions including eye diagrams were computed by Kuramashi et al. (1993) and Fukugita et al. (1995); Bernard and Golterman (1996) showed that the quenched theory modifies Lüscher's relation itself with uncontrolled enhanced corrections; and unquenched, physical-pion-mass determinations of ππ scattering in both isospin channels now exist (Blum et al. 2021; Bruno et al. 2023). The meson-nucleon half remains open to a controlled determination: S-wave πN scattering lengths in both isospin channels have been extracted at unphysically heavy pion mass on a single ensemble (Bulava et al. 2023); both channels have since been reached at the physical pion mass, but only as a by-product of a two-ensemble nucleon-polarizability calculation, without a dedicated analysis of sea-quark, chiral, continuum and finite-volume systematics (Wang et al. 2023); and only the I = 3/2 channel has been reached in a dedicated scattering analysis, on a single ensemble at one lattice spacing and one volume (Alexandrou et al. 2024). A dedicated determination of both lengths at the physical point, with sea-quark, chiral, continuum and finite-volume systematics quantified, therefore remains the target. The closely related nucleon sigma term stands in a live tension between standard direct lattice analyses (about 38-45 MeV) and dispersive Roy-Steiner / two-loop chiral analyses (about 55-60 MeV); that tension motivates the πN amplitudes sought here but is not itself part of the target.

Background

Lattice QCD is the ab initio, non-perturbative regularization of quantum chromodynamics: Euclidean spacetime is replaced by a finite hypercubic grid of spacing aa and spatial extent LL, quark fields live on sites, gauge fields on links, and hadron properties are extracted from Monte Carlo estimates of Euclidean correlation functions over gauge-field configurations ('ensembles'). In a periodic box, a two-hadron state of zero total momentum has an energy EE shifted from the free value m1+m2m_1+m_2; Lüscher proved that, below inelastic threshold and up to exponentially small corrections, this shift determines the elastic S-wave scattering length a0a_0:

E−(m1+m2)=−2π(m1+m2)a0m1m2L3[1+c1a0L+c2(a0L)2]+O(L−6)E-(m_1+m_2) = -\frac{2\pi(m_1+m_2)a_0}{m_1 m_2 L^3}\left[1 + c_1\frac{a_0}{L} + c_2\left(\frac{a_0}{L}\right)^2\right] + O(L^{-6})

with c1=−2.837297c_1 = -2.837297 and c2=6.375183c_2 = 6.375183. The scattering length is the leading coefficient of the effective-range expansion of the S-wave phase shift δ(p)\delta(p), pcot⁡δ(p)=−1/a0+(r0/2)p2+…p\cot\delta(p) = -1/a_0 + (r_0/2)p^2 + \ldots, and sets the strength of low-energy interactions.

The source paper introduced a fluctuation method: on each gauge configuration the effective masses of one- and two-hadron correlators fluctuate around their ensemble means; the second moment of the per-configuration mass distribution encodes the volume-dependent energy shift, which Lüscher's relation converts into a scattering length. The simulations were quenched — the fermion determinant is omitted, so the gauge field does not respond to virtual quark-antiquark pairs — which suppresses all sea-quark effects and leaves only gluon-exchange ('valence') diagrams: quark-annihilation 'eye' diagrams are absent. Two fermion discretizations were used: staggered (Kogut-Susskind) quarks, which retain a remnant chiral symmetry, and Wilson quarks, which break chiral symmetry explicitly so that the massless limit is reached only at a critical value kck_c of the hopping parameter; per-configuration fluctuations of kck_c then act as a fluctuating mass counterterm and can fake chiral behaviour. To avoid annihilation diagrams entirely, the meson-meson channel was built from four different quark flavours (a 'toy' channel). The paper reported a small scattering length (of order 0.02 fm, in physical units, for staggered quarks at β=6.0\beta = 6.0) in the four-flavour channel, a tentative pion-nucleon scattering length of order 0.7 fm in the quark-mass range mπ2/mρ2∈[0.5,0.8]m_\pi^2/m_\rho^2 \in [0.5, 0.8], and only a qualitative indication for the nucleon sigma term.

The nucleon sigma term, σπN=mud ⟨N∣uˉu+dˉd∣N⟩\sigma_{\pi N} = m_{ud}\,\langle N|\bar{u}u + \bar{d}d|N\rangle, measures the light-quark scalar content of the nucleon; by the Feynman-Hellmann theorem it equals the response of the nucleon mass to the light-quark mass, and by the Cheng-Dashen low-energy theorem it is related to the pion-nucleon scattering amplitude evaluated at the subthreshold (unphysical) point. Pion-nucleon scattering and the sigma term are therefore two faces of one problem. σπN\sigma_{\pi N} is a Standard-Model input to the nucleon mass decomposition, to baryon chiral perturbation theory (ChPT, the systematic low-energy effective field theory of pions and nucleons), and — through the nucleon scalar couplings — to dark-matter direct-detection cross sections.

Modern scattering calculations use unquenched ensembles with Nf=2+1N_f = 2+1 or 2+1+12+1+1 sea-quark flavours, improved, domain-wall or twisted-mass discretizations (alternative fermion formulations with better chiral properties, evaluated at several lattice spacings so the continuum limit a→0a \to 0 can be taken), variational (generalized eigenvalue) analyses of two-hadron spectra with large operator bases, and all-to-all propagator techniques for quark-disconnected diagrams. Baryon correlation functions suffer exponential signal-to-noise degradation, making near-threshold NπN\pi levels statistically expensive; the I=1/2I = 1/2 pion-nucleon channel requires quark-disconnected contributions; and the Δ(1232)\Delta(1232) resonance — the lowest NπN\pi resonance, of spin-parity 3/2+3/2^+ and isospin 3/23/2, lying just above threshold — strongly mixes with the S-wave I=3/2I = 3/2 channel. Nucleon two- and three-point functions are contaminated by NπN\pi and NππN\pi\pi excited states — precisely the physics now understood to drive the sigma-term tension.

The source paper closed by acknowledging these limitations: the omission of exchange diagrams (and sea quarks), the insufficient quality of the pion-nucleon signal for a quantitative sigma-term extraction, the Wilson kck_c artefacts, and the restriction to a four-flavour toy world with chiral extrapolation left to future unquenched calculations.

Problem Statement

The open problem: a controlled determination, from unquenched (Nf≥2+1N_f \geq 2+1) lattice QCD, of the S-wave pion-nucleon scattering lengths in both physical isospin channels — I = 1/2, including its quark-disconnected contributions, and I = 3/2, under the nearby Δ(1232) resonance — at the physical pion mass (reached by chiral extrapolation where the simulated pion masses are heavier), with sea-quark, chiral, continuum and finite-volume systematics quantified. The question concerns the physical observables of the source program — the S-wave πN scattering lengths in the physical isospin channels — rather than the particular per-configuration fluctuation technique used in 1990: the group that continued the program pursued the same observables with standard four-point-function methods and identified the blocking deficiency of the 1990 calculation as physical (missing exchange/annihilation diagrams), while Bernard and Golterman (1996) showed that the quenched framework cannot in principle yield controlled scattering lengths, which is why the unquenched condition is required. The ππ part of the source program — annihilation and sea-quark effects, isospin decomposition, and physical-point chiral extrapolation for meson-meson channels — is resolved and is not part of the open target.

The nucleon sigma term σπN\sigma_{\pi N}, its subthreshold-amplitude relation to πN scattering (Cheng-Dashen low-energy theorem), and the present discrepancy between standard direct lattice determinations (about 38-45 MeV) and dispersive Roy-Steiner or two-loop chiral analyses (about 55-60 MeV) are retained here as motivation and context for the πN amplitudes sought; they are not part of the target.

References

  1. M. Guagnelli, E. Marinari, G. Parisi, 'Scattering lengths from fluctuations', Phys. Lett. B 240 (1990) 188. DOI: 10.1016/0370-2693(90)90431-5, https://doi.org/10.1016/0370-2693(90)90431-5
  2. M. Lüscher, 'Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories. 2. Scattering States', Commun. Math. Phys. 105 (1986) 153. DOI: 10.1007/BF01211097, https://doi.org/10.1007/BF01211097
  3. Y. Kuramashi, M. Fukugita, H. Mino, M. Okawa, A. Ukawa, 'Lattice QCD calculation of full pion scattering lengths', Phys. Rev. Lett. 71 (1993) 2387. DOI: 10.1103/PhysRevLett.71.2387, https://doi.org/10.1103/PhysRevLett.71.2387
  4. C. W. Bernard, M. F. L. Golterman, 'Finite-volume two-pion energies and scattering in the quenched approximation', Phys. Rev. D 53 (1996) 476. DOI: 10.1103/PhysRevD.53.476, arXiv:hep-lat/9507004, https://arxiv.org/abs/hep-lat/9507004
  5. M. Fukugita, Y. Kuramashi, M. Okawa, H. Mino, A. Ukawa, 'Hadron scattering lengths in lattice QCD', Phys. Rev. D 52 (1995) 3003. DOI: 10.1103/PhysRevD.52.3003, arXiv:hep-lat/9501024, https://arxiv.org/abs/hep-lat/9501024
  6. M. Fukugita, Y. Kuramashi, M. Okawa, A. Ukawa, 'Pion-nucleon sigma term in lattice QCD', Phys. Rev. D 51 (1995) 5319. DOI: 10.1103/PhysRevD.51.5319, arXiv:hep-lat/9408002, https://arxiv.org/abs/hep-lat/9408002
  7. S. R. Beane, P. F. Bedaque, K. Orginos, M. J. Savage, 'I = 2 ππ scattering from fully-dynamical mixed-action lattice QCD', Phys. Rev. D 73 (2006) 054503. DOI: 10.1103/PhysRevD.73.054503, arXiv:hep-lat/0506013, https://arxiv.org/abs/hep-lat/0506013
  8. X. Feng, K. Jansen, D. B. Renner, 'The π+ π+ scattering length from maximally twisted mass lattice QCD', Phys. Lett. B 684 (2010) 268. DOI: 10.1016/j.physletb.2010.01.018, arXiv:0909.3255, https://arxiv.org/abs/0909.3255
  9. T. Blum, P. A. Boyle, M. Bruno, N. H. Christ, D. Hoying, C. Kelly, C. Lehner, R. D. Mawhinney, A. S. Meyer, D. J. Murphy, C. T. Sachrajda, A. Soni, T. Wang, 'Lattice determination of I=0 and 2 ππ scattering phase shifts with a physical pion mass', Phys. Rev. D 104 (2021) 114506. DOI: 10.1103/PhysRevD.104.114506, arXiv:2103.15131, https://arxiv.org/abs/2103.15131
  10. M. Bruno, D. Hoying, T. Izubuchi, C. Lehner, A. S. Meyer, M. Tomii, 'Isospin 0 and 2 two-pion scattering at physical pion mass using distillation with periodic boundary conditions in lattice QCD', preprint (2023). arXiv:2304.03313, https://arxiv.org/abs/2304.03313
  11. A. Torok, S. R. Beane, W. Detmold, T. C. Luu, K. Orginos, A. Parreño, M. J. Savage, A. Walker-Loud, 'Meson-Baryon Scattering Lengths from Mixed-Action Lattice QCD', Phys. Rev. D 81 (2010) 074506. DOI: 10.1103/PhysRevD.81.074506, arXiv:0907.1913, https://arxiv.org/abs/0907.1913
  12. J. Bulava, 'Meson-Nucleon Scattering Amplitudes from Lattice QCD', AIP Conf. Proc. 2249 (2020) 020006. DOI: 10.1063/5.0008643, arXiv:1909.13097, https://arxiv.org/abs/1909.13097
  13. J. Bulava, A. D. Hanlon, B. Hörz, C. Morningstar, A. Nicholson, F. Romero-López, S. Skinner, P. Vranas, A. Walker-Loud, 'Elastic nucleon-pion scattering at mπ = 200 MeV from lattice QCD', Nucl. Phys. B 987 (2023) 116105. DOI: 10.1016/j.nuclphysb.2023.116105, arXiv:2208.03867, https://arxiv.org/abs/2208.03867
  14. R. Gupta, S. Park, M. Hoferichter, E. Mereghetti, B. Yoon, T. Bhattacharya, 'Pion-Nucleon Sigma Term from Lattice QCD', Phys. Rev. Lett. 127 (2021) 242002. DOI: 10.1103/PhysRevLett.127.242002, arXiv:2105.12095, https://arxiv.org/abs/2105.12095
  15. M. Hoferichter, J. Ruiz de Elvira, B. Kubis, U.-G. Meißner, 'Roy-Steiner-equation analysis of pion-nucleon scattering', Phys. Rep. 625 (2016) 1. DOI: 10.1016/j.physrep.2016.02.002, arXiv:1510.06039, https://arxiv.org/abs/1510.06039
  16. C. Iona, C. Alexandrou, S. Bacchio, J. Finkenrath, G. Koutsou, Y. Li, G. Spanoudes, 'Nucleon axial, tensor, and scalar charges and sigma-terms from lattice QCD', PoS LATTICE2024 (2025) 316. DOI: 10.22323/1.466.0316, arXiv:2502.05274, https://arxiv.org/abs/2502.05274
  17. Z.-R. Liang, H.-X. Chen, F.-K. Guo, Z.-H. Guo, D.-L. Yao, 'Two-Loop Extraction of the Pion-Nucleon Sigma Term', preprint (2025). arXiv:2508.11435, https://arxiv.org/abs/2508.11435
  18. C. Alexandrou, S. Bacchio, G. Koutsou, T. Leontiou, S. Paul, M. Petschlies, F. Pittler, 'Elastic Nucleon-Pion scattering amplitudes in the Δ channel at physical pion mass from Lattice QCD', Phys. Rev. D 109 (2024) 034509. DOI: 10.1103/PhysRevD.109.034509, arXiv:2307.12846, https://arxiv.org/abs/2307.12846
  19. X.-H. Wang, C.-L. Fan, X. Feng, L.-C. Jin, Z.-L. Zhang, 'Nucleon electric polarizabilities and nucleon-pion scattering at physical pion mass', preprint (2023). arXiv:2310.01168, https://arxiv.org/abs/2310.01168