Synchrotron Radiation

Angular Distribution of Radiation from Accelerated Charges closed with a promissory note: a 1/γ1/\gamma searchlight sweeping past a fixed observer compresses the received flash to Δtobs∼ρ/γ3c\Delta t_{\rm obs} \sim \rho/\gamma^3 c, so the spectrum must reach ωc∼γ3c/ρ\omega_c \sim \gamma^3 c/\rho. This article cashes that note. The result is the most consequential spectrum in classical electrodynamics: a particle circulating at one modest frequency — megahertz, in a storage ring — radiating a smooth continuum that stretches ten decades up to hard X-rays, with a universal shape, a universal ω1/3\omega^{1/3} low-frequency rise, an exponential cutoff, and a polarization signature so distinctive it identifies the mechanism across the universe. Every synchrotron light source on Earth — the machines that solve protein structures and image battery electrodes — is a factory built around the critical frequency ωc=32γ3c/ρ\omega_c = \tfrac{3}{2}\gamma^3 c/\rho, and every circular collider design must first survive the energy bill that Relativistic Generalization of Larmor's Formula computes. But the reach is longer: replace the bending magnet by a microgauss interstellar field and the same formulas explain why supernova remnants, radio galaxies, and quasar jets shine with featureless power-law spectra whose polarization maps invisible magnetic fields. One derivation, run at ρ=7\rho = 7 m and ρ=1015\rho = 10^{15} m. The article works three layers honestly: the time-domain argument that fixes ωc\omega_c (executed, not estimated); the universal spectral function F(ω/ωc)F(\omega/\omega_c) with its limits derived and its numbers audited; and the ensemble theorem — a power-law electron population radiates a power-law spectrum with index (p−1)/2(p-1)/2 — that turns radio astronomy into particle spectroscopy.

Key takeaways
  • A circulating charge sweeps its 1/γ1/\gamma cone (owned by Angular Distribution of Radiation from Accelerated Charges) past the observer once per turn; the received flash lasts Δtobs≈ρ/2γ3c\Delta t_{\rm obs} \approx \rho/2\gamma^3 c — one power of γ\gamma from the sweep, two from the head-on subtraction 1−β1 - \beta.
  • The flash width sets the critical frequency ωc=32γ3c/ρ=32γ2eB/m\omega_c = \tfrac{3}{2}\gamma^3 c/\rho = \tfrac{3}{2}\gamma^2 eB/m: the spectrum is smooth up to ∼ωc\sim\omega_c and dies exponentially beyond. Exactly half the radiated power sits above ωc\omega_c, half below.
  • The single-particle spectrum is universal: dP/dω=(93/8π)(Ptot/ωc) F(ω/ωc)dP/d\omega = (9\sqrt{3}/8\pi)(P_{\rm tot}/\omega_c)\,F(\omega/\omega_c) with F(x)=x∫x∞K5/3 dx′F(x) = x\int_x^\infty K_{5/3}\,dx' — peak at 0.29 ωc0.29\,\omega_c, an ω1/3\omega^{1/3} rise at low frequency (independent of particle energy), x e−x\sqrt{x}\,e^{-x} beyond the knee.
  • Polarization is the fingerprint: 100% linear in the orbital plane, elliptical off-plane, and Pσ:Pπ=7:1P_\sigma : P_\pi = 7:1 integrated — 75% linear polarization overall.
  • An electron population N(E)∝E−pN(E) \propto E^{-p} radiates I(ω)∝ω−(p−1)/2I(\omega) \propto \omega^{-(p-1)/2}: the power-law sky of radio astronomy is a particle spectrometer.
  • The total power and the per-turn energy bill are owned by Relativistic Generalization of Larmor's Formula and recalled by name; this article owns where that energy lands in frequency.
  • Key idea: relativistic beaming turns one orbital frequency into ten decades of spectrum — the γ3\gamma^3 pulse compression fixes ωc\omega_c, the Bessel machinery fixes the universal shape, and a power-law ensemble turns both into astronomy.
The core insight

For an ultrarelativistic charge on a circular arc, the 1/γ1/\gamma beaming cone compresses the received pulse by γ3\gamma^3, so the spectrum extends to ωc=32γ3c/ρ\omega_c = \tfrac{3}{2}\gamma^3 c/\rho — a universal continuum, ω1/3\omega^{1/3} below the peak and exponentially cut off above, strongly polarized, with total power ∝γ4\propto \gamma^4.

Where this lives today

Fourth-generation storage rings (MAX IV, ESRF-EBS, APS-U, Sirius) are precision applications of equation (4): multi-bend achromat lattices shrink the electron beam until photon brightness is limited by the radiation's own 1/γ1/\gamma opening angle, and every insertion device — wiggler or undulator — is a scheme for making B(s)B(s) do something cleverer than a single bend. Free-electron lasers push the same electrons into self-organized microbunches whose emission scales as N2N^2 — the coherence theorem of Larmor Formula from the Liénard–Wiechert Fields deployed at X-ray wavelengths. Meanwhile the E4E^4 energy bill still writes accelerator history: it is why the 100 km FCC-ee design pays gigaelectronvolts per turn at top energy (Relativistic Generalization of Larmor's Formula owns those numbers) and why no one proposes a circular 10 TeV electron machine. In the sky, synchrotron emission is the default explanation for any smooth, polarized, power-law continuum: supernova remnants, pulsar wind nebulae, radio-galaxy lobes, blazar jets, burst afterglows. Polarization maps trace magnetic-field geometry kiloparsecs across; spectral breaks date electron populations via the cooling-time clock of worked example 2; and the extreme end — the Crab's gamma-ray flares exceeding the classical radiation-reaction ceiling on synchrotron photon energy (∼\sim160 MeV) — currently stress-tests acceleration theory itself, a frontier this chapter's self-force articles illuminate from the other side.

