Relativistic Generalization of Larmor's Formula

The strong-UG article Relativistic Larmor Formula (Liénard's Result) earned the exact any-speed power twice — once through the κ5\kappa^5 retardation bookkeeping, once through the rest-frame invariance argument — and read off the doctrine: γ6\gamma^6 for linear geometry, γ4\gamma^4 for circular, LEP's ruinous bill versus the linac loophole. This Jackson-register companion completes that result in the three directions graduate work actually uses. First, the identity. That article verified numerically that the four-acceleration's invariant square reproduces Liénard's bracket; here we derive it symbolically — components of aμa^\mu from first principles, the contraction, and the Lagrange-identity step that converts the (β⋅β˙)(\boldsymbol\beta\cdot\dot{\boldsymbol\beta}) form into the (β×β˙)(\boldsymbol\beta\times\dot{\boldsymbol\beta}) form exactly, at every speed. The two spellings of aμaμa_\mu a^\mu are a theorem, not a large-β\beta accident. Second, the fields. Substituting the covariant Lorentz force into the covariant power formula eliminates kinematics altogether: the radiated power of a charge in given external E\mathbf{E} and B\mathbf{B} fields, at given velocity, in one line — no trajectory solved, no acceleration computed. This is the form accelerator and plasma physicists actually evaluate, and its two limits reproduce the linear and circular verdicts as corollaries. Third, the engineering. The per-turn loss δE=q2β3γ4/3ε0ρ\delta E=q^2\beta^3\gamma^4/3\varepsilon_0\rho becomes the working formula δE[MeV]≈0.0885 E4[GeV]/ρ[m]\delta E[\mathrm{MeV}]\approx0.0885\,E^4[\mathrm{GeV}]/\rho[\mathrm{m}] — derived, not quoted — and priced on the machines being designed now: the FCC-ee Z-pole ring and a synchrotron light source, with the radiation-damping timescale read off as the bonus that makes electron rings self-cooling.

Key takeaways
  • Four-acceleration components, derived: a0=γ4(β⋅β˙)ca^0=\gamma^4(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})c, a(3)=γ2c β˙+γ4c (β⋅β˙)β\quad\mathbf{a}_{(3)}=\gamma^2c\,\dot{\boldsymbol\beta}+\gamma^4c\,(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})\boldsymbol\beta.
  • The invariant square in two provably equal spellings: aμaμ=−γ4c2[β˙2+γ2(β⋅β˙)2]=−γ6c2[β˙2−(β×β˙)2]a_\mu a^\mu=-\gamma^4c^2\big[\dot{\boldsymbol\beta}^2+\gamma^2(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2\big]=-\gamma^6c^2\big[\dot{\boldsymbol\beta}^2-(\boldsymbol\beta\times\dot{\boldsymbol\beta})^2\big] — the bridge is the Lagrange identity (β×β˙)2=β2β˙2−(β⋅β˙)2(\boldsymbol\beta\times\dot{\boldsymbol\beta})^2=\beta^2\dot{\boldsymbol\beta}^2-(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2, exact at every β\beta.
  • External-fields form: substituting dpμ/dτ=qFμνuνdp^\mu/d\tau=qF^{\mu\nu}u_\nu gives P=q4γ26πε0m2c3[(E+v×B)2−(β⋅E)2]P=\dfrac{q^4\gamma^2}{6\pi\varepsilon_0m^2c^3}\Big[(\mathbf{E}+\mathbf{v}\times\mathbf{B})^2-(\boldsymbol\beta\cdot\mathbf{E})^2\Big] — power from the fields at the particle, no trajectory required.
  • Limits audit themselves: pure B\mathbf{B} reproduces the circular γ4/ρ2\gamma^4/\rho^2 law; pure E∥v\mathbf{E}\parallel\mathbf{v} gives a γ\gamma-free power at fixed force — the linac loophole, now one substitution deep.
  • Engineering form, derived: δE[MeV]≈0.0885 E4[GeV]/ρ[m]\delta E[\mathrm{MeV}]\approx0.0885\,E^4[\mathrm{GeV}]/\rho[\mathrm{m}] per turn; FCC-ee at the Z pole (45.645.6 GeV, ρ≈10.8\rho\approx10.8 km): ≈36\approx36 MeV/turn; a 6 GeV light source (ρ=23.4\rho=23.4 m): 4.94.9 MeV/turn.
  • Radiation damping: τ∼E/P\tau\sim E/P gives ∼0.4\sim0.4 s (FCC-ee Z) and ∼3.5\sim3.5 ms (light source) — the self-cooling that makes electron storage rings possible.
  • Key idea: the covariant Larmor formula becomes computational three ways — the symbolic aμaμa_\mu a^\mu identity, the field-substitution form that needs no trajectory, and the 0.0885 E4/ρ0.0885\,E^4/\rho engineering spelling that designs real machines.
The core insight

Radiated power is the invariant −μ0q2(aμaμ)/6πc-\mu_0q^2(a_\mu a^\mu)/6\pi c made computational: derive aμaμa_\mu a^\mu symbolically in both spellings, substitute the covariant Lorentz force to get power directly from the external fields, and reduce the circular case to the 0.0885 E4/ρ0.0885\,E^4/\rho per-turn formula that prices accelerators.

Where this lives today

The field form (2) is the radiation module of modern plasma-accelerator and strong-field codes: in a laser wakefield the fields at each macro-particle are known on the grid, and (2) prices the betatron X-ray emission that has become these machines' signature diagnostic — the same formula, evaluated 10910^9 times per step, no trajectories differentiated. Its next-order refinement is a live research frontier: at today's petawatt intensities the radiated momentum begins to react back on the electron within one laser cycle, and "radiation reaction observed in laser–electron collisions" experiments are precisely the regime where (2) must be upgraded to the self-consistent equations of the radiation-reaction articles ahead. The engineering form runs the other industry. Every storage-ring design report opens with 0.0885E4/ρ0.0885E^4/\rho: it sets FCC-ee's RF budget (nine-plus GeV per turn at ttˉt\bar t), fixed the diameter-versus-cost tradeoff of every light source since the first dedicated rings, and — run for protons with the (me/mp)4(m_e/m_p)^4 suppression — explains why the 100 TeV successor FCC-hh can reuse the same tunnel with radiation as a cryogenic nuisance (synchrotron light absorbed at 50 K) rather than an energy wall. Between the two forms, this article's algebra is arguably the most expensive formula in experimental physics: machines are sized, sited, and cancelled by its fourth power.

Origins · 1 minJulian Schwinger and how this was discovered
JS
Julian Schwinger
1918–1994 · who turned Liénard's invariant into the accelerator designer's handbook

Schwinger came to synchrotron radiation sideways, through wartime radar: at the MIT Radiation Laboratory he worked out how electrons in the new betatrons and synchrotrons would radiate, circulating the results in a famous 1945 manuscript years before the polished paper, "On the Classical Radiation of Accelerated Electrons," appeared in 1949. That paper made the covariant power formula useful: the per-turn energy loss in practical units, the proof that the loss grows as the fourth power of energy — sharpening the Iwanenko–Pomeranchuk warning that had seemed to doom circular machines — and the complete spectral and angular structure of the emission (whose critical-frequency story this track defers to the synchrotron article). General Electric's 1947 observation of visible light from their 70 MeV synchrotron confirmed the analysis while the manuscript still circulated by mail. Schwinger's later fame is quantum — QED's Green's functions, the anomalous magnetic moment, the 1965 Nobel shared with Feynman and Tomonaga — but accelerator physicists still work from the classical paper: every "0.0885 E4/ρ0.0885\,E^4/\rho" in a design report is a line item he priced first.

History
1898
— Liénard publishes the arbitrary-velocity power formula, γ factors intact, seven years before relativity names them
1908
— Minkowski's spacetime formalism supplies the four-acceleration that makes Liénard's result a one-line invariant
1944
— Iwanenko and Pomeranchuk predict the fourth-power energy ceiling for circular electron machines
1946
— Blewett measures the predicted orbit shrinkage from radiation loss at the General Electric betatron
1949
— Schwinger's classical-radiation paper: the covariant treatment and the practical per-turn formulas
1961
— AdA at Frascati, the first electron–positron storage ring: the γ⁴ budget becomes routine engineering
2019
— the FCC-ee conceptual design prices its Z-pole run at 36 MeV per turn — this article's arithmetic, at hundred-megawatt scale

The physical picture

The strong-UG treatment established what is true: total radiated power is a Lorentz scalar, computable in the cheapest frame, exploding as γ6\gamma^6 or γ4\gamma^4 depending on geometry. This article is about how the scalar computes. Three pictures organize it. First, the four-acceleration as a bookkeeping machine. Differentiating uμ=γ(c,v)u^\mu=\gamma(c,\mathbf{v}) by proper time — the only covariant clock, per the chapter's ladder principle in Proper Time and Four-Velocity — produces components in which every power of γ\gamma has a job: one γ\gamma from d/dτ=γ d/dtd/d\tau=\gamma\,d/dt, more from γ˙=γ3(β⋅β˙)\dot\gamma=\gamma^3(\boldsymbol\beta\cdot\dot{\boldsymbol\beta}). The invariant square then carries γ4\gamma^4 on the isotropic piece and γ6\gamma^6 on the longitudinal piece — the geometry dependence of radiated power is nothing but the anatomy of aμa^\mu. Second, the force shortcut. Why solve for a trajectory at all? The covariant power formula wants dpμ/dτdp^\mu/d\tau, and for a charge in external fields the equation of motion hands that over directly: dpμ/dτ=qFμνuνdp^\mu/d\tau=qF^{\mu\nu}u_\nu, the covariant Lorentz force whose construction is The Electromagnetic Field Tensor's destination result. Substituting converts "solve the motion, differentiate, contract" into "read the fields at the particle" — radiation as a local function of (E,B,β)(\mathbf{E},\mathbf{B},\boldsymbol\beta). Third, the ring as a radiation engine. On a circular arc the force is fixed by geometry — bend radius ρ\rho demands transverse acceleration whatever the energy — so the γ4\gamma^4 has nowhere to hide, and per-turn loss grows as the fourth power of beam energy. The same formula read constructively: the beam constantly converts stored energy into light and re-absorbs the loss from RF cavities that push only longitudinally — a damping cycle that shrinks oscillations and cools the beam, on a timescale the formula also prices.

Figure
The fourth-power tax, computed — per-turn loss versus beam energy for an FCC-ee-scale ring, log-log
δE per turn = 0.0885 E⁴/ρ (electrons, ρ = 10.76 km) — slope 4 on log-log axes beam energy E (GeV), log scale 20 45.6 80 120 182.5 1 MeV 100 MeV 10 GeV Z pole 45.6 GeV: 36 MeV/turn WW 80 GeV: 0.34 GeV tt̄ 182.5 GeV: 9.1 GeV/turn double the energy → sixteen times the bill computed from the engineering formula derived below; the tt̄ point is why even a 91 km tunnel is the last circular electron machine
Figure
The anatomy of the four-acceleration — where each power of gamma is born and which invariant it feeds
a^μ = du^μ/dτ, unpacked time component a⁰ = γ⁴(β·β̇)c exists only when speed changes: γ from d/dτ, γ³ from γ̇ = γ³(β·β̇) space components a⃗ = γ²c β̇ + γ⁴c (β·β̇)β a transverse piece every motion has, plus a longitudinal piece riding on β the contraction (metric +,−,−,−; a·a spacelike, negative) a·a = −γ⁴c²[ β̇² + γ²(β·β̇)² ] ≡ −γ⁶c²[ β̇² − (β×β̇)² ] bridge: Lagrange identity (β×β̇)² = β²β̇² − (β·β̇)² — exact at every β, not an ultrarelativistic shortcut ⊥ geometry: β·β̇ = 0 → only γ⁴ survives (synchrotron) ∥ geometry: β×β̇ = 0 → full γ⁶ (linear) P = −μ₀q²(a·a)/6πc turns this anatomy directly into the Liénard bracket — Derivation A executes every arrow

The math: three completions of the invariant formula

The starting point is quoted, not re-derived: the strong-UG Relativistic Larmor Formula (Liénard's Result) proved — by rest-frame isotropy and separately by the κ5\kappa^5 integration — that

P=−μ0q26πc aμaμ=−μ0q26πm2c dpμdτdpμdτ,(1)P=-\frac{\mu_0q^2}{6\pi c}\,a_\mu a^\mu=-\frac{\mu_0q^2}{6\pi m^2c}\,\frac{dp^\mu}{d\tau}\frac{dp_\mu}{d\tau}, \qquad (1)

where: aμ=duμ/dτa^\mu=du^\mu/d\tau is the four-acceleration, pμ=muμp^\mu=mu^\mu, and the metric is (+,−,−,−)(+,-,-,-), locked across this track and its mechanics companions. Since aμa^\mu is orthogonal to the timelike uμu^\mu it is spacelike, aμaμ<0a_\mu a^\mu<0, and the minus sign makes PP positive. Everything in this article is (1) being made computational.

Derivation AThe invariant square of the four-acceleration, symbolically

The goal is the exact identity behind Liénard's bracket: compute aμaμa_\mu a^\mu in lab variables from first principles and prove the two standard spellings equal — the (β⋅β˙)(\boldsymbol\beta\cdot\dot{\boldsymbol\beta}) form that falls out of the algebra and the (β×β˙)(\boldsymbol\beta\times\dot{\boldsymbol\beta}) form Liénard wrote — at every speed. The strong-UG twin checked this numerically for random inputs; upgrading the check to a symbolic theorem is worth doing once in full, because the intermediate objects (the components of aμa^\mu, the γ˙\dot\gamma lemma) are reused by every covariant radiation calculation downstream. Strategy: differentiate uμu^\mu by proper time via the chain rule, contract, then apply the Lagrange identity. Watch one thing: keep the factors of cc explicit — the commonest error in this computation is a cc dropped between β˙\dot{\boldsymbol\beta} (units 1/s) and acceleration (m/s²).

Step 1. The γ˙\dot\gamma lemma. From γ=(1−β⋅β)−1/2\gamma=(1-\boldsymbol\beta\cdot\boldsymbol\beta)^{-1/2},

dγdt=γ3 (β⋅β˙).\frac{d\gamma}{dt}=\gamma^3\,(\boldsymbol\beta\cdot\dot{\boldsymbol\beta}) .

Why this lemma is the whole story: every γ\gamma beyond the first two in the final answer descends from this cubic — speeding up costs three powers of γ\gamma per unit of β⋅β˙\boldsymbol\beta\cdot\dot{\boldsymbol\beta}, which is why longitudinal dynamics is the expensive direction.

Step 2. Components of aμa^\mu. With uμ=γ(c,cβ)u^\mu=\gamma(c,c\boldsymbol\beta) and d/dτ=γ d/dtd/d\tau=\gamma\,d/dt:

a0=γd(γc)dt=γ4(β⋅β˙) c,a(3)=γd(γcβ)dt=γ2c β˙+γ4c (β⋅β˙) β.a^0=\gamma\frac{d(\gamma c)}{dt}=\gamma^4(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})\,c,\qquad \mathbf{a}_{(3)}=\gamma\frac{d(\gamma c\boldsymbol\beta)}{dt}=\gamma^2c\,\dot{\boldsymbol\beta}+\gamma^4c\,(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})\,\boldsymbol\beta .

Why the space part has two pieces: the first is the honest acceleration, boosted twice; the second rides along β\boldsymbol\beta and exists only when the speed changes — for circular motion it vanishes identically, which is where the γ6→γ4\gamma^6\to\gamma^4 demotion will come from.

Step 3. Contract. With the metric (+,−,−,−)(+,-,-,-), aμaμ=(a0)2−∣a(3)∣2a_\mu a^\mu=(a^0)^2-|\mathbf{a}_{(3)}|^2. Expanding ∣a(3)∣2|\mathbf{a}_{(3)}|^2 gives γ4c2β˙2+2γ6c2(β⋅β˙)2+γ8c2(β⋅β˙)2β2\gamma^4c^2\dot{\boldsymbol\beta}^2+2\gamma^6c^2(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2+\gamma^8c^2(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2\beta^2; subtracting from (a0)2=γ8c2(β⋅β˙)2(a^0)^2=\gamma^8c^2(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2 and using γ8(1−β2)=γ6\gamma^8(1-\beta^2)=\gamma^6:

aμaμ=−γ4c2[β˙2+γ2(β⋅β˙)2].a_\mu a^\mu=-\gamma^4c^2\Big[\dot{\boldsymbol\beta}^2+\gamma^2(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2\Big] .

Step 4. The Lagrange bridge. For any two vectors, (β×β˙)2=β2β˙2−(β⋅β˙)2(\boldsymbol\beta\times\dot{\boldsymbol\beta})^2=\beta^2\dot{\boldsymbol\beta}^2-(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2. Substituting (β⋅β˙)2=β2β˙2−(β×β˙)2(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2=\beta^2\dot{\boldsymbol\beta}^2-(\boldsymbol\beta\times\dot{\boldsymbol\beta})^2 into Step 3 and collecting with 1+γ2β2=γ21+\gamma^2\beta^2=\gamma^2:

aμaμ=−γ6c2[β˙2−(β×β˙)2].a_\mu a^\mu=-\gamma^6c^2\Big[\dot{\boldsymbol\beta}^2-(\boldsymbol\beta\times\dot{\boldsymbol\beta})^2\Big] .

The two spellings are identically equal at every β\beta — the base article's caveat that the conversion "holds only when β2=1\beta^2=1" is false, and the correction is logged in the provenance. Substituting into (1) gives Liénard's lab form

P=μ0q2γ66πc[β˙2−(β×β˙)2].P=\frac{\mu_0q^2\gamma^6}{6\pi c}\Big[\dot{\boldsymbol\beta}^2-(\boldsymbol\beta\times\dot{\boldsymbol\beta})^2\Big].

Step 5. Audit the geometry limits from the first spelling. Circular (β⋅β˙=0\boldsymbol\beta\cdot\dot{\boldsymbol\beta}=0): aμaμ=−γ4c2β˙2a_\mu a^\mu=-\gamma^4c^2\dot{\boldsymbol\beta}^2 — the γ4\gamma^4 law with no Lagrange gymnastics needed. Linear (β×β˙=0\boldsymbol\beta\times\dot{\boldsymbol\beta}=0, so (β⋅β˙)2=β2β˙2(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})^2=\beta^2\dot\beta^2): the bracket is β˙2(1+γ2β2)=γ2β˙2\dot\beta^2(1+\gamma^2\beta^2)=\gamma^2\dot\beta^2, giving −γ6c2β˙2-\gamma^6c^2\dot\beta^2 — the γ6\gamma^6 law. The doctrine of the strong-UG twin (what is held fixed, which exponent wins) applies unchanged and is not repeated here. Verified: both component expressions and both spellings of aμaμa_\mu a^\mu checked symbolically with a computer-algebra system for this edition (exact zero difference), and the equality of the two spellings confirmed — matching Jackson Eqs. 14.25–14.26 and the mechanics companions' four-vector machinery with the same (+,−,−,−)(+,-,-,-) signature.

Figure
Power straight from the fields — the no-trajectory machine, with its two limit audits
Radiated power without solving the motion inputs, at the particle E, B fields · velocity β no trajectory, no history covariant Lorentz force dp^μ/dτ = qF^μν u_ν the field tensor's destination result output P(E, B, β) one line, eq. (2) audit 1: pure B, v ⊥ B P = q⁴γ²v²B²/6πε₀m²c³ ≡ the γ⁴/ρ² synchrotron law with ρ = γmv/qB (two powers of γ hide inside ρ) audit 2: pure E ∥ v bracket: E² − β²E² = E²/γ² → P = q⁴E²/6πε₀m²c³ γ²(1−β²) = 1: energy-independent at fixed gradient — the linac loophole, one substitution deep same force, two geometries: the bend radiates γ² more than the push — worked example 2 prices it at γ = 10
Go deeper · Derivation B · 2 minRadiated power directly from the external fields
Derivation BRadiated power directly from the external fields

The goal is the working form of the covariant Larmor formula: eliminate the kinematics entirely by feeding the equation of motion into (1), obtaining PP as a pointwise function of the external fields and the particle's velocity. This is the version evaluated in accelerator lattices and plasma codes, where fields are known everywhere and trajectories are expensive. Strategy: write the covariant Lorentz force, contract it with itself, and translate the invariant into (E,B,β)(\mathbf{E},\mathbf{B},\boldsymbol\beta) language; then audit the two classic limits. The one subtlety: the contraction produces a difference of squares — the longitudinal electric force partially escapes the radiation bill, and that minus sign is physics, not bookkeeping.

Step 1. The force. For a charge in external fields, the covariant equation of motion is dpμ/dτ=qFμνuνdp^\mu/d\tau=qF^{\mu\nu}u_\nu — quoted by name from The Electromagnetic Field Tensor, whose construction it crowns. Its components, using uμ=γ(c,v)u^\mu=\gamma(c,\mathbf{v}):

dp0dτ=γqc E⋅v,dpdτ=γq(E+v×B).\frac{dp^0}{d\tau}=\frac{\gamma q}{c}\,\mathbf{E}\cdot\mathbf{v},\qquad \frac{d\mathbf{p}}{d\tau}=\gamma q\big(\mathbf{E}+\mathbf{v}\times\mathbf{B}\big) .

where: the time component is the work rate (energy gain per proper time, over cc) and the space part is the familiar Lorentz force, each carrying one γ\gamma from the proper-time clock.

Step 2. Contract. With (+,−,−,−)(+,-,-,-):

dpμdτdpμdτ=γ2q2[(β⋅E)2−(E+v×B)2].\frac{dp^\mu}{d\tau}\frac{dp_\mu}{d\tau}=\gamma^2q^2\Big[(\boldsymbol\beta\cdot\mathbf{E})^2-\big(\mathbf{E}+\mathbf{v}\times\mathbf{B}\big)^2\Big] .

Step 3. Substitute into (1):

P=q4γ26πε0m2c3[(E+v×B)2−(β⋅E)2](2)P=\frac{q^4\gamma^2}{6\pi\varepsilon_0m^2c^3}\Big[\big(\mathbf{E}+\mathbf{v}\times\mathbf{B}\big)^2-\big(\boldsymbol\beta\cdot\mathbf{E}\big)^2\Big] \qquad (2)

where: every quantity is evaluated at the particle's instantaneous position and velocity — no trajectory has been solved, no acceleration computed. The first term is the full Lorentz force per unit charge, squared; the subtracted (β⋅E)2(\boldsymbol\beta\cdot\mathbf{E})^2 is the longitudinal-work discount: the part of the electric force that changes energy rather than direction radiates less, by exactly the amount that will produce the linac loophole below. Why a γ2\gamma^2 and not γ4\gamma^4 or γ6\gamma^6 sits outside: the bracket holds the force (bounded by the hardware), and the exponent bookkeeping of Derivation A has been absorbed into using τ\tau-derivatives of pμp^\mu — geometry now decides how much of the bracket survives, not how many γ\gamma's multiply it.

Step 4. Audit: pure magnetic field, v⊥B\mathbf{v}\perp\mathbf{B}, E=0\mathbf{E}=0. The bracket is v2B2v^2B^2, so P=q4γ2v2B2/6πε0m2c3P=q^4\gamma^2v^2B^2/6\pi\varepsilon_0m^2c^3. Eliminating BB via the bend radius ρ=γmv/qB\rho=\gamma mv/qB reproduces the circular law P=μ0q2c3β4γ4/6πρ2P=\mu_0q^2c^3\beta^4\gamma^4/6\pi\rho^2 exactly — two of the four γ\gamma's were hiding inside ρ\rho. (Numerically confirmed for this edition at γ=10\gamma=10, B=1B=1 T: both routes give 1.571×10−121.571\times10^{-12} W to twelve digits.)

Step 5. Audit: pure electric field along the motion. The bracket is E2−β2E2=E2/γ2E^2-\beta^2E^2=E^2/\gamma^2, so P=q4E2/6πε0m2c3P=q^4E^2/6\pi\varepsilon_0m^2c^3 — all γ\gamma dependence cancels: at fixed gradient the radiated power is the same at 1 MeV and 1 TeV. This is the linac loophole of the strong-UG twin (its LEP-versus-linac verdicts recalled by name, not recomputed), now visible as a one-line consequence of the minus sign in (2). Verified: equation (2) matches Jackson Eq. 14.13-form (Gaussian → SI translated); the general (2) was checked numerically against the kinematic Liénard form for a random (E,B,β)(\mathbf{E},\mathbf{B},\boldsymbol\beta) configuration by integrating the equation of motion for the acceleration — agreement to machine precision (ratio 1.000000000000) — and the pure-B/pure-E audits reduce to the standard circular and linear results.

Go deeper · Derivation C · 2 minAccelerator systematics — the per-turn loss and its engineering form
Derivation CAccelerator systematics — the per-turn loss and its engineering form

The goal is the formula that designs machines: the energy an electron loses per revolution in a ring of bend radius ρ\rho, first exactly, then in the practical units every accelerator report uses, δE[MeV]=0.0885 E4[GeV]/ρ[m]\delta E[\mathrm{MeV}]=0.0885\,E^4[\mathrm{GeV}]/\rho[\mathrm{m}] — with the 0.08850.0885 derived, not folklore. Strategy: specialize the circular law, multiply by the time spent bending, then translate to engineering units, keeping honest about what "per turn" means when the ring is not all dipoles. New anchors throughout — the strong-UG twin owns the LEP and LHC numbers.

Step 1. Power on the arc. For ultrarelativistic circular motion (β→1\beta\to1), from Derivation B Step 4 or Derivation A Step 5:

P⊥=μ0q2c3β4γ46πρ2=q2c β4γ46πε0ρ2.P_\perp=\frac{\mu_0q^2c^3\beta^4\gamma^4}{6\pi\rho^2}=\frac{q^2c\,\beta^4\gamma^4}{6\pi\varepsilon_0\rho^2} .

Step 2. Loss per turn. The electron bends only while in dipoles; the total time spent on arcs of radius ρ\rho per revolution is 2πρ/βc2\pi\rho/\beta c regardless of how much straight section the ring also has (the straights radiate nothing). So

δE=P⊥⋅2πρβc=q2β3γ43ε0ρ.\delta E=P_\perp\cdot\frac{2\pi\rho}{\beta c}=\frac{q^2\beta^3\gamma^4}{3\varepsilon_0\rho} .

Why this bookkeeping survives real lattices: the loss depends only on the bend radius and the total bend angle (2π2\pi), so distributing the dipoles around a much larger circumference changes the revolution time but not the per-turn bill — which is why the formula needs ρ\rho, the dipole radius, and not the ring's circumference-derived "average radius." Using the latter is the standard way to underestimate a machine's radiation budget.

Step 3. Engineering units. Insert γ=E/mec2\gamma=E/m_ec^2 with mec2=0.511m_ec^2=0.511 MeV, and express δE\delta E in MeV, EE in GeV, ρ\rho in meters:

δE [MeV]=e23ε0 (mec2)4 E4ρ=0.0885 E4 [GeV4]ρ [m]\delta E\ [\mathrm{MeV}]=\frac{e^2}{3\varepsilon_0\,(m_ec^2)^4}\,\frac{E^4}{\rho}=0.0885\,\frac{E^4\,[\mathrm{GeV}^4]}{\rho\,[\mathrm{m}]}

— the coefficient evaluates to 0.088460.08846 MeV·m·GeV−4^{-4} (recomputed for this edition from CODATA constants; the mass to the fourth power in the denominator is why this number is 1.1×10131.1\times10^{13} times smaller for protons).

Step 4. Price two real machines. FCC-ee at the Z pole: E=45.6E=45.6 GeV, effective dipole radius ρ≈10.76\rho\approx10.76 km in the ∼91\sim91 km tunnel — δE=0.0885×45.64/10760=35.5\delta E=0.0885\times45.6^4/10760=35.5 MeV per turn (≈36\approx36 MeV, matching the design figure), rising to 0.340.34 GeV at the WW point (80 GeV) and 9.19.1 GeV per turn at the ttˉt\bar t point (182.5 GeV): the fourth power in action, and the reason FCC-ee is the last circular electron collider anyone plans. A 6 GeV light source with ρ=23.4\rho=23.4 m dipoles: δE=0.0885×1296/23.4=4.9\delta E=0.0885\times1296/23.4=4.9 MeV per turn — three orders of magnitude more fractional loss per meter than the collider, by design: here the radiation is the product.

Step 5. The damping bonus. The beam loses energy in the bends (momentum kick opposite v\mathbf{v}, shrinking all components) and regains it in RF cavities (kick purely longitudinal): transverse oscillations damp on the timescale τ∼E/⟨P⟩=E T0/δE\tau\sim E/\langle P\rangle=E\,T_0/\delta E. FCC-ee Z: T0=326 μT_0=326\ \mus, τ≈45.6 GeV×T0/35.5 MeV≈0.42\tau\approx45.6\ \mathrm{GeV}\times T_0/35.5\ \mathrm{MeV}\approx0.42 s (∼\sim1300 turns). The light source: T0=2.8 μT_0=2.8\ \mus, τ≈3.5\tau\approx3.5 ms. This self-cooling is why electron rings have exquisitely small beams — and why proton rings, with 101310^{13} times weaker radiation, must be cooled by other means or not at all. (The exact partition of damping among the three planes is damping-partition theory — accelerator-physics territory beyond this track.) Verified: the 0.08850.0885 coefficient recomputed from constants (0.0884630.088463); FCC-ee numbers cross-checked at three energies against the E4E^4 scaling (ratios 16.016.0 for EE-doubling); the light-source 4.904.90 MeV agrees with the classic 88.5 E4/ρ88.5\,E^4/\rho keV rule; damping estimates E T0/δEE\,T_0/\delta E evaluate to 0.4180.418 s and 3.53.5 ms from the quoted inputs.

Limiting cases

caseresultmeaning
β→0\beta\to0aμaμ→−c2β˙2=−a2a_\mu a^\mu\to-c^2\dot{\boldsymbol\beta}^2=-a^2, P→P\to Larmorthe rest-frame anchor; bookkeeping article's momentum form
circular (β⋅β˙=0\boldsymbol\beta\cdot\dot{\boldsymbol\beta}=0)a⋅a=−γ4c2β˙2a\cdot a=-\gamma^4c^2\dot\beta^2γ4\gamma^4 read off the first spelling, no identity needed
linear (β×β˙=0\boldsymbol\beta\times\dot{\boldsymbol\beta}=0)a⋅a=−γ6c2β˙2a\cdot a=-\gamma^6c^2\dot\beta^2the full γ6\gamma^6; tame in force variables
pure B\mathbf{B}, v⊥B\mathbf{v}\perp\mathbf{B}P=q4γ2v2B26πε0m2c3P=\dfrac{q^4\gamma^2v^2B^2}{6\pi\varepsilon_0m^2c^3}field form; ≡γ4/ρ2\equiv\gamma^4/\rho^2 law with ρ=γmv/qB\rho=\gamma mv/qB
pure E∥v\mathbf{E}\parallel\mathbf{v}P=q4E26πε0m2c3P=\dfrac{q^4E^2}{6\pi\varepsilon_0m^2c^3}, γ\gamma-freelinac loophole as a one-line corollary of (2)
crossed fields, E+v×B=0\mathbf{E}+\mathbf{v}\times\mathbf{B}=0P=0P=0 exactly (balance forces β⋅E=0\boldsymbol\beta\cdot\mathbf{E}=0 too)a drift-balanced charge does not radiate: force, not fields, is what's priced
per-turn, electronsδE[MeV]=0.0885 E4[GeV]/ρ[m]\delta E[\mathrm{MeV}]=0.0885\,E^4[\mathrm{GeV}]/\rho[\mathrm{m}]the machine-design formula, derived in full
same EE, protonscoefficient ×(me/mp)4\times(m_e/m_p)^41.1×10131.1\times10^{13} suppression — the strong-UG twin's strategic verdict

Worked examples

1 · Pricing the FCC-ee Z-pole run end to end (confidence builder)

Setup. FCC-ee at the Z pole: beam energy E=45.6E=45.6 GeV, dipole bend radius ρ=10.76\rho=10.76 km, circumference C=97.75C=97.75 km. Find the per-turn loss, the average radiated power per electron, and the radiation-damping timescale.

Step 1. Lorentz factor: γ=45.6 GeV/0.511 MeV=8.92×104\gamma=45.6\ \mathrm{GeV}/0.511\ \mathrm{MeV}=8.92\times10^{4}, so γ4=6.34×1019\gamma^4=6.34\times10^{19} and β=1\beta=1 to ten digits.

Step 2. Per-turn loss from the engineering formula:

δE=0.0885 45.6410760 MeV=0.0885×4.32×1061.076×104=35.5 MeV per turn\delta E=0.0885\,\frac{45.6^4}{10760}\ \mathrm{MeV}=0.0885\times\frac{4.32\times10^{6}}{1.076\times10^{4}}=\mathbf{35.5\ MeV\ per\ turn}

— about 0.08%0.08\% of the beam energy every turn, replaced by the RF system; at the design current this is a ∼\sim100 MW wall-plug item, the single largest operating cost of the machine.

Step 3. Average power per electron: T0=C/c=326 μT_0=C/c=326\ \mus, so ⟨P⟩=δE/T0=35.5 MeV×1.602×10−13 J/MeV/3.26×10−4 s=1.7×10−8 W\langle P\rangle=\delta E/T_0=35.5\ \mathrm{MeV}\times1.602\times10^{-13}\ \mathrm{J/MeV}/3.26\times10^{-4}\ \mathrm{s}=\mathbf{1.7\times10^{-8}\ W} — seventeen nanowatts, delivered as keV-scale photons into the 1/γ≈11 μ1/\gamma\approx11\ \murad cone the strong-UG twin's beaming ledger describes.

Step 4. Damping: τ∼E/⟨P⟩=E T0/δE=45.6×103 MeV×326 μs/35.5 MeV=0.42 s\tau\sim E/\langle P\rangle=E\,T_0/\delta E=45.6\times10^{3}\ \mathrm{MeV}\times326\ \mu\mathrm{s}/35.5\ \mathrm{MeV}=\mathbf{0.42\ s}, about 1300 turns — fast enough that injected oscillations shrink within a second, the self-cooling that gives FCC-ee its luminosity. At the ttˉt\bar t point the loss is 9.19.1 GeV/turn (E4E^4: (182.5/45.6)4=256(182.5/45.6)^4=256 times more) and damping takes just a handful of turns — but the RF must then restore ten percent of the beam energy every revolution, which is the wall the fourth power finally builds.

The lesson: one derived coefficient, 0.08850.0885, prices a ten-billion-euro machine's energy budget, its cooling time, and its energy ceiling.

2 · Same force, two geometries — a γ = 10 electron in a bend versus a linac (exam level)

Setup. An electron with γ=10\gamma=10 (β=0.99499\beta=0.99499, kinetic energy 4.6 MeV) crosses a B=1B=1 T field perpendicular to its velocity. Compute the radiated power by both routes — the field formula (2) and the γ4/ρ2\gamma^4/\rho^2 law — then compare with a linac exerting the same force.

Step 1. Field route. The bracket in (2) is (vB)2(vB)^2 with v=0.99499c=2.983×108v=0.99499c=2.983\times10^{8} m/s:

P=e4γ2(vB)26πε0me2c3=(1.602×10−19)4(100)(2.983×108)26π(8.854×10−12)(9.109×10−31)2(2.998×108)3=1.57×10−12 W.P=\frac{e^4\gamma^2(vB)^2}{6\pi\varepsilon_0m_e^2c^3}=\frac{(1.602\times10^{-19})^4(100)(2.983\times10^{8})^2}{6\pi(8.854\times10^{-12})(9.109\times10^{-31})^2(2.998\times10^{8})^3}=\mathbf{1.57\times10^{-12}\ W} .

Step 2. Radius route. ρ=γmev/eB=1.70\rho=\gamma m_ev/eB=1.70 cm, and P=μ0e2c3γ4β4/6πρ2=1.57×10−12P=\mu_0e^2c^3\gamma^4\beta^4/6\pi\rho^2=1.57\times10^{-12} W — identical to twelve digits (this edition's numerical audit). Two powers of γ\gamma moved from the prefactor into ρ2\rho^2; no physics moved anywhere.

Step 3. Per-turn sanity. Revolution time T=2πρ/v=0.36T=2\pi\rho/v=0.36 ns, so the electron loses PT/e≈3.5PT/e\approx3.5 meV per turn — utterly negligible against 4.6 MeV, as expected at γ=10\gamma=10: the fourth power needs thousands of γ\gamma's before it bites.

Step 4. The linac comparison. The magnetic force here is evB=1.60×10−19×2.98×108×1=4.8×10−11evB=1.60\times10^{-19}\times2.98\times10^{8}\times1=4.8\times10^{-11} N — equivalent to an electric gradient Eeq=vB=298E_{\text{eq}}=vB=298 MV/m, three times today's best linac gradients. Applying (2) with that field parallel to v\mathbf{v}: P∥=e4Eeq2/6πε0me2c3=1.57×10−14P_\parallel=e^4E_{\text{eq}}^2/6\pi\varepsilon_0m_e^2c^3=1.57\times10^{-14} W — exactly γ2=100\gamma^2=100 times less than the bend, at identical force. The bend pays the full transverse bracket; the push collects the 1/γ21/\gamma^2 longitudinal discount.

The lesson: equation (2) makes "bend versus push" a substitution exercise — same force, γ2\gamma^2 different bills, and the whole comparison ran without solving a single trajectory.

Going further

References. Jackson, Classical Electrodynamics, 3rd ed., §14.1–14.3 · Griffiths, Introduction to Electrodynamics, 4th ed., §11.2.2 · Zangwill, Modern Electrodynamics, §23.4 · Landau & Lifshitz, Classical Theory of Fields, §73–74 · Schwinger, Phys. Rev. 75, 1912 (1949)

Provenance. Equations are standard textbook results. A verifier should: (1) reproduce the covariant formula P=−μ0q2c aμaμ/6πP=-\mu_0 q^2 c\,a_\mu a^\mu/6\pi from Jackson §14.2 (Eq. 14.26) and Griffiths §11.2.2 (Eq. 11.73); (2) confirm the Liénard form (Jackson Eq. 14.26, Griffiths Eq. 11.74); (3) check the γ6\gamma^6 collinear and γ4\gamma^4 transverse limits (Jackson Eqs. 14.27–14.31); (4) verify the four-acceleration components and γ˙=γ3β⋅β˙\dot\gamma=\gamma^3\boldsymbol\beta\cdot\dot{\boldsymbol\beta} by direct differentiation. No numerical constants beyond μ0\mu_0, cc are quoted. Cross-references: Zangwill §23.4, Purcell Ch. 11. Addendum for this edition (division of labor and corrections): the strong-UG twin Relativistic Larmor Formula (Liénard's Result) owns both derivation routes to the power formula (the κ5\kappa^5 source-clock integration and the rest-frame-isotropy invariance theorem), the γ6\gamma^6-versus-γ4\gamma^4 held-fixed doctrine, and the LEP/LHC/linac anchor numbers — all recalled by name, none re-derived. This edition owns the symbolic aμaμa_\mu a^\mu derivation with the Lagrange-identity bridge (Derivation A, upgrading the twin's numerical check to a theorem, computer-algebra verified: both spellings identically equal), the external-fields form P=q4γ26πε0m2c3[(E+v×B)2−(β⋅E)2]P=\frac{q^4\gamma^2}{6\pi\varepsilon_0m^2c^3}[(\mathbf{E}+\mathbf{v}\times\mathbf{B})^2-(\boldsymbol\beta\cdot\mathbf{E})^2] with its pure-B/pure-E audits and machine-precision numerical cross-check against the kinematic form (Derivation B), and the derived engineering formula δE[MeV]=0.0885 E4[GeV]/ρ[m]\delta E[\mathrm{MeV}]=0.0885\,E^4[\mathrm{GeV}]/\rho[\mathrm{m}] with FCC-ee (35.5 MeV/turn at Z, 9.1 GeV at ttˉt\bar t, ρ=10.76\rho=10.76 km) and 6 GeV light-source (4.90 MeV/turn, ρ=23.4\rho=23.4 m) anchors plus damping estimates E T0/δEE\,T_0/\delta E (0.42 s; 3.5 ms) (Derivation C). Three base defects corrected with this addendum: the base's covariant eqCard (4) and key-idea wrote P=−μ0q2c aμaμ/6πP=-\mu_0q^2c\,a_\mu a^\mu/6\pi, dimensionally short by c2c^2 — the correct coefficient is −μ0q2/6πc-\mu_0q^2/6\pi c (consistent with the rest-frame reduction and Jackson 14.26); the base's Deriv-B remark that the dot-to-cross conversion "holds only when β2=1\beta^2=1" is false (the Lagrange identity is exact at all β\beta — Derivation A Step 4); and the base's space-component expression γ2β˙+γ4(β⋅β˙)βc\gamma^2\dot{\boldsymbol\beta}+\gamma^4(\boldsymbol\beta\cdot\dot{\boldsymbol\beta})\boldsymbol\beta c dropped the factor cc on its first term (both terms carry cc — Derivation A Step 2). Synchrotron spectrum and critical frequency deferred by name to Synchrotron Radiation; angular lobes to Angular Distribution of Radiation from Accelerated Charges; self-consistent back-reaction to the radiation-reaction set.

Builds on: Larmor Formula from the Liénard–Wiechert Fields · Liénard–Wiechert Potentials: Covariant Formulation · Fields of a Relativistic Moving Charge Leads to: Angular Distribution of Radiation from Accelerated Charges · Synchrotron Radiation · Bremsstrahlung

How does this article divide the subject with the strong-UG Relativistic Larmor Formula article?

Cleanly, by route and by role. The strong-UG twin derives the power formula twice — the κ5\kappa^5 retarded-flux integration and the rest-frame invariance theorem — and owns the conceptual doctrine (emitted versus received clocks, γ6\gamma^6 versus γ4\gamma^4 with the held-fixed disclaimers, the LEP/LHC/linac verdicts). This edition assumes the formula and makes it computational: the symbolic four-acceleration algebra (their numeric check, promoted to a theorem), the external-fields substitution that skips trajectories entirely, and the engineering per-turn formula with current-generation machine numbers. Read theirs for why the formula is true; read this one when you need to evaluate it.

Is the radiated power really a pointwise function of the fields — doesn't radiation depend on the history of the motion?

The power is pointwise; the fields at the detector are not. The covariant formula needs only aμa^\mu now, and the Lorentz force ties aμa^\mu to the external fields and velocity now — so the energy-loss rate is local, which is precisely what makes equation (2) code-friendly. What depends on history is where and when that energy shows up far away: the radiation arrives on the light cone of each emission event, carrying retarded-time structure, beaming, and spectral content (the angular-distribution and synchrotron articles). The division of labor is exact: local formula for the budget, retarded machinery for the delivery.

The field form subtracts the beta-dot-E term. What is that discount, physically?

It is the price difference between changing a particle's direction and changing its energy. The longitudinal part of the electric force does work — it feeds the particle's energy account — and work-done shows up in the time component of dpμ/dτdp^\mu/d\tau, which enters the invariant with the opposite sign to the space part. Force spent purely on speeding up therefore radiates less, by exactly 1/γ21/\gamma^2 at fixed force magnitude, than force spent on bending. In the extreme case the two books nearly cancel: a linac's gradient radiates as if the particle were nonrelativistic, forever. This single minus sign is the entire strategic difference between linear and circular colliders, and it was sitting in the metric signature all along.

Why do plasma-wakefield and strong-laser codes use equation (2) instead of integrating trajectories and applying Liénard?

Cost and stability. A particle-in-cell code already knows (E,B)(\mathbf{E},\mathbf{B}) on the grid and β\boldsymbol\beta for each macro-particle; evaluating (2) is a handful of multiplications, while extracting β˙\dot{\boldsymbol\beta} numerically means differentiating noisy trajectories — an operation that amplifies grid noise catastrophically at γ≫1\gamma\gg1, where the bracket involves delicate cancellations. The field form has the cancellation done analytically. The caveat: (2) prices the emission but does not by itself remove the radiated momentum from the particle; once the per-cycle loss stops being small (petawatt lasers, extreme betatron motion), codes must add the back-reaction — the Landau–Lifshitz machinery whose classical foundations are the radiation-reaction articles of the strong-UG chapter and their covariant E11 successors.

Where does the 0.0885 in the accelerator formula actually come from?

From three constants and a unit conversion, nothing else: δE=q2β3γ4/3ε0ρ\delta E=q^2\beta^3\gamma^4/3\varepsilon_0\rho with γ=E/mec2\gamma=E/m_ec^2 gives δE=e23ε0(mec2)4 E4/ρ\delta E=\frac{e^2}{3\varepsilon_0(m_ec^2)^4}\,E^4/\rho; evaluating e2/3ε0(mec2)4e^2/3\varepsilon_0(m_ec^2)^4 with CODATA values and expressing the result in MeV, GeV⁴ and meters yields 0.088460.08846. The fourth power of the electron rest energy in the denominator is worth staring at: it is why the same formula for protons carries (me/mp)4≈9×10−14(m_e/m_p)^4\approx9\times10^{-14}, and why muon-collider studies — mμ=207mem_\mu=207m_e, coefficient smaller by 1.8×1091.8\times10^{9} — keep returning to circular geometry for leptons after all.

The damping estimate divides beam energy by radiated power. Why does radiation damp oscillations instead of just shrinking the whole beam energy?

Because the loss and the replacement act along different directions — that asymmetry is the cooling mechanism. Synchrotron emission kicks opposite the particle's instantaneous velocity, shaving all three momentum components proportionally; the RF cavities restore energy purely longitudinally. Net effect per turn: transverse momenta shrink by ∼δE/E\sim\delta E/E while the longitudinal component is topped back up — an exponential damping of betatron oscillations with time constant ∼ET0/δE\sim E T_0/\delta E (0.42 s at FCC-ee's Z pole, milliseconds in light sources). The exact distribution of damping among horizontal, vertical, and energy oscillations — including the lattice-dependent partition numbers and the quantum-excitation floor set by photon shot noise — is damping-partition theory, deferred to accelerator physics proper.

Does the invariance of P survive when the fields themselves are transformed to another frame?

Yes, and checking it is a nice exercise in the chapter's machinery. Boost to any frame: E\mathbf{E} and B\mathbf{B} mix per Field Transformations, γ\gamma and β\boldsymbol\beta change, and every factor in equation (2) is frame-dependent — yet the combination is not, because (2) is just the scalar −μ0q2(dpμ/dτ)(dpμ/dτ)/6πm2c-\mu_0q^2(dp_\mu/d\tau)(dp^\mu/d\tau)/6\pi m^2c written in one frame's variables. A drift-balanced charge is the instructive case: in the lab it sits in large crossed fields with zero net force (and, since E=−v×B\mathbf{E}=-\mathbf{v}\times\mathbf{B} is perpendicular to v\mathbf{v}, zero work term too); in its rest frame the transverse electric field itself vanishes by the transformation rules — two descriptions, one invariant P=0P=0. When a calculation of radiated power appears to change under a boost, a factor of γ\gamma has been misplaced between dtdt and dτd\tau — the two-clocks bookkeeping of the companion Larmor article, again.

Before this

Where this leads