Partial-wave analysis

Partial-wave analysis, in the context of quantum mechanics, refers to a technique for solving scattering problems by decomposing each wave into its constituent angular-momentum components and solving using boundary conditions. Partial wave analysis is typically useful for low energy scattering where only a few angular momentum components dominate. At high energy were scattering is weak, an alternative called the Born approximation is used.^[1]: 507

Preliminary scattering theory

A steady beam of particles scatters off a spherically symmetric potential ${\displaystyle V(r)}$, which is short-ranged, so that for large distances ${\displaystyle r\to \infty }$, the particles behave like free particles. The incoming beam is assumed to be a collimated plane wave ${\displaystyle \exp(ikz)}$ traveling along the z axis. Because the beam is switched on for times long compared to the time of interaction of the particles with the scattering potential, a steady state is assumed. This means that the stationary Schrödinger equation for the wave function ${\displaystyle \Psi (\mathbf {r} )}$ representing the particle beam should be solved:

${\displaystyle \left[-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}+V(r)\right]\Psi (\mathbf {r} )=E\Psi (\mathbf {r} ).}$

We make the following ansatz:

${\displaystyle \Psi (\mathbf {r} )=\Psi _{0}(\mathbf {r} )+\Psi _{\text{s}}(\mathbf {r} ),}$

where ${\displaystyle \Psi _{0}(\mathbf {r} )\propto \exp(ikz)}$ is the incoming plane wave, and ${\displaystyle \Psi _{\text{s}}(\mathbf {r} )}$ is a scattered part perturbing the original wave function.

It is the asymptotic form of ${\displaystyle \Psi _{\text{s}}(\mathbf {r} )}$ that is of interest, because observations near the scattering center (e.g. an atomic nucleus) are mostly not feasible, and detection of particles takes place far away from the origin. At large distances, the particles should behave like free particles, and ${\displaystyle \Psi _{\text{s}}(\mathbf {r} )}$ should therefore be a solution to the free Schrödinger equation. For a spherically symmetric potential, these solutions should be outgoing spherical waves,${\displaystyle \Psi _{\text{s}}(\mathbf {r} )\propto \exp(ikr)/r}$ at large distances. Thus the asymptotic form of the scattered wave is chosen as^[2]: 371

${\displaystyle \Psi _{\text{s}}(\mathbf {r} )\to f(\theta ,k){\frac {\exp(ikr)}{r}},}$

where ${\displaystyle f(\theta ,k)}$ is the so-called scattering amplitude, which is in this case only dependent on the elevation angle ${\displaystyle \theta }$ and the energy. This gives the following asymptotic expression for the entire wave function:

${\displaystyle \Psi (\mathbf {r} )\to \Psi ^{(+)}(\mathbf {r} )=\exp(ikz)+f(\theta ,k){\frac {\exp(ikr)}{r}}.}$

Partial-wave expansion

In case of a spherically symmetric potential ${\displaystyle V(\mathbf {r} )=V(r)}$, the scattering wave function may be expanded in spherical harmonics, which reduce to Legendre polynomials because of azimuthal symmetry (no dependence on ${\displaystyle \phi }$):

${\displaystyle \Psi (\mathbf {r} )=\sum _{\ell =0}^{\infty }{\frac {u_{\ell }(r)}{r}}P_{\ell }(\cos \theta ).}$

In the standard scattering problem, the incoming beam is assumed to take the form of a plane wave of wave number k, which can be decomposed into partial waves using the plane-wave expansion in terms of spherical Bessel functions and Legendre polynomials:

${\displaystyle \psi _{\text{in}}(\mathbf {r} )=e^{ikz}=\sum _{\ell =0}^{\infty }(2\ell +1)i^{\ell }j_{\ell }(kr)P_{\ell }(\cos \theta ).}$

Here we have assumed a spherical coordinate system in which the z axis is aligned with the beam direction. The radial part of this wave function consists solely of the spherical Bessel function, which can be rewritten as a sum of two spherical Hankel functions:

${\displaystyle j_{\ell }(kr)={\frac {1}{2}}\left(h_{\ell }^{(1)}(kr)+h_{\ell }^{(2)}(kr)\right).}$

This has physical significance: h_ℓ⁽²⁾ asymptotically (i.e. for large r) behaves as i^−(ℓ+1)e^ikr/(kr) and is thus an outgoing wave, whereas h_ℓ⁽¹⁾ asymptotically behaves as i^ℓ+1e^−ikr/(kr) and is thus an incoming wave. The incoming wave is unaffected by the scattering, while the outgoing wave is modified by a factor known as the partial-wave S-matrix element S_ℓ:

${\displaystyle {\frac {u_{\ell }(r)}{r}}{\stackrel {r\to \infty }{\longrightarrow }}{\frac {i^{\ell }k}{\sqrt {2\pi }}}\left(h_{\ell }^{(1)}(kr)+S_{\ell }h_{\ell }^{(2)}(kr)\right),}$

where u_ℓ(r)/r is the radial component of the actual wave function. The scattering phase shift δ_ℓ is defined as half of the phase of S_ℓ:

${\displaystyle S_{\ell }=e^{2i\delta _{\ell }}.}$

If flux is not lost, then |S_ℓ| = 1, and thus the phase shift is real. This is typically the case, unless the potential has an imaginary absorptive component, which is often used in phenomenological models to simulate loss due to other reaction channels.

Therefore, the full asymptotic wave function is

${\displaystyle \psi (\mathbf {r} ){\stackrel {r\to \infty }{\longrightarrow }}\sum _{\ell =0}^{\infty }(2\ell +1)i^{\ell }{\frac {h_{\ell }^{(1)}(kr)+S_{\ell }h_{\ell }^{(2)}(kr)}{2}}P_{\ell }(\cos \theta ).}$

Subtracting ψ_in yields the asymptotic outgoing wave function:

${\displaystyle \psi _{\text{out}}(\mathbf {r} ){\stackrel {r\to \infty }{\longrightarrow }}\sum _{\ell =0}^{\infty }(2\ell +1)i^{\ell }{\frac {S_{\ell }-1}{2}}h_{\ell }^{(2)}(kr)P_{\ell }(\cos \theta ).}$

Making use of the asymptotic behavior of the spherical Hankel functions, one obtains

${\displaystyle \psi _{\text{out}}(\mathbf {r} ){\stackrel {r\to \infty }{\longrightarrow }}{\frac {e^{ikr}}{r}}\sum _{\ell =0}^{\infty }(2\ell +1){\frac {S_{\ell }-1}{2ik}}P_{\ell }(\cos \theta ).}$

Since the scattering amplitude f(θ, k) is defined from

${\displaystyle \psi _{\text{out}}(\mathbf {r} ){\stackrel {r\to \infty }{\longrightarrow }}{\frac {e^{ikr}}{r}}f(\theta ,k),}$

it follows that^[2]: 386

${\displaystyle f(\theta ,k)=\sum _{\ell =0}^{\infty }(2\ell +1){\frac {S_{\ell }-1}{2ik}}P_{\ell }(\cos \theta )=\sum _{\ell =0}^{\infty }(2\ell +1){\frac {e^{i\delta _{\ell }}\sin \delta _{\ell }}{k}}P_{\ell }(\cos \theta ),}$

and thus the differential cross section is given by

${\displaystyle {\frac {d\sigma }{d\Omega }}=|f(\theta ,k)|^{2}={\frac {1}{k^{2}}}\left|\sum _{\ell =0}^{\infty }(2\ell +1)e^{i\delta _{\ell }}\sin \delta _{\ell }P_{\ell }(\cos \theta )\right|^{2}.}$

This works for any short-ranged interaction. For long-ranged interactions (such as the Coulomb interaction), the summation over ℓ may not converge. The general approach for such problems consist in treating the Coulomb interaction separately from the short-ranged interaction, as the Coulomb problem can be solved exactly in terms of Coulomb functions, which take on the role of the Hankel functions in this problem.

See also

• Levinson's theorem