Origins · 1 minJulian Schwinger and how this was discovered
JS
Julian Schwinger
1918–1994 · Harvard — whose 1949 paper turned the synchrotron spectrum into engineering

Schwinger enters this track twice, and the split mirrors the chapter's own division of labor: Relativistic Generalization of Larmor's Formula honors the accelerator-systematics side of his 1949 masterwork "On the Classical Radiation of Accelerated Electrons"; this article honors the spectrum. Schwinger computed the full angle- and frequency-resolved intensity of a charge on a circular arc — the stationary-phase analysis behind equation (3) — organized it into the universal function F(ω/ωc)F(\omega/\omega_c), and fixed the critical frequency's coefficient and the polarization split, all while simultaneously building the QED renormalization program that would bring him the 1965 Nobel Prize. He had derived much of it years earlier for wartime waveguide problems and published only when the GE observation made the question urgent — the theorist's spectrum arriving just in time to explain the technician's blue-white spot. The warm-up exercise now underwrites some seventy light sources and most of high-energy astrophysics.

History
1898
— Liénard computes the velocity-dependent radiation of a fast charge; the angular machinery is in place
1907
— Schott's Adams Prize essay works out orbital radiation harmonics, seeking (wrongly) a classical atom
1944
— Ivanenko and Pomeranchuk show betatron energies are capped by the γ⁴ radiation loss
1946
— Blewett measures the predicted orbit shrinkage at General Electric — the loss is real, the light unseen
1947
— a GE technician looks through a viewing port at the 70 MeV synchrotron and sees the light
1949
— Schwinger's paper delivers the exact spectrum, critical frequency, and polarization
1953
— Shklovsky attributes the Crab Nebula's continuum to synchrotron radiation; polarization confirms it within a year
1956
— Tomboulian and Hartman perform the first synchrotron-light spectroscopy at Cornell — the light source era begins
2016
— MAX IV opens the fourth generation: storage rings brilliant enough to be limited by the radiation's own opening angle

The physical picture

Three clocks run in a synchrotron, and the subject is the arithmetic between them. The orbit clock ticks at ω0=eB/γm\omega_0 = eB/\gamma m — megahertz in a storage ring, slower for more energetic particles at fixed field. The sweep clock measures how long the lighthouse beam illuminates you: the cone of Angular Distribution of Radiation from Accelerated Charges has half-angle 1/γ1/\gamma, so the velocity points your way for only ∼1/γ\sim 1/\gamma of the orbit. The reception clock is the cruel one: during the flash the charge chases its own light head-on, and wavefronts emitted over Δtr\Delta t_r arrive squeezed by 1−β≈1/2γ21 - \beta \approx 1/2\gamma^2. Multiply the three and a megahertz orbit delivers femtosecond flashes. A periodic train of narrow flashes is formally a comb of orbit harmonics reaching up to the inverse flash width — but the interesting harmonics are the ∼γ3\sim\gamma^3-th ones, and any real spread of energies or fields smears them into a continuum. That is the difference between cyclotron radiation — a nonrelativistic charge emitting one spectral line at ω0\omega_0 — and synchrotron radiation: same orbit, same Maxwell, but γ3\gamma^3 compression converts the line into ten decades of continuum. The picture also fixes the polarization before any Bessel function appears. Watched from within the orbital plane, the acceleration projects onto the sky as a purely horizontal wiggle — 100% linear polarization. From slightly above, the circular motion projects to a shallow ellipse traversed in a definite sense — elliptical polarization, opposite handedness above and below the plane. Every beamline and radio polarimeter reads this geometry.

Figure
The lighthouse geometry and the three clocks — the one-per-gamma cone sweeps past the observer once per turn and delivers a pulse train whose flashes are compressed by gamma cubed relative to the orbital period
The synchrotron lighthouse: orbit clock × sweep clock × reception clock orbit radius ρ period T = 2πρ/c charge, γ ≫ 1 beaming cone, half-angle 1/γ rigidly attached to the velocity — sweeps like a lighthouse observer illuminating arc ~ ρ/γ sweep time Δt_r ~ ρ/γc received signal period T (one turn) flash width Δt_obs ≈ ρ/2γ³c — the reception clock pays (1−β) ≈ 1/2γ² spectrum: harmonics of ω₀ filling the comb up to ω_c ~ 1/Δt_obs ~ γ³ω₀ — a continuum in practice

From one orbital frequency to the critical frequency

The orbit supplies exactly one input frequency. A charge qq of energy γmc2\gamma mc^2 in field BB circles at

ω0=eBγm=βcρ,ρ=γmβceB.(1)\omega_0 = \frac{eB}{\gamma m} = \frac{\beta c}{\rho}, \qquad \rho = \frac{\gamma m \beta c}{eB}. \qquad (1)

where: ω0\omega_0 — the relativistic gyration (orbit) frequency; ρ\rho — the bending radius. Note the 1/γ1/\gamma: at fixed BB, faster particles orbit more slowly, because the mass grows while the force does not. Everything else in this article is about how ruthlessly the radiation escapes this single frequency. The escape factor is fixed by the pulse-compression argument — the teaser Angular Distribution of Radiation from Accelerated Charges deferred here by name — and the careful version of that argument defines the one number every light source prints on its brochure:

  ωc≡32 γ3 cρ=32 γ2 eBm  (2)\boxed{\;\omega_c \equiv \frac{3}{2}\,\gamma^3\,\frac{c}{\rho} = \frac{3}{2}\,\gamma^2\,\frac{eB}{m}\;} \qquad (2)

