← Part IV: Relativistic Electrodynamics
Jackson Chapter 16

Radiation Reaction

The Abraham-Lorentz force, the Lorentz-Dirac equation, runaway solutions, preacceleration, and the limits of classical electrodynamics.

3.1 The Abraham-Lorentz Force

A radiating charge loses energy. By Newton's third law, there must be a reaction force on the charge due to its own radiation. The Abraham-Lorentz formula gives this self-force:

$$\boxed{\mathbf{F}_{\text{rad}} = \frac{\mu_0 q^2}{6\pi c}\dot{\mathbf{a}} = m\tau_0\dot{\mathbf{a}}}$$

where $\tau_0 = \mu_0 q^2/(6\pi mc) = q^2/(6\pi\epsilon_0 mc^3) \approx 6.26\times 10^{-24}$ s for an electron.

Derivation 1: Abraham-Lorentz Force from Energy Conservation

We derive the radiation reaction force by demanding that the work done against it equals the radiated power over a complete cycle.

Step 1: Energy balance

The power radiated by a non-relativistic charge is $P = \mu_0 q^2 a^2/(6\pi c)$. The work done by the radiation reaction force is $\mathbf{F}_{\text{rad}}\cdot\mathbf{v}$. Energy conservation requires:

$$\int_{t_1}^{t_2}\mathbf{F}_{\text{rad}}\cdot\mathbf{v}\,dt = -\int_{t_1}^{t_2}\frac{\mu_0 q^2}{6\pi c}a^2\,dt$$

Step 2: Integrate by parts

$\int a^2\,dt = \int \dot{v}\cdot\dot{v}\,dt = [\dot{v}\cdot v]_{t_1}^{t_2} - \int \ddot{v}\cdot v\,dt$. If the motion is periodic or the boundary terms vanish:

$$\int a^2\,dt = -\int\dot{\mathbf{a}}\cdot\mathbf{v}\,dt$$

Step 3: Identify the force

$$\int\mathbf{F}_{\text{rad}}\cdot\mathbf{v}\,dt = \frac{\mu_0 q^2}{6\pi c}\int\dot{\mathbf{a}}\cdot\mathbf{v}\,dt$$

Step 4: Read off the force

$$\boxed{\mathbf{F}_{\text{rad}} = \frac{\mu_0 q^2}{6\pi c}\dot{\mathbf{a}} = m\tau_0\dot{\mathbf{a}}}$$

Step 5: Equation of motion

The full equation becomes: $m\mathbf{a} = \mathbf{F}_{\text{ext}} + m\tau_0\dot{\mathbf{a}}$. This is a third-order ODE, requiring three initial conditions. The extra degree of freedom leads to the notorious runaway and preacceleration problems.

3.2 Runaway Solutions

The Abraham-Lorentz equation admits unphysical "runaway" solutions where the acceleration grows exponentially even in the absence of external forces:

$$m\mathbf{a} = m\tau_0\dot{\mathbf{a}} \implies \mathbf{a}(t) = \mathbf{a}_0 e^{t/\tau_0}$$

Derivation 2: Eliminating Runaways - The Integro-Differential Equation

We replace the Abraham-Lorentz equation with an integro-differential form that eliminates runaways but introduces preacceleration.

Step 1: Rewrite the equation

$m(\mathbf{a} - \tau_0\dot{\mathbf{a}}) = \mathbf{F}_{\text{ext}}$. Multiply by the integrating factor $e^{-t/\tau_0}$:

$$\frac{d}{dt}(\mathbf{a}\,e^{-t/\tau_0}) = -\frac{1}{m\tau_0}\mathbf{F}_{\text{ext}}\,e^{-t/\tau_0}$$

Step 2: Integrate from t to infinity

Requiring $\mathbf{a}(t)\to 0$ as $t\to\infty$ (no runaways):

$$\mathbf{a}(t) = \frac{1}{m\tau_0}\int_t^\infty \mathbf{F}_{\text{ext}}(t')\,e^{-(t'-t)/\tau_0}\,dt'$$

Step 3: Preacceleration

The acceleration at time $t$ depends on the force at future times $t' > t$. For a force that turns on at $t = 0$, the particle begins accelerating at $t < 0$.

Step 4: Timescale of acausality

The preacceleration decays with time constant $\tau_0 \approx 6\times 10^{-24}$ s, far smaller than any classical measurement timescale. Quantum effects (pair creation, Compton scattering) become important at this scale.

