Relativistic Particles & Radiation
Synchrotron radiation, the relativistic Larmor formula, betatron motion, angular distributions of relativistic radiators, and radiation energy loss.
2.1 Relativistic Equation of Motion
The relativistic equation of motion for a charged particle in an electromagnetic field is:
In 3-vector form: $\frac{d(\gamma m\mathbf{v})}{dt} = q(\mathbf{E}+\mathbf{v}\times\mathbf{B})$. The magnetic force does no work; only $\mathbf{E}$ changes the particle energy.
Derivation 1: Relativistic Larmor Formula
We generalize the Larmor formula to fully relativistic motion using 4-vector notation.
Step 1: Covariant form of radiated power
The radiated power must be a Lorentz scalar. The unique scalar quadratic in the 4-acceleration $a^\mu = du^\mu/d\tau$ is:
$$P = -\frac{\mu_0 q^2 c}{6\pi}a_\mu a^\mu$$
Step 2: Express $a_\mu a^\mu$ in terms of lab quantities
The 4-acceleration has components $a^\mu = \gamma^2(\gamma^2\mathbf{v}\cdot\mathbf{a}/c, \gamma^2(\mathbf{v}\cdot\mathbf{a})\mathbf{v}/c^2 + \mathbf{a})$. Computing the invariant:
$$-a_\mu a^\mu = \gamma^4\left[a^2 + \gamma^2\frac{(\mathbf{v}\cdot\mathbf{a})^2}{c^2}\right] = \gamma^6\left[a^2 - \frac{|\mathbf{v}\times\mathbf{a}|^2}{c^2}\right]$$
Step 3: Relativistic Larmor formula
$$\boxed{P = \frac{\mu_0 q^2 c}{6\pi}\gamma^6\left(a^2 - \frac{|\mathbf{v}\times\mathbf{a}|^2}{c^2}\right)}$$
Step 4: Acceleration parallel to velocity
When $\mathbf{a} \parallel \mathbf{v}$: $P = \mu_0 q^2 c\gamma^6 a^2/(6\pi)$. This applies to linear accelerators. The $\gamma^6$ factor makes relativistic linear acceleration extremely radiatively costly.
Step 5: Acceleration perpendicular to velocity
When $\mathbf{a} \perp \mathbf{v}$: $P = \mu_0 q^2 c\gamma^4 a^2/(6\pi)$. This applies to synchrotrons and storage rings. The power scales as $\gamma^4$, which is why electron synchrotrons are limited to ~100 GeV while proton synchrotrons reach TeV energies.
2.2 Synchrotron Radiation
A relativistic electron in a magnetic field emits synchrotron radiation, characterized by a broad spectrum peaking near the critical frequency:
where $R = \gamma mc/(eB)$ is the orbit radius and $\omega_0 = eB/(\gamma m)$ is the cyclotron frequency.
Derivation 2: Critical Frequency of Synchrotron Radiation
We derive the characteristic frequency of synchrotron radiation using the forward-beaming argument.
Step 1: Opening angle of the radiation cone
Relativistic beaming concentrates the radiation into a forward cone of half-angle $\Delta\theta \sim 1/\gamma$.
Step 2: Duration of the pulse
An observer sees radiation only while the beam sweeps past, covering an arc $\Delta\phi \sim 2/\gamma$. The retarded time interval is $\Delta t_{\text{ret}} = 2R/(c\gamma)$.
Step 3: Time compression by the Doppler effect
Due to the relativistic approach, the observed pulse duration is compressed: $\Delta t_{\text{obs}} = \Delta t_{\text{ret}}(1-\beta) \approx \Delta t_{\text{ret}}/(2\gamma^2)$.
Step 4: Characteristic frequency
$$\omega_c \sim \frac{1}{\Delta t_{\text{obs}}} \sim \frac{2\gamma^2}{2R/(c\gamma)} = \frac{\gamma^3 c}{R}$$
The exact result includes a factor of $3/2$: $\omega_c = \frac{3}{2}\gamma^3(c/R)$.
Step 5: Spectrum
The spectral distribution involves the Airy integral: $dP/d\omega \propto (\omega/\omega_c)\int_{\omega/\omega_c}^\infty K_{5/3}(x)\,dx$ where $K_{5/3}$ is a modified Bessel function. It rises as $\omega^{1/3}$ for $\omega \ll \omega_c$ and falls exponentially for $\omega \gg \omega_c$.
2.3 Angular Distribution of Relativistic Radiation
The radiation from a relativistic particle is concentrated in a narrow forward cone due to the Lorentz transformation. The angular power distribution for acceleration perpendicular to velocity is:
Derivation 3: Forward Beaming from Relativistic Aberration
We derive the beaming effect from the relativistic aberration formula.
Step 1: Aberration formula
If a photon is emitted at angle $\theta'$ in the rest frame, in the lab frame: $\cos\theta = (\cos\theta' + \beta)/(1+\beta\cos\theta')$.
Step 2: Maximum emission angle
In the rest frame, maximum emission is at $\theta' = \pi/2$. In the lab: $\cos\theta_{\max} = \beta$, i.e., $\sin\theta_{\max} = 1/\gamma$. For $\gamma \gg 1$: $\theta_{\max} \approx 1/\gamma$.
Step 3: Solid angle compression
The isotropic rest-frame radiation ($4\pi$ sr) is compressed into a forward cone of solid angle $\Delta\Omega \sim \pi/\gamma^2$. The intensity is enhanced by a factor of $\sim\gamma^2$.
Step 4: Doppler boosting of power
The emitted power transforms as $dP/d\Omega = \mathcal{D}^4(dP'/d\Omega')$ where $\mathcal{D} = [\gamma(1-\beta\cos\theta)]^{-1}$ is the Doppler factor. For $\theta \approx 0$: $\mathcal{D} \approx 2\gamma$.
Step 5: Observational consequences
A synchrotron source moving toward the observer has apparent luminosity boosted by $\sim\gamma^4$, while one moving away is suppressed. This explains the one-sidedness of relativistic jets in AGN and the extreme brightness of blazars (jets aimed at Earth).
2.4 Betatron & Cyclotron Dynamics
The betatron uses a time-varying magnetic field to accelerate electrons in a circular orbit. The betatron condition requires the field at the orbit to be half the average field enclosed:
Derivation 4: The Betatron Condition
We derive the condition for a stable circular orbit of constant radius in a betatron.
Step 1: Circular orbit condition
At radius $R_0$: $p = eBR_0$ (relativistically exact). For the orbit to remain at $R_0$, both $p$ and $B$ must change together.
Step 2: Faraday's law gives the accelerating field
The induced electric field at $R_0$: $E \cdot 2\pi R_0 = -d\Phi/dt$ where $\Phi = \int_0^{R_0}B(r)2\pi r\,dr$.
Step 3: Rate of momentum change
$dp/dt = eE = \frac{e}{2\pi R_0}\frac{d\Phi}{dt}$. From the orbit condition: $dp/dt = eR_0\,dB(R_0)/dt$.
Step 4: Equate and integrate
$eR_0\,dB(R_0)/dt = \frac{e}{2\pi R_0}\frac{d\Phi}{dt}$. Integrating: $\Phi = 2\pi R_0^2 B(R_0)$.
Step 5: The betatron condition
$$\boxed{B(R_0) = \frac{\langle B\rangle}{2} = \frac{\Phi}{2\pi R_0^2}}$$
The field at the orbit must equal half the average field within the orbit. This is achieved by shaping the pole pieces. The first betatron (Kerst, 1940) accelerated electrons to 2.3 MeV; modern betatrons reach 300 MeV.
2.5 Radiation Energy Loss in Accelerators
The energy radiated per revolution in a circular accelerator is:
For electrons at LEP (E = 100 GeV, R = 4.3 km): $\Delta E \approx 3$ GeV/turn. For protons at the same energy: $\Delta E \sim 10^{-6}$ eV/turn.
Derivation 5: Energy Loss per Revolution
We compute the energy radiated in one complete orbit of a circular accelerator.
Step 1: Power for circular motion
For perpendicular acceleration in a circular orbit: $a = v^2/R \approx c^2/R$ and $P = \mu_0 q^2 c\gamma^4 a^2/(6\pi)$.
Step 2: Substitute $a = c^2/R$
$$P = \frac{\mu_0 q^2 c^5\gamma^4}{6\pi R^2} = \frac{2r_0 c}{3}\frac{E^4}{(mc^2)^4 R^2}$$
where $r_0 = q^2/(4\pi\epsilon_0 mc^2)$ is the classical particle radius.
Step 3: Revolution period
$T = 2\pi R/v \approx 2\pi R/c$ for ultrarelativistic particles.
Step 4: Energy loss per revolution
$$\Delta E = PT = \frac{4}{3}\frac{r_0}{(mc^2)^3}\frac{E^4}{R}$$
Step 5: Mass dependence
$$\boxed{\Delta E \propto \frac{E^4}{m^4 R}}$$
Since $m_p/m_e \approx 1836$, protons lose $(1836)^4 \approx 10^{13}$ times less energy than electrons at the same energy. This is why the largest circular electron collider (LEP) reached 209 GeV while the LHC reaches 14 TeV with protons.
Historical Notes
Schott (1912): George Adolphus Schott first computed the radiation from a charged particle moving in a circle, laying the theoretical foundation for synchrotron radiation decades before it was observed.
Kerst (1940): Donald Kerst built the first working betatron at the University of Illinois, accelerating electrons to 2.3 MeV using the time-varying magnetic flux principle.
Elder et al. (1947): Synchrotron radiation was first observed visually as a bright arc of light at the General Electric 70 MeV synchrotron, confirming the theoretical predictions.
Schwinger (1949): Julian Schwinger provided the definitive theoretical treatment of synchrotron radiation, computing the complete angular and spectral distributions that are used to this day.
Applications
Synchrotron Light Sources
Third-generation light sources (APS, ESRF, SPring-8) use bending magnets and undulators to produce intense X-ray beams. The critical energy $\epsilon_c = \frac{3}{2}\hbar\gamma^3 c/R$ is tuned from UV to hard X-rays.
Free Electron Lasers
FELs exploit the coherent emission from relativistic electrons in undulators. LCLS and European XFEL produce femtosecond X-ray pulses with peak brilliance $10^{10}$ times greater than synchrotrons.
Astrophysical Synchrotron
Synchrotron radiation from relativistic electrons in magnetic fields produces the radio emission of supernova remnants, pulsar wind nebulae, and AGN jets. The power-law spectrum reveals the electron energy distribution.
Circular Collider Design
Synchrotron radiation loss sets the energy limit of circular electron colliders. The proposed FCC-ee at CERN (91 km) would reach 365 GeV, while the FCC-hh with protons could reach 100 TeV.
Simulation: Relativistic Particles
This simulation shows the synchrotron radiation spectrum, relativistic forward beaming of radiation patterns, a betatron orbit with radiation energy loss, and synchrotron power and critical energy vs. Lorentz factor.
Relativistic Particles: Synchrotron, Betatron & Radiation Loss
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server