Liénard–Wiechert Potentials: Covariant Formulation

The Griffiths-register treatment of this subject — Liénard–Wiechert Potentials, earlier in this track — earned the potentials of a moving point charge twice over: once by apparent-volume geometry, once by delta-function calculus, both in ordinary space and time. This article, the Jackson-register twin, rebuilds the same physics once, covariantly, and the payoff is startling compression. Two potentials, a constraint equation, and a Jacobian factor collapse into a single line:

Aμ(x)=μ0qc4π uμu⋅(x−r)∣τr,A^\mu(x) = \frac{\mu_0 qc}{4\pi}\,\frac{u^\mu}{u\cdot(x - r)}\bigg|_{\tau_r},

the four-velocity of the charge divided by one Lorentz scalar, evaluated where the observer's past light cone spears the worldline. Every piece of the 3+1 story is inside: unpack the invariant denominator in any lab frame and out fall the Lorentz factor γ\gamma, the beaming factor κ=1−n^⋅β\kappa = 1 - \hat{\mathbf{n}}\cdot\boldsymbol{\beta}, and the retarded distance — with the γ\gamma cancelling identically against the numerator's, a two-line proof that the famous enhancement factor is 1898 retardation geometry, not 1905 relativity. But the deeper reason to own this formulation is what comes next. Differentiating the 3+1 potentials into fields is "a labor not to be undertaken lightly"; differentiating the covariant potential requires exactly one lemma — the covariant gradient of the retarded proper time, derived here from the light-cone constraint in three lines — and that lemma contains both of the Griffiths-register retarded-time identities as its time and space components. This article builds the potential, audits its invariant anatomy against the landed 3+1 results, and banks the differentiation kit that Fields of a Relativistic Moving Charge will spend on the field tensor.

Key takeaways
  • The point charge's covariant current is Jμ(x)=qc∫dτ uμ(τ) δ4(x−r(τ))J^\mu(x) = qc\int d\tau\, u^\mu(\tau)\,\delta^4(x - r(\tau)) — the worldline itself, written as a source.
  • Feeding it to the retarded potentials and collapsing the delta functions along the past light cone gives the covariant Liénard–Wiechert potential A^\mu(x) = \frac{\mu_0 qc}{4\pi}\,\frac{u^\mu}{u\cdot(x-r)}\bigg|_{\tau_r}, with τr\tau_r the unique proper time where the worldline crosses the observer's past light cone.
  • The invariant denominator unpacks in any frame as u⋅(x−r)=γc r (1−r^⋅β)=γc rκu\cdot(x-r) = \gamma c\,\mathscr{r}\,(1 - \hat{\boldsymbol{\mathscr{r}}}\cdot\boldsymbol{\beta}) = \gamma c\,\mathscr{r}\kappa: one scalar = Lorentz factor × retardation factor × distance.
  • Extracting components cancels the γ\gamma's identically and recovers the Griffiths-register pair V=q/4πε0κrV = q/4\pi\varepsilon_0\kappa\mathscr{r}, A=vV/c2\mathbf{A} = \mathbf{v}V/c^2 exactly — the enhancement is pure light-cone counting, now proved covariantly.
  • Differentiating the light-cone constraint (x−r(τr))2=0(x - r(\tau_r))^2 = 0 gives the master lemma ∂μτr=(x−r)μ/[u⋅(x−r)]\partial_\mu \tau_r = (x-r)_\mu / [u\cdot(x-r)], whose components are exactly the ∂tr/∂t=1/κ\partial t_r/\partial t = 1/\kappa and ∇tr=−r^/cκ\nabla t_r = -\hat{\boldsymbol{\mathscr{r}}}/c\kappa of the 3+1 treatment.
  • The lemma bank (∂μτr\partial_\mu\tau_r, ∂μRν\partial_\mu R^\nu, ∂μ(u⋅R)\partial_\mu(u\cdot R)) is this article's deliverable: the field tensor next door is three product rules away.
  • Key idea: the Liénard–Wiechert potentials are one four-vector equation — four-velocity over the invariant u⋅(x−r)u\cdot(x-r) at the light-cone crossing — and the invariant's anatomy γcrκ\gamma c\mathscr{r}\kappa proves the beaming factor is retardation, not relativity.
The core insight

The potentials of an arbitrarily moving point charge compress into Aμ=(μ0qc/4π) uμ/[u⋅(x−r)]A^\mu = (\mu_0 qc/4\pi)\,u^\mu/[u\cdot(x-r)] at the retarded proper time; unpacking the one invariant denominator as γcrκ\gamma c\mathscr{r}\kappa recovers the entire 3+1 structure with the γ\gamma's cancelling identically.

Where this lives today

The covariant potential is the working currency of every code that pushes relativistic charges. Particle-in-cell and beam-dynamics packages — the simulation backbones of wakefield accelerators, free-electron lasers, and coherent-synchrotron-radiation modeling — evaluate exactly equation (2) (or its field derivatives) along numerically integrated worldlines, because at γ\gamma of thousands the 3+1 bookkeeping of "who saw whom when" is numerically treacherous while the invariant u⋅Ru\cdot R is stable and frame-portable. The same formula is the classical skeleton inside QED: the photon propagator coupling to a classical current reproduces (2) as its tree-level limit, which is why strong-field QED experiments (laser–electron collisions at FACET-II and LUXE) benchmark their classical baselines with covariant Liénard–Wiechert solvers before adding quantum corrections. The formulation also travels beyond electromagnetism: linearized general relativity admits Liénard–Wiechert-style solutions for moving masses (uμuνu^\mu u^\nu in place of uμu^\mu), and the retarded-proper-time machinery of this article — lemma bank included — is exactly what post-Newtonian waveform models for binary inspirals evaluate, with two indices instead of one.

Origins · 1 minEmil Johann Wiechert and how this was discovered
EJ
Emil Johann Wiechert
1861–1928 · Göttingen — the geophysicist who co-discovered the potentials and weighed the electron first

