Low-Energy Scattering and the Scattering Length

Cool a collision far enough and the projectile goes blind. A thermal neutron carries a de Broglie wavelength of 1.801.80 Å — four orders of magnitude larger than the nucleus it strikes. At that resolution the potential's depth, shape, and wiggles are invisible; the only question the wave can still ask is where the target pushes its phase. The answer developed here: one real number — the scattering length asa_s — replaces the entire potential. This article owns the threshold machinery of the chapter. We derive the Wigner threshold law tan⁡δℓ∝(kR)2ℓ+1\tan\delta_\ell\propto(kR)^{2\ell+1} from the matching formula of Phase Shifts — cashing the statement that Effective Potential and the Centrifugal Barrier explicitly banked to this chapter — which is why only the s-wave survives as k→0k\to0. We derive the zero-energy trio δ0→−kas\delta_0\to-ka_s, f→−asf\to-a_s, σ→4πas2\sigma\to4\pi a_s^2, resolving the factor-of-four surprise that Differential and Total Cross Section stated and deeded here. We derive the effective range expansion kcot⁡δ0=−1/as+12r0k2k\cot\delta_0=-1/a_s+\tfrac12 r_0k^2 by Bethe's two-energy comparison and verify it numerically on an exactly solvable well. And we compute the universal curve as(V0)a_s(V_0) for the square well — divergences where bound states are born, zeros where the target turns transparent. The payoff is the most famous number in low-energy nuclear physics: assembling the spin-weighted neutron–proton cross section from two measured lengths, π(3at2+as2)=20.5\pi(3a_t^2+a_s^2)=20.5 b — the anchor the chapter opener quoted, finally computed from first principles.

Key takeaways
  • At ka≪1ka\ll1 the centrifugal barrier silences every ℓ≥1\ell\ge1 channel: tan⁡δℓ∝(ka)2ℓ+1\tan\delta_\ell\propto(ka)^{2\ell+1} (the Wigner threshold law, derived here from the matching formula). Low-energy scattering is an s-wave monopoly.
  • The s-wave phase shift vanishes linearly, δ0→−kas\delta_0\to-k a_s, defining the scattering length as=−lim⁡k→0tan⁡δ0/ka_s=-\lim_{k\to0}\tan\delta_0/k — geometrically, the rr-intercept of the straight-line zero-energy wavefunction u0∝(r−as)u_0\propto(r-a_s).
  • Everything at threshold follows from asa_s: amplitude f→−asf\to-a_s (isotropic), cross section σ→4πas2\sigma\to4\pi a_s^2 — four times the classical silhouette, a pure wave effect.
  • The next correction is one more length: kcot⁡δ0=−1/as+12r0k2k\cot\delta_0=-1/a_s+\tfrac12 r_0k^2, the effective range expansion — verified here on the square well to 0.1%0.1\% (analytic as=11.13a_s=11.13 fm, r0=2.27r_0=2.27 fm against a numerical fit).
  • The square-well master curve as=R[1−tan⁡(K0R)/(K0R)]a_s=R[1-\tan(K_0R)/(K_0R)] diverges each time a new bound state arrives (K0R=π/2, 3π/2,…K_0R=\pi/2,\,3\pi/2,\dots) and crosses zero at tan⁡(K0R)=K0R\tan(K_0R)=K_0R — near-threshold poles and transparencies, the doorway to Bound States and Zero-Energy Resonances.
  • Flagship: at=+5.42a_t=+5.42 fm and as=−23.7a_s=-23.7 fm assemble into σnp=π(3at2+as2)=20.5\sigma_{np}=\pi(3a_t^2+a_s^2)=20.5 b — the thermal neutron–proton cross section, spin bookkeeping included.
  • Key idea: at low energy a short-range potential scatters only s-waves, and a single length — the scattering length asa_s — replaces the entire potential, giving an isotropic amplitude f→−asf\to-a_s and a constant cross section σ→4πas2\sigma\to4\pi a_s^2.

Where this lives today

