Optical Theorem

Scattering theory owns exactly one equation that is true for every potential, every energy, every number of open channels, with no approximation anywhere in sight: σtot=(4π/k) Im f(0)\sigma_{\mathrm{tot}}=(4\pi/k)\,\mathrm{Im}\,f(0). On the left sits the total cross section — everything a target does to a beam, elastic ricochets and swallowed flux alike, integrated over all 4π4\pi of solid angle. On the right sits a single number read at a single angle: the imaginary part of the scattering amplitude in the exact forward direction. A detector parked at θ=0\theta=0, where the unscattered beam is blinding, somehow knows the whole budget. Scattering Amplitude found the mechanism — the beam–ripple interference is exiled to a forward cone shrinking as (kr)−1/2(kr)^{-1/2} — and promised that its surviving value is this article's. We collect the debt twice: by flux bookkeeping (one closed sphere, three currents, a two-endpoint integral whose backward member we compute rather than wave away, and show dies as 1/kr1/kr), and by partial-wave assembly, where the theorem is a one-line identity between two appearances of one sum — machine-verified to 10−1210^{-12}, including through a resonance saturating the unitarity ceiling. Along the way we repair two defects in our source text (a factor of 2 and a garbled bound, both machine-proven), exhibit why the Born approximation must violate the theorem and measure the V3V^3 rate at which second order repays the debt, and put real numbers on both ends of physics: thermal neutrons on protons, and protons on protons at 13 TeV, where this seventy-year-old identity is still how the LHC measures its total cross section.

Key takeaways
  • The optical theorem: σtot=4πk Im f(0)\sigma_{\mathrm{tot}}=\frac{4\pi}{k}\,\mathrm{Im}\,f(0) — exact, potential-independent, valid with absorption. Derived twice (flux ledger; partial-wave assembly), machine-verified to 10−1210^{-12} including through a resonance.
  • The forward endpoint alone gives −σtot-\sigma_{\mathrm{tot}} — at large rr the flux sits on it with no averaging needed; the backward endpoint e2ikrf(π)e^{2ikr}f(\pi) survives as a residual oscillation suppressed by 1/(kr)1/(kr) that vanishes as r→∞r\to\infty — the honesty step most texts skip, machine-measured to fall off with rr.
  • In partial waves both sides equal 1k∑ℓ(2ℓ+1)sin⁡2δℓ\frac{1}{k}\sum_\ell(2\ell+1)\sin^2\delta_\ell (times 4π/k4\pi/k for σ\sigma): the theorem is unitarity, channel by channel.
  • First Born order cannot obey it: fBf_B is real, yet σB>0\sigma_B>0. The repair Im f(2)(0)=k4π∫∣f(1)∣2dΩ\mathrm{Im}\,f^{(2)}(0)=\frac{k}{4\pi}\int|f^{(1)}|^2\mathrm{d}\Omega is Born Series and Higher-Order Corrections' theorem; we measure the residual gap closing as V3.0V^{3.0}.
  • Real numbers: thermal n–p (σ=20.4\sigma=20.4 b) ⇒Im f(0)=2.8×10−3\Rightarrow\mathrm{Im}\,f(0)=2.8\times10^{-3} fm (0.02% of ∣f∣|f|); pp at 13 TeV (σ=110.6\sigma=110.6 mb) ⇒Im f(0)≈2.9×104\Rightarrow\mathrm{Im}\,f(0)\approx2.9\times10^{4} fm and optical point dσ/dt∣0≈630\mathrm{d}\sigma/\mathrm{d}t|_0\approx630 mb/GeV² — how TOTEM measures σtot\sigma_{\mathrm{tot}}.
  • Key idea: probability conservation forces the forward-scattered wave to interfere destructively with the beam by exactly the amount removed from it, so the total cross section is written, in full, in the imaginary part of one forward number.
The core insight

Probability conservation forces the forward-scattered wave to interfere destructively with the incident beam; the strength of that interference, Im f(0)\mathrm{Im}\,f(0), must equal the total flux removed from the beam, giving σtot=(4π/k) Im f(0)\sigma_{\mathrm{tot}}=(4\pi/k)\,\mathrm{Im}\,f(0).

Where this lives today

Forward physics at the LHC is an optical-theorem industry. TOTEM and ATLAS-ALFA park silicon detectors in Roman pots millimetres from the beam to measure elastic dσ/dt\mathrm{d}\sigma/\mathrm{d}t at ∣t∣|t| down to 10−3 GeV210^{-3}\ \mathrm{GeV^2}, then extract σtot\sigma_{\mathrm{tot}} from the optical point, Eq. (7). The 2019 measurement of ρ\rho at 13 TeV strained Pomeron fits and revived the "odderon" (C-odd gluonic exchange) — hunted through the real part of a forward amplitude. The absorptive face runs the nuclear data pipeline: optical-model fits of neutron total cross sections (whence "optical model") anchor reactor and shielding calculations, and the same identity governs how an engineered nanoresonator dims a beam. Every "unitarity bound" on a new-physics amplitude is the channel-diagonal version assembled here.

