Partial Wave Expansion of Plane Waves

This article closes the chapter by writing down its Rosetta stone. Two complete descriptions of a free particle now coexist: the plane wave eikze^{ikz} — the beam every accelerator and every scattering theory starts from — and the spherical channels jℓ(kr)Pℓ(cos⁡θ)j_\ell(kr)P_\ell(\cos\theta) built over the last four serials. The translation between them is a single exact identity, the Rayleigh expansion:

eikz=∑ℓ=0∞(2ℓ+1) iℓ jℓ(kr) Pℓ(cos⁡θ)e^{ikz}=\sum_{\ell=0}^{\infty}(2\ell+1)\,i^\ell\,j_\ell(kr)\,P_\ell(\cos\theta)

Every partial-wave computation in quantum mechanics — every phase shift, every cross section, every "the s-wave dominates" argument — begins by quoting this line. Here it is earned: derived twice, by two genuinely independent routes (orthogonality projection, and leading-power matching), with the coefficients (2ℓ+1)iℓ(2\ell+1)i^\ell verified numerically eighteen ways and the general-direction form eik⋅r=4π∑ℓmiℓjℓ(kr)Yℓm∗(k^)Yℓm(r^)e^{i\mathbf{k}\cdot\mathbf{r}}=4\pi\sum_{\ell m}i^\ell j_\ell(kr)Y_\ell^{m*}(\hat{k})Y_\ell^m(\hat{r}) checked to machine precision at random orientations. Beyond the identity itself, the article owns its two readings. The budget reading: how much of a beam lives in each ℓ-channel — the weights (2ℓ+1)jℓ2(kr)(2\ell+1)j_\ell^2(kr) sum exactly to one, 90% of a kr=10kr=10 beam sits below ℓ≈10\ell\approx10, and the semiclassical dictionary ℓ+12=kb\ell+\tfrac12=kb turns channel number into impact parameter. The ledger reading: split each channel's far field into incoming and outgoing spherical waves and the plane wave is revealed as a perfectly balanced account — equal arrivals and departures in every channel. That balance is the zero point of scattering theory: a potential can only rephase the outgoing column, and the machinery that exploits this — amplitudes, phase shifts, the S-matrix — is named here and owned by the scattering chapter. With this dictionary banked, the central-potentials chapter is complete: twenty serials, from separation to translation.

Key takeaways
  • THE identity: eikz=∑ℓ(2ℓ+1)iℓjℓ(kr)Pℓ(cos⁡θ)e^{ikz}=\sum_\ell(2\ell+1)i^\ell j_\ell(kr)P_\ell(\cos\theta) — exact, all rr, all θ\theta; coefficients derived twice and verified numerically for ℓ≤5\ell\le5 at three radii (18 checks).
  • Route A (projection): Legendre orthogonality peels off channel ℓ; the radial integral is (2ℓ+1)iℓjℓ(2\ell+1)i^\ell j_\ell — no asymptotics anywhere.
  • Route B (leading powers): match coefficients of xℓμℓx^\ell\mu^\ell using jℓ→xℓ/(2ℓ+1)!!j_\ell\to x^\ell/(2\ell+1)!! and PℓP_\ell's leading coefficient — the same cℓc_\ell falls out (symbolic, ℓ≤4\ell\le4).
  • General direction: eik⋅r=4π∑ℓmiℓjℓ(kr)Yℓm∗(k^)Yℓm(r^)e^{i\mathbf{k}\cdot\mathbf{r}}=4\pi\sum_{\ell m}i^\ell j_\ell(kr)Y_\ell^{m*}(\hat{k})Y_\ell^m(\hat{r}) — via the addition theorem owned by Associated Legendre Functions (verified at random orientations to 10−1510^{-15}).
  • The ℓ-budget: channel weights (2ℓ+1)jℓ2(kr)(2\ell+1)j_\ell^2(kr) sum to exactly 1; at kr=10kr=10, 90% of the beam sits in ℓ≤10\ell\le10; semiclassically ℓ+12=kb\ell+\tfrac12=kb (impact parameter).
  • Convergence is geometric-then-abrupt: partial sums need L≈kr+L\approx kr+ a few for 10−610^{-6} (measured: 7, 15, 36 at kr=1,5,20kr=1,5,20) — the turnstile again: channels beyond krkr barely exist.
  • The ledger: asymptotically each channel splits as [eikr−(−1)ℓe−ikr]/2ikr[e^{ikr}-(-1)^\ell e^{-ikr}]/2ikr — equal in/out magnitudes: a plane wave is a standing wave in every channel; scattering can only twist outgoing phases (owners: the scattering chapter, named).
  • Key idea: one exact identity translates beam language into channel language — and its in/out balance is the zero point from which all scattering is measured.