The scattering length may be the most used single number in quantum physics. Every cold-atom paper begins by quoting asa_s: the Gross–Pitaevskii coupling g=4πℏ2as/mg=4\pi\hbar^2a_s/m, healing length, loss rates, and the very sign of stability all flow from it, and near a Feshbach resonance a magnetic field sweeps asa_s through ±∞\pm\infty — deliberately driving the system across the divergences of the universal curve to make molecules, unitary Fermi gases, and Efimov trimers to order. In neutron science the coherent spin-channel combination is the material constant of reflectometry and crystallography. And pionless effective field theory expands in r0/asr_0/a_s, treating the deuteron and the np virtual state as two faces of one near-threshold pole — the unification the next article completes.

The physical picture

Picture a slow wave, wavelength meters long, washing over a pebble. It cannot resolve the pebble's shape — but it does respond. Outside the pebble the zero-energy s-wave is exact and trivial: u0′′=0u_0''=0, whose complete solution is the straight line u0(r)∝(r−as)u_0(r)\propto(r-a_s). Whatever violence the potential does inside, the exterior remembers one number: where that line crosses zero. That intercept is the scattering length. The sign tells a story. A repulsive core pushes the wavefunction out: the line crosses at positive rr, and for a hard sphere the crossing sits exactly at the surface, as=Ra_s=R. A weak attraction pulls the wavefunction in: the line now rises from a negative intercept — as<0a_s<0. Strengthen the attraction toward the point where a bound state is about to form and the exterior line tilts toward horizontal: its intercept races off to infinity. A diverging scattering length is a bound state announcing itself at threshold — the pole story that Bound States and Zero-Energy Resonances owns.

Figure
The zero-energy wavefunction is a straight line and its intercept is the scattering length
u₀(r) ∝ (r − aₛ): three potentials, three intercepts r u₀ potential region r ≲ R aₛ positive (repulsive: pushed out) aₛ negative weak attraction: pulled in, intercept behind the origin near binding threshold: line almost flat — aₛ → ∞ outside the potential u₀″ = 0 exactly: a straight line. The intercept is the ONLY thing a k→0 wave can measure.

Why is the cross section 4πas24\pi a_s^2 and not πas2\pi a_s^2? Because a k→0k\to0 wave does not cast a shadow — it breathes. The scattered wave at threshold is a pure isotropic s-wave emitted coherently by the whole target; its strength is set by the intercept asa_s, and integrating the isotropic ∣f∣2=as2|f|^2=a_s^2 over all 4π4\pi steradians gives the full sphere's worth, not the disk's. Differential and Total Cross Section flagged exactly this — "four times the classical shadow πR2\pi R^2" — and named this article the owner of the derivation, delivered below. One more face, the experimentalist's: plot kcot⁡δ0k\cot\delta_0 against k2k^2 and low-energy data collapse onto a straight line. The intercept reads off −1/as-1/a_s; the slope reads off half the effective range r0r_0. Two numbers, measured from a line on graph paper, summarize any short-range potential — that is the shape independence that built low-energy nuclear physics.

The threshold laws

Everything rests on the matching formula of Phase Shifts: the interior fixes a logarithmic derivative at the potential's edge; comparison with free spherical waves yields tan⁡δℓ\tan\delta_\ell. We need only its low-kk limit.

Derivation AThe Wigner threshold law from one matching formula

We prove that tan⁡δℓ∝k2ℓ+1\tan\delta_\ell\propto k^{2\ell+1} as k→0k\to0 for any potential of finite range RR — the law whose statement Effective Potential and the Centrifugal Barrier explicitly banked to this chapter ("owned with its derivation by the scattering chapter"). The strategy: take the master matching formula of Phase Shifts as given, insert the small-argument behavior of the free spherical waves, and watch the powers of kk organize themselves. The only subtlety worth guarding is that the interior logarithmic derivative tends to a finite, kk-independent constant as k→0k\to0 — true for finite-range potentials because the interior equation at E=0E=0 is perfectly regular.

Step 1. Start from the matching formula (recalled by name from Phase Shifts, not re-derived): with γℓ\gamma_\ell the logarithmic derivative of the interior radial function at r=Rr=R,

tan⁡δℓ=k jℓ′(kR)−γℓ jℓ(kR)k nℓ′(kR)−γℓ nℓ(kR)(1)\tan\delta_\ell=\frac{k\,j_\ell'(kR)-\gamma_\ell\, j_\ell(kR)}{k\,n_\ell'(kR)-\gamma_\ell\, n_\ell(kR)} \qquad (1)