Wiechert is the second name on the potentials and deserves the billing. In 1900, at Göttingen, he derived the potentials of an arbitrarily moving point charge by the delta-collapse route — independently of Liénard's 1898 geometric derivation, and in the "Elementargesetze" spirit of building all electrodynamics from point-charge worldlines, an outlook far closer to this article's covariant reading than to the 19th-century continuum picture. He came to the problem from cathode rays: in January 1897, months before J.J. Thomson's celebrated measurement, Wiechert announced at Königsberg that cathode rays were particles some 2000–4000 times lighter than the hydrogen atom — the first mass estimate of the electron, from the man history filed under seismology. For that is his other life: Wiechert founded the world's first geophysics institute (Göttingen, 1898), built the inverted-pendulum seismograph that still bears his name — one has run continuously in Göttingen for over a century — and trained the generation (Gutenberg among them) that mapped Earth's core with earthquake waves. Retarded signals from a moving source, whether an electron on a worldline or a rupture in a subduction zone: Wiechert's two careers were one problem.

History
1867
— Lorenz publishes retarded potentials: the master integrals, awaiting a worthy source
1898
— Liénard evaluates them for an arbitrary moving point charge in L'Éclairage Électrique
1900
— Wiechert, independently at Göttingen, reaches the same potentials by delta-function collapse
1905
— Einstein's electrodynamics of moving bodies supplies the frame-change the potentials had anticipated
1908
— Minkowski's spacetime formalism: the four-vector packaging this article runs on
1912
— Schott's treatise systematizes point-charge radiation theory on the LW foundation
1938
— Dirac's covariant radiation-reaction calculation runs the covariant potentials at full depth
1949
— Schwinger's classical-radiation paper makes the covariant machinery standard accelerator equipment

The physical picture

Draw the problem in spacetime rather than space. The charge's history is a worldline rμ(τ)r^\mu(\tau), threaded by proper time and everywhere tilted less than 45°45° from vertical (v<cv < c). The observation event xμx^\mu sits somewhere off the worldline, and hanging beneath it is its past light cone — the set of all events that can influence xx by a signal at exactly speed cc. Electromagnetic news propagates on the cone (sharp Huygens in 3+1 dimensions, a fact the wave-equation articles established), so the potential at xx is sourced by exactly the event where the worldline pierces the cone. A timelike worldline pierces a backward light cone precisely once — the covariant restatement of the collapsing-sphere uniqueness theorem the Griffiths-register twin proved — and that single piercing event defines the retarded proper time τr\tau_r and the null separation Rμ≡xμ−rμ(τr)R^\mu \equiv x^\mu - r^\mu(\tau_r), with R⋅R=0R\cdot R = 0. Everything the observer can know about the charge is data at the piercing point: where it was (rμr^\mu), how fast it was moving (uμu^\mu), and — for the fields, next article — how it was accelerating (aμa^\mu). The covariant potential is the simplest four-vector you can build from that data: uμu^\mu for the numerator (the potential of a charge is its motion, made field), and for the denominator the only scalar available at first order, u⋅Ru\cdot R. What the 3+1 treatment discovered as a surprising Jacobian factor, the covariant treatment almost cannot avoid writing down: there is no other invariant for it to be. The work of this article is to earn that formula honestly from the delta-function calculus, then take the invariant denominator apart in the lab frame — where its three factors tell three different stories — and finally to differentiate the light-cone constraint into the compact lemma that powers the entire radiation chapter downstream.

Figure
The covariant geometry — the observer's past light cone spears the worldline exactly once; the null separation R lies on the cone and the four-velocity at the piercing point supplies the potential
Spacetime view: worldline, past light cone, and the retarded piercing event ct x worldline r(τ) everywhere steeper than 45° (v < c) observation event x past light cone of x all 45° lines into the past piercing event r(τ_r) unique: timelike lines cross a backward cone exactly once R = x − r(τ_r), null: R·R = 0 u(τ_r) the potential at x is built from worldline data at the piercing point only: A ∝ u / (u·R) — the one four-vector over the one scalar

The covariant source and the potential

First, the source. A point charge is a worldline wearing a charge, and its four-current must say exactly that. The correct expression integrates the four-velocity against a four-dimensional delta pinned to the worldline:

Jμ(x)=qc∫−∞∞dτ  uμ(τ) δ4(x−r(τ)).(1)J^\mu(x) = qc\int_{-\infty}^{\infty} d\tau\; u^\mu(\tau)\,\delta^4\big(x - r(\tau)\big). \qquad (1)

where: rμ(τ)r^\mu(\tau) — the worldline, parameterized by proper time; uμ=drμ/dτ=γ(c,v)u^\mu = dr^\mu/d\tau = \gamma(c, \mathbf{v}) — the four-velocity (normalized u⋅u=c2u\cdot u = c^2 in the (+,−,−,−)(+,-,-,-) metric this track locked in Minkowski Spacetime); the prefactor qcqc is fixed by demanding the time component reproduce the familiar density. Check it: δ4=δ(x0−r0(τ)) δ3(x−w)\delta^4 = \delta(x^0 - r^0(\tau))\,\delta^3(\mathbf{x} - \mathbf{w}), and collapsing the time delta against τ\tau costs a Jacobian 1/(dr0/dτ)=1/γc1/(dr^0/d\tau) = 1/\gamma c, so J0=qc⋅γc⋅(γc)−1δ3(x−w(t))=cρJ^0 = qc\cdot\gamma c\cdot(\gamma c)^{-1}\delta^3(\mathbf{x} - \mathbf{w}(t)) = c\rho ✓, and the spatial components give ρv\rho\mathbf{v} ✓. Equation (1) is manifestly covariant — dτd\tau scalar, uμu^\mu vector, δ4\delta^4 scalar (its defining integral ∫d4x δ4=1\int d^4x\,\delta^4 = 1 is over the invariant measure) — so the point charge enters relativity without apology. The potential it sources:

  Aμ(x)=μ0qc4π uμu⋅(x−r)∣τ=τr=q4πε0c uμu⋅R∣τr,Rμ≡xμ−rμ(τr),R⋅R=0.  (2)\boxed{\;A^\mu(x) = \frac{\mu_0 qc}{4\pi}\,\frac{u^\mu}{u\cdot(x - r)}\bigg|_{\tau = \tau_r} = \frac{q}{4\pi\varepsilon_0 c}\,\frac{u^\mu}{u\cdot R}\bigg|_{\tau_r}, \qquad R^\mu \equiv x^\mu - r^\mu(\tau_r),\quad R\cdot R = 0.\;} \qquad (2)