The core insight

The plane wave decomposes exactly into free spherical channels with universal coefficients (2ℓ+1)iℓ(2\ell+1)i^\ell: an infinite, perfectly balanced superposition — incoming equals outgoing in every ℓ — whose channel weights die abruptly beyond ℓ≈kr\ell\approx kr, so that a beam meeting a target of size RR engages only ∼kR\sim kR partial waves.

Where this lives today

Partial-wave analysis is the daily bread of hadron spectroscopy: amplitude analyses at facilities like Jefferson Lab and BESIII decompose measured angular distributions into exactly this article's channels to hunt resonances — modern "PWA" collaborations are named after the acronym. Ultracold-atom physics lives at the opposite end of the budget: kr≪1kr\ll1 collapses the expansion to one channel, making the s-wave scattering length the universal currency of quantum gases (with deliberate p-wave experiments tuning past the turnstile). And electromagnetic simulation reuses the identity verbatim — Mie codes and T-matrix methods for nanoparticle optics expand incident laser fields in vector versions of equation (4), so every plasmonics paper's "multipole decomposition" is Rayleigh's expansion, still on duty a century and a half later.

The physical picture

A plane wave has no origin; a channel decomposition needs one. Plant the origin (the future target's seat) anywhere, and ask: from this vantage, how much of the passing beam carries angular momentum ℓ\ell? Classical intuition answers first: a particle passing at impact parameter bb carries L=pb=ℏkbL=pb=\hbar kb, so angular momentum ℓℏ\ell\hbar corresponds to a ring of impact parameters around b=(ℓ+12)/kb=(\ell+\tfrac12)/k. A uniform beam illuminates all rings — every ℓ, with weight growing like the ring's circumference (2ℓ+1)(2\ell+1) — until the ring radius exceeds the beam's coherent reach at radius rr, where the weights collapse. That picture is the expansion: coefficients (2ℓ+1)(2\ell+1), weights cut off at ℓ≈kr\ell\approx kr, verified below to the digit.

Figure
Impact-parameter rings: the classical skeleton of the partial-wave budget
beam → rings of impact parameter b = (ℓ+½)/k → channels ℓ plane wave: all impact parameters at once ℓ=0: b ≈ ½k⁻¹ ℓ=1 ℓ=2 ℓ ≳ kR: misses a target of size R weights (2ℓ+1)j_ℓ²(kr): sum = 1 exactly (verified); 90% of a kr=10 beam sits in ℓ ≤ 10 5 MeV α on a 7 fm nucleus: kR = 6.85 (computed) — about seven channels do all the work

The quantum content the classical picture misses is phase. The coefficients are not (2ℓ+1)(2\ell+1) but (2ℓ+1)iℓ(2\ell+1)i^\ell: each channel enters a quarter-turn rotated from its neighbor. That iℓi^\ell is bookkeeping with teeth — it is exactly what cancels the −ℓπ/2-\ell\pi/2 phase lags of the far-field jℓj_\ell's (Spherical Bessel and Neumann Functions' ladder), so that all channels rebuild a flat wavefront when summed. Drop it, or misplace one power, and the partial sums reconstruct garbage; the verifier's convergence tables double as an iℓi^\ell audit.

The expansion, derived twice