where: jℓ,nℓj_\ell,n_\ell — spherical Bessel and Neumann functions; primes are derivatives with respect to the argument. As k→0k\to0, γℓ→\gamma_\ell\to const.

Step 2. Insert the small-argument forms jℓ(x)→xℓ/(2ℓ+1)!!j_\ell(x)\to x^\ell/(2\ell+1)!! and nℓ(x)→−(2ℓ−1)!!/xℓ+1n_\ell(x)\to-(2\ell-1)!!/x^{\ell+1}. The numerator scales as the regular family: every term carries (kR)ℓ⋅k(kR)^\ell\cdot k or (kR)ℓγℓ(kR)^\ell\gamma_\ell — overall kℓk^\ell. The denominator is dominated by the irregular family: nℓ∼k−(ℓ+1)n_\ell\sim k^{-(\ell+1)}. Why the powers cannot conspire away. The numerator's two terms could in principle cancel at leading order — that is precisely what happens at a zero-energy resonance, and it is non-generic. Generically the ratio scales as kℓ/k−(ℓ+1)=k2ℓ+1k^{\ell}/k^{-(\ell+1)}=k^{2\ell+1}.

Step 3. Hence

tan⁡δℓ  →k→0  cℓ (kR)2ℓ+1,cℓ=O(1)(2)\tan\delta_\ell \;\xrightarrow{k\to0}\; c_\ell\,(kR)^{2\ell+1},\qquad c_\ell=O(1) \qquad (2)

where: cℓc_\ell depends on the interior physics through γℓ\gamma_\ell but the power is universal. Each unit of ℓ\ell costs two powers of kaka — the centrifugal barrier's tunneling toll, now quantitative.

Step 4. Read the hierarchy: δ0∼kR\delta_0\sim kR survives; δ1∼(kR)3\delta_1\sim(kR)^3 and δ2∼(kR)5\delta_2\sim(kR)^5 die. For a thermal neutron on a nucleus kR∼10−4kR\sim10^{-4}: the p-wave is down eight orders of magnitude. One channel, hence — next derivation — one number. Verified: the sympy expansion of formula (1) with the exact jℓ,nℓj_\ell,n_\ell for ℓ=0,1,2\ell=0,1,2 returns lowest powers k1,k3,k5k^1,k^3,k^5 exactly; the exponent for ℓ=1\ell=1 was also measured on a numerically integrated square well (log–log slope 3.0003.000). Article verifier, 26 checks.

The law explains the monopoly; the next card cashes it into observables.

Go deeper · Derivation B · 1 minOne number at threshold: the intercept, the amplitude, and four pi a squared
Derivation BOne number at threshold: the intercept, the amplitude, and four pi a squared

We define the scattering length three equivalent ways — slope of the phase shift, intercept of the zero-energy wavefunction, zero-energy amplitude — and prove the threshold cross section σ→4πas2\sigma\to4\pi a_s^2, resolving the hard-sphere factor-of-four that Differential and Total Cross Section deeded to this article. The route: linearize the s-wave phase shift, then let the exterior wavefunction's geometry speak.

Step 1. Define, from the threshold law's ℓ=0\ell=0 rung,

as≡−lim⁡k→0tan⁡δ0(k)k(3)a_s\equiv-\lim_{k\to0}\frac{\tan\delta_0(k)}{k} \qquad (3)

where: the sign is the Sakurai–Shankar convention, chosen so a hard sphere gives as=+Ra_s=+R (some nuclear texts flip it — see the pitfall box). Dimensionally asa_s is a length: the only length the k→0k\to0 limit can produce.

Step 2. Geometry. Outside the potential the s-wave is u0(r)∝sin⁡(kr+δ0)u_0(r)\propto\sin(kr+\delta_0); as k→0k\to0 the sine linearizes, u0∝krcos⁡δ0+sin⁡δ0∝r+tan⁡δ0/k→r−asu_0\propto k r\cos\delta_0+\sin\delta_0\propto r+\tan\delta_0/k\to r-a_s. The zero-energy exterior is the straight line of the hero figure, and its rr-intercept is asa_s. Why this is the honest definition. Equation (3) references scattering data; the straight line references the wavefunction. They agree identically, which is what lets bound-state physics (wavefunctions) and scattering physics (phase shifts) trade information — the bridge Bound States and Zero-Energy Resonances will drive across.