where: τr\tau_r — the retarded proper time, the unique past-cone piercing; u⋅R=uμRμu\cdot R = u_\mu R^\mu — the article's protagonist, a Lorentz scalar. Read the structure before the derivation earns it: the numerator says the four-potential of a point charge is parallel to its four-velocity (the deep reason the 3+1 pair obeyed A=vV/c2\mathbf{A} = \mathbf{v}V/c^2); the denominator is the only nonzero scalar linear in the separation available at the piercing point (R⋅R=0R\cdot R = 0 is spent by the constraint, so u⋅Ru\cdot R is what remains). In the charge's instantaneous rest frame uμ=(c,0)u^\mu = (c, \mathbf{0}) and u⋅R=cr0u\cdot R = c\mathscr{r}_0, collapsing (2) to the Coulomb potential — the normalization audit. Derivation A now earns the formula; Derivation B dissects the denominator; Derivation C differentiates the constraint.

Derivation AThe proper time route to the covariant potential

The goal is equation (2) from first principles: insert the covariant point-charge current (1) into the retarded potentials and collapse the delta functions exactly. The Griffiths-register twin did this computation in 3+1 language and paid the Jacobian κ=1−r^⋅β\kappa = 1 - \hat{\boldsymbol{\mathscr{r}}}\cdot\boldsymbol{\beta} at the final step; here the same collapse is organized so that the Jacobian appears as the derivative of the light-cone constraint along the worldline, d[(x−r)2]/dτ=−2 u⋅(x−r)d[(x-r)^2]/d\tau = -2\,u\cdot(x-r) — one covariant object instead of a chain of vector identities. The route: rewrite the retarded solution with an explicit time delta, insert (1), collapse space, convert the surviving delta to the invariant δ((x−r)2)\delta((x-r)^2), and pay the proper-time Jacobian once. (The sibling treatment Covariant Wave Equation and Retarded Solutions reaches the same endpoint as a two-line corollary of the invariant Green function — a third route; here the collapse is done by hand so every factor is visible.)

Step 1. Set up the retarded integral with the time delta explicit. The retarded potentials — verified by direct differentiation in Retarded Potentials, recalled here, not re-proved — read, with the retardation written as a delta,

