Scattering by a Screened Coulomb Potential

Born Approximation built the machine: to first order, the scattering amplitude is a Fourier transform of the potential, read at the momentum transfer. This article feeds the machine its most important meal. The screened Coulomb — or Yukawa — potential V(r)=β e−μr/rV(r)=\beta\,e^{-\mu r}/r is a Coulomb force wearing an exponential coat: full strength inside the screening length 1/μ1/\mu, essentially gone beyond it. It is what a charge actually looks like in nature — a nucleus dressed by its electron cloud, an ion in a plasma wrapped in its Debye sheath, a nucleon feeling another through pion exchange. Three things get done here, each computed rather than asserted. First, the Yukawa evaluation: the Fourier transform is an elementary Lorentzian, V~∝1/(q2+μ2)\tilde V\propto1/(q^2+\mu^2), and the whole angular distribution follows in four lines. Second, the Rutherford recovery, this article's centerpiece: switching the screening off (μ→0\mu\to0) reproduces Rutherford's 1911 formula exactly — the same law the classical orbit analysis gives (Rutherford Scattering & the Nuclear Atom, recalled, not rebuilt), and, by a celebrated coincidence, the same as the exact quantum Coulomb answer. Third, the divergence ledger: the total cross section comes out in closed form, σ∝1/[μ2(4k2+μ2)]\sigma\propto1/[\mu^2(4k^2+\mu^2)], and its 1/μ21/\mu^2 blow-up as screening is removed is the quantum edition of the infinite-range scandal that Differential and Total Cross Section flagged and deeded here. We then audit Geiger and Marsden's own experiment: Thomas–Fermi screening confines all deviations from Rutherford's law below θmin⁡≈10−4 rad\theta_{\min}\approx10^{-4}\,\mathrm{rad} — five hundred times smaller than any angle they measured. Every number is machine-verified (36 checks).

Key takeaways
  • The Yukawa potential V(r)=β e−μr/rV(r)=\beta\,e^{-\mu r}/r (β=Z1Z2e2/4πε0\beta=Z_1Z_2e^2/4\pi\varepsilon_0 for charges) is the prototype short-range potential: Coulomb inside the screening length 1/μ1/\mu, dead beyond it.
  • Its first Born amplitude is a pure Lorentzian, f(1)=−2mβℏ2(q2+μ2)f^{(1)}=-\dfrac{2m\beta}{\hbar^2(q^2+\mu^2)} with q=2ksin⁡(θ/2)q=2k\sin(\theta/2) — the momentum-space propagator of a mediator of mass ℏμ/c\hbar\mu/c.
  • Removing the screening at fixed angle gives Rutherford exactly: dσ/dΩ=(Z1Z2e2/16πε0E)2csc⁡4(θ/2)\mathrm{d}\sigma/\mathrm{d}\Omega=(Z_1Z_2e^2/16\pi\varepsilon_0E)^2\csc^4(\theta/2) — classical, first-Born, and exact quantum answers all coincide for 1/r1/r (the Coulomb miracle).
  • The total cross section is finite only because of screening: σ=16πm2β2ℏ4μ2(4k2+μ2)→∞\sigma=\dfrac{16\pi m^2\beta^2}{\hbar^4\mu^2(4k^2+\mu^2)}\to\infty as μ→0\mu\to0, at the rate 1/μ21/\mu^2 — every impact parameter out to 1/μ1/\mu scatters.
  • Screening audit of the real 1909–1913 experiment: for 5.5 MeV α on gold, θmin⁡=1/(kaTF)≈9×10−5\theta_{\min}=1/(ka_{\mathrm{TF}})\approx9\times10^{-5} rad; in the measured window (θ≥5°\theta\ge5°) the screening correction is below 2×10−62\times10^{-6} — Rutherford's unscreened formula was safe.
  • Honesty: at η=Z1Z2αc/v≈21\eta=Z_1Z_2\alpha c/v\approx21 the Born derivation fails for α–gold; the angular law survives anyway (miracle), but Born magnitudes like σtot\sigma_{\mathrm{tot}} do not — quoted only where η≪1\eta\ll1.
  • Key idea: the first Born amplitude is the potential's Fourier content at the momentum transfer; for a screened Coulomb force that content is a Lorentzian whose μ→0\mu\to0 limit is Rutherford's law, and whose finite μ\mu is the only thing keeping the total cross section finite.
The core insight

In the first Born approximation the scattering amplitude is the Fourier transform of the potential at the momentum transfer q=2ksin⁡(θ/2)q=2k\sin(\theta/2). For the Yukawa potential this is a simple Lorentzian, and letting the screening mass μ→0\mu\to0 analytically continues the result to the Rutherford formula.

Where this lives today

Every collider physicist reads 1/(q2+μ2)1/(q^2+\mu^2) as a propagator: Yukawa scattering is tree-level exchange of a massive scalar, the μ→0\mu\to0 limit is the photon pole, and "Rutherford behavior at small qq" remains the baseline against which parton distributions and form factors are measured. In materials analysis, Rutherford backscattering spectrometry (RBS) uses Eq. (4) daily to profile thin films — with screening corrections applied at exactly the θmin⁡\theta_{\min} level this article computes. And plasma physics runs on the regulated divergence: the Coulomb logarithm in every collision operator (Spitzer resistivity, stopping power, fusion ignition margins) is the μ→0\mu\to0 blow-up of Derivation C tamed by the Debye screen, the same ln⁡Λ≈ln⁡(λD/bmin⁡)\ln\Lambda\approx\ln(\lambda_D/b_{\min}) bookkeeping as Example 2. When dark-matter theorists posit a "dark photon", the fit function under the data is still Eq. (3)'s squared Lorentzian — with μ\mu promoted to the new particle's mass.

The physical picture

Think of the potential as a warehouse of Fourier components. The Born rule (Born Approximation's theorem, recalled): the amplitude to deflect the beam through momentum transfer ℏq\hbar\mathbf q is proportional to how much eiq⋅re^{i\mathbf q\cdot\mathbf r}-content the potential stocks. A broad potential — small μ\mu, long range — stocks mostly small-qq components: scattering hugs the forward direction. A sharp one — large μ\mu — spreads its content out to large qq: wide-angle scattering. The Yukawa warehouse inventory is exact and simple: 1/(q2+μ2)1/(q^2+\mu^2), flat for q≲μq\lesssim\mu, falling as 1/q21/q^2 beyond. That crossover at q≈μq\approx\mu is the screening length announcing itself in angle: features of the potential at range 1/μ1/\mu show up at angles θ≈μ/k\theta\approx\mu/k. The screening length is also a maximum impact parameter. A projectile passing farther than 1/μ1/\mu from the target feels essentially nothing; closer in, it sees a bare Coulomb charge. Remove the screen and arbitrarily distant passes still deflect — that is why the bare Coulomb total cross section is infinite, classically (Total Cross-Section and Scattering Limits proved the power-law-tail divergence) and quantum mechanically (this article's Derivation C, the version Differential and Total Cross Section named as emerging here). The differential cross section at any fixed angle survives the limit untouched, because a fixed angle probes a fixed qq, and q2+μ2→q2q^2+\mu^2\to q^2 smoothly. There is a third reading, and it points at quantum field theory. The Lorentzian 1/(q2+μ2)1/(q^2+\mu^2) is the static propagator of a mediating field of mass mmed=ℏμ/cm_{\mathrm{med}}=\hbar\mu/c. Yukawa ran this dictionary in reverse in 1934: given a nuclear force of range ∼1.4\sim1.4 fm, the mediator must weigh ℏc/(1.4 fm)≈141\hbar c/(1.4\,\mathrm{fm})\approx141 MeV — the pion, found in 1947. Spherical Hankel Functions and Outgoing Waves carries the position-space half of that dictionary (the e−κr/re^{-\kappa r}/r tail as an evanescent outgoing wave); here we meet its momentum-space twin, sitting inside a measurable cross section.

Figure
One potential, two costumes: real-space coat and momentum-space Lorentzian
V(r) = β e^(−μr)/r — screened in r, Lorentzian in q r bare 1/r Yukawa: e^(−μr)/r 1/μ = range q = 2k sin(θ/2) Ṽ(q) ∝ 1/(q²+μ²): flat plateau, then 1/q² q ≈ μ plateau ⇒ isotropic low-E scattering · 1/q² tail ⇒ Rutherford at q ≫ μ · propagator of mass ℏμ/c

The Yukawa amplitude in four lines

Born Approximation proved the working formula for any central potential — its owned reduction, quoted here as the starting object:

f(1)(q)=−2mℏ2q∫0∞r V(r) sin⁡(qr) dr,q=∣k′−k∣=2ksin⁡θ2(1)f^{(1)}(q)=-\frac{2m}{\hbar^2 q}\int_0^\infty r\,V(r)\,\sin(qr)\,\mathrm{d}r, \qquad q=|\mathbf k'-\mathbf k|=2k\sin\frac{\theta}{2} \qquad (1)

where: mm is the reduced mass, ℏk=2mE\hbar k=\sqrt{2mE} the incident momentum, and qq the elastic momentum transfer running from 00 (forward) to 2k2k (backscattering) — Scattering Amplitude's kinematic dictionary. Everything below is this one integral fed with V=βe−μr/rV=\beta e^{-\mu r}/r.

Derivation AThe Yukawa evaluation and its Lorentzian

We compute the first Born amplitude for the screened Coulomb potential in closed form. The strategy is the shortest in the whole scattering chapter: the 1/r1/r of the potential cancels the rr of the radial measure, leaving a textbook exponential–sine integral. What to watch: convergence of that integral requires μ>0\mu>0 — the screening is not decoration but the very thing that makes the transform exist — and the answer must depend on angle only through qq.

Step 1. Insert V=βe−μr/rV=\beta e^{-\mu r}/r into Eq. (1). The integrand becomes βe−μrsin⁡(qr)\beta e^{-\mu r}\sin(qr): the potential's 1/r1/r has been eaten by the measure's rr,

f(1)=−2mβℏ2q∫0∞e−μrsin⁡(qr) dr.f^{(1)}=-\frac{2m\beta}{\hbar^2 q}\int_0^\infty e^{-\mu r}\sin(qr)\,\mathrm{d}r.

Why the cancellation matters: it is special to 1/r1/r-type potentials and is the algebraic reason the Yukawa/Coulomb family is the one exactly Fourier-solvable case.

Step 2. Evaluate the integral: ∫0∞e−μrsin⁡(qr) dr=q/(q2+μ2)\int_0^\infty e^{-\mu r}\sin(qr)\,\mathrm{d}r=q/(q^2+\mu^2), convergent only for μ>0\mu>0. The stray qq cancels:

f(1)(θ)=−2mβℏ2 (q2+μ2)(2)f^{(1)}(\theta)=-\frac{2m\beta}{\hbar^2\,(q^2+\mu^2)} \qquad (2)

where: the denominator is the Lorentzian — the static propagator of a mediator of mass ℏμ/c\hbar\mu/c. The amplitude is real (a first-order signature; see the FAQ on the optical theorem) and negative for repulsion (β>0\beta>0).