Step 3. Amplitude. The exact s-wave amplitude from the partial-wave ladder of Partial-Wave Expansion is f0=1/(kcot⁡δ0−ik)f_0=1/(k\cot\delta_0-ik). As k→0k\to0, kcot⁡δ0→−1/ask\cot\delta_0\to-1/a_s while ik→0ik\to0:

f  →k→0  −as(4)f\;\xrightarrow{k\to0}\;-a_s \qquad (4)

where: the amplitude is real, isotropic, and energy-independent — a point target with a single voice.

Step 4. Cross section. With only the s-wave alive, σ=(4π/k2)sin⁡2δ0→(4π/k2)(kas)2\sigma=(4\pi/k^2)\sin^2\delta_0\to(4\pi/k^2)(ka_s)^2:

σ  →k→0  4πas2(5)\sigma\;\xrightarrow{k\to0}\;4\pi a_s^2 \qquad (5)

where: for the hard sphere (as=Ra_s=R, from δ0=−kR\delta_0=-kR exactly) this is 4πR24\pi R^2 — the entire surface area, four times the geometric shadow πR2\pi R^2. The wave scatters coherently off the whole sphere, not off a silhouette. Verified: f0→−asf_0\to-a_s and σ→4πas2\sigma\to4\pi a_s^2 symbolic; hard-sphere tan⁡δ0=−tan⁡kR\tan\delta_0=-\tan kR hence as=Ra_s=R symbolic; hard-sphere tan⁡δ1→−(kR)3/3\tan\delta_1\to-(kR)^3/3 (negative, as repulsion demands). Article verifier.

The effective range expansion

The scattering length is the k=0k=0 snapshot. Real experiments run at small but finite kk, and the first correction turns out to be another length — not a new function. This robustness is why two numbers per channel organized all of low-energy nuclear physics.

Go deeper · Derivation C · 2 minThe effective range from a two-energy comparison
Derivation CThe effective range from a two-energy comparison

We derive kcot⁡δ0=−1/as+12r0k2+O(k4)k\cot\delta_0=-1/a_s+\tfrac12r_0k^2+O(k^4) and identify r0r_0 as an integral over the region where the true wavefunction differs from its own asymptote — Bethe's 1949 comparison argument, the one he called the shape-independent approximation. The strategy: write the s-wave problem at two energies, cross-multiply the two Schrödinger equations, integrate, and let the boundary terms do the work. What to watch: the comparison functions must be normalized to the same asymptotic line, or the boundary terms lie.

Step 1. Let uk(r)u_k(r) solve the s-wave equation at energy ℏ2k2/2m\hbar^2k^2/2m with uk(0)=0u_k(0)=0, and let vk(r)v_k(r) be the free solution matching uku_k exactly outside the range, normalized so vk(0)=1v_k(0)=1 at every energy. Write the radial equations at k1k_1 and k2k_2, cross-multiply, subtract, integrate 0→∞0\to\infty: the potential terms cancel identically, leaving boundary terms.

Step 2. The identical manipulation applied to the free pair (vk1,vk2)(v_{k_1},v_{k_2}) gives boundary terms at both ends; for the true pair (uk1,uk2)(u_{k_1},u_{k_2}) the origin contributes nothing since u(0)=0u(0)=0. Subtracting the two identities:

k2cot⁡δ0(k2)−k1cot⁡δ0(k1)=(k22−k12)∫0∞ ⁣(vk1vk2−uk1uk2) dr(6)k_2\cot\delta_0(k_2)-k_1\cot\delta_0(k_1)=(k_2^2-k_1^2)\int_0^\infty\!\big(v_{k_1}v_{k_2}-u_{k_1}u_{k_2}\big)\,\mathrm{d}r \qquad (6)

where: the integrand lives only inside the potential — outside, u=vu=v by construction and the integrand vanishes. All shape dependence is confined to a short-range integral.

Step 3. Set k1=0k_1=0, k2=kk_2=k, and expand the integral at small kk (its kk-dependence is analytic because the integrand has finite support):

kcot⁡δ0(k)=−1as+12r0k2+O(k4),r0=2∫0∞ ⁣(v02−u02) dr(7)k\cot\delta_0(k)=-\frac{1}{a_s}+\frac{1}{2}r_0k^2+O(k^4),\qquad r_0=2\int_0^\infty\!\big(v_0^2-u_0^2\big)\,\mathrm{d}r \qquad (7)

