Skip to content

Scattering Theory

12.1 Partial Wave Analysis

For a spherically symmetric potential V(r)V(r)The scattering amplitude can be expanded in partial waves:

f(θ)=12ikl=0(2l+1)(e2iδl1)Pl(cosθ)=1kl=0(2l+1)eiδlsinδlPl(cosθ)f(\theta) = \frac{1}{2ik}\sum_{l=0}^{\infty}(2l + 1)(e^{2i\delta_l} - 1)P_l(\cos\theta) = \frac{1}{k}\sum_{l=0}^{\infty}(2l + 1)e^{i\delta_l}\sin\delta_l\,P_l(\cos\theta)

Where δl\delta_l is the phase shift for partial wave ll.

Optical theorem:

σtotal=4πkImf(0)=4πk2l=0(2l+1)sin2δl\sigma_{\text{total} = \frac{4\pi}{k}\,\text{Im}\,f(0) = \frac{4\pi}{k^2}\sum_{l=0}^{\infty}(2l+1)\sin^2\delta_l}

Partial wave unitarity bound: sin2δl1\sin^2\delta_l \leq 1So the maximum contribution of partial wave ll to the cross section is:

σlmax=4πk2(2l+1)\sigma_l^{\max} = \frac{4\pi}{k^2}(2l + 1)

12.2 The Born Approximation

For a weak potential, the scattering amplitude to first order is:

f(θ,ϕ)m2π2eik"rV(r)d3rf(\theta, \phi) \approx -\frac{m}{2\pi\hbar^2}\int e^{-i\mathbf{k}"\cdot\mathbf{r}'}V(\mathbf{r}')\,d^3r'

Where k\mathbf{k}' is the scattered wave vector and q=kk\mathbf{q} = \mathbf{k}' - \mathbf{k} is the momentum transfer.

For the Yukawa potential V(r)=(V0/r)eμrV(r) = (V_0/r)e^{-\mu r}:

f(θ)=2mV02(μ2+q2),q=2ksin(θ/2)f(\theta) = -\frac{2m V_0}{\hbar^2(\mu^2 + q^2)}, \quad q = 2k\sin(\theta/2)

Setting μ=0\mu = 0 (Coulomb potential), this reproduces the Rutherford scattering formula.

12.3 Resonance Scattering

When the scattering energy is near a quasi-bound state, the phase shift passes through π/2\pi/2 (resonance):

δl(E)δbg+arctan ⁣(Γ/2ERE)\delta_l(E) \approx \delta_{\text{bg} + \arctan\!\left(\frac{\Gamma/2}{E_R - E}\right)}

Where ERE_R is the resonance energy and Γ\Gamma is the width. The cross section has the Breit—Wigner form:

σl(E)=4πk2(2l+1)(Γ/2)2(EER)2+(Γ/2)2\sigma_l(E) = \frac{4\pi}{k^2}(2l+1)\frac{(\Gamma/2)^2}{(E - E_R)^2 + (\Gamma/2)^2}

At resonance (E=ERE = E_R): σlmax=4πk2(2l+1)\sigma_l^{\max} = \frac{4\pi}{k^2}(2l+1) (unitarity limit).

Worked Example 12.1: Low-Energy Scattering and Scattering Length

For ss-wave scattering (l=0l = 0) at low energy (ka1ka \ll 1), only the l=0l = 0 phase shift contributes:

σ4πk2sin2δ04πas2\sigma \approx \frac{4\pi}{k^2}\sin^2\delta_0 \approx 4\pi a_s^2

Where the scattering length asa_s is defined by kcotδ01/ask\cot\delta_0 \to -1/a_s as k0k \to 0.

For a hard sphere of radius aa: δ0=ka\delta_0 = -ka (exact), giving as=aa_s = a and σ=4πa2\sigma = 4\pi a^2 (four times the geometric cross section πa2\pi a^2 --- a purely quantum result).

For the 3^3He—4^4He system: as1.4a_s \approx 1.4 nm (positive, indicating a repulsive effective potential). For the neutron—proton system (triplet): as5.4a_s \approx 5.4 fm (positive, with a bound state --- the deuteron). For singlet: as23.7a_s \approx -23.7 fm (negative, indicating a virtual state).