Set up on the beam axis: eikz=eikrμe^{ikz}=e^{ikr\mu} with μ=cos⁡θ\mu=\cos\theta. The left side is a free-particle energy eigenstate; the free spherical waves at the same energy, regular at the origin, are jℓ(kr)Pℓ(μ)j_\ell(kr)P_\ell(\mu) (the m=0m=0 sector suffices — the beam is axially symmetric, so no ϕ\phi-dependence, no m≠0m\ne0; that observation kills infinitely many coefficients before any computation). Completeness of the channel basis (The Free Particle in Three Dimensions' resolution of the identity, consumed) licenses

eikrμ=∑ℓ=0∞cℓ jℓ(kr) Pℓ(μ)(1)e^{ikr\mu}=\sum_{\ell=0}^{\infty}c_\ell\,j_\ell(kr)\,P_\ell(\mu) \qquad (1)

where: the only unknowns are the constants cℓc_\ell — constants, not functions of rr, precisely because both sides solve the same equation at the same energy and the regular radial solution per channel is unique up to normalization (the Wronskian independence certificate of Spherical Bessel and Neumann Functions, cited). Everything now hangs on extracting numbers.

Derivation ARoute A — projection: Legendre orthogonality reads off each channel

We extract cℓc_\ell by projecting equation (1) against PℓP_\ell, using the orthogonality relation owned by Associated Legendre Functions (consumed, never re-proved): ∫−11PℓPℓ′ dμ=22ℓ+1δℓℓ′\int_{-1}^1P_\ell P_{\ell'}\,d\mu=\tfrac{2}{2\ell+1}\delta_{\ell\ell'}. The derivation is exact at every radius — no asymptotic step anywhere — which is its advantage over the folklore versions.

Step 1. Multiply (1) by Pℓ(μ)P_\ell(\mu) and integrate over μ\mu:

∫−11Pℓ(μ) eikrμ dμ=2 cℓ2ℓ+1 jℓ(kr)(2)\int_{-1}^{1}P_\ell(\mu)\,e^{ikr\mu}\,d\mu=\frac{2\,c_\ell}{2\ell+1}\,j_\ell(kr) \qquad (2)

The left side is a classical integral (a Legendre–Fourier transform); the statement to establish is that it equals 2 iℓjℓ(kr)2\,i^\ell j_\ell(kr) — from which cℓ=(2ℓ+1)iℓc_\ell=(2\ell+1)i^\ell immediately.

Step 2. Establish the integral's identity honestly. Both sides of "12∫Pℓeixμdμ=iℓjℓ(x)\tfrac12\int P_\ell e^{ix\mu}d\mu=i^\ell j_\ell(x)" satisfy the same spherical Bessel equation in xx (differentiate under the integral, integrate by parts twice using Legendre's equation — the standard manipulation, consumed from the Legendre article's toolkit) and both are regular at x=0x=0; by the uniqueness-of-regular-solution argument they are proportional. The constant is fixed at leading order: expand eixμe^{ix\mu}, the first surviving term against PℓP_\ell is (ix)ℓℓ!∫Pℓμℓdμ=(ix)ℓℓ!⋅2ℓ+1ℓ!2(2ℓ+1)!\tfrac{(ix)^\ell}{\ell!}\int P_\ell\mu^\ell d\mu=\tfrac{(ix)^\ell}{\ell!}\cdot\tfrac{2^{\ell+1}\ell!^2}{(2\ell+1)!}, and matching against iℓxℓ/(2ℓ+1)!!i^\ell x^\ell/(2\ell+1)!! closes exactly (the double-factorial arithmetic verified symbolically).

Step 3. Numerical audit, because coefficient derivations are exactly where sign and factor errors hide: the projection integral computed by quadrature equals (2ℓ+1)iℓjℓ(x)/(2ℓ+1)⋅2(2\ell+1)i^\ell j_\ell(x)/(2\ell+1)\cdot2... — concretely, 2ℓ+12∫Pℓeixμdμ\tfrac{2\ell+1}{2}\int P_\ell e^{ix\mu}d\mu against (2ℓ+1)iℓjℓ(x)(2\ell+1)i^\ell j_\ell(x): agreement to 10−1010^{-10} for all ℓ≤5\ell\le5 at x=0.7, 3.3, 11.0x=0.7,\,3.3,\,11.0 (eighteen checks). The identity is exact and the constant is (2ℓ+1)iℓ(2\ell+1)i^\ell — no surviving convention freedom. Why the iℓi^\ell is inevitable: eikrμe^{ikr\mu} is complex and even/odd in μ\mu alternately by channel; jℓj_\ell and PℓP_\ell are real. The alternating parity forces alternating real/imaginary projections — a factor iℓi^\ell, not a choice. Real expansions (e.g. of cos⁡(krμ)\cos(kr\mu)) would split into even channels only. *Verified: 18 quadrature checks; the leading-order double-factorial match symbolic for ℓ≤4\ell\le4 (article verifier).*

Go deeper · Derivation B · 2 minRoute B — leading powers, and the balanced far-field ledger
Derivation BRoute B — leading powers, and the balanced far-field ledger

We re-derive cℓc_\ell by an independent route — matching the lowest surviving power of rr channel by channel — then push the expansion to large rr and split it into the incoming/outgoing ledger that scattering theory will inherit. (This derivation absorbs the coefficient-matching argument formerly resident in the free-particle article's base — relocated here, where the expansion lives.)

Step 1. Expand both sides of (1) in powers of x=krx=kr at fixed μ\mu and collect the coefficient of xℓμℓx^\ell\mu^\ell. Left: from eixμ=∑(ixμ)n/n!e^{ix\mu}=\sum(ix\mu)^n/n!, the term iℓ/ℓ!i^\ell/\ell!. Right: only channels ℓ′≤ℓ\ell'\le\ell reach xℓx^\ell (the turnstile jℓ′∼xℓ′j_{\ell'}\sim x^{\ell'}), and among them only ℓ′=ℓ\ell'=\ell supplies μℓ\mu^\ell (Pℓ′P_{\ell'} has degree ℓ′\ell') — one term: cℓ⋅1(2ℓ+1)!!⋅(2ℓ)!2ℓℓ!2c_\ell\cdot\tfrac{1}{(2\ell+1)!!}\cdot\tfrac{(2\ell)!}{2^\ell\ell!^2}, using PℓP_\ell's leading coefficient. Equate: cℓ=(2ℓ+1)iℓc_\ell=(2\ell+1)i^\ell (the factorial identity verified symbolically for ℓ≤4\ell\le4). Two routes, one answer — the agreement is itself a nontrivial audit of the whole function toolkit.

Step 2. The far-field ledger. Insert the verified asymptotics jℓ→sin⁡(kr−ℓπ/2)/krj_\ell\to\sin(kr-\ell\pi/2)/kr and use the exact identity iℓsin⁡(x−ℓπ/2)=12i[eix−(−1)ℓe−ix]i^\ell\sin(x-\ell\pi/2)=\tfrac{1}{2i}[e^{ix}-(-1)^\ell e^{-ix}] (verified symbolically, ℓ≤3\ell\le3):

eikz  →kr≫ℓ2  12ikr∑ℓ(2ℓ+1)[eikr−(−1)ℓe−ikr]Pℓ(cos⁡θ)(3)e^{ikz}\;\xrightarrow{kr\gg\ell^2}\;\frac{1}{2ikr}\sum_\ell(2\ell+1)\Big[e^{ikr}-(-1)^\ell e^{-ikr}\Big]P_\ell(\cos\theta) \qquad (3)

Step 3. Read the ledger. In every channel, the outgoing amplitude (eikre^{ikr} column) and the incoming amplitude (e−ikre^{-ikr} column) have equal magnitude (2ℓ+1)/2kr(2\ell+1)/2k r — the plane wave is a standing wave in each channel, transporting no net radial probability (the flux calculus of Spherical Hankel Functions and Outgoing Waves, cited: this is jℓ=12(h(1)+h(2))j_\ell=\tfrac12(h^{(1)}+h^{(2)}) channel by channel). The relative sign −(−1)ℓ-(-1)^\ell is the reflection phase off the centrifugal core. Why this is the scattering chapter's birth certificate: place any central target at the origin. It cannot touch the incoming column (those waves haven't arrived yet — causality) and, if elastic, cannot change outgoing magnitudes (the exterior doctrine's flux balance) — it can only multiply each outgoing term by a phase e2iδℓe^{2i\delta_\ell}. Everything measurable then lives in those phases. That machinery — amplitude, f(θ)f(\theta), cross sections — belongs BY NAME to Scattering Amplitude, Partial-Wave Expansion (the scattering chapter's own serial of that name) and Phase Shifts; this article stops at the balanced ledger they will unbalance. *Verified: split identity symbolic; convergence of partial sums to eikze^{ikz} at 10−610^{-6} with L=7,15,36L=7,15,36 for kr=1,5,20kr=1,5,20 (article verifier).*

The axially symmetric identity generalizes to any beam direction via THE addition theorem of Associated Legendre Functions (its proof, cited not re-run): Pℓ(k^⋅r^)=4π2ℓ+1∑mYℓm∗(k^)Yℓm(r^)P_\ell(\hat{k}\cdot\hat{r})=\tfrac{4\pi}{2\ell+1}\sum_mY_\ell^{m*}(\hat{k})Y_\ell^m(\hat{r}), giving

eik⋅r=4π∑ℓ=0∞∑m=−ℓℓiℓ jℓ(kr) Yℓm∗(k^) Yℓm(r^)(4)e^{i\mathbf{k}\cdot\mathbf{r}}=4\pi\sum_{\ell=0}^{\infty}\sum_{m=-\ell}^{\ell}i^\ell\,j_\ell(kr)\,Y_\ell^{m*}(\hat{k})\,Y_\ell^m(\hat{r}) \qquad (4)

where: the YℓmY_\ell^m catalogue is Spherical Harmonics' (consumed), and the formula's symmetry between k^\hat{k} and r^\hat{r} is worth a pause — the beam direction and the observation direction enter as conjugate twins, which is why the same identity serves scattering (fixed k^\hat{k}, varying r^\hat{r}) and momentum-space analysis (the reverse). Verified at randomly drawn k^,r^\hat{k},\hat{r}: sum through ℓ=40\ell=40 matches eik⋅re^{i\mathbf{k}\cdot\mathbf{r}} to 2×10−152\times10^{-15} (computed). Finally, the budget quantified. Weight channel ℓ by its share of the local probability: wℓ(kr)=(2ℓ+1)jℓ2(kr)w_\ell(kr)=(2\ell+1)j_\ell^2(kr). The completeness of the expansion forces ∑ℓwℓ=1\sum_\ell w_\ell=1 exactly (verified numerically to 10−610^{-6} at kr=10kr=10 — an identity known classically as the "addition theorem for spherical Bessel functions"), so the wℓw_\ell are honest probabilities: at kr=10kr=10, 90% of the beam sits in ℓ≤10\ell\le10 (computed), with the maximum near ℓ≈kr\ell\approx kr and a collapse beyond — the semiclassical rings, weighed.

Figure
The balanced ledger: every channel of a plane wave holds equal incoming and outgoing columns
far field, channel by channel: (2ℓ+1)/2ikr × [ e^{ikr} − (−1)^ℓ e^{−ikr} ] channel outgoing e^{ikr} incoming e^{−ikr} ℓ = 0 sign − ℓ = 1 ℓ = 2 sign + equal magnitudes (2ℓ+1)/2kr in every row — zero net flux: a standing account an elastic target may only rotate the red column's phase: ×e^{2iδ_ℓ} — the scattering chapter's variable bar lengths ∝ (2ℓ+1) at fixed kr; alternating relative sign −(−1)^ℓ = reflection phase off the centrifugal core verified: split identity i^ℓ sin(x−ℓπ/2) = [e^{ix}−(−1)^ℓ e^{−ix}]/2i, symbolic ℓ ≤ 3 (article verifier)

Limiting cases

caseresultmeaning
r→0r\to0only ℓ=0\ell=0 survives: eikz→1=j0(0)e^{ikz}\to1=j_0(0)the origin hears one channel — contact physics is s-wave
kr≪1kr\ll1wℓ∼(kr)2ℓw_\ell\sim(kr)^{2\ell}: budget collapses onto ℓ=0\ell=0thermal-neutron regime (the free-particle article's census)
kr≫1kr\gg1weights flat up to ℓ≈kr\ell\approx kr, then cliffclassical rings recovered; Lneeded≈kr+L_{\text{needed}}\approx kr+ few (measured 7/15/36)
θ=0\theta=0Pℓ(1)=1P_\ell(1)=1: all channels in phaseforward direction = perfect reconstruction; seed of the optical theorem (owner named)
real part onlyeven channels of cos⁡(krμ)\cos(kr\mu)the iℓi^\ell as parity bookkeeping
V≠0V\neq0 addedoutgoing column rephased: e2iδℓe^{2i\delta_\ell}the scattering chapter begins (its serials named)

Worked examples

Example 1 — how many channels is a plane wave, operationally? Reconstruct eikze^{ikz} at kr=1, 5, 20kr=1,\,5,\,20 and count.

Step 1. Sum the expansion to increasing LL and measure the worst-angle error against the exact exponential. Results (computed): error below 10−610^{-6} at L=7L=7 for kr=1kr=1; L=15L=15 for kr=5kr=5; L=36L=36 for kr=20kr=20.

Step 2. Read the pattern: L≈kr+L\approx kr+ a nearly constant margin (the margin grows only like (kr)1/3(kr)^{1/3} — the turnstile's soft edge). The identity is an infinite sum; the physics at radius rr is a finite one with ∼kr\sim kr terms.

Step 3. Budget check at kr=10kr=10: the weights (2ℓ+1)jℓ2(2\ell+1)j_\ell^2 sum to 1.0000001.000000 (computed) and accumulate to 90% by ℓ=10\ell=10. Answer: "a plane wave at radius rr" is, to experimental accuracy, kr+O((kr)1/3)kr+\mathcal{O}\big((kr)^{1/3}\big) channels — 7, 15, 36 at kr=1,5,20kr=1,5,20 for 10−610^{-6}.

Example 2 — the ℓ-budget of an α-particle on gold. A 5 MeV α approaches a nucleus of radius R=7R=7 fm. Which channels can matter?

Step 1. Wavenumber: k=2mαE/ℏ=0.978 fm−1k=\sqrt{2m_\alpha E}/\hbar=0.978\,\text{fm}^{-1} (computed, mαc2=3727.4m_\alpha c^2=3727.4 MeV).

Step 2. The target engages channels whose impact-parameter rings pass within it: ℓ+12≲kR=6.85\ell+\tfrac12\lesssim kR=6.85 (computed) — about seven channels (ℓ=0…6\ell=0\dots6), no more. Channels beyond sail past (their turnstile keeps them outside RR), contributing to the beam but not to the collision.

Step 3. Consequence for the experiment's analysis: a partial-wave fit of α–gold elastic scattering at this energy needs seven phase shifts, not seventy — the expansion turns an infinite-dimensional unknown (the potential) into a handful of numbers at each energy. (For the honest record: the α is charged, so the long Coulomb tail modifies the far field logarithmically — the Coulomb-function caveat of the free-particle article, one line, scattering chapter's machinery.) Answer: kR=6.85kR=6.85 ⇒ roughly seven active channels — why partial-wave analysis is practical, not just exact.

Origins · 1 minPartial Wave Expansion and how this was discovered
PW
Partial Wave Expansion
1872–1933 · Terling Place, Kristiania, Bristol

The identity at this article's heart was born acoustic, went quantum in one paper, and became an industry in a textbook. Rayleigh — whose perturbation-theory side lives in another portrait of this track — derived the expansion of a plane sound wave in spherical harmonics in his Theory of Sound era (1872 onward), needing exactly what quantum mechanics would need: to know how a flat wavefront addresses a round obstacle, channel by channel. For half a century it remained a tool of acoustics and optics. Then in 1927, in Kristiania (Oslo)'s orbit, Hilding Faxén and Johan Holtsmark ported the whole apparatus into the brand-new wave mechanics to attack a concrete puzzle — the Ramsauer–Townsend transparency of noble gases to slow electrons (a story whose physics belongs to the scattering chapter) — and in doing so wrote down quantum partial-wave analysis essentially in its modern form: Rayleigh's expansion for the incident beam, phase shifts per channel for the target's reply. The method's canonization came with Mott and Massey's The Theory of Atomic Collisions (1933, from Bristol and beyond), the treatise that made "expand in partial waves, fit the phase shifts" the default grammar of collision physics — the grammar every nuclear, atomic, and particle scattering paper still speaks. A century-spanning relay — Rayleigh's identity, Faxén–Holtsmark's port, Mott–Massey's codification — with this article's eighteen quadrature checks as the modern audit of the baton.

History
1872
— Rayleigh expands plane waves in spherical harmonics for acoustic scattering
1927
— Faxén and Holtsmark port the expansion into quantum mechanics for electron-atom collisions
1928
— Holtsmark's phase-shift computations explain the Ramsauer–Townsend transparency
1933
— Mott and Massey's treatise codifies partial-wave analysis as the language of collision theory
1950
s — Nuclear and particle physics adopt phase-shift analysis wholesale (pion-nucleon, nucleon-nucleon)

Going further

  • Sakurai & Napolitano §6.4 — the expansion in its scattering habitat.
  • Griffiths & Schroeter Ch. 10 — the identity and its cross-section payoff at the source register.
  • Shankar Ch. 19 — the free-particle dictionary developed in full.
  • Mott & Massey, The Theory of Atomic Collisions — the codifying treatise (historical depth).

Provenance note: Rayleigh expansion and coefficient derivation verifiable in Arfken–Weber §15 and Sakurai §6.4; verifier should reproduce cℓ=iℓ(2ℓ+1)c_\ell=i^\ell(2\ell+1) via the integral ∫−11uℓPℓ(u) du=2ℓ+1(ℓ!)2/(2ℓ+1)!\int_{-1}^1 u^\ell P_\ell(u)\,\mathrm{d}u=2^{\ell+1}(\ell!)^2/(2\ell+1)! combined with jℓ(x)→xℓ/(2ℓ+1)!!j_\ell(x)\to x^\ell/(2\ell+1)!!. Asymptotic form jℓ→sin⁡(x−ℓπ/2)/xj_\ell\to\sin(x-\ell\pi/2)/x and the scattering amplitude f=12ik∑(2ℓ+1)(e2iδℓ−1)Pℓf=\frac{1}{2ik}\sum(2\ell+1)(e^{2i\delta_\ell}-1)P_\ell are standard (Griffiths Ch. 11, Sakurai §6.4). Lowest jℓj_\ell explicit forms checkable by direct substitution into the spherical Bessel ODE. — The coefficient reproduction is executed in this article's verify script BOTH ways (the quoted moment-integral route = Route B's leading-power match, symbolic ℓ≤4\ell\le4; plus 18 independent quadrature projections at ℓ≤5\ell\le5), together with the split identity (ℓ≤3\ell\le3 symbolic), measured convergence L=7/15/36L=7/15/36 at kr=1/5/20kr=1/5/20, and the general-direction form at random orientations (10−1510^{-15}), the unit-sum of channel weights, the 90% budget line, and the α-on-gold numbers. The provenance's scattering-amplitude formula f=12ik∑(2ℓ+1)(e2iδℓ−1)Pℓf=\frac{1}{2ik}\sum(2\ell+1)(e^{2i\delta_\ell}-1)P_\ell is quoted for the record and NOT derived or verified here — it is the property of the scattering chapter's serials named below. Edition note: the coefficient derivation formerly in the base of The Free Particle in Three Dimensions is owned HERE (relocation noted at both ends). The scattering payoff — amplitudes f(θ)f(\theta), σ=(4π/k2)∑(2ℓ+1)sin⁡2δℓ\sigma=(4\pi/k^2)\sum(2\ell+1)\sin^2\delta_\ell, phase-shift extraction — belongs BY NAME to the scattering chapter's Scattering Amplitude, Partial-Wave Expansion, and Phase Shifts (whose slug, note, differs from this article's: that serial expands the scattering amplitude; this one expands the free beam). The optical theorem's forward-direction seed is likewise its owner's (Optical Theorem — a cutover triple-claimant, re-flagged here for the record).

Is the expansion of any use away from scattering problems?

Everywhere a spherical boundary meets a flat wave. Antenna theory expands incident fields on spherical near-field scan surfaces; acoustics expands engine noise on spherical microphone arrays; quantum chemistry expands plane-wave basis sets around atomic centers (the "projector augmented wave" machinery of modern DFT codes contains exactly these integrals). The identity is a coordinate-translation tool; scattering is merely its most famous customer.

The sum over ℓ is infinite — does the expansion converge at every point?

Yes, absolutely and uniformly on any bounded region (the terms decay super-geometrically once ℓ>ekr/2\ell>ekr/2, by the turnstile estimate). What fails is uniformity in rr: for fixed truncation the error creeps outward — the warning box's trap. At any fixed point the series converges fast; there is no Gibbs phenomenon because the function being expanded is smooth.

What exactly does a potential do to the expansion — one sentence per column?

Incoming column: untouched — causality forbids the target from editing waves still arriving. Outgoing column: each term multiplied by Sℓ=e2iδℓS_\ell=e^{2i\delta_\ell} (elastic) or by a complex number of modulus ≤1\le1 (absorptive). The difference between the modified and free outgoing columns, collected over channels, is the scattered wave — one sentence that is the entire scattering chapter's program, executed serial by serial there.

Why is the weight (2ℓ+1)j_ℓ² a probability — probabilities of what, exactly?

Of finding the particle's angular momentum to be ℓ\ell, conditioned on it being at radius rr. Project the plane wave onto the shell of radius rr: the angular content decomposes over the orthonormal Yℓ0Y_\ell^0, and the squared amplitude of each component is (2ℓ+1)jℓ2(kr)(2\ell+1)j_\ell^2(kr) (after the shell normalization that makes them sum to 1 — the identity verified numerically). It is an honest measurement distribution: an L^2\hat{L}^2 measurement on a beam particle known to be at radius rr returns ℓ(ℓ+1)ℏ2\ell(\ell+1)\hbar^2 with exactly this weight.

Does the identity survive for wave packets, which are the physical beams?

Linearly, yes: superpose plane waves of different kk and each expands independently — the packet's channel content is the kk-weighted average of (2ℓ+1)iℓjℓ(kr)(2\ell+1)i^\ell j_\ell(kr) coefficients. For packets narrow in kk (every real beam), the monochromatic budget is an excellent proxy, which is why cross sections computed from single-kk partial waves describe real accelerator data. Wide packets need the kk-integral kept — a bookkeeping, not a conceptual, complication (the wave-packet articles of the earlier chapters own that grammar).

Is there a version for particles with spin?

Yes — the orbital identity is unchanged, but the useful basis couples YℓmY_\ell^m with spin states into total-angular-momentum channels (the Clebsch–Gordan machinery of the previous chapter). A spin-½ beam expands into channels labeled (ℓ,j=ℓ±12)(\ell,j=\ell\pm\tfrac12), and targets with spin-orbit forces shift the two jj-partners differently — the origin of polarization phenomena in nuclear scattering. The pure-orbital expansion here is the s=0s=0 case and the scaffold for all others.

Before this

Where this leads