where: r0r_0 — the effective range — measures the spatial extent over which the true zero-energy wavefunction fails to be its own asymptote. It is a range, not a strength: potentials of wildly different depth but similar reach share r0r_0.

Step 4. Instant payoff — the hard sphere. There u0=0u_0=0 inside (r<Rr<R) and u0=v0=1−r/Ru_0=v_0=1-r/R outside, so r0=2∫0R(1−r/R)2dr=2R/3r_0=2\int_0^R(1-r/R)^2\mathrm{d}r=2R/3. A hard sphere of radius 22 fm has as=2a_s=2 fm, r0=1.33r_0=1.33 fm — both lengths of order the geometry, as they should be.

Step 5. Machine check on a solvable well (V0=20V_0=20 MeV, R=2.5R=2.5 fm, np reduced mass). Analytic: as=R−tan⁡(K0R)/K0=11.13a_s=R-\tan(K_0R)/K_0=11.13 fm, r0=R−1/(K02as)−R3/(3as2)=2.272r_0=R-1/(K_0^2a_s)-R^3/(3a_s^2)=2.272 fm. Numerical: exact δ0(k)\delta_0(k) at twelve momenta, fit kcot⁡δ0k\cot\delta_0 vs k2k^2 — a line to 0.1%0.1\%, returning as=11.13a_s=11.13 fm, r0=2.274r_0=2.274 fm. Two lengths exhaust the low-energy content. Verified: hard-sphere r0=2R/3r_0=2R/3 symbolic (Bethe integral); square-well asa_s closed form symbolic; the kcot⁡δ0k\cot\delta_0 linearity, fitted (as,r0)(a_s,r_0) vs analytic values, all in the article verifier.

Figure
The experimentalist's straight line: k·cot δ₀ versus k² for the solvable well
k cot δ₀ = −1/aₛ + ½ r₀ k² (square well, V₀ = 20 MeV, R = 2.5 fm) k² (fm⁻²) k cot δ₀ (fm⁻¹) intercept = −1/aₛ = −0.0898 fm⁻¹ ⇒ aₛ = 11.13 fm slope = r₀/2 = 1.136 fm ⇒ r₀ = 2.27 fm exact phase shifts (dots) vs the two-parameter line: residual under 0.1% — “shape independence”: only aₛ and r₀ survive at low k.

The universal curve: where scattering lengths come from

For the attractive square well the closed form is worth staring at. With K0=2mV0/ℏK_0=\sqrt{2mV_0}/\hbar the zero-energy interior wavenumber,

as=R[1−tan⁡(K0R)K0R](8)a_s=R\left[1-\frac{\tan(K_0R)}{K_0R}\right] \qquad (8)

where: the whole personality of the well enters through the single dial K0RK_0R. Weak well (K0R≪1K_0R\ll1): as≈−R(K0R)2/3<0a_s\approx-R(K_0R)^2/3<0 — the generic weak-attraction signature. As K0R→π/2K_0R\to\pi/2 the tangent blows up: as→−∞a_s\to-\infty, reappearing at +∞+\infty past threshold. Each divergence marks a bound state's birth; each zero (at tan⁡K0R=K0R=4.493,…\tan K_0R=K_0R=4.493,\dots) makes the well invisible to zero-energy projectiles.