Step 3. Square for the observable (Scattering Amplitude's current theorem supplies dσ/dΩ=∣f∣2\mathrm d\sigma/\mathrm d\Omega=|f|^2):

dσdΩ=(2mβℏ2)21[4k2sin⁡2(θ/2)+μ2]2(3)\frac{\mathrm{d}\sigma}{\mathrm{d}\Omega}=\left(\frac{2m\beta}{\hbar^2}\right)^2\frac{1}{\left[4k^2\sin^2(\theta/2)+\mu^2\right]^2} \qquad (3)

where: the angular distribution is a plateau inside θ≲μ/k\theta\lesssim\mu/k and a csc⁡4(θ/2)\csc^4(\theta/2) cliff beyond — the warehouse inventory made visible. At k≪μk\ll\mu the whole sphere sits on the plateau: isotropic scattering with f→−aBf\to-a_B, aB≡2mβ/ℏ2μ2a_B\equiv2m\beta/\hbar^2\mu^2 (the Born scattering length; the genuine article, with its bound-state subtleties, is Low-Energy Scattering and the Scattering Length's). Why dσ/dΩ\mathrm d\sigma/\mathrm d\Omega is blind to the sign of β\beta: it is even in β\beta at this order — attraction and repulsion scatter identically until second order interferes. Machine-checked. *Verified: the exponential–sine integral, Eq. (2) against the 3D-transform route −(m/2πℏ2)V~(q)-(m/2\pi\hbar^2)\tilde V(q) with V~=4πβ/(q2+μ2)\tilde V=4\pi\beta/(q^2+\mu^2), the sign-blindness, and the f(0)=−aBf(0)=-a_B plateau are all symbolic checks in this article's verifier.*

The Rutherford recovery — one limit, three theories

Now the payoff. Send μ→0\mu\to0 at fixed scattering angle, so q=2ksin⁡(θ/2)q=2k\sin(\theta/2) stays finite and simply swallows μ\mu:

Go deeper · Derivation B · 1 minSwitching off the screen: Rutherford from the Born amplitude
Derivation BSwitching off the screen: Rutherford from the Born amplitude

We take the unscreened limit of Eq. (3) and land, exactly, on the 1911 formula. The route is two substitutions; the physics is in what the limit does and does not preserve. Watch two things: the limit is taken at fixed θ>0\theta>0 — the forward direction q→0q\to0 is precisely where it fails — and the final form must expose the celebrated E−2E^{-2} and csc⁡4(θ/2)\csc^4(\theta/2) scalings.

Step 1. Drop μ2\mu^2 against q2=4k2sin⁡2(θ/2)q^2=4k^2\sin^2(\theta/2) in Eq. (3):

dσdΩ→ μ→0 (2mβℏ2q2)2=(mβ2ℏ2k2)21sin⁡4(θ/2).\frac{\mathrm{d}\sigma}{\mathrm{d}\Omega}\xrightarrow{\ \mu\to0\ }\left(\frac{2m\beta}{\hbar^2q^2}\right)^2=\left(\frac{m\beta}{2\hbar^2k^2}\right)^2\frac{1}{\sin^4(\theta/2)}.

Why fixed angle is essential: at any θ>0\theta>0, screening corrections enter as [q2/(q2+μ2)]2[q^2/(q^2+\mu^2)]^2, which →1\to1; at θ≲μ/k\theta\lesssim\mu/k they never become small — the divergence has been squeezed into a forward sliver, not abolished.

Step 2. Trade kk for the kinetic energy, E=ℏ2k2/2mE=\hbar^2k^2/2m, so 2ℏ2k2=4mE2\hbar^2k^2=4mE, and insert β=Z1Z2e2/4πε0\beta=Z_1Z_2e^2/4\pi\varepsilon_0:

dσdΩ=(β4E)21sin⁡4(θ/2)=(Z1Z2e216πε0E)2csc⁡4θ2(4)\frac{\mathrm{d}\sigma}{\mathrm{d}\Omega}=\left(\frac{\beta}{4E}\right)^2\frac{1}{\sin^4(\theta/2)}=\left(\frac{Z_1Z_2e^2}{16\pi\varepsilon_0 E}\right)^2\csc^4\frac{\theta}{2} \qquad (4)

where: every advertised scaling is on display — Z12Z22Z_1^2Z_2^2 (charge), E−2E^{-2} (energy), csc⁡4(θ/2)\csc^4(\theta/2) (angle) — and ℏ\hbar has cancelled out: the formula carries no quantum fingerprint.

Step 3. Read the three-way coincidence. The classical orbit calculation — Rutherford Scattering & the Nuclear Atom's owned result, dσ/dΩ=(κ/4E)2sin⁡−4(Θ/2)\mathrm d\sigma/\mathrm d\Omega=(\kappa/4E)^2\sin^{-4}(\Theta/2) from the hyperbolic deflection function — gives Eq. (4) verbatim. And the exact quantum Coulomb problem, solved in parabolic coordinates (Gordon, 1928; stated with source, not derived — the wavefunctions live outside this track), gives the same ∣f∣2|f|^2 once more: the exact Coulomb amplitude differs from the Born one only by a phase eiα(θ)e^{i\alpha(\theta)}, invisible in ∣f∣2|f|^2. Classical = first Born = exact. The disappearance of ℏ\hbar in Step 2 is the trace of this conspiracy, and the Kepler/Coulomb hidden symmetry (the conserved Runge–Lenz vector of the 1/r1/r problem) is its structural root. For any other power law the three answers disagree. Verified: Eq. (4) from Eq. (3) symbolically; the classical formula's identity is quoted against the mechanics article's served statement (gate table); the phase-only discrepancy of the exact amplitude is stated with source.

The coincidence has a sharp practical edge. For 5.5 MeV α particles on gold the Born expansion parameter is the Sommerfeld number η=Z1Z2α c/v≈21\eta=Z_1Z_2\alpha\,c/v\approx21 — the "weak potential" premise fails by a factor of twenty, and honest error analysis says the first Born term should be garbage. It is not, for this one potential: the angular law is protected by the miracle. The protection does not extend to quantities the miracle doesn't cover — Born magnitudes like the screened total cross section below are trustworthy only when η≪1\eta\ll1 (Born Approximation owns the criteria; here we only apply them). That double status — right answer, illegitimate derivation — is why Rutherford's classical 1911 analysis survived quantum mechanics unscathed.

The divergence ledger — what screening buys

Go deeper · Derivation C · 1 minThe total cross section and its inverse-square blow-up
Derivation CThe total cross section and its inverse-square blow-up

We integrate Eq. (3) over the sphere in closed form, then watch the result diverge as the screen is removed. The goal is a ledger: exactly how fast σ\sigma blows up, where the divergent contribution lives in angle, and what physical screening scale regulates real experiments. The trick is the standard substitution u=q2u=q^2, which turns the angular integral into a rational one.

Step 1. With u=q2=4k2sin⁡2(θ/2)u=q^2=4k^2\sin^2(\theta/2), du=2k2sin⁡θ dθ\mathrm du=2k^2\sin\theta\,\mathrm d\theta (jacobian machine-checked), and uu running 0→4k20\to4k^2:

σ=2π(2mβℏ2)212k2∫04k2du(u+μ2)2.\sigma=2\pi\left(\frac{2m\beta}{\hbar^2}\right)^2\frac{1}{2k^2}\int_0^{4k^2}\frac{\mathrm{d}u}{(u+\mu^2)^2}.

Step 2. The integral is elementary, ∫04k2(u+μ2)−2du=4k2/[μ2(4k2+μ2)]\int_0^{4k^2}(u+\mu^2)^{-2}\mathrm du=4k^2/[\mu^2(4k^2+\mu^2)], giving the closed form

σ=16πm2β2ℏ4 μ2 (4k2+μ2)(5)\sigma=\frac{16\pi m^2\beta^2}{\hbar^4\,\mu^2\,(4k^2+\mu^2)} \qquad (5)

where: at k→0k\to0 this is 4πaB24\pi a_B^2 (isotropic plateau over the whole sphere); at k≫μk\gg\mu it falls as 1/k21/k^2 — faster beams see a smaller effective disk.

Step 3. The blow-up. As μ→0\mu\to0 at fixed kk, Eq. (5) diverges as σ→4πm2β2/(ℏ4k2μ2)\sigma\to4\pi m^2\beta^2/(\hbar^4k^2\mu^2) — an inverse-square divergence in the screening mass (coefficient machine-verified). The geometry behind the algebra: the integrand's divergent weight sits at u≲μ2u\lesssim\mu^2, i.e. angles θ≲θmin⁡\theta\lesssim\theta_{\min} with

θmin⁡≃μk=1k ascr(6)\theta_{\min}\simeq\frac{\mu}{k}=\frac{1}{k\,a_{\mathrm{scr}}} \qquad (6)

where: ascr=1/μa_{\mathrm{scr}}=1/\mu is the screening length. This is the semiclassical statement "impact parameters beyond ascra_{\mathrm{scr}} don't scatter" translated into angle — the quantum edition of the classical regularization Total Cross-Section and Scattering Limits built with impact-parameter cutoffs. At θ=θmin⁡\theta=\theta_{\min} the screening suppression factor [q2/(q2+μ2)]2[q^2/(q^2+\mu^2)]^2 is exactly 1/41/4 (machine-checked): the angle where Rutherford's law has lost half its amplitude. Why the divergence is physical, not a Born artifact: the exact Coulomb problem has no finite σ\sigma either — infinite range scatters every impact parameter. The ledger's message is that "total cross section of a bare charge" is not a measurable; what nature measures is always screened. *Verified: Eq. (5) symbolically by the u-substitution AND numerically against direct solid-angle quadrature to 10−1010^{-10}; the k→0k\to0 and μ→0\mu\to0 coefficients; the factor-1/4 at θmin⁡\theta_{\min}.*

The two real screens. In matter, the nuclear charge is screened by atomic electrons at the Thomas–Fermi radius aTF≈0.885 a0Z−1/3a_{\mathrm{TF}}\approx0.885\,a_0Z^{-1/3} — for gold, 0.110.11 Å. In a plasma or electrolyte, mobile charges build the Debye cloud, ascr=λDa_{\mathrm{scr}}=\lambda_D — fractions of a millimetre in fusion-class plasmas. Both plug straight into Eqs. (5)–(6); the punchline figure runs the first one against history.

Figure
Punchline — auditing Geiger and Marsden: screening correction R(θ) vs angle, with the 1909–13 data window
R(θ) = dσ_screened / dσ_Rutherford = [q²/(q²+μ²)]² — 5.5 MeV α on gold, μ = 1/a_TF θ (rad, log scale: 10⁻⁵ → π) R(θ) 10⁻⁵ 10⁻³ 5° π R = 1 (pure Rutherford) θ_min = 1/(k·a_TF) = 9.1×10⁻⁵ rad — R = 1/4 here Geiger–Marsden window: 5° → 150° screening correction below 2×10⁻⁶ here all of σ's divergence lives below θ_min — invisible to any detector the 1911 analysis never needed screening: its data sat 500× above θ_min, where the bare-Coulomb law is exact to parts per million

Gaussian vs SI: one symbol, two dialects

Scattering literature runs on Gaussian units: there β=Z1Z2e2\beta=Z_1Z_2e^2 and the Rutherford prefactor reads (Z1Z2e2/4E)2(Z_1Z_2e^2/4E)^2; in SI, β=Z1Z2e2/4πε0\beta=Z_1Z_2e^2/4\pi\varepsilon_0 and the prefactor is (Z1Z2e2/16πε0E)2(Z_1Z_2e^2/16\pi\varepsilon_0E)^2. Both are the same (β/4E)2(\beta/4E)^2. The practical bridge is e2/4πε0=αℏc=1.440e^2/4\pi\varepsilon_0=\alpha\hbar c=1.440 MeV·fm, which turns every formula in this article into calculator arithmetic — it is how both worked examples below run. The Sommerfeld parameter is η=Z1Z2 α c/v\eta=Z_1Z_2\,\alpha\,c/v in either dialect.

Limiting cases

caseresultmeaning
μ→0\mu\to0, fixed θ>0\theta>0dσ/dΩ→(β/4E)2csc⁡4(θ/2)\mathrm d\sigma/\mathrm d\Omega\to(\beta/4E)^2\csc^4(\theta/2)Rutherford recovered exactly; classical = Born = exact (1/r1/r miracle)
μ→0\mu\to0, totalσ≃4πm2β2/(ℏ4k2μ2)→∞\sigma\simeq4\pi m^2\beta^2/(\hbar^4k^2\mu^2)\to\inftyinverse-square blow-up; the quantum face of the infinite-range scandal
θ→0\theta\to0, μ>0\mu>0f(0)=−aB=−2mβ/ℏ2μ2f(0)=-a_B=-2m\beta/\hbar^2\mu^2screening caps the forward peak at the Born scattering length
k≪μk\ll\muf→−aBf\to-a_B, σ→4πaB2\sigma\to4\pi a_B^2isotropic s-wave-like plateau; the genuine low-energy theory is Low-Energy Scattering and the Scattering Length's
q≫μq\gg\mudσ/dΩ→(2mβ/ℏ2q2)2\mathrm d\sigma/\mathrm d\Omega\to(2m\beta/\hbar^2q^2)^2large-angle probes see the bare core — screening invisible past θmin⁡\theta_{\min}
η=Z1Z2αc/v≳1\eta=Z_1Z_2\alpha c/v\gtrsim1angular law survives, magnitudes don'tthe miracle protects csc⁡4\csc^4, not Born normalizations like Eq. (5)

Worked examples

Example 1 — Rutherford's actual experiment, end to end. A 10610^6 α/s beam of 5.5 MeV α particles hits a 0.5 μm gold foil; a 1 cm² detector sits 10 cm away at θ=60°\theta=60°. Find the count rate, and audit whether screening could have distorted the 1909–13 data.

Step 1. Kinematics. Reduced mass m=mαMAu/(mα+MAu)=3653m=m_\alpha M_{\mathrm{Au}}/(m_\alpha+M_{\mathrm{Au}})=3653 MeV/c2/c^2; CM energy E=5.5×MAu/(MAu+mα)=5.39E=5.5\times M_{\mathrm{Au}}/(M_{\mathrm{Au}}+m_\alpha)=5.39 MeV; k=2mE/ℏc=1.006k=\sqrt{2mE}/\hbar c=1.006 fm⁻¹.

Step 2. The prefactor. β=Z1Z2 αℏc=2×79×1.440=227.5\beta=Z_1Z_2\,\alpha\hbar c=2\times79\times1.440=227.5 MeV·fm, so (β/4E)2=(227.5/21.56)2=111.3(\beta/4E)^2=(227.5/21.56)^2=111.3 fm². At 60°60°: csc⁡4(30°)=16\csc^4(30°)=16, so dσ/dΩ=1781 fm2/sr=\mathrm d\sigma/\mathrm d\Omega=1781\ \mathrm{fm^2/sr}= 17.8 b/sr.

Step 3. Count rate. Gold: n=ρNA/A=5.91×1022n=\rho N_A/A=5.91\times10^{22} cm⁻³, so the foil presents nt=2.95×1018nt=2.95\times10^{18} nuclei/cm². Detector solid angle ΔΩ=1/102=10−2\Delta\Omega=1/10^2=10^{-2} sr. Rate =106×2.95×1018×17.8×10−24×10−2==10^6\times2.95\times10^{18}\times17.8\times10^{-24}\times10^{-2}= 0.53 counts/s — patient scintillation-counting territory, which is exactly what Geiger's dark-adapted eyes endured.

Step 4. The screening audit. Thomas–Fermi radius of gold aTF=0.885 a0×79−1/3=1.09×104a_{\mathrm{TF}}=0.885\,a_0\times79^{-1/3}=1.09\times10^4 fm, so θmin⁡=1/(kaTF)=9.1×10−5\theta_{\min}=1/(ka_{\mathrm{TF}})=9.1\times10^{-5} rad. The smallest angle in the data, 5°=0.0875°=0.087 rad, sits ∼103 θmin⁡\sim10^3\,\theta_{\min} out: the correction factor is 1−R<2×10−61-R<2\times10^{-6}. Integrated over the window θ>5°\theta>5°, σ>5°=4π(β/4E)2[csc⁡2(2.5°)−1]=\sigma_{>5°}=4\pi(\beta/4E)^2[\csc^2(2.5°)-1]= 7.3 kb — finite, screening-free Rutherford physics. The unscreened formula was safe to parts per million.

Example 2 — where the Born magnitude is honest: a Debye-screened plasma. A 1 keV electron scatters off a proton in a fusion-class plasma (Te=1T_e=1 keV, ne=1014n_e=10^{14} cm⁻³). Here v/c=0.0626v/c=0.0626 and η=αc/v=0.117≪1\eta=\alpha c/v=0.117\ll1: first Born is legitimate, magnitudes included.

Step 1. The screen. λD=7430TeV/ncm−3 cm=0.235\lambda_D=7430\sqrt{T_{\mathrm{eV}}/n_{\mathrm{cm^{-3}}}}\ \mathrm{cm}=0.235 mm — enormous compared with atomic scales, so μ=1/λD\mu=1/\lambda_D is tiny and 4k2≫μ24k^2\gg\mu^2 by 30 orders of magnitude.

Step 2. Total cross section, Eq. (5) in its 4k2≫μ24k^2\gg\mu^2 form: with mβ/ℏ2=1/a0m\beta/\hbar^2=1/a_0 (hydrogen coupling), μ=1/λD\mu=1/\lambda_D and k=1.62×1011k=1.62\times10^{11} m⁻¹, it collapses to σ=4πλD2/(a0k)2\sigma=4\pi\lambda_D^2/(a_0k)^2. Numbers: σ=\sigma= 9.4×10−99.4\times10^{-9} m² =0.17 λD2=0.17\,\lambda_D^2 (exact closed form matches to 10−1510^{-15}).

Step 3. Read it honestly. The cross section is a sixth of the entire Debye disk — because every electron crossing the sheath is deflected a little, and σ\sigma counts any deflection. The forward peak is astonishing: f(0)=2λD2/a0≈2f(0)=2\lambda_D^2/a_0\approx2 km(!), confined to a cone θc=1/(kλD)=2.6×10−8\theta_c=1/(k\lambda_D)=2.6\times10^{-8} rad; the narrow-cone bookkeeping ∣f(0)∣2πθc2|f(0)|^2\pi\theta_c^2 reproduces σ\sigma (machine-checked). Transport theory therefore never uses σ\sigma raw: it weights by (1−cos⁡θ)(1-\cos\theta), converting the 1/μ21/\mu^2 sensitivity into the tame Coulomb logarithm ln⁡(λD/bmin⁡)\ln(\lambda_D/b_{\min}) — the reason plasma resistivity is finite. σ = 9.4×10⁻⁹ m², dominated by micro-deflections; momentum transport is what's physical.

Origins · 1 minYukawa Potential and how this was discovered
YP
Yukawa Potential
1923–1949 · Zurich, Osaka

The functional form came first, the particle later. In 1923 Peter Debye and Erich Hückel, computing how an ion's field survives inside an electrolyte, found the exponentially-coated Coulomb potential e−r/λD/re^{-r/\lambda_D}/r — pure classical statistical mechanics, no quanta anywhere. Llewellyn Thomas (1926) and Enrico Fermi (1927) rediscovered the coat inside the atom: the electron cloud screens the nucleus at the radius ∼a0Z−1/3\sim a_0Z^{-1/3} that this article's Geiger–Marsden audit uses. Hideki Yukawa's 1934 leap (published February 1935) was to run the dictionary backwards: if the nuclear force has range 1.4 fm not because something screens it but because the exchanged quantum itself is massive, then the mass must be ℏμc≈\hbar\mu c\approx 100–200 MeV. The pion's discovery in 1947 made that reading literal, and the 1949 Nobel Prize — Japan's first — sealed it. The same Lorentzian 1/(q2+μ2)1/(q^2+\mu^2) this article extracts in the Born approximation is, in quantum field theory, the propagator of the exchanged boson: this scattering problem is where students first touch a Feynman-diagram denominator. (Yukawa's person-portrait, with the position-space h0(1)(iκr)h_0^{(1)}(i\kappa r) tail dictionary, hangs in Spherical Hankel Functions and Outgoing Waves; this card is the potential's own biography.)

History
1911
— Rutherford derives the csc⁴(θ/2) law classically and infers the nuclear atom from Geiger–Marsden data
1923
— Debye and Hückel derive the screened Coulomb form for electrolytes: the potential exists before the meson
1926
— Born's collision papers make the amplitude a Fourier transform; Thomas (1926) and Fermi (1927) find atomic screening
1928
— Gordon solves the exact quantum Coulomb problem: same cross section as classical and Born — the miracle certified
1935
— Yukawa publishes the massive-mediator reading of the short-range nuclear force, predicting the ~100–200 MeV meson
1947
— the pion is discovered in cosmic-ray emulsions; 1949 — Yukawa's Nobel Prize

Going further

References. Griffiths, Introduction to Quantum Mechanics, 3rd ed., Ch. 11 (Scattering) · Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., Ch. 6 · Shankar, Principles of Quantum Mechanics, 2nd ed., Ch. 19 · Cohen-Tannoudji, Diu & Laloë, Quantum Mechanics, Vol. 2, Ch. VIII · Townsend, A Modern Approach to Quantum Mechanics, 2nd ed., Ch. 13.

Provenance. Equations derived from the standard central-potential Born formula; the exponential–sine integral ∫0∞e−μrsin⁡(qr)dr=q/(q2+μ2)\int_0^\infty e^{-\mu r}\sin(qr)\mathrm{d}r=q/(q^2+\mu^2) is a table result. A verifier should reproduce f(1)=−2mβ/[ℏ2(q2+μ2)]f^{(1)}=-2m\beta/[\hbar^2(q^2+\mu^2)] (Griffiths Ch.11, Sakurai §6), the total cross section integral σ=16πm2β2/[ℏ4μ2(4k2+μ2)]\sigma=16\pi m^2\beta^2/[\hbar^4\mu^2(4k^2+\mu^2)], and the μ→0\mu\to0 Rutherford limit ( Z1Z2e2/16πε0E)2csc⁡4(θ/2)(\,Z_1Z_2 e^2/16\pi\varepsilon_0 E)^2\csc^4(\theta/2). Cross-check angular kinematics q=2ksin⁡(θ/2)q=2k\sin(\theta/2) and the equality of Born and exact Rutherford differential cross sections against Sakurai, Shankar Ch.19, and Cohen-Tannoudji Vol.2 Ch.VIII. — All three reproductions, both integral routes, and the kinematics are executed in this article's verify script (36 checks); the classical-formula identity is quoted against the mechanics article's served text.

Builds on: Born Approximation · Scattering Amplitude · Differential and Total Cross Section · Lippmann-Schwinger Equation Leads to: Born Series and Higher-Order Corrections · Optical Theorem · Partial-Wave Expansion · Phase Shifts

Is the screened potential just a mathematical regulator, or physically real?

Physically real, in three independent costumes. In a plasma or electrolyte the Debye cloud is a measurable charge rearrangement, and λD\lambda_D enters wave dispersion and transport data. In an atom the Thomas–Fermi cloud is the actual electron density; electron-microscopy form factors measure exactly the screened field. In nuclear physics the "screen" is not a medium at all but the mediator's own mass — the pion's 140 MeV. The bare Coulomb potential is the idealization; nature always hands you some μ>0\mu>0.

If the Born series fails for Coulomb-strength coupling, why do we trust anything on this page for α on gold?

Only one output is trusted there: the angular law, Eq. (4), which the Coulomb miracle protects — the exact Coulomb ∣f∣2|f|^2 happens to equal the Born ∣f∣2|f|^2. Everything else (the total cross section's magnitude, the forward amplitude, the phase of ff) is genuinely first-order-only and needs η=Z1Z2αc/v≪1\eta=Z_1Z_2\alpha c/v\ll1. The article's warning box quantifies the failure: Eq. (5) at η=21\eta=21 overshoots the geometric Thomas–Fermi disk by a factor ~1800. Where the series and its convergence actually stand is Born Series and Higher-Order Corrections' business.

Why does the amplitude come out real, and doesn't that violate the optical theorem?

Real central potentials give real first-Born amplitudes — Eq. (2) is manifestly real. The optical theorem demands Imf(0)=kσtot/4π>0\mathrm{Im}f(0)=k\sigma_{\mathrm{tot}}/4\pi>0, so first order violates it at order β2\beta^2: the deficit is repaired order by order (the second Born term's forward imaginary part supplies exactly the missing kσ(1)/4πk\sigma^{(1)}/4\pi — Born Series and Higher-Order Corrections proves the ladder, Optical Theorem owns the exact statement). A perturbative truncation conserves probability only to the order kept.

Where exactly did ℏ go in the Rutherford formula?

It cancelled between the amplitude's m/ℏ2m/\hbar^2 and the momentum transfer's k2=2mE/ℏ2k^2=2mE/\hbar^2. That cancellation is the algebraic shadow of the classical–quantum coincidence: a formula with no ℏ\hbar can agree with classical mechanics, and this one does. Don't over-read it — the underlying amplitudes differ by phases (measurable in interference, e.g. identical-particle Mott scattering), and for any non-1/r1/r potential the ℏ\hbar's survive and the answers split.

Is the Born scattering length in Eq. (3) the same object as the real scattering length?

Same definition (f→−af\to-a as k→0k\to0), different reliability. aB=2mβ/ℏ2μ2a_B=2m\beta/\hbar^2\mu^2 is linear in the potential: it can never show the sign flips, divergences, or zero crossings that the true aa exhibits as a well deepens toward binding a state — nonperturbative physics owned by Low-Energy Scattering and the Scattering Length and Bound States and Zero-Energy Resonances. For a weakly-coupled Yukawa (2m∣β∣/ℏ2μ≪12m|\beta|/\hbar^2\mu\ll1) they agree; for the deuteron-strength nuclear well they do not even share a sign convention's worth of accuracy.

Why does the divergent σ not make plasmas opaque?

Because σ counts deflections of any size, and almost all of the Debye-disk-sized cross section of Example 2 is micro-deflections of order 10−810^{-8} rad. Transport — resistivity, stopping, thermalization — weights each collision by its momentum transfer, (1−cos⁡θ)(1-\cos\theta), and that integral is only logarithmically sensitive to the screen: the Coulomb logarithm ln⁡Λ∼ln⁡(λD/bmin⁡)≈10 ⁣− ⁣20\ln\Lambda\sim\ln(\lambda_D/b_{\min})\approx10\!-\!20 in fusion plasmas. A plasma is nearly transparent to any single large-angle event yet collisional in aggregate.

The exact Coulomb amplitude differs from Born only by a phase — is that phase ever measurable?

Yes, twice over. Scatter identical charged particles and the direct and exchange amplitudes f(θ)f(\theta), f(π−θ)f(\pi-\theta) interfere: the Coulomb phase produces the Mott oscillations that distinguish quantum from classical Rutherford scattering (owned by the identical-particle article of the next chapter). And in any multi-channel situation (Coulomb + nuclear force at the same vertex) the relative phase between the two amplitudes shifts interference minima — precisely how low-energy nuclear phase shifts are extracted beneath a Coulomb foreground, Phase Shifts' territory.

Before this

Where this leads