Step 5: Interpretation

The preacceleration is not a true violation of causality but an indication that classical electrodynamics is incomplete at the scale $\tau_0 \sim r_e/c$. Quantum electrodynamics resolves these issues consistently, with radiation reaction arising naturally from the interaction of the electron with its own quantized field.

3.3 The Lorentz-Dirac Equation

The covariant generalization of the Abraham-Lorentz equation is the Lorentz-Dirac equation:

$$\boxed{m a^\mu = F_{\text{ext}}^\mu + m\tau_0\left(\dot{a}^\mu - \frac{a_\nu a^\nu}{c^2}u^\mu\right)}$$

The last term ensures $u_\mu a^\mu = 0$ is preserved. Here $\dot{a}^\mu = da^\mu/d\tau$.

Derivation 3: Covariant Radiation Reaction from the Larmor Invariant

We derive the Lorentz-Dirac equation by requiring the reaction force to balance the covariant Larmor power.

Step 1: Covariant power balance

The 4-vector radiation reaction force $\Gamma^\mu$ must satisfy $\Gamma^\mu u_\mu = 0$ (orthogonal to 4-velocity) and reproduce the Larmor power: $\Gamma^\mu u_\mu/c = 0$.

Step 2: Most general form

The most general 4-vector quadratic in the motion is $\Gamma^\mu = \alpha\dot{a}^\mu + \beta a_\nu a^\nu u^\mu$. Requiring $u_\mu\Gamma^\mu = 0$ and using $u_\mu\dot{a}^\mu = -a_\mu a^\mu$ fixes $\beta = \alpha/c^2$.

Step 3: Fix the coefficient

Integrating the energy component $\Gamma^0 v$ over a period and comparing with the integrated Larmor power fixes $\alpha = m\tau_0$.

Step 4: Full equation

$$m a^\mu = F_{\text{ext}}^\mu + m\tau_0\left(\dot{a}^\mu - \frac{a_\nu a^\nu}{c^2}u^\mu\right)$$

The term $-(a_\nu a^\nu/c^2)u^\mu$ is the Schott term, responsible for the relativistic correction to radiation reaction.

Step 5: Non-relativistic limit

For $v \ll c$: $a^\mu a_\mu \approx -a^2$ and $\dot{a}^i \approx \dot{a}_i$. The spatial part reduces to $m\mathbf{a} = \mathbf{F}_{\text{ext}} + m\tau_0\dot{\mathbf{a}}$, recovering the Abraham-Lorentz equation.

3.4 Classical Models of the Electron

The radiation reaction problem is intimately connected to the self-energy of a classical charged particle. A point charge has infinite electromagnetic self-energy:

$$W_{\text{self}} = \frac{1}{2}\int\frac{|\mathbf{E}|^2}{\mu_0 c^2}\,d^3r = \frac{e^2}{8\pi\epsilon_0}\int_0^\infty\frac{dr}{r^2} \to \infty$$

Derivation 4: Electromagnetic Mass and the 4/3 Problem

We derive the electromagnetic contribution to the electron mass and encounter the famous factor-of-4/3 discrepancy.

Step 1: Self-energy of a charged sphere

A uniformly charged sphere of radius $a$ has electrostatic self-energy $W = e^2/(8\pi\epsilon_0 a) \times 3/5$ (for a uniform volume distribution) or $W = e^2/(8\pi\epsilon_0 a)$ (surface charge).

Step 2: Electromagnetic mass

Identifying $W = m_{\text{em}}c^2$: $m_{\text{em}} = e^2/(8\pi\epsilon_0 ac^2)$. If all the electron mass is electromagnetic, then $a = r_e/2$ where $r_e = e^2/(4\pi\epsilon_0 m_e c^2)$.

Step 3: Momentum of the field

The electromagnetic momentum of a moving charged sphere is $\mathbf{p}_{\text{em}} = (4/3)m_{\text{em}}\gamma\mathbf{v}$. The factor 4/3 (not 1) means $E \neq pc$ for the field alone.

Step 4: Poincare stresses

Poincare showed that non-electromagnetic cohesive forces (Poincare stresses) are needed to hold the electron together. These stresses contribute $-W/3$ to the total momentum, restoring the correct relativistic relation $p = \gamma mv$.

Step 5: Resolution in QED