Figure
The universal curve aₛ against well depth: divergences at bound-state thresholds, zeros at transparency
aₛ/R = 1 − tan(K₀R)/(K₀R): the scattering length as the well deepens K₀R aₛ/R K₀R = π/2 K₀R = 3π/2 aₛ = 0 at tan(K₀R) = K₀R = 4.493: transparency (the Ramsauer–Townsend Effect is this zero in an atom) deuteron-like: one shallow bound state, aₛ large, positive np singlet: almost bound, aₛ = −23.7 fm (virtual) no bound state one bound state two bound states each divergence = a bound state born at threshold (pole crossing k = 0 — the next article's story); weak wells start small and NEGATIVE.

The two marked points are the whole of low-energy nuclear physics. The neutron–proton triplet channel sits just past the first divergence: one shallow bound state — the deuteron — and a large positive at=+5.42a_t=+5.42 fm. The singlet channel sits just before it: no bound state, but so nearly one that as=−23.7a_s=-23.7 fm — enormous and negative. What "almost bound" means precisely — an S-matrix pole on the wrong axis, a virtual state — is owned by Bound States and Zero-Energy Resonances. The deuteron connection is quantitative with tools already on the table. Spherical Hankel Functions and Outgoing Waves computed the deuteron's exterior decay length 1/κ=4.321/\kappa=4.32 fm from EB=2.2246E_B=2.2246 MeV. That is not the scattering length — the effective range bridges the 20% gap: continuing kcot⁡δ0k\cot\delta_0 to the pole gives 1/at=κ−12rtκ21/a_t=\kappa-\tfrac12 r_t\kappa^2, and rt=1.76r_t=1.76 fm predicts at=5.42a_t=5.42 fm from the binding energy alone — matching the measured value. Two lengths, one channel, closed books.

Limiting cases

caseresultmeaning
hard sphere, k→0k\to0as=Ra_s=R, r0=2R/3r_0=2R/3, σ=4πR2\sigma=4\pi R^2intercept at the surface; full surface area, not the shadow
weak attraction (K0R≪1K_0R\ll1)as≈−R(K0R)2/3<0a_s\approx-R(K_0R)^2/3<0small negative asa_s: generic signature of an unbinding attraction
binding threshold K0R→π/2K_0R\to\pi/2as→±∞a_s\to\pm\infty, σ→\sigma\to hugezero-energy resonance; pole at k=0k=0 (owner: Bound States and Zero-Energy Resonances)
transparency tan⁡K0R=K0R\tan K_0R=K_0Ras=0a_s=0: no k→0k\to0 scatteringthe well is invisible — the Ramsauer–Townsend mechanism (its article, by name)
kas≪1ka_s\ll1, asa_s fixedσ=4πas2\sigma=4\pi a_s^2 constantthreshold plateau; energy-independence is the s-wave fingerprint
k∣as∣≫1k\vert a_s\vert\gg1 (huge asa_s)σ→4π/k2\sigma\to4\pi/k^2the unitarity ceiling of The S-Matrix and Unitarity in Scattering: asa_s drops out entirely
ℓ≥1\ell\ge1 at thresholdtan⁡δℓ∝(kR)2ℓ+1\tan\delta_\ell\propto(kR)^{2\ell+1}Wigner law (derived here): centrifugal barrier silences all but s-waves

Worked examples

Example 1 — the 20.5 barns of the thermal neutron. Differential and Total Cross Section anchored its benchmark ladder with "n+p elastic: 20.4 b" while its geometric estimate was ∼1\sim1 b. Assemble that number from the two measured scattering lengths: at=+5.424a_t=+5.424 fm (spin-triplet), as=−23.749a_s=-23.749 fm (spin-singlet).

Step 1. Spin bookkeeping. An unpolarized neutron meets an unpolarized proton in a random relative spin state: of the four two-spin states, three are triplet and one is singlet, so the statistical weights are 3/43/4 and 1/41/4.

Step 2. Each channel scatters with its own length: σt=4πat2=4π(5.424)2=369.7 fm2=3.70\sigma_t=4\pi a_t^2=4\pi(5.424)^2=369.7\ \mathrm{fm}^2=3.70 b and σs=4πas2=4π(23.749)2=7087 fm2=70.9\sigma_s=4\pi a_s^2=4\pi(23.749)^2=7087\ \mathrm{fm}^2=70.9 b. The nearly-bound singlet channel scatters twenty times harder than the channel with the actual bound state — proximity to threshold, not binding, is what matters.

Step 3. Weight and add (cross sections add because triplet and singlet final states are distinguishable in principle):

σnp=34σt+14σs=π(3at2+as2)=2049 fm2\sigma_{np}=\tfrac34\sigma_t+\tfrac14\sigma_s=\pi\big(3a_t^2+a_s^2\big)=2049\ \mathrm{fm}^2

where: 11 b =100 fm2=100\ \mathrm{fm}^2, so σnp=20.5\sigma_{np}=20.5 b — the chapter opener's anchor, recovered from two intercepts. Answer: σnp=20.5\sigma_{np}=20.5 b.

Step 4. The instructive wrong answer: averaging lengths first — aˉ=(3at+as)/4=−1.87\bar a=(3a_t+a_s)/4=-1.87 fm, 4πaˉ2=0.444\pi\bar a^2=0.44 b — misses by 47×. Probabilities average, amplitudes don't. (That coherent aˉ\bar a is physical in neutron optics as the coherent scattering length — a different, interference, question.)

Example 2 — an ultracold gas is a scattering length with a thermometer. ⁸⁷Rb atoms have as≈100 a0=5.29a_s\approx100\,a_0=5.29 nm. Compute the collision cross section at T=100T=100 nK.

Step 1. Check the monopoly: relative momentum k=2μkBT/ℏk=\sqrt{2\mu k_BT}/\hbar, μ=mRb/2\mu=m_{\mathrm{Rb}}/2, gives kas=0.022≪1ka_s=0.022\ll1 — deep threshold, p-wave suppression ∼(ka)3≈10−5\sim(ka)^3\approx10^{-5}.

Step 2. Identical bosons in one internal state interfere constructively in the s-wave, doubling the cross section: σ=8πas2\sigma=8\pi a_s^2 (symmetrization belongs to Identical-Particle Scattering; we import the factor).

Step 3. Numbers: σ=8π(5.29×10−9 m)2=7.0×10−16 m2\sigma=8\pi(5.29\times10^{-9}\ \mathrm{m})^2=7.0\times10^{-16}\ \mathrm{m}^2 — about 25002500 times the atom's geometric size. Answer: σ≈7.0×10−16 m2\sigma\approx7.0\times10^{-16}\ \mathrm{m}^2, with kas=0.022ka_s=0.022 certifying pure threshold physics. Every BEC mean-field number — healing length, sound speed, stability — enters through this one length via g=4πℏ2as/mg=4\pi\hbar^2a_s/m.

Origins · 1 minScattering Length and how this was discovered
SL
Scattering Length
1936–1949 · Rome, Ithaca

The scattering length was born twice in one decade. In Rome, Enrico Fermi's 1936 slow-neutron work introduced the pseudopotential — a point interaction tuned to reproduce one measured length — so neutron waves in hydrogenous materials could be computed without knowing the nuclear force; its regularized descendant g δ3(r)g\,\delta^3(\mathbf r), g=4πℏ2as/mg=4\pi\hbar^2a_s/m, still runs cold-atom theory. After the war the systematics crystallized: Julian Schwinger's variational treatment of np scattering circulated in the late 1940s, and in 1949 Hans Bethe published the transparent two-energy derivation of the effective range expansion — the "shape-independent approximation" — with Blatt and Jackson fitting the np data the same year. Two lengths per channel carried nuclear physics for decades, then found a second life when ultracold-atom experiments learned to tune asa_s through Feshbach resonances, making the scattering length the most engineered number in many-body physics.

History
1936
— Fermi introduces the pseudopotential for slow neutrons in hydrogenous media: one length stands in for the nuclear force
1947
— Schwinger's variational effective-range results for np scattering circulate (Harvard notes; published abstracts)
1949
— Bethe's two-energy derivation of the effective range expansion (Phys. Rev. 76, 38); Blatt and Jackson fit the np lengths the same year
1957
— neutron interferometry and reactor physics standardize coherent scattering lengths as tabulated material constants
1998
— a magnetically tuned Feshbach resonance in a sodium BEC (Ketterle group) turns the scattering length into an experimental dial

Going further

  • Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., §6.4–6.5 — low-energy scattering, the scattering length, and the square-well curve.
  • Shankar, Principles of Quantum Mechanics, 2nd ed., Ch. 19 — phase shifts and threshold behavior.
  • Griffiths, Introduction to Quantum Mechanics, 3rd ed., §11.3 — the hard sphere and partial waves at low energy.
  • Cohen-Tannoudji, Diu & Laloë, Quantum Mechanics, Vol. II, Ch. VIII — the square-well scattering length and its divergences.
  • Taylor, Scattering Theory, Ch. 11; Newton, Scattering Theory of Waves and Particles, Ch. 11 — effective-range theory done rigorously.
  • Bethe, Phys. Rev. 76, 38 (1949) — the two-energy derivation; Blatt & Jackson, Phys. Rev. 76, 18 (1949) — the np fits.

Provenance. All standard results from low-energy scattering theory. A verifier should: (1) confirm f0→−asf_0\to -a_s and σ→4πas2\sigma\to4\pi a_s^2 from the partial-wave amplitude (Sakurai §6.5, Shankar Ch.19); (2) reproduce the hard-sphere values δ0=−kR\delta_0=-kR, as=Ra_s=R, r0=2R/3r_0=2R/3 (Griffiths §11.3); (3) verify the square-well matching as=R[1−tan⁡(K0R)/(K0R)]a_s=R[1-\tan(K_0R)/(K_0R)] and its divergence at K0R=(n+12)πK_0R=(n+\tfrac12)\pi (Cohen-Tannoudji Vol. II); (4) check the Bethe effective-range derivation and integral for r0r_0 (Newton, Taylor). Quoted np data (at≈+5.4a_t\approx+5.4 fm, r0t≈1.75r_{0t}\approx1.75 fm, as≈−23.7a_s\approx-23.7 fm, deuteron EB≈2.22E_B\approx2.22 MeV) are standard tabulated nuclear values (e.g. Krane, Introductory Nuclear Physics). Contact coupling g=4πℏ2as/mg=4\pi\hbar^2a_s/m is standard cold-atom theory (Pethick & Smith). — All four verifier items are executed in this article's verification script (26 checks), including the square-well fit and the np assembly; the threshold-law and 4πas24\pi a_s^2 derivations live in this article by chapter ownership, with the matching formula recalled from Phase Shifts and the (kR)2ℓ+1(kR)^{2\ell+1} statement cashed from Effective Potential and the Centrifugal Barrier.

Builds on: Partial-Wave Expansion · Phase Shifts · Asymptotic Form of the Scattering Wave Function · Scattering Amplitude Leads to: Optical Theorem · The S-Matrix and Unitarity in Scattering · Differential and Total Cross Section

Is the scattering length something you can measure directly, or only infer?

Both. In scattering you infer it from the low-energy plateau σ=4πas2\sigma=4\pi a_s^2 or the intercept of the kcot⁡δ0k\cot\delta_0 line. In neutron optics it is almost directly visible: the coherent scattering length sets a material's neutron refractive index, so total-reflection angles and interferometer phases read out scattering lengths to four significant figures — the np and nn lengths are known this precisely.

Why does one number suffice — where did all the information about the potential go?

Still there, but unresolvable. The exterior wavefunction at k→0k\to0 is a straight line no matter what; the interior can only fix its intercept. Every other detail enters through terms that vanish with powers of kRkR — the first survivor being the effective range at order k2k^2. It is multipole logic: at long distance you see the monopole first, and long wavelength is the momentum-space version of far away.

Can the scattering length be complex?

For elastic scattering off a real potential, never — δ0\delta_0 is real, so asa_s is real. With open inelastic channels (absorption, reactions) the s-wave "length" gains a negative imaginary part encoding loss, parallel to ∣Sℓ∣<1|S_\ell|<1 in Inelastic Scattering and Multichannel Theory — why tabulated lengths for cadmium, gadolinium, boron are complex.

If the singlet np channel has the bigger length, why is there no singlet deuteron?

Because as=−23.7a_s=-23.7 fm is negative: the near-threshold pole sits on the wrong side of k=0k=0 — a virtual state. Scattering magnitude feels only ∣as∣|a_s|, so the singlet channel scatters ferociously, but binding requires the pole on the positive imaginary axis. The pole taxonomy, and how experiment told the cases apart, is the business of Bound States and Zero-Energy Resonances.

Does the effective range expansion always converge?

It is an expansion of kcot⁡δ0k\cot\delta_0 in k2k^2, analytic for strictly finite-range potentials; convergence is limited by the nearest singularity — the potential's range scale, or for a Yukawa tail e−μr/re^{-\mu r}/r the branch point at k=iμ/2k=i\mu/2. Power-law tails degrade it: the van der Waals −C6/r6-C_6/r^6 tail of atomic physics allows only a modified expansion. Ask what tail your potential has before trusting the straight line far from k=0k=0.

Is a large scattering length fine-tuned or generic?

Fine-tuned: on the universal curve, ∣as∣≫R|a_s|\gg R occupies narrow windows around the divergences. Its natural occurrences are accidents — np sits near threshold by grace of QCD, helium-4 dimers by van der Waals luck. Feshbach resonances removed the luck: park the system anywhere on the curve, which is why large-asa_s universality went from curiosity to research field.

Before this

Where this leads