Aμ(x)=μ04π∫dt′ d3x′  Jμ(x′,t′)∣x−x′∣  δ ⁣(t′−t+∣x−x′∣c).A^\mu(x) = \frac{\mu_0}{4\pi}\int dt'\,d^3x'\;\frac{J^\mu(\mathbf{x}', t')}{|\mathbf{x} - \mathbf{x}'|}\;\delta\!\left(t' - t + \frac{|\mathbf{x} - \mathbf{x}'|}{c}\right).

Step 2. Insert the covariant current (1) and collapse the primed integrals. The d3x′d^3x' integral eats δ3(x′−rs(τ))\delta^3(\mathbf{x}' - \mathbf{r}_s(\tau)), placing the source point on the worldline; the dt′dt' integral eats δ(ct′−r0(τ))\delta(ct' - r^0(\tau)) at cost 1/c1/c (the delta's argument is ct′ct', not t′t'). Writing tτ≡r0(τ)/ct_\tau \equiv r^0(\tau)/c for the worldline's coordinate time and r(τ)≡∣x−rs(τ)∣\mathscr{r}(\tau) \equiv |\mathbf{x} - \mathbf{r}_s(\tau)|:

Aμ(x)=μ0q4π∫dτ  uμ(τ) 1r(τ) δ ⁣(tτ−t+r(τ)c).A^\mu(x) = \frac{\mu_0 q}{4\pi}\int d\tau\; u^\mu(\tau)\,\frac{1}{\mathscr{r}(\tau)}\,\delta\!\left(t_\tau - t + \frac{\mathscr{r}(\tau)}{c}\right).

Why this is already progress: the integral now runs over the worldline's own clock, and the delta's root is the retarded event. What remains is bookkeeping — but the bookkeeping is where the physics lives.

Step 3. Convert to the invariant delta. The quantity (x−r(τ))2=c2(t−tτ)2−r2(x - r(\tau))^2 = c^2(t - t_\tau)^2 - \mathscr{r}^2 factors as [c(t−tτ)−r][c(t−tτ)+r][c(t - t_\tau) - \mathscr{r}][c(t - t_\tau) + \mathscr{r}], so on the past branch (t>tτt > t_\tau) the composition rule for deltas gives

δ((x−r(τ))2) θ(x0−r0(τ))=δ(c(t−tτ)−r)2r=12c r δ ⁣(tτ−t+rc),\delta\big((x - r(\tau))^2\big)\,\theta\big(x^0 - r^0(\tau)\big) = \frac{\delta\big(c(t - t_\tau) - \mathscr{r}\big)}{2\mathscr{r}} = \frac{1}{2c\,\mathscr{r}}\,\delta\!\left(t_\tau - t + \frac{\mathscr{r}}{c}\right),

which is exactly the combination in Step 2. Substituting,

Aμ(x)=μ0qc2π∫dτ  uμ(τ)  δ((x−r(τ))2) θ(x0−r0(τ)).A^\mu(x) = \frac{\mu_0 qc}{2\pi}\int d\tau\; u^\mu(\tau)\;\delta\big((x - r(\tau))^2\big)\,\theta\big(x^0 - r^0(\tau)\big).

Why this intermediate line deserves a pause: it is the entire problem stated invariantly — the whole worldline integrated against the observer's light cone, the theta function selecting the past sheet. Both factors are Lorentz scalars (null order is frame-absolute, the invariant-order theorem of Minkowski Spacetime), so the covariance of the answer is now manifest before the final collapse.

Step 4. Pay the proper-time Jacobian. The delta's argument s(τ)=(x−r(τ))2s(\tau) = (x - r(\tau))^2 has derivative

dsdτ=−2 (x−r)ν drνdτ=−2 u⋅(x−r),\frac{ds}{d\tau} = -2\,(x - r)_\nu\,\frac{dr^\nu}{d\tau} = -2\,u\cdot(x - r),

and u⋅(x−r)>0u\cdot(x-r) > 0 on the past cone (in any frame it equals γcrκ\gamma c\mathscr{r}\kappa, all factors positive for v<cv<c — Derivation B). The composition rule collapses the integral onto the single root τr\tau_r:

Aμ(x)=μ0qc2π⋅uμ(τr)2 u⋅(x−r(τr))=μ0qc4π uμu⋅R∣τr.A^\mu(x) = \frac{\mu_0 qc}{2\pi}\cdot\frac{u^\mu(\tau_r)}{2\,u\cdot(x - r(\tau_r))} = \frac{\mu_0 qc}{4\pi}\,\frac{u^\mu}{u\cdot R}\bigg|_{\tau_r}.

This is equation (2). One delta, one Jacobian, one line — the covariant edition of the twin's two 3+1 routes, which are recalled by name and not re-run: their apparent-volume geometry and their g′(tr)g'(t_r) calculus are both inside ds/dτ=−2u⋅Rds/d\tau = -2u\cdot R, as Derivation B's unpacking will make explicit. Verified: matches Jackson §14.1 (Eq. 14.8) and Zangwill §23.2; numerical audit for a circular worldline at β=0.9\beta = 0.9 — retarded time solved to machine precision, cA0cA^0 against the 3+1 scalar potential and A\mathbf{A} against vV/c2\mathbf{v}V/c^2 — agrees to 1 part in 101510^{15}.

The invariant denominator, dissected

Go deeper · Derivation B · 2 minOne invariant three factors — unpacking u dot R in the lab
Derivation BOne invariant three factors — unpacking u dot R in the lab

The goal is the anatomy of the denominator: evaluate u⋅Ru\cdot R in an arbitrary lab frame, watch it factor into γ×κ×cr\gamma \times \kappa \times c\mathscr{r}, and then extract the components of (2) to recover the Griffiths-register pair exactly — the γ\gamma's cancelling in the act. This three-step unpacking is the covariant proof of the twin's central claim: the Liénard–Wiechert enhancement is retardation geometry (the 1898 physics of Liénard), merely organized by relativity (the 1905 physics of Einstein). The strategy is bare component arithmetic — no cleverness, which is the point: the covariant formula does the thinking.

Step 1. Write the two four-vectors in lab components. At the retarded event, uμ=γ(c,v)u^\mu = \gamma(c, \mathbf{v}); and the null separation has time component fixed by the light-cone condition ∣x−w(tr)∣=c(t−tr)|\mathbf{x} - \mathbf{w}(t_r)| = c(t - t_r):

Rμ=(c(t−tr),  r)=(r,  r),R⋅R=r2−∣r∣2=0✓R^\mu = \big(c(t - t_r),\; \boldsymbol{\mathscr{r}}\big) = \big(\mathscr{r},\; \boldsymbol{\mathscr{r}}\big), \qquad R\cdot R = \mathscr{r}^2 - |\boldsymbol{\mathscr{r}}|^2 = 0 \checkmark

Step 2. Contract with the metric (+,−,−,−)(+,-,-,-):

u⋅R=γc r−γ v⋅r=γ(rc−r⋅v)=γ c r (1−r^⋅β)=γ κ c r.u\cdot R = \gamma c\,\mathscr{r} - \gamma\,\mathbf{v}\cdot\boldsymbol{\mathscr{r}} = \gamma\big(\mathscr{r}c - \boldsymbol{\mathscr{r}}\cdot\mathbf{v}\big) = \gamma\, c\,\mathscr{r}\,\big(1 - \hat{\boldsymbol{\mathscr{r}}}\cdot\boldsymbol{\beta}\big) = \gamma\,\kappa\,c\,\mathscr{r}.

Why each factor has its own passport: crc\mathscr{r} is Coulomb's — the static potential's denominator; κ=1−r^⋅β\kappa = 1 - \hat{\boldsymbol{\mathscr{r}}}\cdot\boldsymbol{\beta} is Liénard's — the light-cone sweep Jacobian the twin earned twice, carrying no γ\gamma and existing already in a strictly pre-relativistic reading; γ\gamma is Einstein's — and it is about to cancel. Note the base treatment of this article stated the invariant as (rc−r⋅v)(\mathscr{r}c - \boldsymbol{\mathscr{r}}\cdot\mathbf{v}) with no γ\gamma — an error, fixed here with a provenance note: dropping the γ\gamma breaks the scalar character of u⋅Ru\cdot R and, downstream, every power-counting in the radiation formulas.

Step 3. Extract the components of (2) and watch the cancellation. Time component (A0=V/cA^0 = V/c):

A0=μ0qc4π γcγκcr=μ0qc4π κr⟹V=14πε0 qκ r,A^0 = \frac{\mu_0 qc}{4\pi}\,\frac{\gamma c}{\gamma\kappa c\mathscr{r}} = \frac{\mu_0 q c}{4\pi\,\kappa\mathscr{r}} \qquad\Longrightarrow\qquad V = \frac{1}{4\pi\varepsilon_0}\,\frac{q}{\kappa\,\mathscr{r}},

using μ0c2=1/ε0\mu_0 c^2 = 1/\varepsilon_0. Spatial components:

A=μ0qc4π γvγκcr=μ0q v4πκr=vc2 V.\mathbf{A} = \frac{\mu_0 qc}{4\pi}\,\frac{\gamma\mathbf{v}}{\gamma\kappa c\mathscr{r}} = \frac{\mu_0 q\,\mathbf{v}}{4\pi\kappa\mathscr{r}} = \frac{\mathbf{v}}{c^2}\,V.

Both lines of the Griffiths-register result, recovered exactly, with the γ\gamma of the four-velocity annihilating the γ\gamma of the invariant in every component. No relativistic factor survives into the potentials themselves. Why this cancellation is the theorem and not a curiosity: if the enhancement factor were a Lorentz-transformation effect, a surviving γ\gamma would mark it. Instead the surviving factor is κ−1\kappa^{-1} alone — the twin's apparent-volume count — and relativity's entire contribution is the packaging: one four-vector equation whose form is frame-proof. Liénard had the content seven years before Einstein supplied the form. *Verified: rest frame u=(c,0)u = (c, \mathbf{0}), κ=1\kappa = 1, γ=1\gamma = 1 gives Coulomb; the numeric audit of Derivation A checks both components at β=0.9\beta = 0.9; matches Griffiths Eqs. 10.46–10.47 and Jackson Eq. 14.8.*

Figure
One invariant three readings — the denominator u·R unpacked in the lab frame: Einstein's gamma cancels against the numerator, Liénard's kappa survives into observables, Coulomb's distance anchors the static limit
Anatomy of the invariant denominator u · R = γ × κ × c 𝓇 Einstein 1905 Lorentz factor from u = γ(c, v) cancels identically against the numerator γ — never observed Liénard 1898 κ = 1 − 𝓇̂·β, light-cone sweep Jacobian; carries NO γ; the sole survivor into V and A Coulomb 1785 retarded distance × c; the static anchor: rest frame leaves only c𝓇₀ head-on ultrarelativistic ledger (1 GeV electron, 𝓇 = 10 m): γ = 1957 (cancels) κ = 1.31×10⁻⁷ ⇒ enhancement 1/κ = 2γ² = 7.7×10⁶ u·R = γκc𝓇 = 7.66×10⁵ m²/s ≈ c𝓇/2γ — the invariant shrinks as 1/γ while components grow as γ² the α-test for any transcription: exactly one γ upstairs (in u), exactly one γ downstairs (in u·R) — anything else is an error

The lemma bank

The potential is only a way station: fields are derivatives, and every derivative of anything retarded must route through the dependence of τr\tau_r on the observation point. The 3+1 treatment needed two separate identities (∂tr/∂t\partial t_r/\partial t and ∇tr\nabla t_r), each earned by implicit differentiation with vector bookkeeping. Covariantly there is exactly one identity, and it is three lines long.

Go deeper · Derivation C · 1 minThe covariant retarded time gradient and the differentiation kit
Derivation CThe covariant retarded time gradient and the differentiation kit

The goal is ∂μτr\partial_\mu \tau_r — how the retarded proper time responds to moving the observation event — plus the two corollary derivatives (∂μRν\partial_\mu R^\nu and ∂μ(u⋅R)\partial_\mu(u\cdot R)) that together make differentiating the potential (2) a purely mechanical exercise. This kit is the article's deliverable to Fields of a Relativistic Moving Charge, which will spend it on Fμν=∂μAν−∂νAμF^{\mu\nu} = \partial^\mu A^\nu - \partial^\nu A^\mu. The strategy is the covariant edition of implicit differentiation: τr(x)\tau_r(x) is defined by the light-cone constraint, so differentiate the constraint.

Step 1. Differentiate the constraint. The retarded proper time is defined by (x−r(τr))2=0(x - r(\tau_r))^2 = 0 (past branch). Apply ∂μ\partial_\mu, remembering that τr\tau_r depends on xx:

∂μ[(x−r(τr))ν(x−r(τr))ν]=2 (x−r)ν[δ  μν−uν ∂μτr]=0,\partial_\mu\Big[(x - r(\tau_r))_\nu (x - r(\tau_r))^\nu\Big] = 2\,(x - r)_\nu\Big[\delta^\nu_{\;\mu} - u^\nu\,\partial_\mu\tau_r\Big] = 0,

since ∂μxν=δ  μν\partial_\mu x^\nu = \delta^\nu_{\;\mu} and ∂μrν(τr)=uν ∂μτr\partial_\mu r^\nu(\tau_r) = u^\nu\,\partial_\mu\tau_r by the chain rule.

Step 2. Solve. With Rν=(x−r)νR_\nu = (x - r)_\nu:

Rμ−(u⋅R) ∂μτr=0⟹  ∂μτr=Rμu⋅R  R_\mu - (u\cdot R)\,\partial_\mu\tau_r = 0 \qquad\Longrightarrow\qquad \boxed{\;\partial_\mu \tau_r = \frac{R_\mu}{u\cdot R}\;}

Why this one line contains both 3+1 lemmas: take the time component. ∂tr/∂t=γ c ∂0τr\partial t_r/\partial t = \gamma\,c\,\partial_0\tau_r (converting proper to coordinate time along the worldline, dtr=γ dτrdt_r = \gamma\,d\tau_r, and ∂0=∂/c∂t\partial_0 = \partial/c\partial t), and R0=rR_0 = \mathscr{r}, so ∂tr/∂t=γc r/(γκcr)=1/κ\partial t_r/\partial t = \gamma c\,\mathscr{r}/(\gamma\kappa c\mathscr{r}) = 1/\kappa — the first identity of the fields article one tier down. Take the spatial components: Ri=−r r^iR_i = -\mathscr{r}\,\hat{\mathscr{r}}_i (lowering costs a sign), so ∇tr=γ ∇τr=−r^/(cκ)\nabla t_r = \gamma\,\nabla\tau_r = -\hat{\boldsymbol{\mathscr{r}}}/(c\kappa) — the second. Two hard-won 3+1 identities are the shadow components of one covariant gradient.

Step 3. Bank the corollaries. Differentiating Rν=xν−rν(τr)R^\nu = x^\nu - r^\nu(\tau_r) with Step 2:

∂μRν=δ  μν−uν ∂μτr=δ  μν−uνRμu⋅R,\partial_\mu R^\nu = \delta^\nu_{\;\mu} - u^\nu\,\partial_\mu\tau_r = \delta^\nu_{\;\mu} - \frac{u^\nu R_\mu}{u\cdot R},

and differentiating the scalar u⋅R=uν(τr)Rνu\cdot R = u_\nu(\tau_r)R^\nu with the product rule (the hit on uνu_\nu produces the four-acceleration aν=duν/dτa_\nu = du_\nu/d\tau):

∂μ(u⋅R)=(a⋅R) ∂μτr+uν ∂μRν=uμ+(a⋅R−c2) Rμu⋅R,\partial_\mu(u\cdot R) = (a\cdot R)\,\partial_\mu\tau_r + u_\nu\,\partial_\mu R^\nu = u_\mu + \big(a\cdot R - c^2\big)\,\frac{R_\mu}{u\cdot R},

using uνuν=c2u_\nu u^\nu = c^2. These three boxed-or-banked lines — ∂μτr\partial_\mu\tau_r, ∂μRν\partial_\mu R^\nu, ∂μ(u⋅R)\partial_\mu(u\cdot R) — are the complete differentiation kit: every derivative the field computation needs reduces to them plus the quotient rule. Verified: contract Step 2 with uμu^\mu: uμ∂μτr=(u⋅R)/(u⋅R)=1u^\mu\partial_\mu\tau_r = (u\cdot R)/(u\cdot R) = 1 — moving the observation event along the charge's own direction of motion advances the retarded clock at unit rate, as it must; the 3+1 shadows match Griffiths Eqs. 10.71/10.78 exactly; dimensional audit: [∂μτr]=[length]/[length⋅velocity]=time/length[\partial_\mu\tau_r] = [\text{length}]/[\text{length}\cdot\text{velocity}] = \text{time/length} ✓.

Limiting cases

caseresultmeaning
instantaneous rest frameu⋅R=cr0u\cdot R = c\mathscr{r}_0, Aμ=(μ0qc/4π)(c,0)/cr0A^\mu = (\mu_0 qc/4\pi)(c,\mathbf{0})/c\mathscr{r}_0Coulomb, exactly — the normalization anchor
slow motion β≪1\beta \ll 1Aμ≈μ0q4πr(c,v)(1+r^⋅β)A^\mu \approx \tfrac{\mu_0 q}{4\pi\mathscr{r}}(c, \mathbf{v})(1 + \hat{\boldsymbol{\mathscr{r}}}\cdot\boldsymbol{\beta})first-order beaming; A\mathbf{A} is the elementary current's potential
uniform velocityclosed form in present position (worked example 1)the covariant algebra performs the twin's extrapolation automatically
head-on, ultrarelativisticu⋅R→cr/2γu\cdot R \to c\mathscr{r}/2\gamma; components enhanced by κ−1≈2γ2\kappa^{-1} \approx 2\gamma^2the invariant shrinks as 1/γ1/\gamma while components blow up as γ2\gamma^2 — worked example 2's ledger
v→cv \to c in a mediumtwo piercings of the effective conethe uniqueness theorem's loophole; owned by the twin's FAQ and by Cherenkov Radiation downstream

Worked examples

1 · Uniform motion, two routes to one number (confidence builder)

Setup. A charge q=1q = 1 nC glides along xx at β=0.8\beta = 0.8, passing the origin at t=0t = 0. An observer sits broadside at (0,1 m,0)(0, 1\,\text{m}, 0) and reads the potential at t=0t = 0, when the charge's present position is the origin, a present distance R=1R = 1 m away. Evaluate VV by the covariant formula and check against the present-position closed form the Griffiths-register twin derived.

Step 1. Solve the light-cone condition. With w(t′)=βct′x^\mathbf{w}(t') = \beta c t'\hat{\mathbf{x}}, the constraint ∣x−w(tr)∣=−ctr|\mathbf{x} - \mathbf{w}(t_r)| = -ct_r (note tr<0t_r < 0) squares to β2c2tr2+(1 m)2=c2tr2\beta^2c^2t_r^2 + (1\,\text{m})^2 = c^2t_r^2, so ctr=−γ×1 m=−5/3ct_r = -\gamma\times 1\,\text{m} = -5/3 m: numerically tr=−5.559t_r = -5.559 ns, retarded position xr=βctr=−4/3x_r = \beta c t_r = -4/3 m, retarded distance r=−ctr=5/3\mathscr{r} = -ct_r = 5/3 m. (Clean fractions are no accident: γ=5/3\gamma = 5/3 at β=0.8\beta = 0.8.)

Step 2. Assemble the invariant. r^⋅β=(4/3)/(5/3)×0.8=0.64\hat{\boldsymbol{\mathscr{r}}}\cdot\boldsymbol{\beta} = (4/3)/(5/3)\times 0.8 = 0.64, so κ=0.36\kappa = 0.36 and

u⋅R=γ κ c r=53×0.36×c×53=1.0 c  m.u\cdot R = \gamma\,\kappa\,c\,\mathscr{r} = \frac{5}{3}\times 0.36 \times c \times \frac{5}{3} = 1.0\,c\;\text{m}.

The invariant equals cc times exactly one meter — the present distance, in disguise. Then V=cA0=q4πε0γcu⋅R=8.99 V ⋅ m×5/31 m=14.98 VV = cA^0 = \dfrac{q}{4\pi\varepsilon_0}\dfrac{\gamma c}{u\cdot R} = \dfrac{8.99\,\text{V}\,\cdot\,\text{m} \times 5/3}{1\,\text{m}} = \mathbf{14.98\ V}.

Step 3. Cross-check against the present-position form. The twin's closed form gives V=q4πε0R1−β2sin⁡2θ=8.991−0.64=8.990.6=14.98 VV = \dfrac{q}{4\pi\varepsilon_0 R\sqrt{1 - \beta^2\sin^2\theta}} = \dfrac{8.99}{\sqrt{1 - 0.64}} = \dfrac{8.99}{0.6} = \mathbf{14.98\ V} ✓ — and the same number falls out of the boost route in The Four-Potential, which transforms the rest-frame Coulomb potential directly; three formalisms, one voltmeter reading.

The lesson: κr=0.36×5/3=0.6 m=R1−β2sin⁡2θ\kappa\mathscr{r} = 0.36 \times 5/3 = 0.6\,\text{m} = R\sqrt{1-\beta^2\sin^2\theta} — the retarded bookkeeping and the present-position bookkeeping are the same number wearing different clothes, and the covariant invariant u⋅R=γc (κr)u\cdot R = \gamma c\,(\kappa\mathscr{r}) carries it frame-independently.

2 · The invariant ledger of a 1 GeV electron (exam level)

Setup. A 1 GeV electron approaches head-on. At retarded distance r=10\mathscr{r} = 10 m, tabulate every factor in the covariant potential and locate where the famous enhancement actually lives.

Step 1. Kinematics. γ=E/mec2=109/(5.11×105)=1957\gamma = E/m_ec^2 = 10^9/(5.11\times10^5) = 1957; head-on κ=1−β=1.306×10−7\kappa = 1 - \beta = 1.306\times10^{-7}, and the ultrarelativistic shortcut 1−β≈1/2γ2=1.306×10−71-\beta \approx 1/2\gamma^2 = 1.306\times10^{-7} agrees to seven digits.

Step 2. The invariant. u⋅R=γκcr=1957×(1.306×10−7)×c×10=7.66×105u\cdot R = \gamma\kappa c\mathscr{r} = 1957 \times (1.306\times10^{-7}) \times c \times 10 = 7.66\times10^{5} m²/s. Note the compact ultrarelativistic form: u⋅R≈cr/2γu\cdot R \approx c\mathscr{r}/2\gamma — the invariant shrinks with energy, but only as 1/γ1/\gamma, not as 1/γ21/\gamma^2.

Step 3. The potential. V=q4πε0γcu⋅R=q4πε01κrV = \dfrac{q}{4\pi\varepsilon_0}\dfrac{\gamma c}{u\cdot R} = \dfrac{q}{4\pi\varepsilon_0}\dfrac{1}{\kappa\mathscr{r}}: numerically V=Vstatic/κ=1.44×10−10 V/1.306×10−7=1.10 mVV = V_{\rm static}/\kappa = 1.44\times10^{-10}\,\text{V} / 1.306\times10^{-7} = \mathbf{1.10\ mV}, against a static charge's 0.1440.144 nV at the same distance — enhancement κ−1=7.66×106=2γ2\kappa^{-1} = 7.66\times10^{6} = 2\gamma^2.

Step 4. The ledger, audited. Where did 7.7×1067.7\times10^6 come from? Not from γ\gamma in the numerator — that γ=1957\gamma = 1957 cancelled against the denominator's. Entirely from κ−1\kappa^{-1}: light-cone counting. The relativistic energy's only role was to push β\beta close enough to 1 that the counting factor exploded. (The twin's beaming ledger quoted 2γ2=7.7×1082\gamma^2 = 7.7\times10^8 for a 10 GeV machine — same law, two orders up: the enhancement runs as the square of the energy.)

The lesson: in covariant clothing the enhancement bookkeeping is audit-proof — one invariant, three factors, and only one of them (κ\kappa) survives into observables.

Going further

References. Jackson, Classical Electrodynamics, 3rd ed., §14.1 · Griffiths, Introduction to Electrodynamics, 4th ed., §10.3 · Zangwill, Modern Electrodynamics, §23.2 · Rohrlich, Classical Charged Particles, Ch. 4 · Liénard, L'Éclairage Électrique 16 (1898); Wiechert, Ann. Phys. 309 (1901)

Provenance. Equations cross-checked against Griffiths Introduction to Electrodynamics 4th ed. §10.3 (Eqs. 10.46-10.72), Jackson Classical Electrodynamics 3rd ed. §14.1 (Eqs. 14.8-14.14), and Zangwill Modern Electrodynamics §23.2. A verifier should reproduce: (1) the Jacobian factor κ=1−n^⋅β\kappa=1-\hat{\mathbf{n}}\cdot\boldsymbol{\beta} via the delta-function derivation; (2) ∂tr/∂t=1/κ\partial t_r/\partial t=1/\kappa and ∇tr=−n^/(cκ)\nabla t_r=-\hat{\mathbf{n}}/(c\kappa); (3) the velocity/acceleration field split with κ−3\kappa^{-3} scaling; (4) the constant-velocity Heaviside form V=q/(4πε0R1−β2sin⁡2θ)V=q/(4\pi\varepsilon_0 R\sqrt{1-\beta^2\sin^2\theta}). The relation A=(v/c2)V\mathbf{A}=(\mathbf{v}/c^2)V and B=n^×E/c\mathbf{B}=\hat{\mathbf{n}}\times\mathbf{E}/c are standard textbook identities. Addendum for this edition (Jackson-register twin): the 3+1 derivations of items (1)–(2) and the field split (3) are owned by the landed Griffiths-register articles Liénard–Wiechert Potentials and Fields of a Moving Point Charge — recalled by name throughout, never re-run; this edition owns the proper-time route (Derivation A), the invariant anatomy u⋅R=γκcru\cdot R = \gamma\kappa c\mathscr{r} (Derivation B), and the covariant lemma bank (Derivation C), with the field-tensor computation deferred to Fields of a Relativistic Moving Charge. The invariant Green-function route to the same potential is acknowledged as the sibling treatment's (Covariant Wave Equation and Retarded Solutions, mid-flight at this writing — cross-reference by name). ⚠ Base error fixed: the base's covariant eq-card glossed the invariant as u⋅(x−w)=(rc−r⋅v)u\cdot(x-w) = (\mathscr{r}c - \boldsymbol{\mathscr{r}}\cdot\mathbf{v}), omitting the γ\gamma; the correct evaluation is γ(rc−r⋅v)\gamma(\mathscr{r}c - \boldsymbol{\mathscr{r}}\cdot\mathbf{v}), and the omission would corrupt every downstream power of γ\gamma. Title retitled from the base's "Liénard-Wiechert Potentials" to break the within-track duplicate against the E9 article (track recon; trap-6 policy). Worked numbers (tr=−5.559t_r = -5.559 ns, r=5/3\mathscr{r} = 5/3 m, κ=0.36\kappa = 0.36, V=14.98V = 14.98 V dual-route; γ=1957\gamma = 1957, κ=1.306×10−7\kappa = 1.306\times10^{-7}, u⋅R=7.66×105u\cdot R = 7.66\times10^5 m²/s, V=1.10V = 1.10 mV vs 0.1440.144 nV static) recomputed for this edition, including a machine-precision numerical audit of (2) against the 3+1 pair on a circular β=0.9\beta = 0.9 worldline.

Builds on: Retarded Potentials · Lorenz Gauge and the Wave Equation · Special Relativity Kinematics · Dirac Delta Function Leads to: Fields of a Relativistic Moving Charge · Larmor Formula from the Liénard–Wiechert Fields · Relativistic Generalization of Larmor's Formula · Angular Distribution of Radiation from Accelerated Charges · Synchrotron Radiation · Bremsstrahlung

Why bother with the covariant formulation if the 3+1 potentials already contain the same physics?

Three reasons, in ascending order. Economy: two potentials plus a constraint plus two differentiation identities become one four-vector, one scalar, one lemma. Auditability: the γ\gamma-vs-κ\kappa bookkeeping — which the 3+1 treatment must argue about (is the factor relativistic or not?) — becomes a two-line cancellation you can watch happen; that is how this edition caught the missing γ\gamma in its own base text. Leverage: everything downstream (the field tensor, the covariant Larmor formula, the null-field theorem, synchrotron theory) is built by differentiating this object, and the covariant lemma bank makes those computations mechanical where the 3+1 versions are heroic. The Griffiths-register twin is where the physics is discovered; this formulation is where it becomes equipment.

Is the delta-function in proper time the same trick as the 3+1 auxiliary-time integral?

It is the same trick executed one level deeper. The twin's Derivation B introduced an auxiliary time t′t' to linearize the composed delta and then paid the Jacobian g′(tr)=κg'(t_r) = \kappa. Here the worldline parameter τ\tau plays the auxiliary role from the start, and the Jacobian is ∣ds/dτ∣=2u⋅R|ds/d\tau| = 2u\cdot R with s=(x−r)2s = (x-r)^2 — which, unpacked, is 2γκcr2\gamma\kappa c\mathscr{r}: the twin's κ\kappa times the factors (γ\gamma from the clock conversion, 2cr2c\mathscr{r} from the quadratic constraint) that the covariant packaging keeps attached. Neither route is more correct; the covariant one simply refuses to let you separate the Jacobian from its invariant meaning, which is why transcription errors (like a dropped γ\gamma) are easier to catch.

The formula says the four-potential is parallel to the four-velocity. Is that a gauge statement or physics?

Both, carefully separated. In the Lorenz gauge — which the retarded potentials inhabit, and which is the gauge this whole chapter runs in — the statement is exact physics-of-this-gauge: the only four-vector data at the piercing point, at first order, is uμu^\mu, so Aμ∝uμA^\mu \propto u^\mu is forced by covariance. Under a gauge transformation Aμ→Aμ+∂μλA^\mu \to A^\mu + \partial^\mu\lambda the parallelism can be destroyed, so it is not gauge-invariant content by itself. What is invariant is everything built antisymmetrically from derivatives: the field tensor of the next article, whose velocity part remembers the parallelism as the boosted-Coulomb structure. The honest summary: Aμ∥uμA^\mu \parallel u^\mu is the Lorenz-gauge fingerprint of a structureless point charge — a composite or polarized source adds terms.

Where exactly did the derivation choose retarded over advanced boundary conditions?

In one theta function. Step 3 of Derivation A inserted θ(x0−r0(τ))\theta(x^0 - r^0(\tau)) to select the past sheet of the double cone δ((x−r)2)\delta((x-r)^2) — drop it and the delta happily fires on the future sheet too, producing the advanced potential with denominator u⋅Radvu\cdot R_{\rm adv} (which unpacks to γcr(1+r^a⋅βa)\gamma c\mathscr{r}(1 + \hat{\boldsymbol{\mathscr{r}}}_a\cdot\boldsymbol{\beta}_a), the sign-flipped Jacobian the twin's FAQ derived). The selection is data, not law — the doctrine settled in Advanced Potentials and Causality: quiescent-past boundary conditions pick retarded. The covariant packaging adds one guarantee worth stating: because null ordering is frame-invariant (the invariant-order theorem), the split into past and future sheets is Lorentz-proof — no boost can turn a retarded solution into an advanced one.

Does the covariant potential apply when the charge accelerates violently, or is it a constant-velocity result?

Fully general — that is its quiet strength. Nothing in Derivation A restricted the worldline: rμ(τ)r^\mu(\tau) may wiggle, spiral, or slam to a stop, and (2) remains exact, with all data read at the piercing point τr\tau_r. The acceleration simply does not appear in the potential — it first enters observables through differentiation, because ∂μ\partial_\mu moves the piercing point along the worldline and thereby samples aμ=duμ/dτa^\mu = du^\mu/d\tau (the a⋅Ra\cdot R term in the banked derivative of u⋅Ru\cdot R). This is the covariant restatement of the twin's punchline that all news is news of acceleration: the potential of a swerving charge looks exactly like the potential of a gliding one until you take a derivative. The fields article cashes precisely this observation into the velocity/acceleration split.

How do I actually solve for the retarded proper time in practice?

Analytically only in privileged cases: uniform velocity (a quadratic — worked example 1), uniform proper acceleration (hyperbolic motion, still closed-form), and circular motion (transcendental but one-dimensional). Everything else is numerical, and the standard recipe is exactly what this article's own verification used: parameterize the worldline, define f(τ)=(x−r(τ))2f(\tau) = (x - r(\tau))^2 restricted to r0<x0r^0 < x^0, and root-find — ff is monotonic in the relevant bracket for timelike worldlines (the uniqueness theorem guarantees a single sign change), so bisection or Brent's method converges unconditionally. Production beam codes tabulate worldline history and interpolate the crossing. The one practical warning: at ultrarelativistic head-on geometry the root becomes ill-conditioned in coordinate time (the light cone and worldline nearly parallel — κ∼10−7\kappa \sim 10^{-7} in worked example 2), and solving in proper time, where the Jacobian is the well-behaved u⋅Ru\cdot R, is numerically kinder. The covariant formulation is not just prettier; it conditions better.

Before this

Where this leads