Worked Example 12.2: Born Approximation for a Gaussian Potential

Consider V(r)=V0er2/(2a2)V(r) = V_0\,e^{-r^2/(2a^2)}.

f(θ)=m2π2eiqrV0er2/(2a2)d3rf(\theta) = -\frac{m}{2\pi\hbar^2}\int e^{-i\mathbf{q}\cdot\mathbf{r}}V_0 e^{-r^2/(2a^2)}\,d^3r

=mV02π2(2πa2)3/2eq2a2/2=mV02(2π)1/2a3e2k2a2sin2(θ/2)= -\frac{m V_0}{2\pi\hbar^2}(2\pi a^2)^{3/2}e^{-q^2 a^2/2} = -\frac{m V_0}{\hbar^2}(2\pi)^{1/2}a^3\,e^{-2k^2a^2\sin^2(\theta/2)}

The total cross section:

σ=f2dΩ=2π0πf2sinθdθ\sigma = \int |f|^2\,d\Omega = 2\pi\int_0^\pi |f|^2\sin\theta\,d\theta

At low energy (ka1ka \ll 1): fmV02(2π)1/2a3f \approx -\frac{mV_0}{\hbar^2}(2\pi)^{1/2}a^3 (independent of θ\theta), giving:

σ4π(mV02)22πa6=8π2m2V02a64\sigma \approx 4\pi\left(\frac{mV_0}{\hbar^2}\right)^2 2\pi\,a^6 = \frac{8\pi^2 m^2 V_0^2 a^6}{\hbar^4}

The Born approximation is valid when V02/(ma2)|V_0| \ll \hbar^2/(ma^2)I.e., the potential is weak compared to the kinetic energy associated with the length scale aa.

Common Pitfalls (Additional)

  1. Symmetrisation applies to the full wavefunction: For fermions, the overall wavefunction (spatial \otimes spin \otimes any other degrees of freedom) must be antisymmetric. A symmetric spatial part requires an antisymmetric spin part (singlet), and vice versa. Do not apply (anti)symmetrisation to spatial and spin parts separately without ensuring the correct combined symmetry.

  2. The variational principle gives an upper bound: EtrialE0E_{\text{trial} \geq E_0} always. If you obtain a variational energy lower than the known exact ground state energy, you have made an error in the calculation (wrong normalisation, incorrect matrix element, or the trial function is not in the correct Hilbert space).

  3. Fermi’s Golden Rule applies to transitions to a continuum: For transitions to discrete states, use the Rabi formula instead. The density of states ρ(Ef)\rho(E_f) is essential --- if it is zero, the transition rate is zero regardless of the matrix element.

  4. The Born approximation assumes a weak potential: The condition is mV0a/21m|V_0|a/\hbar^2 \ll 1 where aa is the range of the potential. For strong potentials (like the nuclear potential or hard spheres), the Born approximation gives qualitatively wrong results. Use partial wave analysis instead.

  5. Resonances require careful treatment: Near a resonance, perturbation theory breaks down. The Breit—Wigner formula is non-perturbative in the width Γ\Gamma. The scattering length can be much larger than the range of the potential near a resonance (the unitarity limit).

Problems (Additional)

Problem 19: Exchange Energy in Lithium

Lithium (Z=3Z = 3) has the electron configuration 1s22s11s^2 2s^1. Using the variational method with ZeffZ_{\text{eff}} for the 1s1s electrons:

(a) Calculate ZeffZ_{\text{eff}} for the 1s1s electrons, treating the 2s2s electron as a perturbation.

(b) Calculate the ionisation energy (removing the 2s2s electron) and compare with the experimental value of 5.39 eV.

(c) Explain why the 2s2s electron is effectively screened by Zeff1.26Z_{\text{eff} \approx 1.26}.

Solution:

(a) For the 1s1s electrons, the effective charge is reduced from Z=3Z = 3 by screening from the other 1s1s electron and partial penetration of the 2s2s electron. The 1s1s electrons screen each other partially: using the helium result, Zeff(1s)Z5/16=30.3125=2.69Z_{\text{eff}(1s) \approx Z - 5/16 = 3 - 0.3125 = 2.69}.