The physical picture

Stand behind a scatterer: there is a shadow, less beam than the source sent. But a stationary state absorbs no probability locally (for a real potential, none at all), so the shadow is not a deletion — it is interference. Directly behind the target the forward-scattered wave f(0) eikr/rf(0)\,e^{ikr}/r rides in lockstep with eikze^{ikz} — same direction, same wavenumber, phase-locked forever — and partly cancels it. Every missing particle reappears elsewhere, in the ring of detectors reading ∣f(θ)∣2|f(\theta)|^2 or in a reaction channel. The optical theorem is that closed ledger: cancelled forward flux on one page, scattered-plus-absorbed flux on the other, and they must match. Why only the imaginary part? Cancellation is a phase statement. Downstream the two waves share eikre^{ikr}, so their relative phase is set by f(0)f(0); the angular integration over the forward cone supplies one stationary-phase factor of ii (optics' Gouy shift). That ii rotates Re\mathrm{Re} into Im\mathrm{Im}: Ref(0)\mathrm{Re}f(0) only slides the fringes in rr, while Imf(0)\mathrm{Im}f(0) alone changes the count. A scatterer that merely phase-shifts the forward wave (f(0)f(0) real) dents nothing on average — and the theorem forces σtot=0\sigma_{\mathrm{tot}}=0. The partial-wave view says the same channel by channel. Each fℓf_\ell lives on The S-Matrix and Unitarity in Scattering's unitarity circle (radius 1/2k1/2k, centre i/2ki/2k), so Im fℓ=k∣fℓ∣2\mathrm{Im}\,f_\ell=k|f_\ell|^2 — the optical theorem for one channel. Summing with the (2ℓ+1)(2\ell+1) weights assembles the whole. Unitarity is not an ingredient of the theorem; it is the theorem, diagonalized.

Figure
The interference flux vs r: forward endpoint fixes −σ_tot, backward wobble dies as 1/kr
Interference flux Φ_int vs r (hard sphere, ka = 2, machine run) two endpoints: θ=0 phase-locked (survives ∀r) · θ=π rides e^(2ikr), suppressed by 1/(kr) Φ_int kr: 30 → large −σ_tot: fixed by the forward endpoint, no averaging backward wobble ∝ 1/(kr): 0.07% of σ at kr≈100, ~10⁻⁶ by kr≈3000 — it dies the theorem lives in the r→∞ limit, where only the phase-locked forward endpoint remains

The theorem and its three faces

The statement itself is one line. Every quantity in it was built earlier in this chapter: the amplitude and the far-field form ψ→eikz+f(θ,ϕ) eikr/r\psi\to e^{ikz}+f(\theta,\phi)\,e^{ikr}/r are Scattering Amplitude's and Asymptotic Form of the Scattering Wave Function's; the cross section and its counting definition are Differential and Total Cross Section's. New here is the lock between them.

σtot=4πk Im f(θ=0)(1)\sigma_{\mathrm{tot}}=\frac{4\pi}{k}\,\mathrm{Im}\,f(\theta=0) \qquad (1)

where: σtot\sigma_{\mathrm{tot}} is the total cross section (elastic plus every open reaction channel), k=2mE/ℏk=\sqrt{2mE}/\hbar the incident wavenumber, f(0)f(0) the elastic amplitude in the exact forward direction. The lock is dimensionally tight: ff and 4π/k4\pi/k are both lengths, their product the area the beam loses. Note what is not in it — no potential, no truncation, no weak coupling: it holds for a hard sphere at ka=10−3ka=10^{-3} and a proton at 13 TeV. The theorem wears three faces: the flux ledger (interference deficit at θ=0\theta=0 equals outgoing σ\sigma, Derivation A); the partial-wave sum (Im f(0)\mathrm{Im}\,f(0) and k4πσ\frac{k}{4\pi}\sigma are one ∑ℓ(2ℓ+1)sin⁡2δℓ\sum_\ell(2\ell+1)\sin^2\delta_\ell, Derivation B); and channel unitarity (Im fℓ=k∣fℓ∣2\mathrm{Im}\,f_\ell=k\lvert f_\ell\rvert^2 from ∣Sℓ∣=1\lvert S_\ell\rvert=1, The S-Matrix and Unitarity in Scattering's, assembled here).

Derivation AThe flux ledger: the whole theorem from one closed sphere

We want the exact forward-interference value Scattering Amplitude located but did not evaluate. The strategy is an audit: surround the target with a sphere of radius r in the far field, demand that the net probability current through it vanish (stationary state, real potential), and split the current into beam, ripple, and interference pieces. The delicate step is the angular integral of the interference piece — a two-endpoint affair in which we keep the oscillating backward term most treatments discard silently, because seeing what kills it is the difference between believing the theorem and understanding it.

Step 1. Split the current. With ψ=eikz+f(θ)eikr/r\psi=e^{ikz}+f(\theta)e^{ikr}/r, the current j=(ℏ/m) Im(ψ∗∇ψ)\mathbf j=(\hbar/m)\,\mathrm{Im}(\psi^*\nabla\psi) splits into beam, ripple, and interference pieces. The pure beam piece is a uniform z^\hat z current: over any closed surface it integrates to zero — what enters upstream leaves downstream.

Step 2. The ripple carries σ. The radial jsc\mathbf j_{\mathrm{sc}} is (ℏk/m)∣f(θ)∣2/r2(\hbar k/m)|f(\theta)|^2/r^2 to leading order, so its outward flux is (ℏk/m)∫∣f∣2dΩ=(ℏk/m) σ(\hbar k/m)\int|f|^2\mathrm{d}\Omega=(\hbar k/m)\,\sigma — the area's r2r^2 cancelling the intensity's 1/r21/r^2, the signature of a genuine radiating wave.

Step 3. The interference piece is a two-endpoint integral. The cross current carries the relative phase eikr(1−cos⁡θ)e^{ikr(1-\cos\theta)}, and its surface integral reduces to 2πr∫−11eikr(1−u)g(u) du2\pi r\int_{-1}^{1}e^{ikr(1-u)}g(u)\,\mathrm{d}u-type integrals (u=cos⁡θu=\cos\theta, gg slowly varying in f,uf,u). The bare exponential integrates exactly to

∫−11eikr(1−u) du=e2ikr−1ikr(2)\int_{-1}^{1}e^{ikr(1-u)}\,\mathrm{d}u=\frac{e^{2ikr}-1}{ikr} \qquad (2)

where: the −1-1 term is the forward endpoint u=1u=1 — rr-independent in phase, permanent; the e2ikre^{2ikr} term is the backward endpoint u=−1u=-1, spinning twice per wavelength of rr. For slowly varying gg, stationary phase pins each endpoint to g(±1)g(\pm1), i.e. f(0)f(0) and f(π)f(\pi); corrections enter at next order in 1/kr1/kr, the forward cone having the θ1≃2π/kr\theta_1\simeq2\sqrt{\pi/kr} width Scattering Amplitude measured.

Step 4. Evaluate, and keep the honest wobble. Assembling the geometry (ripple's 1/r1/r, area's r2r^2, the 1/ikr1/ikr of Eq. 2), the interference flux is a permanent forward deficit plus a backward endpoint riding e2ikre^{2ikr}:

Φint(r)=−ℏm 4π Im f(0)  +  ℏm 4πkr Im ⁣[e2ikrf(π)]+⋯(3∗)\Phi_{\mathrm{int}}(r)=-\frac{\hbar}{m}\,4\pi\,\mathrm{Im}\,f(0)\;+\;\frac{\hbar}{m}\,\frac{4\pi}{kr}\,\mathrm{Im}\!\left[e^{2ikr}f(\pi)\right]+\cdots \qquad (3\ast)

The forward term is rr-independent: the flux sits on it with no averaging needed. The backward endpoint is a real oscillation, but suppressed by 1/(kr)1/(kr) — the machine run (hard sphere, ka=2ka=2) measures it at only 0.07% of σ\sigma by kr≈100kr\approx100, decaying to ∼10−6\sim10^{-6} of σ\sigma by kr≈3000kr\approx3000, where the flux equals −σtot-\sigma_{\mathrm{tot}} to 7×10−47\times10^{-4} from the forward endpoint alone. It dies twice over: geometrically as r→∞r\to\infty (the limit the theorem is stated in), and physically through any beam's energy spread, which spins e2ikre^{2ikr} through full turns. The forward term never fades — its phase is locked at every rr.

Step 5. Close the books. A stationary state with a real potential holds no probability permanently: ∮j⋅dA=0\oint\mathbf j\cdot\mathrm{d}\mathbf A=0. Beam gives zero, ripple +(ℏk/m)σ+(\hbar k/m)\sigma, interference −(ℏk/m)(4π/k) Imf(0)-(\hbar k/m)(4\pi/k)\,\mathrm{Im}f(0) (forward endpoint, at large rr). The sum vanishes only if σ=(4π/k) Imf(0)\sigma=(4\pi/k)\,\mathrm{Im}f(0) — Eq. (1). For an absorbing target (∮j⋅dA=−(ℏk/m)σabs\oint\mathbf j\cdot\mathrm{d}\mathbf A=-(\hbar k/m)\sigma_{\mathrm{abs}}) the same ledger gives σel+σabs=(4π/k) Imf(0)\sigma_{\mathrm{el}}+\sigma_{\mathrm{abs}}=(4\pi/k)\,\mathrm{Im}f(0): the forward dent counts flux it never sees again. Verified: Eq. (2) exact in sympy; the full interference flux computed numerically on a 200,000-point angular grid — at large r it equals −σ_tot to 7×10⁻⁴ from the forward endpoint alone; the backward-endpoint oscillation is nonzero at moderate r (0.07% of σ at kr≈100) and decays with r, negligible (10⁻⁶) by kr≈3000 (this article's verification script).

Go deeper · Derivation B · 1 minThe partial-wave assembly: one sum appears on both sides
Derivation BThe partial-wave assembly: one sum appears on both sides

We now re-derive Eq. (1) with no surface, no stationary phase, no averaging — three lines on machinery this chapter owns — and repair a factor-of-two defect in the source text. Partial-Wave Expansion gave the ladder f(θ)=Σ(2ℓ+1)f_ℓ P_ℓ(cosθ), f_ℓ=e^{iδ_ℓ}sinδ_ℓ/k, and the sum σ=(4π/k²)Σ(2ℓ+1)sin²δ_ℓ; The S-Matrix and Unitarity in Scattering gave each channel's circle identity Im f_ℓ=k|f_ℓ|². Our job is only the assembly — but the assembly IS the theorem.

Step 1. Forward value. At θ=0\theta=0 every Legendre polynomial equals one, Pℓ(1)=1P_\ell(1)=1, so

Im f(0)=1k∑ℓ(2ℓ+1) Im ⁣[eiδℓsin⁡δℓ]=1k∑ℓ(2ℓ+1)sin⁡2δℓ(4)\mathrm{Im}\,f(0)=\frac{1}{k}\sum_\ell(2\ell+1)\,\mathrm{Im}\!\left[e^{i\delta_\ell}\sin\delta_\ell\right]=\frac{1}{k}\sum_\ell(2\ell+1)\sin^2\delta_\ell \qquad (4)

where: the only trigonometry is Im(eiδsin⁡δ)=sin⁡2δ\mathrm{Im}(e^{i\delta}\sin\delta)=\sin^2\delta — each channel contributes its sin⁡2δℓ\sin^2\delta_\ell, the same quantity measuring how strongly it scatters. This is the circle identity: Imfℓ=sin⁡2δℓ/k=k∣fℓ∣2\mathrm{Im}f_\ell=\sin^2\delta_\ell/k=k|f_\ell|^2.

Step 2. Compare with the cross-section sum. Partial-Wave Expansion's σ=(4π/k2)∑ℓ(2ℓ+1)sin⁡2δℓ\sigma=(4\pi/k^2)\sum_\ell(2\ell+1)\sin^2\delta_\ell is the same sum times 4π/k24\pi/k^2. Divide: σ/Imf(0)=4π/k\sigma/\mathrm{Im}f(0)=4\pi/k term by term — Eq. (1), now a channel-by-channel conservation statement, not a global accident. The machine check runs both sides independently (quadrature of ∫∣f∣2dΩ\int|f|^2\mathrm{d}\Omega vs the forward imaginary part) for hard spheres at ka=0.5,2,5ka=0.5,2,5 and a δ-shell resonance model swept through its peak, where σ0\sigma_0 saturates the 4π/k24\pi/k^2 ceiling: agreement to 10−1210^{-12} everywhere.

Step 3. Repair the source's factor of two. Our base text states the unitarity route as "2 Imf(0)=(k/4π)∫∣f∣2dΩ2\,\mathrm{Im}f(0)=(k/4\pi)\int|f|^2\mathrm{d}\Omega" — off by exactly 2 (the machine evaluates both sides for generic phase shifts; the ratio is 2, identically). The slip is instructive: the operator relation −i(T−T†)=T†T-i(T-T^\dagger)=T^\dagger T has left side 2 Im T2\,\mathrm{Im}\,T, and it is tempting to carry that 2 into the amplitude statement. But the normalization connecting ⟨k′∣T∣k⟩\langle\mathbf k'|T|\mathbf k\rangle to ff absorbs it: with dσ/dΩ=∣f∣2\mathrm{d}\sigma/\mathrm{d}\Omega=|f|^2, the on-shell relation is Imf(0)=(k/4π)∫∣f∣2dΩ\mathrm{Im}f(0)=(k/4\pi)\int|f|^2\mathrm{d}\Omega, no 2. Eq. (4) — three lines of trigonometry — is the arbiter. Verified: symbolic two-route equality for generic δ₀…δ₃ (sympy, exact); the base's factor-2 line reproduced and shown ≡ 2× the true relation; numeric two-route equality ≤10⁻¹² for hard sphere ka=0.5/2/5 and through the δ-shell resonance; ceiling saturation at the resonance confirmed to 0.2% (verification script).

What first order owes, and how the ledger collects

Here is the theorem's most useful day job: policing approximations. Born Approximation gives a real fBf_B for any real VV (a symmetric Fourier transform), so ImfB(0)=0\mathrm{Im}f_B(0)=0 while σB=∫∣fB∣2dΩ>0\sigma_B=\int|f_B|^2\mathrm{d}\Omega>0: first order violates the theorem by construction, scattering flux it never debits from the forward beam. Nothing is wrong — the theorem relates order nn amplitudes to order 2n2n cross sections, so the O(V2)O(V^2) deficit in σ\sigma is repaid by the O(V2)O(V^2) amplitude, second Born order. That repayment is exact, Born Series and Higher-Order Corrections' structural theorem: Imf(2)(0)=k4π∫∣f(1)∣2dΩ\mathrm{Im}f^{(2)}(0)=\frac{k}{4\pi}\int|f^{(1)}|^2\mathrm{d}\Omega, the on-shell (iεi\varepsilon) piece of the second-order propagator supplying the missing imaginary part. This article adds the measurement. For an attractive square well (ℏ=m=1\hbar=m=1, a=1a=1, k=1k=1) we tracked the deficit ∣Imfexact(0)−k4πσB∣\big|\mathrm{Im}f_{\mathrm{exact}}(0)-\tfrac{k}{4\pi}\sigma_B\big| as the depth was halved: it fell by 23.02^{3.0} — the residual closes as V3V^{3}, exactly as an identity satisfied through O(V2)O(V^2) demands. Symbolically too: expanding δℓ=λeℓ+O(λ2)\delta_\ell=\lambda e_\ell+O(\lambda^2), the O(λ2)O(\lambda^2) imaginary part of f(0)f(0) equals (k/4π)∫∣f(1)∣2dΩ(k/4\pi)\int|f^{(1)}|^2\mathrm{d}\Omega term for term, for arbitrary channel coefficients.

Figure
One number, both routes — and where Born sits
Im f(0) vs kσ_tot/4π — exact theories sit on the diagonal; first Born cannot kσ_tot/4π Im f(0) the theorem: slope exactly 1 hard sphere ka=0.5 ka=2 ka=5 δ-shell ON resonance: σ₀ = 4π/k² ceiling — still on the line (10⁻¹²) first Born: Im f_B(0) = 0 — off the line by the whole σ_B second Born repays: gap closes as V^3.0 (measured) every exact point machine-computed by two independent routes (quadrature vs forward value); Born points from the closed-form f_B

The forward floor: what the theorem forces on measurement

Because ∣f(0)∣2=[Ref(0)]2+[Imf(0)]2|f(0)|^2=[\mathrm{Re}f(0)]^2+[\mathrm{Im}f(0)]^2, the theorem sets a hard floor under the forward differential cross section. Our source text shipped this bound garbled (a literal "?" for the inequality's left side — an export corruption, repaired here and flagged in the cheat sheet); the correct statement is one line, worth deriving because a billion-dollar measurement program rests on it.

Go deeper · Derivation C · 1 minThe forward floor and the optical point
Derivation CThe forward floor and the optical point

We want the minimum forward peak a scatterer must show given its total cross section, and the dictionary turning this floor into the collider observable dσ/dt at t=0 — the "optical point". The payoff is the actual method by which total hadronic cross sections are measured: σ_tot is not directly countable at a collider (you cannot instrument 4π at millidegree angles), so experiments measure the forward elastic peak and let this floor-turned-equality do the rest.

Step 1. The floor. Split f(0)f(0) into real and imaginary parts and use Eq. (1):

dσdΩ∣0=∣f(0)∣2  ≥  [Im f(0)]2=(k σtot4π) ⁣2(5)\frac{\mathrm{d}\sigma}{\mathrm{d}\Omega}\bigg|_{0}=|f(0)|^2\;\ge\;[\mathrm{Im}\,f(0)]^2=\left(\frac{k\,\sigma_{\mathrm{tot}}}{4\pi}\right)^{\!2} \qquad (5)

where: equality holds when Ref(0)=0\mathrm{Re}f(0)=0 (purely absorptive). The floor is quadratic in σtot\sigma_{\mathrm{tot}} and grows with k2k^2: fast, strongly interacting systems are forced to be forward-peaked. A measured peak below this floor is a normalization error, not new physics.

Step 2. Promote to an equality with one parameter. Define ρ≡Ref(0)/Imf(0)\rho\equiv\mathrm{Re}f(0)/\mathrm{Im}f(0). Then exactly

dσdΩ∣0=(1+ρ2)(k σtot4π) ⁣2(6)\frac{\mathrm{d}\sigma}{\mathrm{d}\Omega}\bigg|_{0}=(1+\rho^2)\left(\frac{k\,\sigma_{\mathrm{tot}}}{4\pi}\right)^{\!2} \qquad (6)

where: ρ\rho is small when absorption dominates (pp at LHC: ρ≈0.10\rho\approx0.10, a 1% ρ2\rho^2 correction). Forward peak plus ρ\rho gives σtot\sigma_{\mathrm{tot}} — no 4π4\pi detector needed.

Step 3. The collider dictionary. Experiments bin in t=−q2=−2k2(1−cos⁡θ)t=-q^2=-2k^2(1-\cos\theta), so dσ/dt=(π/k2) dσ/dΩ\mathrm{d}\sigma/\mathrm{d}t=(\pi/k^2)\,\mathrm{d}\sigma/\mathrm{d}\Omega. At t=0t=0, Eq. (6) becomes the optical-point formula

dσdt∣0=(1+ρ2) σtot216π(7)\frac{\mathrm{d}\sigma}{\mathrm{d}t}\bigg|_{0}=\frac{(1+\rho^2)\,\sigma_{\mathrm{tot}}^{2}}{16\pi} \qquad (7)

where: every factor of kk cancels — the optical point depends only on σtot\sigma_{\mathrm{tot}} and ρ\rho, a clean experimental target. (Units: with σ\sigma in GeV⁻² via (ℏc)2=0.3894(\hbar c)^2=0.3894 mb·GeV², Eq. (7) returns GeV⁻⁴; ×(ℏc)2(\hbar c)^2 again for mb/GeV².) Verified: Eqs. (5)–(7) checked numerically on the np and pp numbers of the worked examples; the np floor is satisfied by 7 orders of magnitude of headroom, and the pp optical point evaluates to 631 mb/GeV² (verification script).

Limiting cases

caseresultmeaning
real potential, purely elasticσel=(4π/k) Imf(0)\sigma_{\mathrm{el}}=(4\pi/k)\,\mathrm{Im}f(0)every dented beam particle reappears in ∣f∣2\lvert f\rvert^2
absorbing target (∣Sℓ∣<1\lvert S_\ell\rvert<1)σel+σabs=(4π/k) Imf(0)\sigma_{\mathrm{el}}+\sigma_{\mathrm{abs}}=(4\pi/k)\,\mathrm{Im}f(0)forward dent counts undetected channels; anatomy in Inelastic Scattering and Multichannel Theory
on a resonance, δℓ=π/2\delta_\ell=\pi/2both sides saturate the 4π(2ℓ+1)/k24\pi(2\ell+1)/k^2 ceilingverified here to 10−1210^{-12} at the δ-shell peak; lineshape is Resonances and the Breit-Wigner Formula's
k→0k\to0 (s-wave world)Imf(0)=ka2/(1+k2a2)=kσ/4π\mathrm{Im}f(0)=ka^2/(1+k^2a^2)=k\sigma/4\pi exactlyscattering-length regime: both sides vanish linearly in kk; the amplitude turns real at threshold, owned by Low-Energy Scattering and the Scattering Length
first Born orderImfB(0)=0\mathrm{Im}f_B(0)=0 while σB>0\sigma_B>0the theorem fails at O(V2)O(V^2) and is repaid exactly at second order — Born Series and Higher-Order Corrections
black disk at high energyσel=σabs=σtot/2\sigma_{\mathrm{el}}=\sigma_{\mathrm{abs}}=\sigma_{\mathrm{tot}}/2maximal absorption still satisfies the theorem — half the "total" is the shadow's own diffraction

Worked examples

Example 1 — thermal neutrons on protons: how small an exact statement can be. The thermal n–p total cross section is σ=20.4\sigma=20.4 b — Differential and Total Cross Section's ladder anchor. What does the theorem force on the forward amplitude?

Step 1. Kinematics in the CM. A neutron at Elab=25.3E_{\mathrm{lab}}=25.3 meV on a nearly free proton: klab=2mnE/ℏc=3.49×10−5 fm−1k_{\mathrm{lab}}=\sqrt{2m_nE}/\hbar c=3.49\times10^{-5}\ \mathrm{fm}^{-1}, and with μ≈mn/2\mu\approx m_n/2 the relative wavenumber is k=1.75×10−5 fm−1k=1.75\times10^{-5}\ \mathrm{fm}^{-1}.

Step 2. Apply the theorem. Im f(0)=kσ/4π=(1.75×10−5×2044 fm)/4π\mathrm{Im}\,f(0)=k\sigma/4\pi=(1.75\times10^{-5}\times2044\ \mathrm{fm})/4\pi =2.8×10−3 fm=2.8\times10^{-3}\ \mathrm{fm}.

Step 3. Compare sizes. S-wave dominance gives ∣f∣≃σ/4π=12.8|f|\simeq\sqrt{\sigma/4\pi}=12.8 fm, so Imf(0)/∣f(0)∣≈2.2×10−4\mathrm{Im}f(0)/|f(0)|\approx2.2\times10^{-4}: the forward amplitude is 99.98% real. No contradiction — the theorem never says Imf\mathrm{Im}f is large, only exactly kσ/4πk\sigma/4\pi, and at threshold both sides vanish linearly in kk. Indeed f=−a/(1+ika)f=-a/(1+ika) gives Imf(0)=ka2/(1+k2a2)=kσ/4π\mathrm{Im}f(0)=ka^2/(1+k^2a^2)=k\sigma/4\pi identically (sympy-verified). Answer: Imf(0)=2.8×10−3\mathrm{Im}f(0)=2.8\times10^{-3} fm on a 12.8 fm nearly-real amplitude.

Example 2 — pp at 13 TeV: measuring σ_tot with the theorem. TOTEM's measured σtot=110.6\sigma_{\mathrm{tot}}=110.6 mb at s=13\sqrt s=13 TeV. Run the machine forward, then explain how the experiment actually ran it backward.

Step 1. Convert: σ=110.6/0.3894=284 GeV−2\sigma=110.6/0.3894=284\ \mathrm{GeV^{-2}}; CM momentum k≃s/2=6500k\simeq\sqrt s/2=6500 GeV.

Step 2. Forward amplitude: Imf(0)=kσ/4π=6500×284/4π=1.47×105 GeV−1\mathrm{Im}f(0)=k\sigma/4\pi=6500\times284/4\pi=1.47\times10^{5}\ \mathrm{GeV^{-1}}, i.e. × ℏc=0.1973 GeV fm\times\,\hbar c=0.1973\ \mathrm{GeV\,fm}: 2.9×104\mathbf{2.9\times10^{4}} fm — thirty nanometers of forward amplitude from a femtometer-sized proton. Forward peaking at high energy is not an option; it is Eq. (5).

Step 3. Optical point, Eq. (7): dσ/dt∣0=(1+0.102)×2842/16π=631 mb/GeV2\mathrm{d}\sigma/\mathrm{d}t|_0=(1+0.10^2)\times284^2/16\pi=\mathbf{631\ mb/GeV^2}. In practice the logic inverts: TOTEM measures the elastic peak, extrapolates to t=0t=0, takes ρ\rho from Coulomb–nuclear interference, and solves Eq. (7) for σtot\sigma_{\mathrm{tot}} — the optical theorem is the measurement, not a check. Answer: Imf(0)≈1.5×105\mathrm{Im}f(0)\approx1.5\times10^{5} GeV⁻¹; optical point ≈630\approx630 mb/GeV².

Origins · 1 minOptical Theorem and how this was discovered
OT
Optical Theorem
1871–1955 · Cambridge, Copenhagen

The theorem is older than quantum mechanics. In the 1870s optics already knew — through Sellmeier's dispersion work and Rayleigh's scattered-light analysis — that a beam's attenuation through a medium is fixed by the forward response of its scatterers, the optics "extinction theorem." Quantum mechanics inherited the ledger: in 1932 Eugene Feenberg proved the relation for matter waves (electron–atom scattering). Its full generality — total cross section including all reaction channels, from the forward elastic amplitude alone — crystallized in a 1939 nuclear-dispersion analysis by Bohr, Peierls and Placzek that was never published (the war intervened; physics cites the manuscript anyway). Postwar, van de Hulst made the extinction form standard and the "optical" name migrated into particle physics, where the theorem now converts LHC forward-peak measurements into total cross sections no detector could integrate directly.

History
1871
— Rayleigh's scattering analysis and Sellmeier's dispersion theory tie a medium's attenuation to the forward response of its scatterers — the optics ancestry of the theorem
1932
— Eugene Feenberg proves the quantum relation between total cross section and forward amplitude for electron scattering (Phys. Rev. 40, 40)
1939
— Bohr, Peierls and Placzek generalize it to all reaction channels in a nuclear-dispersion manuscript famously never published, yet cited for decades
1957
— van de Hulst's "Light Scattering by Small Particles" canonizes the extinction theorem in optics; the "optical theorem" name settles into scattering theory
2019
— TOTEM at the LHC measures ρ = Re f(0)/Im f(0) ≈ 0.09–0.10 at 13 TeV via Coulomb-nuclear interference — precision forward physics built on Eq. (1)

Going further

References. Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., §6.3 · Griffiths, Introduction to Quantum Mechanics, 3rd ed., Ch. 11 · Shankar, Principles of Quantum Mechanics, 2nd ed., Ch. 19 · Cohen-Tannoudji, Diu & Laloë, Quantum Mechanics, Vol. II, Ch. VIII · van de Hulst, Light Scattering by Small Particles (1957) · R. G. Newton, "Optical theorem and beyond," Am. J. Phys. 44, 639 (1976)

Provenance. Verifier should recompute: (1) the flux-integral proof, confirming the stationary-phase cross term gives −(4π/k)Im f(0)-(4\pi/k)\mathrm{Im}\,f(0) (Sakurai §6.3, Cohen-Tannoudji Vol. II Ch. VIII); (2) the partial-wave identity that both Im f(0)\mathrm{Im}\,f(0) and σtot\sigma_{\mathrm{tot}} equal 4πk2∑(2ℓ+1)sin⁡2δℓ\frac{4\pi}{k^2}\sum(2\ell+1)\sin^2\delta_\ell up to the 4π/k4\pi/k factor, using Pℓ(1)=1P_\ell(1)=1 and Legendre orthogonality; (3) the unitarity-limit value 4πk2(2ℓ+1)\frac{4\pi}{k^2}(2\ell+1) at δℓ=π/2\delta_\ell=\pi/2; (4) the second-Born relation Im f(2)(0)=k4π∫∣f(1)∣2dΩ\mathrm{Im}\,f^{(2)}(0)=\frac{k}{4\pi}\int|f^{(1)}|^2\mathrm{d}\Omega. All formulas are standard (Sakurai, Griffiths Ch. 11, Shankar Ch. 19). (Lane note: items (1)–(3) are executed in this article's verification script, item (4) symbolically at order λ²; the systematic second-Born derivation is chapter property of Born Series and Higher-Order Corrections. Two base display defects repaired, machine-proven: the "generalized optical theorem" line carried a spurious factor 2, and the forward-bound display shipped with a literal "?" — both corrected in Derivation B and Derivation C respectively.)

Builds on: Scattering Amplitude · Asymptotic Form of the Scattering Wave Function · Differential and Total Cross Section · Partial-Wave Expansion

Leads to: The S-Matrix and Unitarity in Scattering · Phase Shifts · Resonances and the Breit-Wigner Formula · Born Series and Higher-Order Corrections

If the interference term only survives at θ = 0, why does anything else matter to the theorem?

Because the value the forward direction retains is an integral over everything else. Stationary phase pins the surviving interference to f(0)f(0), but conservation equates it to ∫∣f∣2dΩ\int|f|^2\mathrm{d}\Omega plus all absorption — that is the theorem. The forward direction is not where the physics happens; it is where the accounting is filed. Which is why one measured point (plus ρ\rho) fixes σtot\sigma_{\mathrm{tot}} at a collider.

Is Im f(0) > 0 guaranteed?

Imf(0)=kσtot/4π≥0\mathrm{Im}f(0)=k\sigma_{\mathrm{tot}}/4\pi\ge0 always, with equality only for no scattering. That fixed sign and size powers dispersion relations: the energy dependence of Imf(0)∝σtot(E)\mathrm{Im}f(0)\propto\sigma_{\mathrm{tot}}(E) constrains Ref(0,E)\mathrm{Re}f(0,E) at each energy — the LHC's measured ρ\rho parameter is exactly such a dispersion-relation output, which is why its 2019 value carried theoretical weight.

Does the theorem hold at second order in Born, or only in the exact theory?

Order by order, interlocking: Imf(2)(0)=(k/4π)∫∣f(1)∣2dΩ\mathrm{Im}f^{(2)}(0)=(k/4\pi)\int|f^{(1)}|^2\mathrm{d}\Omega, then Imf(3)\mathrm{Im}f^{(3)} pairs with the f(1)∗f(2)f^{(1)*}f^{(2)} cross terms, and so on — each order's imaginary part pays the previous orders' bill. The exact amplitude satisfies it identically; truncate and you owe the highest order's debt, which is why our square-well deficit closes as V3V^3.

The Coulomb potential scatters plenty — why does the theorem break there?

Because σtot\sigma_{\mathrm{tot}} diverges (the 1/sin⁡4(θ/2)1/\sin^4(\theta/2) forward divergence integrates to infinity) and the asymptotic state never becomes a clean plane wave — logarithmic phase distortions persist at all distances, as Asymptotic Form of the Scattering Wave Function details. Both sides of Eq. (1) are infinite; the finite statement lives one article over: screen the potential (Scattering by a Screened Coulomb Potential) and the theorem holds with the screening length setting the forward scale.

How can the LHC measure a total cross section it cannot see? Detectors cover almost no forward angle at all.

That is precisely the theorem's job. Elastic protons at milliradian angles are tracked in Roman-pot detectors far downstream; the measured dσ/dt\mathrm{d}\sigma/\mathrm{d}t is extrapolated to t=0t=0; Coulomb–nuclear interference at ultra-small ∣t∣|t| fixes ρ\rho; and Eq. (7) converts the optical point into σtot\sigma_{\mathrm{tot}}. The luminosity-independent variant compares elastic and inelastic rates so that even the beam normalization drops out. Every quoted "total pp cross section" above ISR energies is an optical-theorem number.

Is there a version for identical particles or for classical waves?

Yes to both. For identical particles the theorem applies to the properly (anti)symmetrized f(0)f(0) — the bookkeeping shifts but conservation does not; Identical-Particle Scattering handles it. And the identity predates quantum mechanics: any linear wave with a conserved quadratic flux — light, sound, water waves around a piling — obeys an extinction theorem of identical structure, which is where the "optical" name came from and why antenna engineers use Eq. (1) with ff a classical field amplitude.

Before this

Where this leads