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 — the beam every accelerator and every scattering theory starts from — and the spherical channels built over the last four serials. The translation between them is a single exact identity, the Rayleigh expansion:
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 verified numerically eighteen ways and the general-direction form 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 sum exactly to one, 90% of a beam sits below , and the semiclassical dictionary 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.
- THE identity: — exact, all , all ; coefficients derived twice and verified numerically for at three radii (18 checks).
- Route A (projection): Legendre orthogonality peels off channel ℓ; the radial integral is — no asymptotics anywhere.
- Route B (leading powers): match coefficients of using and 's leading coefficient — the same falls out (symbolic, ).
- General direction: — via the addition theorem owned by Associated Legendre Functions (verified at random orientations to ).
- The ℓ-budget: channel weights sum to exactly 1; at , 90% of the beam sits in ; semiclassically (impact parameter).
- Convergence is geometric-then-abrupt: partial sums need a few for (measured: 7, 15, 36 at ) — the turnstile again: channels beyond barely exist.
- The ledger: asymptotically each channel splits as — 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 plane wave decomposes exactly into free spherical channels with universal coefficients : an infinite, perfectly balanced superposition — incoming equals outgoing in every ℓ — whose channel weights die abruptly beyond , so that a beam meeting a target of size engages only 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: 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 ? Classical intuition answers first: a particle passing at impact parameter carries , so angular momentum corresponds to a ring of impact parameters around . A uniform beam illuminates all rings — every ℓ, with weight growing like the ring's circumference — until the ring radius exceeds the beam's coherent reach at radius , where the weights collapse. That picture is the expansion: coefficients , weights cut off at , verified below to the digit.
The quantum content the classical picture misses is phase. The coefficients are not but : each channel enters a quarter-turn rotated from its neighbor. That is bookkeeping with teeth — it is exactly what cancels the phase lags of the far-field '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 audit.
The expansion, derived twice
Set up on the beam axis: with . The left side is a free-particle energy eigenstate; the free spherical waves at the same energy, regular at the origin, are (the sector suffices — the beam is axially symmetric, so no -dependence, no ; 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
where: the only unknowns are the constants — constants, not functions of , 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.
We extract by projecting equation (1) against , using the orthogonality relation owned by Associated Legendre Functions (consumed, never re-proved): . The derivation is exact at every radius — no asymptotic step anywhere — which is its advantage over the folklore versions.
Step 1. Multiply (1) by and integrate over :
The left side is a classical integral (a Legendre–Fourier transform); the statement to establish is that it equals — from which immediately.
Step 2. Establish the integral's identity honestly. Both sides of "" satisfy the same spherical Bessel equation in (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 ; by the uniqueness-of-regular-solution argument they are proportional. The constant is fixed at leading order: expand , the first surviving term against is , and matching against 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 ... — concretely, against : agreement to for all at (eighteen checks). The identity is exact and the constant is — no surviving convention freedom. Why the is inevitable: is complex and even/odd in alternately by channel; and are real. The alternating parity forces alternating real/imaginary projections — a factor , not a choice. Real expansions (e.g. of ) would split into even channels only. *Verified: 18 quadrature checks; the leading-order double-factorial match symbolic for (article verifier).*
Go deeper · Derivation B · 2 minRoute B — leading powers, and the balanced far-field ledger
We re-derive by an independent route — matching the lowest surviving power of channel by channel — then push the expansion to large 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 at fixed and collect the coefficient of . Left: from , the term . Right: only channels reach (the turnstile ), and among them only supplies ( has degree ) — one term: , using 's leading coefficient. Equate: (the factorial identity verified symbolically for ). 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 and use the exact identity (verified symbolically, ):
Step 3. Read the ledger. In every channel, the outgoing amplitude ( column) and the incoming amplitude ( column) have equal magnitude — 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 channel by channel). The relative sign 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 . Everything measurable then lives in those phases. That machinery — amplitude, , 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 at with for (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): , giving
where: the catalogue is Spherical Harmonics' (consumed), and the formula's symmetry between and is worth a pause — the beam direction and the observation direction enter as conjugate twins, which is why the same identity serves scattering (fixed , varying ) and momentum-space analysis (the reverse). Verified at randomly drawn : sum through matches to (computed). Finally, the budget quantified. Weight channel ℓ by its share of the local probability: . The completeness of the expansion forces exactly (verified numerically to at — an identity known classically as the "addition theorem for spherical Bessel functions"), so the are honest probabilities: at , 90% of the beam sits in (computed), with the maximum near and a collapse beyond — the semiclassical rings, weighed.
Limiting cases
| case | result | meaning |
|---|---|---|
| only survives: | the origin hears one channel — contact physics is s-wave | |
| : budget collapses onto | thermal-neutron regime (the free-particle article's census) | |
| weights flat up to , then cliff | classical rings recovered; few (measured 7/15/36) | |
| : all channels in phase | forward direction = perfect reconstruction; seed of the optical theorem (owner named) | |
| real part only | even channels of | the as parity bookkeeping |
| added | outgoing column rephased: | the scattering chapter begins (its serials named) |
Worked examples
Example 1 — how many channels is a plane wave, operationally? Reconstruct at and count.
Step 1. Sum the expansion to increasing and measure the worst-angle error against the exact exponential. Results (computed): error below at for ; for ; for .
Step 2. Read the pattern: a nearly constant margin (the margin grows only like — the turnstile's soft edge). The identity is an infinite sum; the physics at radius is a finite one with terms.
Step 3. Budget check at : the weights sum to (computed) and accumulate to 90% by . Answer: "a plane wave at radius " is, to experimental accuracy, channels — 7, 15, 36 at for .
Example 2 — the ℓ-budget of an α-particle on gold. A 5 MeV α approaches a nucleus of radius fm. Which channels can matter?
Step 1. Wavenumber: (computed, MeV).
Step 2. The target engages channels whose impact-parameter rings pass within it: (computed) — about seven channels (), no more. Channels beyond sail past (their turnstile keeps them outside ), 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: ⇒ roughly seven active channels — why partial-wave analysis is practical, not just exact.
Origins · 1 minPartial Wave Expansion and how this was discovered
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.
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 via the integral combined with . Asymptotic form and the scattering amplitude are standard (Griffiths Ch. 11, Sakurai §6.4). Lowest 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 ; plus 18 independent quadrature projections at ), together with the split identity ( symbolic), measured convergence at , and the general-direction form at random orientations (), the unit-sum of channel weights, the 90% budget line, and the α-on-gold numbers. The provenance's scattering-amplitude formula 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 , , 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 , by the turnstile estimate). What fails is uniformity in : 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 (elastic) or by a complex number of modulus (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 , conditioned on it being at radius . Project the plane wave onto the shell of radius : the angular content decomposes over the orthonormal , and the squared amplitude of each component is (after the shell normalization that makes them sum to 1 — the identity verified numerically). It is an honest measurement distribution: an measurement on a beam particle known to be at radius returns with exactly this weight.
Does the identity survive for wave packets, which are the physical beams?
Linearly, yes: superpose plane waves of different and each expands independently — the packet's channel content is the -weighted average of coefficients. For packets narrow in (every real beam), the monochromatic budget is an excellent proxy, which is why cross sections computed from single- partial waves describe real accelerator data. Wide packets need the -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 with spin states into total-angular-momentum channels (the Clebsch–Gordan machinery of the previous chapter). A spin-½ beam expands into channels labeled , and targets with spin-orbit forces shift the two -partners differently — the origin of polarization phenomena in nuclear scattering. The pure-orbital expansion here is the case and the scaffold for all others.