(b) The 2s2s electron sees an effective nuclear charge of Zeff(2s)32×0.85=1.3Z_{\text{eff}(2s) \approx 3 - 2 \times 0.85 = 1.3} (Slater’s rules). The energy:

E_{2s} = -\frac{Z_{\text{eff}^2}{n^2}\times 13.6\ \text{eV} = -\frac{1.3^2}{4}\times 13.6 = -\frac{1.69}{4}\times 13.6 = -5.75\ \text{eV}}

The ionisation energy is E2s=5.75|E_{2s}| = 5.75 eV, close to the experimental 5.39 eV. The discrepancy reflects the crudeness of the Slater screening constants.

(c) The 2s2s electron has significant radial extent beyond the 1s1s core, so it sees a nearly bare nuclear charge at small rr but is screened by both 1s1s electrons at large rr. The effective charge Zeff1.26Z_{\text{eff} \approx 1.26} (using Hartree—Fock) represents this average screening.

Problem 20: Partial Wave Analysis for Square Well

Consider scattering from the attractive square well V(r)=V0V(r) = -V_0 for r<ar < a and V(r)=0V(r) = 0 for r>ar > a.

(a) Show that the ss-wave phase shift satisfies:

δ0=ka+arctan ⁣(kκtan(κa))\delta_0 = -ka + \arctan\!\left(\frac{k}{\kappa}\tan(\kappa a)\right)

Where κ=2m(V0+E)/\kappa = \sqrt{2m(V_0 + E)}/\hbar and k=2mE/k = \sqrt{2mE}/\hbar.

(b) Show that a bound state exists at energy E=EE = -|E| when κ0a=π/2\kappa_0 a = \pi/2 where κ0=2m(V0E)/\kappa_0 = \sqrt{2m(V_0 - |E|)}/\hbar.

(c) Show that the scattering length diverges as a new bound state appears.

Solution:

(a) Inside the well (r<ar < a), the radial wavefunction for l=0l = 0 is u(r)=Asin(κr)u(r) = A\sin(\kappa r). Outside (r>ar > a), u(r)=Bsin(kr+δ0)u(r) = B\sin(kr + \delta_0).

Matching uu and uu' at r=ar = a:

κcos(κa)=kcos(ka+δ0)/sin(ka+δ0)k\kappa\cos(\kappa a) = k\cos(ka + \delta_0)/\sin(ka + \delta_0)\cdot k

Wait: κcot(κa)=kcot(ka+δ0)\kappa\cot(\kappa a) = k\cot(ka + \delta_0).

cot(ka+δ0)=κkcot(κa)\cot(ka + \delta_0) = \frac{\kappa}{k}\cot(\kappa a)

ka+δ0=\arccot ⁣(κkcot(κa))=arctan ⁣(kκtan(κa))ka + \delta_0 = \arccot\!\left(\frac{\kappa}{k}\cot(\kappa a)\right) = \arctan\!\left(\frac{k}{\kappa}\tan(\kappa a)\right)

δ0=ka+arctan ⁣(kκtan(κa))\delta_0 = -ka + \arctan\!\left(\frac{k}{\kappa}\tan(\kappa a)\right)

(b) A bound state has E<0E < 0So k=iκk = i\kappa' where κ=2mE/\kappa' = \sqrt{2m|E|}/\hbar. The bound state condition is that the exterior solution decays exponentially: u(r)=Beκru(r) = Be^{-\kappa' r}. Matching:

κcot(κa)=κ\kappa\cot(\kappa a) = -\kappa'

As E0|E| \to 0: κ0\kappa' \to 0So κcot(κa)0\kappa\cot(\kappa a) \to 0Giving κa=π/2\kappa a = \pi/2 (the threshold for the first bound state).

(c) The scattering length as=limk0(δ0/k)a_s = -\lim_{k \to 0}(\delta_0/k). As κaπ/2\kappa a \to \pi/2, tan(κa)\tan(\kappa a) \to \inftySo:

δ0ka+arctan()=ka+π/2\delta_0 \approx -ka + \arctan(\infty) = -ka + \pi/2

as=1k(ka+π2)=aπ2ka_s = -\frac{1}{k}\left(-ka + \frac{\pi}{2}\right) = a - \frac{\pi}{2k}

As k0k \to 0: as±a_s \to \pm\infty (diverges), changing sign as the bound state appears.