where: ωc\omega_c — the critical frequency, the knee of the spectrum. The two forms differ by equation (1): at fixed radius the reach grows as γ3\gamma^3; at fixed field, as γ2\gamma^2 (one power is spent as the orbit slows). The 3/23/2 is a convention anchored to an exact property — half the power is radiated above ωc\omega_c, half below (audited in Derivation B) — not to the spectral peak, which sits lower, at 0.29 ωc0.29\,\omega_c. In photon-energy units, for electrons, equation (2) becomes the workhorse engineering rule εc[keV]=0.665 E2[GeV2] B[T]\varepsilon_c[\mathrm{keV}] = 0.665\,E^2[\mathrm{GeV}^2]\,B[\mathrm{T}] (recomputed for this edition: a 3 GeV ring with 1.4 T bends gives εc=8.38\varepsilon_c = 8.38 keV by both routes).

Derivation APulse compression executed and the critical frequency

The goal is equation (2) from the time domain: take the sweep geometry of Angular Distribution of Radiation from Accelerated Charges (whose lobes and κ\kappa-doctrine are recalled, not re-derived), execute the emitted-to-received clock conversion honestly — the step where two of the three powers of γ\gamma appear, and where careless treatments lose factors — and read off the spectral reach. The watch-out: the compression is not time dilation; it is the classical chase between a charge and its own wavefronts, and it must be computed with the (1−β)(1-\beta) subtraction, not with 1/γ1/\gamma.

Step 1. The sweep. The beaming cone, half-angle ≈1/γ\approx 1/\gamma, rides the instantaneous velocity. As the charge rounds the orbit at angular rate βc/ρ\beta c/\rho, the observer stays inside the cone while the velocity rotates through ≈2/γ\approx 2/\gamma — an emission time

Δtr≈2ργβc.\Delta t_r \approx \frac{2\rho}{\gamma \beta c}.

Why this is one power of gamma: the cone narrows as 1/γ1/\gamma while the rotation rate is γ\gamma-independent at fixed ρ\rho — the lighthouse spins at a fixed rate but its beam gets thinner.

Step 2. The chase. The flash's first wavefront is emitted at the arc's start; its last at the arc's end, a path length βc Δtr\beta c\,\Delta t_r closer to the observer (head-on geometry). The reception window is therefore not Δtr\Delta t_r but

Δtobs=Δtr−βc Δtrc=(1−β) Δtr≈Δtr2γ2,\Delta t_{\rm obs} = \Delta t_r - \frac{\beta c\,\Delta t_r}{c} = (1 - \beta)\,\Delta t_r \approx \frac{\Delta t_r}{2\gamma^2},

using 1−β=(1−β2)/(1+β)≈1/2γ21 - \beta = (1-\beta^2)/(1+\beta) \approx 1/2\gamma^2. This is the emitted-to-received Jacobian κ\kappa of the E9 doctrine at its head-on floor — a subtraction of two nearly equal travel times, which is why the naive replacement 1−β→1/γ1-\beta \to 1/\gamma (time dilation) is wrong by a full power of γ\gamma.

Step 3. Assemble the width and invert it:

Δtobs≈2ργβc⋅12γ2=ργ3βc⟹ωreach∼1Δtobs∼γ3cρ.\Delta t_{\rm obs} \approx \frac{2\rho}{\gamma\beta c}\cdot\frac{1}{2\gamma^2} = \frac{\rho}{\gamma^3 \beta c} \qquad\Longrightarrow\qquad \omega_{\rm reach} \sim \frac{1}{\Delta t_{\rm obs}} \sim \frac{\gamma^3 c}{\rho}.

A pulse of width Δtobs\Delta t_{\rm obs} has Fourier support up to ∼1/Δtobs\sim 1/\Delta t_{\rm obs} and little beyond — the spectral ledger machinery of Larmor Formula from the Liénard–Wiechert Fields applied to the sharpest pulse in the chapter.

Step 4. Fix the convention. The full spectral computation (quoted below) concentrates its knee at ξ=ωρ3c(1/γ2+ψ2)3/2∼1\xi = \tfrac{\omega\rho}{3c}(1/\gamma^2 + \psi^2)^{3/2} \sim 1 — on the orbital plane, ω∼3γ3c/ρ\omega \sim 3\gamma^3 c/\rho — and the community defines ωc=32γ3c/ρ\omega_c = \tfrac{3}{2}\gamma^3 c/\rho, the choice for which the radiated power splits exactly in half at ωc\omega_c (audited in Derivation B: the half-power point of FF lands at x=1.0000x = 1.0000). Corollary worth memorizing: the harmonic comb extends to nc=ωc/ω0=32γ3n_c = \omega_c/\omega_0 = \tfrac{3}{2}\gamma^3 — at γ≈5900\gamma \approx 5900, some 3×10113\times10^{11} harmonics, far too finely spaced for any real beam to resolve. Verified: matches Jackson §14.6 and Rybicki & Lightman §6.4; Δtobs=0.12\Delta t_{\rm obs} = 0.12 as and nc=3.0×1011n_c = 3.0\times10^{11} recomputed for a 3 GeV, 1.4 T ring; the (1−β)(1-\beta)-vs-1/γ1/\gamma distinction checked numerically at γ=20\gamma = 20 (factor 40 difference in predicted width).

The universal spectrum

The honest boundary of this article runs here. The full angle- and frequency-resolved intensity is obtained by Fourier-transforming the Liénard–Wiechert fields along the orbital arc — a stationary-phase computation that Jackson §14.6 executes in full and this track quotes with anatomy rather than re-deriving:

d2Idω dΩ=q216π3ε0c(2ωρ3c)2(1γ2+ψ2)2[K2/32(ξ)+ψ21/γ2+ψ2 K1/32(ξ)](3)\frac{d^2 I}{d\omega\,d\Omega} = \frac{q^2}{16\pi^3\varepsilon_0 c}\left(\frac{2\omega\rho}{3c}\right)^{2}\left(\frac{1}{\gamma^2}+\psi^2\right)^{2}\left[K_{2/3}^2(\xi) + \frac{\psi^2}{1/\gamma^2+\psi^2}\,K_{1/3}^2(\xi)\right] \qquad (3)

where: ψ\psi — the angle out of the orbital plane; ξ=ωρ3c(1/γ2+ψ2)3/2\xi = \tfrac{\omega\rho}{3c}(1/\gamma^2+\psi^2)^{3/2} — the Bessel argument, which reaches order unity exactly when ω∼ωc\omega \sim \omega_c and ψ≲1/γ\psi \lesssim 1/\gamma; K1/3,K2/3K_{1/3}, K_{2/3} — modified Bessel functions. The two terms are the two polarizations, and the structure is readable without any special-function fluency: the K2/32K_{2/3}^2 term is the field component in the orbital plane (the σ\sigma mode, driven by the projected acceleration); the K1/32K_{1/3}^2 term, carrying an explicit factor ψ2\psi^2, is the out-of-plane (π\pi) mode — it vanishes identically at ψ=0\psi = 0, which is the Bessel-language statement of the geometric fact from the physical picture: in-plane observers see 100% linear polarization. Both KK's decay exponentially once ξ≫1\xi \gg 1, enforcing simultaneously the frequency cutoff at ωc\omega_c and the angular cutoff at ψ∼1/γ\psi \sim 1/\gamma — one exponential guards both borders, and softer photons emerge in a wider cone — ψeff≈(1/γ)(2ωc/ω)1/3\psi_{\rm eff} \approx (1/\gamma)(2\omega_c/\omega)^{1/3} below the knee — so hard X-rays hug the plane tighter than infrared. Integrating (3) over angles collapses everything into one function of one variable — the universal synchrotron spectrum:

dPdω=938π Ptotωc F ⁣(ωωc),F(x)=x∫x∞K5/3(x′) dx′.(4)\frac{dP}{d\omega} = \frac{9\sqrt{3}}{8\pi}\,\frac{P_{\rm tot}}{\omega_c}\,F\!\left(\frac{\omega}{\omega_c}\right), \qquad F(x) = x\int_x^\infty K_{5/3}(x')\,dx'. \qquad (4)

where: PtotP_{\rm tot} — the total synchrotron power, owned by Relativistic Generalization of Larmor's Formula (its external-fields form gives Ptot=q4γ2B2β2/6πε0m2cP_{\rm tot} = q^4\gamma^2B^2\beta^2/6\pi\varepsilon_0 m^2c for pure-B circular motion, recalled by name); FF — the dimensionless shape every electron in the universe shares. The prefactor 93/8π9\sqrt{3}/8\pi is forced by ∫0∞F=8π/93\int_0^\infty F = 8\pi/9\sqrt{3}, so that (4) integrates back to PtotP_{\rm tot} exactly — audited numerically for this edition (and silently omitted in the source text; see the provenance addendum).

Figure
The universal synchrotron function F computed for this edition — omega to the one-third rise, peak at 0.29 omega-c, exponential cutoff, and the exact half-power split at omega-c
F(x) = x ∫ K₅/₃ dx′ — the shape every synchrotron source shares (x = ω/ω_c, log scale) 10⁻³ 10⁻² 10⁻¹ 1 8 0.95 0 ω^(1/3) rise — energy-independent √x e⁻ˣ cutoff peak: x = 0.286, F = 0.918 x = 1 (ω = ω_c): exactly half the power above, half below audits for this edition: ∫F = 8π/9√3 = 1.6123 ✓ · mean photon at 0.308 ω_c ✓ · half-power at x = 1.0000 ✓ (scipy quadrature of K₅/₃)
Go deeper · Derivation B · 2 minThe one-third power law and the shape audits
Derivation BThe one-third power law and the shape audits

The goal is the complete working anatomy of F(x)F(x): the universal ω1/3\omega^{1/3} low-frequency rise derived from the small-argument Bessel limit (with the physical reason it is energy-independent), the exponential cutoff, and the calibration numbers — peak position, half-power point, mean photon energy — computed rather than quoted. The strategy: asymptotics at both ends, quadrature in the middle.

Step 1. Small-argument limit. For ν>0\nu > 0, Kν(ξ)→12Γ(ν) (2/ξ)νK_\nu(\xi) \to \tfrac{1}{2}\Gamma(\nu)\,(2/\xi)^\nu as ξ→0\xi \to 0. The integrand of FF near its lower limit therefore behaves as K5/3(x′)∼x′−5/3K_{5/3}(x') \sim x'^{-5/3}, so the integral diverges at small xx like ∫xx′−5/3dx′∝x−2/3\int_x x'^{-5/3}dx' \propto x^{-2/3}:

F(x)=x∫x∞K5/3(x′) dx′  →x→0  x⋅C x−2/3=4π3 Γ(1/3)(x2)1/3.F(x) = x\int_x^\infty K_{5/3}(x')\,dx' \;\xrightarrow{x\to0}\; x\cdot C\,x^{-2/3} = \frac{4\pi}{\sqrt{3}\,\Gamma(1/3)}\left(\frac{x}{2}\right)^{1/3}.

Why the constant matters less than the exponent: the 1/31/3 is a fixed point of the theory — the asymptote matches full quadrature to 0.2% at x=10−4x = 10^{-4} (checked).

Step 2. Energy independence of the soft tail. Restore units: dP/dω∝(Ptot/ωc)(ω/ωc)1/3dP/d\omega \propto (P_{\rm tot}/\omega_c)(\omega/\omega_c)^{1/3}, and in the fixed-ρ\rho form the γ\gamma's cancel exactly: Ptot∝γ4/ρ2P_{\rm tot} \propto \gamma^4/\rho^2, ωc∝γ3/ρ\omega_c \propto \gamma^3/\rho give dP/dω∝ω1/3ρ−2/3×γ0dP/d\omega \propto \omega^{1/3}\rho^{-2/3} \times \gamma^0. Below its own ωc\omega_c, every electron contributes the same soft spectrum regardless of its energy — why ensembles with messy energy distributions still show a clean ω1/3\omega^{1/3} envelope at the bottom (and why Derivation C's theorem only probes electrons whose ωc\omega_c sits near the observed band).

Step 3. The hard cutoff. For ξ≫1\xi \gg 1, Kν(ξ)→π/2ξ e−ξK_\nu(\xi) \to \sqrt{\pi/2\xi}\,e^{-\xi} independent of ν\nu at leading order; Laplace-evaluating the integral gives

F(x)  →x≫1  π2 x1/2e−x,F(x) \;\xrightarrow{x\gg1}\; \sqrt{\frac{\pi}{2}}\,x^{1/2}e^{-x},

within 7% of exact quadrature already at x=10x = 10 (computed). The knee is gentle on a log axis, but beyond a few ωc\omega_c the source is simply dark: at x=20x = 20, FF has fallen below 10−810^{-8}.

Step 4. Calibration by quadrature (this edition's audit, scipy K5/3K_{5/3} integration): peak at xpeak=0.2858x_{\rm peak} = 0.2858 with Fmax=0.918F_{\rm max} = 0.918; normalization ∫0∞F=1.61227=8π/93\int_0^\infty F = 1.61227 = 8\pi/9\sqrt{3} to six digits; half-power point at x=1.0000x = 1.0000 — the exact property the 3/23/2 in equation (2) was chosen to produce; mean photon energy ⟨ℏω⟩=∫F dx/∫(F/x) dx=0.3079 ℏωc=8153ℏωc\langle\hbar\omega\rangle = \int F\,dx / \int (F/x)\,dx = 0.3079\,\hbar\omega_c = \tfrac{8}{15\sqrt{3}}\hbar\omega_c (the number-weighted mean, distinct from the power-weighted mean 1.32 ωc1.32\,\omega_c — quote the right one when estimating photon counts from power).

Step 5. Polarization bookkeeping from equation (3), integrated over ω\omega and ψ\psi for γ≫1\gamma \gg 1 (double quadrature, this edition): Pσ/Pπ=7.000P_\sigma/P_\pi = 7.000, so the integrated degree of linear polarization is (Pσ−Pπ)/(Pσ+Pπ)=75.0%(P_\sigma - P_\pi)/(P_\sigma + P_\pi) = 75.0\% — the classic 7:17:1 split, with the σ\sigma mode's dominance strongest in the plane (100%) and the handedness of the residual circular component flipping across it. Verified: all five numbers against Rybicki & Lightman §6.4 (their eq. 6.34 low-frequency coefficient reproduced) and Jackson §14.6; quadratures converged to the digits quoted; the 7:1 ratio recovered from equation (3) directly, which simultaneously audits (3)'s two-mode structure.

The ensemble theorem — power-law particles, power-law sky

A single electron radiates the universal bump of Derivation B. Astrophysical sources contain populations, and shock acceleration generically delivers power laws in energy, N(E) dE∝E−p dEN(E)\,dE \propto E^{-p}\,dE between cutoffs. The observable is then a superposition of bumps, each centered at its electron's ωc(E)∝E2B\omega_c(E) \propto E^2 B, weighted by the population — and the superposition has a closed-form slope.

Go deeper · Derivation C · 1 minPower-law electron populations radiate power-law spectra
Derivation CPower-law electron populations radiate power-law spectra

The goal is the astrophysics workhorse: for N(E)∝E−pN(E) \propto E^{-p}, the emitted spectrum obeys I(ω)∝ω−(p−1)/2I(\omega) \propto \omega^{-(p-1)/2}, so a measured radio slope hands you the exponent of an invisible particle accelerator. The route is a scaling substitution — the derivation is short but the bookkeeping of which electrons contribute where is the physics.

Step 1. Set up the superposition. Each electron of energy EE contributes spectrum (4) with Ptot(E)∝E2B2P_{\rm tot}(E) \propto E^2 B^2 and ωc(E)=aE2\omega_c(E) = a E^2 (fixed field, a∝B/m3c5a \propto B/m^3c^5 collecting constants):

I(ω)=∫dE N(E) dPdω(ω;E)∝∫dE E−p E2ωc(E) F ⁣(ωωc(E))∝∫dE E−p F ⁣(ωaE2),I(\omega) = \int dE\,N(E)\,\frac{dP}{d\omega}(\omega; E) \propto \int dE\,E^{-p}\,\frac{E^2}{\omega_c(E)}\,F\!\left(\frac{\omega}{\omega_c(E)}\right) \propto \int dE\,E^{-p}\,F\!\left(\frac{\omega}{aE^2}\right),

the E2/ωc∝E0E^2/\omega_c \propto E^0 cancellation being the first pleasant surprise: the weight is the population alone.

Step 2. Change variable to each electron's scaled frequency x=ω/aE2x = \omega/aE^2, so E=(ω/ax)1/2E = (\omega/ax)^{1/2} and dE=−12(ω/a)1/2x−3/2dxdE = -\tfrac{1}{2}(\omega/a)^{1/2}x^{-3/2}dx:

I(ω)∝ω−(p−1)/2∫0∞x(p−3)/2 F(x) dx.I(\omega) \propto \omega^{-(p-1)/2}\int_0^\infty x^{(p-3)/2}\,F(x)\,dx.

Step 3. Check the integral is a pure number. F(x)∼x1/3F(x) \sim x^{1/3} at small xx makes the integrand ∼x(p−3)/2+1/3\sim x^{(p-3)/2+1/3}, integrable at zero for p>1/3p > 1/3; the exponential cutoff handles infinity. For any physical pp the integral converges (it evaluates to a ratio of Gamma functions — Rybicki & Lightman eq. 6.36), leaving the entire frequency dependence outside:

  I(ω)∝ω−α,α=p−12.  \boxed{\;I(\omega) \propto \omega^{-\alpha}, \qquad \alpha = \frac{p-1}{2}.\;}

Step 4. Read it as a spectrometer, with its fine print. Strong-shock Fermi acceleration predicts p≈2p \approx 2–2.42.4, hence α≈0.5\alpha \approx 0.5–0.70.7 — the range radio surveys find for supernova remnants and radio galaxies. The mapping runs both ways: the Crab's radio slope α≈0.30\alpha \approx 0.30 gives p=1.6p = 1.6. The fine print: the power law holds only where electrons with ωc(E)∼ω\omega_c(E) \sim \omega exist; outside the population's cutoffs the spectrum bends — and radiative aging (the E2E^2 loss rate depleting high energies first, on the cooling timescale of worked example 2) steepens α\alpha by exactly 1/21/2 above a time-dependent break. The bends are the population's biography. Verified: substitution algebra checked symbolically; convergence exponents against the Step-1 asymptotics of Derivation B; α=(p−1)/2\alpha = (p-1)/2 and the aging steepening Δα=1/2\Delta\alpha = 1/2 against Rybicki & Lightman §6.3–6.5.

Figure
Nine electrons on a power-law energy ladder each contribute the universal bump at their own critical frequency — the envelope is a clean power law with slope minus (p minus 1) over 2, computed here for p equal to 2.5
The ensemble theorem: Σ (universal bumps × E⁻ᵖ ladder) = power law (log–log, p = 2.5) log ω (6 decades) log I(ω) (4 decades) low-E electrons: many, soft ω_c high-E: few, hard ω_c slope −(p−1)/2 = −0.75 envelope = observed spectrum measure the radio slope α, recover the electron index p = 2α + 1 — the sky as a particle spectrometer (curves computed from F(x) for this edition)

Limiting cases

caseresultmeaning
γ→1\gamma \to 1single line at ω0=eB/m\omega_0 = eB/m, broad sin⁡2\sin^2 patterncyclotron radiation — same orbit, no compression
mildly relativisticfew harmonics ω0,2ω0,3ω0\omega_0, 2\omega_0, 3\omega_0gyro-synchrotron; line-like spectrum
γ≫1\gamma \gg 1continuum to ωc=32γ3c/ρ\omega_c = \tfrac{3}{2}\gamma^3 c/\rho, cone 1/γ1/\gammatrue synchrotron; harmonics nc=32γ3n_c = \tfrac32\gamma^3 unresolvable
ω≪ωc\omega \ll \omega_cdP/dω∝ω1/3dP/d\omega \propto \omega^{1/3}, γ\gamma-independentuniversal soft tail; survives ensemble averaging
ω≫ωc\omega \gg \omega_cF→π/2 x1/2e−xF \to \sqrt{\pi/2}\,x^{1/2}e^{-x}exponential darkness beyond the knee
ensemble N∝E−pN \propto E^{-p}I∝ω−(p−1)/2I \propto \omega^{-(p-1)/2}the astrophysics spectrometer; aging adds Δα=1/2\Delta\alpha = 1/2

Worked examples

1 · A 3 GeV light source bend, end to end (confidence builder)

Setup. A 3.0 GeV electron storage ring with 1.4 T bending magnets — Diamond-class numbers. Compute every headline observable of the bending-magnet radiation.

Step 1. Kinematics. γ=E/mec2=3×109/0.511×106=5871\gamma = E/m_ec^2 = 3\times10^9/0.511\times10^6 = 5871; bending radius ρ=γmec/eB=7.15\rho = \gamma m_e c/eB = \mathbf{7.15} m.

Step 2. The critical energy, two routes. From (2): ωc=32γ3c/ρ=1.27×1019\omega_c = \tfrac{3}{2}\gamma^3c/\rho = 1.27\times10^{19} rad/s, so εc=ℏωc=8.38\varepsilon_c = \hbar\omega_c = \mathbf{8.38} keV. Engineering rule: 0.665×32×1.4=8.380.665 \times 3^2 \times 1.4 = 8.38 keV — identical, as it must be. Hard X-rays from a machine whose RF runs at hundreds of MHz.

Step 3. The geometry and the clock. Vertical cone 1/γ=1701/\gamma = \mathbf{170} μrad — the beam hugs the orbital plane so tightly that beamline slits sit kilometers of γψ\gamma\psi-space away at a few hundred microradians. Single-electron flash width Δtobs≈ρ/γ3c=0.12\Delta t_{\rm obs} \approx \rho/\gamma^3c = \mathbf{0.12} attoseconds; harmonic count nc=32γ3=3.0×1011n_c = \tfrac32\gamma^3 = 3.0\times10^{11} — the comb is a continuum for any real beam.

Step 4. Where the photons sit. Spectral peak at 0.29 εc=2.40.29\,\varepsilon_c = 2.4 keV; mean photon 0.31 εc=2.60.31\,\varepsilon_c = 2.6 keV; exactly half the wattage above 8.38 keV. The per-turn loss is Relativistic Generalization of Larmor's Formula's business — its 0.0885 E4/ρ0.0885\,E^4/\rho MeV rule gives 1.0 MeV/turn, the bill the RF pays — while this article's ledger says how that megaelectronvolt is distributed: ω1/3\omega^{1/3} rise through the visible, knee in the hard X-ray, darkness beyond ∼\sim25 keV.

The lesson: one γ\gamma and one BB fix the entire photon factory: εc∝γ2B\varepsilon_c \propto \gamma^2 B, cone ∝1/γ\propto 1/\gamma, and the universal FF does the rest.

2 · The Crab Nebula's X-rays date its electrons (exam level)

Setup. The Crab Nebula emits synchrotron X-rays at ε≈1\varepsilon \approx 1 keV in a nebular field B≈100 μG=10−8B \approx 100\ \mu\mathrm{G} = 10^{-8} T. What electrons does that require, how long do they survive, and what does the answer force you to conclude? (Radio slope: α≈0.30\alpha \approx 0.30.)

Step 1. Invert the critical frequency (field form of (2)): γ=2meωc/3eB\gamma = \sqrt{2m_e\omega_c/3eB} with ωc=ε/ℏ=1.52×1018\omega_c = \varepsilon/\hbar = 1.52\times10^{18} rad/s gives γ=2.4×107\gamma = \mathbf{2.4\times10^{7}} — electron energy γmec2≈12\gamma m_ec^2 \approx \mathbf{12} TeV, far beyond any accelerator on Earth, inferred from classical electrodynamics alone.

Step 2. The radiative lifetime. Power per electron (pitch 90°90°, the external-fields form recalled by name): P=e4γ2B2c2/6πε0me2c3=9.1×10−16P = e^4\gamma^2B^2c^2/6\pi\varepsilon_0m_e^2c^3 = 9.1\times10^{-16} W. Energy over power: τ=γmec2/P=2.2×109\tau = \gamma m_ec^2/P = 2.2\times10^9 s ≈70\approx \mathbf{70} years.

Step 3. The forced conclusion. The supernova was seen in 1054 — the nebula is ∼\sim970 years old, fourteen cooling times older than its X-ray electrons. That population cannot be a relic of the explosion; something inside — the pulsar wind's termination shock — must be accelerating electrons to 12 TeV today. Shklovsky's 1950s argument turned a brightness map into the discovery of ongoing cosmic acceleration, decades before the pulsar was found.

Step 4. Cross-check with the ensemble theorem. The radio slope α=0.30\alpha = 0.30 gives injection index p=1.6p = 1.6; radiative aging (Derivation C) steepens the X-ray slope by Δα≈0.5\Delta\alpha \approx 0.5 — and the observed ≈1.1\approx 1.1 obliges. The spectrum's shape independently confirms what the cooling time demanded.

The lesson: ωc\omega_c read backwards is a particle detector, the cooling time a clock — together they do astronomy no telescope could do alone.

Going further

References. Jackson, Classical Electrodynamics, 3rd ed., §14.4–14.6 · Griffiths, 4th ed., §11.2 · Zangwill, §23.6 · Purcell & Morin, Ch. 5 · Rybicki & Lightman, Radiative Processes in Astrophysics, Ch. 6 · Hofmann, The Physics of Synchrotron Radiation

Provenance. Power formulas: recompute Liénard's relativistic Larmor result and specialize to β˙⊥β\dot{\boldsymbol\beta}\perp\boldsymbol\beta — Jackson §14.2–14.4 eqs. (14.26),(14.31). Critical frequency ωc=32γ3c/ρ\omega_c=\tfrac32\gamma^3 c/\rho and Bessel-function spectrum: Jackson §14.6, Rybicki & Lightman Ch. 6. Universal function F(x)F(x), its x1/3x^{1/3} and x1/2e−xx^{1/2}e^{-x} limits, and normalization ∫F=8π/93\int F=8\pi/9\sqrt3: standard (Rybicki & Lightman §6.4). Electron loss constant 88.5 keV/(GeV4^4·m): verify from ΔE=q2γ4/(3ε0ρ)\Delta E=q^2\gamma^4/(3\varepsilon_0\rho). Astrophysical index α=(p−1)/2\alpha=(p-1)/2: Rybicki & Lightman §6.3. All numbers are standard textbook values; verifier should reproduce the γ\gamma-scalings and the constant from the boxed expressions. Addendum for this edition: the Liénard totals, covariant power, and per-turn constant (0.0885 MeV·GeV⁻⁴·m) are owned by the E9 Relativistic Larmor Formula (Liénard's Result) and Relativistic Generalization of Larmor's Formula, recalled by name; the lobes and the pulse-compression teaser are Angular Distribution of Radiation from Accelerated Charges's, cashed here. Two source-text corrections, verified: the source's Liénard-bracket eqCard carried γ4\gamma^4 where the general prefactor is γ6\gamma^6 (its own steps state γ6\gamma^6), and its spectrum card omitted the 93/8π9\sqrt3/8\pi normalization required for dP/dωdP/d\omega to integrate to PtotP_{\rm tot} — equation (4) restores it. Quadrature audits for this edition: peak 0.28580.2858 (Fmax⁡=0.918F_{\max}=0.918), half-power at x=1.0000x=1.0000, mean photon 0.3079 ωc=8/(153) ωc0.3079\,\omega_c = 8/(15\sqrt3)\,\omega_c, ∫F=1.61227\int F = 1.61227, Pσ/Pπ=7.000P_\sigma/P_\pi = 7.000 (75.0%); machine (γ=5871\gamma=5871, ρ=7.15\rho=7.15 m, εc=8.38\varepsilon_c=8.38 keV both routes, 0.12 as, nc=3.0×1011n_c=3.0\times10^{11}) and Crab (2.4×1072.4\times10^7, 12 TeV, τ≈70\tau\approx70 yr) numbers recomputed; SVG curves traced from computed F(x)F(x).

Builds on: Relativistic Generalization of Larmor's Formula · Fields of a Relativistic Moving Charge · Liénard–Wiechert Potentials: Covariant Formulation · Angular Distribution of Radiation from Accelerated Charges Leads to: Bremsstrahlung · Antenna Arrays and Directivity

If each flash is attoseconds long, why do beamlines quote picosecond X-ray pulses?

Because the beam is not one electron. The attosecond figure ρ/γ3c\rho/\gamma^3c is the flash of a single electron; a stored bunch spreads ∼1010\sim10^{10} electrons over millimeters, and their flashes arrive over the bunch's transit time — tens of picoseconds. The single-electron width still matters: it fixes the spectrum (each flash carries the full F(x)F(x) continuum to ωc\omega_c), while the bunch length fixes the timing resolution of pump-probe experiments. The scales decouple because electrons radiate incoherently at X-ray wavelengths. Free-electron lasers close exactly this gap: microbunching at the radiation wavelength makes the bunch radiate as one coherent macro-charge — femtosecond pulses with N2N^2 intensity, the coherence theorem of Larmor Formula from the Liénard–Wiechert Fields put to work.

Does the electron really emit at frequencies it never oscillates at? Where do the high harmonics physically come from?

Nothing oscillates at 101810^{18} Hz in the source — and nothing needs to. Fourier content is a property of the received field's time profile. The orbit modulates the field at ω0\omega_0, but the observer receives that modulation through the retarded-time map t=tr+R(tr)/ct = t_r + R(t_r)/c, whose derivative near the head-on instant is κ=1−β≈1/2γ2\kappa = 1-\beta \approx 1/2\gamma^2: a nonlinear time-warp that crams a smooth emission history into a spike, and spikes have broad spectra. In the frequency picture the received signal is strictly periodic, so its spectrum is a comb of harmonics nω0n\omega_0 — but the pulse shape populates the comb up to n∼32γ3n \sim \tfrac32\gamma^3, and finite energy spread blurs the teeth. High harmonics are the price of a sharp pulse, not evidence of hidden fast oscillators.

Why does the critical frequency scale as gamma cubed at fixed radius but only gamma squared at fixed field?

Because at fixed field the orbit itself grows with energy. The compression argument always delivers ωc∼γ3 (c/ρ)\omega_c \sim \gamma^3\,(c/\rho). At fixed ρ\rho (a bought tunnel), you keep all three powers. At fixed BB, the radius ρ=γmc/eB\rho = \gamma mc/eB stretches linearly with energy, so c/ρ∝1/γc/\rho \propto 1/\gamma eats one: ωc=32γ2eB/m\omega_c = \tfrac32\gamma^2 eB/m. Machine designers live in the fixed-ρ\rho world (BB ramps with energy), astrophysicists in the fixed-BB world (the nebula's field is what it is) — quoting the wrong scaling in the wrong context is a factor-of-γ\gamma error in disguise.

Is synchrotron radiation coherent? The formula says power per electron times N, but lasers get N squared.

Bending-magnet and wiggler radiation is incoherent at optical and X-ray wavelengths: electrons sit at random phases within the bunch, cross-terms average away, and powers add — N×N \times the single-electron spectrum. Coherence switches on when electrons cluster within a wavelength, making fields add before squaring: N2N^2. That happens naturally at wavelengths longer than the bunch (millimeter/THz "coherent synchrotron radiation," which can destabilize short bunches) and artificially in free-electron lasers, which microbunch the beam at the X-ray wavelength, converting 10910^9 electrons into one coherent super-electron. The dividing line is always wavelength versus charge-distribution scale — the N2N^2-versus-NN theorem owned by Larmor Formula from the Liénard–Wiechert Fields.

Why is the in-plane radiation 100 percent polarized but the total only 75 percent?

Geometry, then bookkeeping. An in-plane observer sees the circular motion edge-on: the acceleration projects onto the sky as a purely horizontal line — perfect linear polarization (the π\pi-mode amplitude carries an explicit ψ\psi and dies at ψ=0\psi = 0). An observer above the plane sees a slightly opened ellipse: both components present with a fixed phase — elliptical polarization, handedness set by the side. Integrating over the cone mixes the viewpoints; quadrature of the two Bessel terms gives Pσ:Pπ=7:1P_\sigma : P_\pi = 7 : 1, degree 75%. Radio astronomers rarely see even that: field tangling and Faraday rotation depolarize further, so a measured 40% polarization is itself evidence of a remarkably ordered field.

What happens when the synchrotron photon energy approaches the electron energy itself?

Classical electrodynamics quietly assumes ℏωc≪γmec2\hbar\omega_c \ll \gamma m_ec^2 — the electron barely notices emitting any single photon. The ratio χ∝γB\chi \propto \gamma B grows with energy and field; near unity the emission turns quantum: discrete recoil kicks (quantum excitation — already a design driver for damping rings at χ∼10−3\chi \sim 10^{-3}), a truncated spectrum, eventually pair creation. In parallel, classically, there is the radiation-reaction ceiling: balancing the accelerating force against the drag of Relativistic Generalization of Larmor's Formula caps synchrotron photons near ∼160\sim160 MeV regardless of BB — a bound the Crab's gamma-ray flares violate, which is why they electrify theorists. Both frontiers — quantum synchrotron and self-force — are exactly where this chapter's closing articles and the QED bridge take over, by name.

Before this

Where this leads