In quantum electrodynamics, the self-energy divergence is handled by renormalization: the bare mass $m_0$ and the electromagnetic self-energy $\delta m$ combine to give the observed physical mass $m = m_0 + \delta m$. Only the physical mass is observable, and the theory makes finite, testable predictions (e.g., the anomalous magnetic moment to 10 decimal places).

3.5 Radiation Damping of Oscillating Systems

For a bound electron oscillating at frequency $\omega_0$, the radiation reaction acts as a damping force with rate $\Gamma = \tau_0\omega_0^2$:

$$\ddot{x} + \Gamma\dot{x} + \omega_0^2 x = \frac{F_{\text{ext}}}{m}, \qquad \Gamma = \frac{e^2\omega_0^2}{6\pi\epsilon_0 m_e c^3}$$

Derivation 5: Natural Linewidth from Radiation Damping

We derive the classical linewidth of spectral lines from the radiation damping of an atomic oscillator.

Step 1: Damped oscillator solution

For $\Gamma \ll \omega_0$: $x(t) = x_0 e^{-\Gamma t/2}\cos(\omega_0 t)$. The energy decays as $e^{-\Gamma t}$.

Step 2: Fourier transform of the signal

The Fourier transform of $e^{-\Gamma t/2}e^{-i\omega_0 t}$ for $t > 0$ is a Lorentzian:

$$|\hat{x}(\omega)|^2 \propto \frac{1}{(\omega-\omega_0)^2 + (\Gamma/2)^2}$$

Step 3: Full width at half maximum

$$\Delta\omega = \Gamma = \tau_0\omega_0^2$$

Step 4: Classical lifetime

The $e$-folding time for energy decay is $\tau = 1/\Gamma$. For the Lyman-alpha line ($\lambda = 121.6$ nm): $\Gamma \approx 6.3\times 10^8$ s$^{-1}$, giving $\tau \approx 1.6$ ns.

Step 5: Comparison with quantum result

The classical result gives the correct order of magnitude. The quantum mechanical linewidth for hydrogen Lyman-alpha is $\Gamma_{\text{QM}} = \omega_0^3 r_e/(3c) \times f_{12}$, where $f_{12}$ is the oscillator strength. The classical and quantum widths agree to within factors of order unity, one of the remarkable successes of the classical oscillator model.

Historical Notes

Abraham (1903) and Lorentz (1904): Max Abraham and Hendrik Lorentz independently derived the non-relativistic radiation reaction force. Abraham's rigid electron model led to the famous 4/3 problem.

Poincare (1906): Henri Poincare showed that non-electromagnetic stresses are needed to stabilize the classical electron, resolving the 4/3 discrepancy and anticipating the need for new physics at small scales.

Dirac (1938): Paul Dirac derived the covariant equation of motion including radiation reaction (the Lorentz-Dirac equation) using only the conservation of energy-momentum and the retarded Green's function.

Wheeler and Feynman (1945): John Wheeler and Richard Feynman proposed the absorber theory of radiation, which eliminates self-interaction by using half-retarded plus half-advanced fields, offering an alternative perspective on radiation reaction.

Applications

Laser-Plasma Interactions

At intensities above $10^{23}$ W/cm$^2$, radiation reaction significantly affects electron dynamics in ultra-intense laser fields. This regime is being explored at facilities like ELI and APOLLON.

Spectral Line Broadening

The natural linewidth of atomic transitions arises from radiation damping. This sets the fundamental limit on spectral resolution and is essential for laser physics and precision spectroscopy.

Beam Dynamics in Storage Rings

Radiation damping in electron storage rings reduces beam emittance (radiation cooling). The equilibrium beam size results from the balance between radiation damping and quantum excitation.

Gravitational Wave Emission

The gravitational radiation reaction (Burke-Thorne formula) is the gravitational analog of the Abraham-Lorentz force. The orbital decay of the Hulse-Taylor binary pulsar confirmed this effect, earning the 1993 Nobel Prize.

Simulation: Radiation Reaction

This simulation visualizes radiative damping of an oscillating electron, compares runaway and physical (preacceleration) solutions to the Abraham-Lorentz equation, shows the driven oscillator response with radiation damping, and compares characteristic timescales of classical electrodynamics.

Radiation Reaction: Abraham-Lorentz, Dirac & Damping

Python
script.py129 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Rate this chapter: