Jefimenko's Equations

Two articles ago a warning was posted and left unpaid: you may not retard the fields. The retarded potentials are exact, but differentiating them — the step that produces E\mathbf{E} and B\mathbf{B} — acts on the retarded time hiding inside every integrand, and the naive "Coulomb's law with stale news" misses whole terms. This article does the differentiation honestly, once, in full, and writes down what actually comes out: Jefimenko's equations, the exact, causal, gauge-free field formulas of classical electrodynamics. Every field that has ever existed — around every antenna, atom, and lightning stroke — is an instance of these two integrals. Their anatomy repays study. The electric field carries three terms: the retarded Coulomb field (ρ/r2\rho/\mathscr{r}^2 — statics with stale news, after all), a middle term driven by the rate of change of charge density (ρ˙/cr\dot\rho/c\mathscr{r}), and a term driven by the rate of change of current (J˙/c2r\dot{\mathbf{J}}/c^2\mathscr{r}). The magnetic field carries two: retarded Biot–Savart plus a J˙\dot{\mathbf{J}} correction. The 1/r1/\mathscr{r} terms are the ones the naive guess forgot — and they are the radiation: only they survive at infinity, and only changing currents generate them. Two chapter debts are settled on the spot. The Coulomb-gauge article's causality acquittal gets its receipts: here are the fields, built from past sources only, no gauge anywhere in sight. And a deeper puzzle is resolved — why did quasistatics ever work? The answer is a small miracle of cancellation: for slowly varying currents, the first-order retardation error in Biot–Savart is exactly repaid by the J˙\dot{\mathbf{J}} term, making the static law good to second order. Your intuition survived four chapters of time-dependent physics because Jefimenko's equations were quietly covering for it.

Key takeaways
  • Differentiating the retarded potentials honestly gives the exact fields of any localized source: \mathbf{E}(\mathbf{r},t) = \frac{1}{4\pi\varepsilon_0}\int\left[\frac{\rho(t_r)}{\mathscr{r}^2}\hat{\boldsymbol{\mathscr{r}}} + \frac{\dot\rho(t_r)}{c\,\mathscr{r}}\hat{\boldsymbol{\mathscr{r}}} - \frac{\dot{\mathbf{J}}(t_r)}{c^2\mathscr{r}}\right]d^3r' \mathbf{B}(\mathbf{r},t) = \frac{\mu_0}{4\pi}\int\left[\frac{\mathbf{J}(t_r)}{\mathscr{r}^2} + \frac{\dot{\mathbf{J}}(t_r)}{c\,\mathscr{r}}\right]\times\hat{\boldsymbol{\mathscr{r}}}\;d^3r'.
  • These are the time-dependent generalizations of Coulomb and Biot–Savart: the 1/r21/\mathscr{r}^2 terms are the static laws with retarded sources; the 1/r1/\mathscr{r} terms — sourced by ρ˙\dot\rho and J˙\dot{\mathbf{J}} — are new, and they are the radiation.
  • The equations are manifestly causal and manifestly gauge-free: fields at (r,t)(\mathbf{r},t) are built from charge and current on the past light cone, period — the receipts promised by the Coulomb-gauge causality acquittal.
  • Quasistatic vindication: expanding in the retardation, the first-order error of naive Biot–Savart is exactly cancelled by the J˙\dot{\mathbf{J}} term — the static magnetic law with the present current is correct through first order, which is why circuit-era magnetostatics never failed.
  • The electric analog holds too: Coulomb's law with the present ρ\rho is also good to second order — but the J˙\dot{\mathbf{J}} term has no static partner and survives as genuine radiation.
  • Radiation requires changing currents: steady ρ\rho and J\mathbf{J}, however arranged, produce only 1/r21/\mathscr{r}^2 fields; only J˙≠0\dot{\mathbf{J}} \neq 0 (accelerating charge) reaches infinity.
  • Key idea: Jefimenko's equations are the honest derivatives of the retarded potentials — retarded Coulomb and Biot–Savart plus ρ˙\dot\rho and J˙\dot{\mathbf{J}} terms of 1/r1/\mathscr{r} range that carry all the radiation and enforce causality with no reference to any gauge.
The core insight

Jefimenko's equations — obtained by honestly differentiating the retarded potentials — express the exact fields as retarded Coulomb and Biot–Savart integrals plus ρ˙\dot\rho and J˙\dot{\mathbf{J}} terms of 1/r1/\mathscr{r} range that carry all the radiation.

Where this lives today

Jefimenko's equations are the conscience of computational electromagnetics. Time-domain integral-equation solvers evaluate exactly these convolutions over meshed source histories; near-field-to-far-field transformations in FDTD codes are (1)–(2) applied on a bounding surface; and the "causality audits" of metamaterial and superluminal-media claims consist of checking reported fields against retarded-source integrals of precisely this form. In plasma and beam physics, the equations are the standard for validating particle-in-cell field solvers, because they give the exact fields of a prescribed current history with no gauge or grid artifacts to argue about. They also settled a pedagogical war. The 1991 Griffiths–Heald paper ("Time-dependent generalizations of the Biot–Savart and Coulomb laws") brought the equations into mainstream teaching and pinned down exactly which naive retarded guesses work (potentials; B with present current, to second order) and which fail (retarded Coulomb for E\mathbf{E}) — the structure this article inherited. And the vindication theorem quietly underwrites entire industries: magnetic resonance imaging computes its gradient and RF coils with Biot–Savart at kilohertz-to-megahertz frequencies, transformer design lives at 50–60 Hz, and both are licensed not by "retardation is small" hand-waving but by the exact first-order cancellation — engineering trust, secured by a conspiracy between staleness and change.

Origins · 1 minOliver Heaviside and how this was discovered
OH
Oliver Heaviside
1850–1925 · Camden Town and Devon — the self-taught engineer who put the fields first

Heaviside left school at sixteen, learned electromagnetism as a telegraph operator, and then — jobless, half-deaf, living with his parents — rewrote Maxwell. The twenty quaternion-laden equations of the Treatise became, in Heaviside's hands (1884–85), the four vector equations every student now learns; the curl, divergence, and the modern E\mathbf{E}, B\mathbf{B} notation are substantially his. He championed exactly this article's viewpoint: the fields are the physics, the potentials mere scaffolding (he campaigned, with characteristic venom, to "murder ψ\psi and A\mathbf{A}"), and he computed the fields of moving charges directly — anticipating pieces of the Liénard–Wiechert story, including the deformed field of a fast charge, before either namesake. His operational calculus, invented to solve telegraph-line transients and scorned by mathematicians ("Shall I refuse my dinner because I do not fully understand the process of digestion?"), was vindicated decades later as Laplace-transform theory. The exact causal field equations he never quite assembled were finally organized by Panofsky–Phillips (1962) and named for Jefimenko (1966) — but their fields-first, source-history-in, fields-out architecture is Heaviside's program, completed.

History
1867
— Lorenz's retarded potentials contain the fields implicitly; nobody differentiates them in print
1884
— Heaviside recasts Maxwell into the four modern field equations; begins computing moving-charge fields directly
1888
— Heaviside publishes the distorted field of a moving point charge — fields-first electrodynamics in action
1912
— Schott's treatise systematizes radiation theory, still potential-first
1962
— Panofsky and Phillips's textbook presents the explicit retarded field integrals
1966
— Jefimenko's Electricity and Magnetism makes the equations a pedagogical centerpiece; the name attaches
1991
— Griffiths and Heald's paper fixes the modern teaching canon: which retarded guesses work, which fail
2000
s — time-domain integral-equation solvers make the equations industrial machinery

The physical picture

Think of the field at your location as a report compiled from the past light cone, with three species of correspondent filing to the electric desk. The first correspondent reports where the charge was: the retarded Coulomb term, 1/r21/\mathscr{r}^2, pointing along the stale line of sight — the only reporter statics ever employed. The second reports how fast the charge count was changing: the ρ˙\dot\rho term, also radial but falling only as 1/r1/\mathscr{r} — a longer-range stringer whose job, it will turn out, is mostly to correct the first one's stale bookkeeping. The third reports how the currents were accelerating: the J˙\dot{\mathbf{J}} term, 1/r1/\mathscr{r}, pointing along the (negative) current acceleration rather than the line of sight — the only correspondent whose dispatches survive to infinity, because a 1/r1/\mathscr{r} field carries a 1/r21/\mathscr{r}^2 Poynting flux through spheres of area r2\mathscr{r}^2: finite power forever outward. Radiation is the third correspondent, and nothing else. The magnetic desk runs leaner: retarded Biot–Savart (J×r^/r2\mathbf{J}\times\hat{\boldsymbol{\mathscr{r}}}/\mathscr{r}^2) plus one J˙×r^/cr\dot{\mathbf{J}}\times\hat{\boldsymbol{\mathscr{r}}}/c\mathscr{r} stringer. No charge-density reporter at all — magnetic fields never cared where charge sits, only where it flows.

Figure
Anatomy of Jefimenko's electric field — three terms, three ranges, three jobs; only the current-acceleration term survives to infinity
E(r, t) = (1/4πε₀) ∫ [ term 1 + term 2 + term 3 ] dτ′ — all sampled at t_r ρ(t_r) 𝓇̂ / 𝓇² retarded Coulomb where the charge WAS range 1/𝓇² — near zone statics' only term ρ̇(t_r) 𝓇̂ / c𝓇 charge-rate term how the count was CHANGING range 1/𝓇 — but radial mostly corrects term 1's staleness − J̇(t_r) / c²𝓇 current-acceleration term not radial — along −J̇ range 1/𝓇 — reaches infinity THE radiation 𝓇 1/𝓇² terms: die fast — power through far spheres → 0 1/𝓇 terms: flux 𝓇²·(1/𝓇)² = constant — finite power to infinity B's anatomy is leaner: retarded Biot–Savart (J×𝓇̂/𝓇²) + one J̇×𝓇̂/c𝓇 term — no ρ reporter at the magnetic desk gauge-free and causal by construction: every term reads the past light cone only

Two doctrinal points before the algebra. First, gauge-freedom: these formulas contain ρ\rho, J\mathbf{J}, and geometry — no potentials, no gauge condition, nothing to choose. They are what the Coulomb-gauge causality acquittal appealed to as "the receipts": however weird a gauge's bookkeeping (instantaneous VV, nonlocal transverse currents), the fields it computes must equal these integrals, which manifestly read only the past light cone. Second, historical irony: the exact field equations of electrodynamics were not written in a textbook until 1962 (Panofsky & Phillips) and not named until Oleg Jefimenko's 1966 text made them a teaching centerpiece — nine decades after Maxwell. Nobody needed them: potentials were easier for every practical calculation, and the quasistatic cancellation below meant the static laws quietly kept working far beyond their warranty.

The equations, earned

Everything flows from the retarded potentials of two articles back, plus their differentiation kit: ∇tr=−r^/c\nabla t_r = -\hat{\boldsymbol{\mathscr{r}}}/c, so any retarded source obeys ∇f(tr)=−f˙ r^/c\nabla f(t_r) = -\dot f\,\hat{\boldsymbol{\mathscr{r}}}/c.

E(r,t)=14πε0∫[ρ(r′,tr)r2 r^+ρ˙(r′,tr)c r r^−J˙(r′,tr)c2 r]d3r′,(1)\mathbf{E}(\mathbf{r},t) = \frac{1}{4\pi\varepsilon_0}\int\left[\frac{\rho(\mathbf{r}',t_r)}{\mathscr{r}^2}\,\hat{\boldsymbol{\mathscr{r}}} + \frac{\dot\rho(\mathbf{r}',t_r)}{c\,\mathscr{r}}\,\hat{\boldsymbol{\mathscr{r}}} - \frac{\dot{\mathbf{J}}(\mathbf{r}',t_r)}{c^2\,\mathscr{r}}\right]d^3r', \qquad (1)
B(r,t)=μ04π∫[J(r′,tr)r2+J˙(r′,tr)c r]×r^  d3r′.(2)\mathbf{B}(\mathbf{r},t) = \frac{\mu_0}{4\pi}\int\left[\frac{\mathbf{J}(\mathbf{r}',t_r)}{\mathscr{r}^2} + \frac{\dot{\mathbf{J}}(\mathbf{r}',t_r)}{c\,\mathscr{r}}\right]\times\hat{\boldsymbol{\mathscr{r}}}\;d^3r'. \qquad (2)

where: r=∣r−r′∣\mathscr{r} = |\mathbf{r}-\mathbf{r}'|, tr=t−r/ct_r = t - \mathscr{r}/c, dots are time derivatives evaluated at trt_r. Read the structure before deriving it: set all dots to zero and (1)–(2) collapse to Coulomb and Biot–Savart exactly — statics is the no-news-is-new-news limit. The dotted terms carry the dynamics, and their 1/r1/\mathscr{r} range divides the world: near a source (r≪cτ\mathscr{r} \ll c\tau, with τ\tau the source timescale) the static-shaped terms dominate; far away (r≫cτ\mathscr{r} \gg c\tau) only the J˙\dot{\mathbf{J}} term (and ρ˙\dot\rho term, though it is radial and — after the transversality bookkeeping of the coming radiation articles — drops out of the far-zone transverse field) matters. Note also what is absent: no ρ¨\ddot\rho, no memory integrals, no fields-sourcing-fields. Given the source history, two integrals, done.

Derivation AThe electric field — differentiating the retarded potentials honestly

The goal is equation (1): E=−∇V−∂tA\mathbf{E} = -\nabla V - \partial_t\mathbf{A} evaluated on the retarded potentials with every retarded-time dependence tracked. The strategy is pure bookkeeping discipline, and half the work is already banked: the gradient of the retarded VV was computed in the verification derivation of Retarded Potentials (its Step 1), and is recalled rather than recomputed. What remains is the time derivative of A\mathbf{A} — trivial, since ∂tr/∂t=1\partial t_r/\partial t = 1 — and assembly. Watching which term comes from which parent is the point: the naive "retarded Coulomb" guess corresponds to keeping only the kernel-gradient piece and forgetting that moving the field point also re-dates every dispatch.

Step 1. Recall the gradient of the retarded scalar potential (banked two articles ago):

∇V=14πε0∫[−ρ˙(tr)cr^r−ρ(tr)r^r2]d3r′.\nabla V = \frac{1}{4\pi\varepsilon_0}\int\left[-\frac{\dot\rho(t_r)}{c}\frac{\hat{\boldsymbol{\mathscr{r}}}}{\mathscr{r}} - \rho(t_r)\frac{\hat{\boldsymbol{\mathscr{r}}}}{\mathscr{r}^2}\right]d^3r'.

Why two pieces: the second is the static gradient acting on the 1/r1/\mathscr{r} kernel; the first exists because stepping the field point closer re-dates the news fresher — the derivative of staleness. This first piece is precisely what the naive guess forgets.

Step 2. Differentiate the vector potential in time. Since tr=t−r/ct_r = t - \mathscr{r}/c has ∂tr/∂t=1\partial t_r/\partial t = 1 at fixed positions, the derivative slides through the integral onto the source:

∂A∂t=μ04π∫J˙(tr)r d3r′.\frac{\partial\mathbf{A}}{\partial t} = \frac{\mu_0}{4\pi}\int\frac{\dot{\mathbf{J}}(t_r)}{\mathscr{r}}\,d^3r'.

Step 3. Assemble E=−∇V−∂tA\mathbf{E} = -\nabla V - \partial_t\mathbf{A}, using μ0=1/ε0c2\mu_0 = 1/\varepsilon_0 c^2 to put both contributions over the common constant:

E=14πε0∫[ρ(tr)r2r^+ρ˙(tr)crr^−J˙(tr)c2r]d3r′.\mathbf{E} = \frac{1}{4\pi\varepsilon_0}\int\left[\frac{\rho(t_r)}{\mathscr{r}^2}\hat{\boldsymbol{\mathscr{r}}} + \frac{\dot\rho(t_r)}{c\mathscr{r}}\hat{\boldsymbol{\mathscr{r}}} - \frac{\dot{\mathbf{J}}(t_r)}{c^2\mathscr{r}}\right]d^3r'.

Why each parent matters: the Coulomb term and the ρ˙\dot\rho term both descend from VV — they are radial, tied to the line of sight; the J˙\dot{\mathbf{J}} term descends from A\mathbf{A} and points along the current's acceleration, with no obligation to the line of sight. That directional freedom is what will let it be transverse in the far zone — the geometric prerequisite for radiation.

Step 4. The magnetic companion. B=∇×A\mathbf{B} = \nabla\times\mathbf{A} requires the curl of J(tr)/r\mathbf{J}(t_r)/\mathscr{r}; the product rule gives a kernel piece and a retardation piece, ∇×[J(tr)/r]=J(tr)×r^/r2+J˙(tr)×r^/(cr)\nabla\times[\mathbf{J}(t_r)/\mathscr{r}] = \mathbf{J}(t_r)\times\hat{\boldsymbol{\mathscr{r}}}/\mathscr{r}^2 + \dot{\mathbf{J}}(t_r)\times\hat{\boldsymbol{\mathscr{r}}}/(c\mathscr{r}) — the second term from ∇×J(tr)=−1c∇r×J˙=J˙×r^/c\nabla\times\mathbf{J}(t_r) = -\tfrac{1}{c}\nabla\mathscr{r}\times\dot{\mathbf{J}} = \dot{\mathbf{J}}\times\hat{\boldsymbol{\mathscr{r}}}/c — and integrating yields (2). Verified: matches Griffiths Eqs. 10.36–10.38 and Jackson §6.5 (Eqs. 6.55–6.56); setting all time derivatives to zero recovers Coulomb and Biot–Savart exactly; the pair reproduces the wire fields of the retarded-potentials article when applied to its switch-on problem.

The quasistatic vindication

Here is the puzzle four chapters of magnetostatics quietly posed. Biot–Savart was derived for steady currents; engineers apply it to 50 Hz transformers, audio circuits, even MHz coils, with the present current I(t)I(t) pasted in — and it works superbly. Why? The naive answer — "retardation is small" — is only half the story, because retardation enters at first order in r/cτ\mathscr{r}/c\tau, and first-order errors are not usually invisible. The true answer is a conspiracy.

Figure
The quasistatic conspiracy — the staleness error of retarded Biot–Savart and the J-dot term cancel at first order, leaving the present-current static law accurate to second order
Expanding Jefimenko's B about the present time t J(t_r)/𝓇² term J(t)/𝓇² − (1/c𝓇)·J̇(t) + (𝓇⁰/2c²)·J̈(t) static law + staleness error + curvature error J̇(t_r)/c𝓇 term + (1/c𝓇)·J̇(t) − (𝓇⁰/c²)·J̈(t) the new dynamical term, expanded first-order J̇ pieces cancel EXACTLY B = Biot–Savart with the PRESENT current × [1 + O((𝓇/cτ)²)] 50 Hz at 1 m: error ~ 10⁻¹² — magnetostatics was an identity with invisible corrections linear ramp I = Kt: J̈ = 0, cancellation exact (worked example 2 catches it in the act)
Go deeper · Derivation B · 2 minWhy Biot–Savart never failed you — the first-order cancellation
Derivation BWhy Biot–Savart never failed you — the first-order cancellation

The goal is the quasistatic expansion of Jefimenko's B\mathbf{B}, equation (2), for sources varying slowly on the light-crossing time — and the punch line that the naive Biot–Savart law with the present current is correct not to zeroth but to second order in the retardation. The strategy: Taylor-expand both retarded terms about the present time tt and watch the first-order pieces annihilate. This cancellation is the magnetic sibling of the middle-term cancellation in the retarded-potential verification, and it is the honest license behind every quasistatic magnetic calculation in the earlier chapters.

Step 1. Expand the sources about the present. With tr−t=−r/ct_r - t = -\mathscr{r}/c:

J(tr)=J(t)−rcJ˙(t)+r22c2J¨(t)−⋯ ,J˙(tr)=J˙(t)−rcJ¨(t)+⋯\mathbf{J}(t_r) = \mathbf{J}(t) - \frac{\mathscr{r}}{c}\dot{\mathbf{J}}(t) + \frac{\mathscr{r}^2}{2c^2}\ddot{\mathbf{J}}(t) - \cdots, \qquad \dot{\mathbf{J}}(t_r) = \dot{\mathbf{J}}(t) - \frac{\mathscr{r}}{c}\ddot{\mathbf{J}}(t) + \cdots

Why the expansion is legitimate: it requires the source to change smoothly over the light-crossing time to the field point — the quasistatic condition r≪cτ\mathscr{r} \ll c\tau, i.e. field points well inside a source wavelength.

Step 2. Insert into the integrand of (2) and collect orders:

J(tr)r2+J˙(tr)c r=J(t)r2−J˙(t)c r+J˙(t)c r⏟= 0+J¨(t)2c2−J¨(t)c2⏟− J¨(t)/2c2+⋯\frac{\mathbf{J}(t_r)}{\mathscr{r}^2} + \frac{\dot{\mathbf{J}}(t_r)}{c\,\mathscr{r}} = \frac{\mathbf{J}(t)}{\mathscr{r}^2} \underbrace{- \frac{\dot{\mathbf{J}}(t)}{c\,\mathscr{r}} + \frac{\dot{\mathbf{J}}(t)}{c\,\mathscr{r}}}_{=\ 0} + \underbrace{\frac{\ddot{\mathbf{J}}(t)}{2c^2} - \frac{\ddot{\mathbf{J}}(t)}{c^2}}_{-\,\ddot{\mathbf{J}}(t)/2c^2} + \cdots

The first-order terms cancel exactly: the staleness correction to the Biot–Savart term is repaid, to the penny, by the leading piece of the J˙\dot{\mathbf{J}} term. Why this is a conspiracy and not luck: the two terms have entirely different parentage — one is the Taylor error of sampling the past, the other a genuine new term of the dynamical theory — yet their coefficients are locked by the same structure that made the retarded potentials solve the wave equation. Nature arranged for the static law to be better than its derivation.

Step 3. Read off the theorem:

B(r,t)=μ04π∫J(r′,t)×r^r2 d3r′×[1+O ⁣((rcτ) ⁣2)]\mathbf{B}(\mathbf{r},t) = \frac{\mu_0}{4\pi}\int \frac{\mathbf{J}(\mathbf{r}',t)\times\hat{\boldsymbol{\mathscr{r}}}}{\mathscr{r}^2}\,d^3r' \times\left[1 + \mathcal{O}\!\left(\left(\frac{\mathscr{r}}{c\tau}\right)^{\!2}\right)\right]

— Biot–Savart with the present current, with no first-order error. At 50 Hz and one meter, the neglected terms are of order (2πr/λ)2∼10−12(2\pi \mathscr{r}/\lambda)^2 \sim 10^{-12}: magnetostatics was never really an approximation in the lab; it was an identity with invisible corrections.

Step 4. The electric side, and where the vindication ends. The same expansion applied to (1)'s first two terms gives Coulomb's law with the present ρ\rho, also good through first order — but the third term, −J˙/c2r-\dot{\mathbf{J}}/c^2\mathscr{r}, has no static partner to hide in. In the quasistatic zone it is small (O(v/c)2\mathcal{O}(v/c)^2 relative); in the far zone it is everything. The vindication is a near-zone theorem; radiation is precisely the part of the field that no present-time static law can mimic. Verified: Griffiths §10.2.2 (Eq. 10.39 and the surrounding discussion); the exact wire computation of worked example 2 exhibits the cancellation non-perturbatively — the two integrand pieces merge into the static form exactly.

Limiting cases

caseresultmeaning
statics (all dots zero)(1)–(2) → Coulomb + Biot–Savartthe parent laws, recovered exactly — with retarded sampling now irrelevant
quasistatic zone (r≪cτ\mathscr{r}\ll c\tau)static laws with present sources, errors O(r/cτ)2\mathcal{O}(\mathscr{r}/c\tau)^2the vindication theorem: first-order retardation cancels
far zone (r→∞\mathscr{r}\to\infty)only 1/r1/\mathscr{r} terms survive; transverse part of −J˙-\dot{\mathbf{J}} dominates E\mathbf{E}radiation: finite Poynting flux to infinity, dipole formula three articles ahead
point charge on a trajectorydelta-function sources; integrals collapse with the LW Jacobiannext two articles — Jefimenko is correct but the collapse is subtle, so the potential route is used instead
no time-varying current anywhereE,B∼1/r2\mathbf{E}, \mathbf{B} \sim 1/\mathscr{r}^2 onlyno radiation without J˙\dot{\mathbf{J}}: accelerating charge is the only transmitter nature sells

Worked examples

1 · The correction dial — when do the new terms matter? (confidence builder)

Setup. The dotted terms enter relative to the static ones at the ratio (ρ˙/cr)/(ρ/r2)=rρ˙/cρ∼ωr/c=2πr/λ\left(\dot\rho/c\mathscr{r}\right)/\left(\rho/\mathscr{r}^2\right) = \mathscr{r}\dot\rho/c\rho \sim \omega\mathscr{r}/c = 2\pi\mathscr{r}/\lambda for oscillation at ω\omega. Evaluate the dial across three regimes.

Step 1. Power engineering: f=50f = 50 Hz at r=1\mathscr{r} = 1 m: ωr/c=(314)(1)/(3×108)=1.0×10−6\omega\mathscr{r}/c = (314)(1)/(3\times10^8) = 1.0\times10^{-6}. The Jefimenko corrections are parts per million of the field — and by the vindication theorem, the surviving error is this number squared: 10−1210^{-12}. Circuit theory is safe by twelve orders of magnitude.

Step 2. RF on a bench: f=1f = 1 GHz at r=10\mathscr{r} = 10 cm: ωr/c=2.1\omega\mathscr{r}/c = 2.1. Order unity — the "correction" terms are the field. There is no quasistatic reading of a microwave circuit at arm's length; this is why RF engineering is transmission-line and antenna theory rather than circuit theory.

Step 3. Broadcast reception: f=1f = 1 MHz at r=100\mathscr{r} = 100 m: ωr/c=2.1\omega\mathscr{r}/c = 2.1 again — same dial, same verdict. An AM receiver 100 m from the tower already lives in Jefimenko's 1/r1/\mathscr{r} terms; the 1/r21/\mathscr{r}^2 Coulomb–Biot–Savart pieces of the tower's field died within a wavelength (300 m) of the mast.

The lesson: one dimensionless dial, 2πr/λ2\pi\mathscr{r}/\lambda, decides which terms of (1)–(2) run the show — the same d/λd/\lambda logic as the retardation budget, now applied to field structure rather than source sampling.

2 · The steadily ramping wire — exact fields, and the cancellation caught in the act (exam level)

Setup. An infinite straight wire carries I(t)=KtI(t) = Kt for t>0t > 0 (zero before), with K=100K = 100 A/s. Compute B\mathbf{B} and E\mathbf{E} at distance s=0.5s = 0.5 m exactly from Jefimenko's equations, and evaluate at t=1 μt = 1\ \mus. (Griffiths Prob. 10.12 territory — done here as a Jefimenko showcase.)

Step 1. Set up the B\mathbf{B} integral. For a line current, (2) becomes Bφ=μ04π∫dz s[I(tr)r3+I˙(tr)c r2]B_\varphi = \frac{\mu_0}{4\pi}\int dz\, s\left[\frac{I(t_r)}{\mathscr{r}^3} + \frac{\dot I(t_r)}{c\,\mathscr{r}^2}\right] with r=s2+z2\mathscr{r} = \sqrt{s^2+z^2}, the factor s/rs/\mathscr{r} coming from ∣z^×r^∣|\hat{\mathbf{z}}\times\hat{\boldsymbol{\mathscr{r}}}|. Only elements with tr>0t_r > 0, i.e. ∣z∣<c2t2−s2|z| < \sqrt{c^2t^2 - s^2}, have reported in.

Step 2. Watch the integrand collapse. With I(tr)=K(t−r/c)I(t_r) = K(t - \mathscr{r}/c) and I˙=K\dot I = K:

K(t−r/c)r3+Kc r2=Ktr3−Kc r2+Kc r2=Ktr3.\frac{K(t - \mathscr{r}/c)}{\mathscr{r}^3} + \frac{K}{c\,\mathscr{r}^2} = \frac{Kt}{\mathscr{r}^3} - \frac{K}{c\,\mathscr{r}^2} + \frac{K}{c\,\mathscr{r}^2} = \frac{Kt}{\mathscr{r}^3}.

The retardation of the Biot–Savart term and the entire I˙\dot I term cancel exactly — Derivation B's conspiracy, exhibited non-perturbatively (a linear ramp has I¨=0\ddot I = 0, so the cancellation has no remainder from the expansion; only the finite reporting window survives as a correction).

Step 3. Integrate: Bφ=μ0Kt s4π∫−z0z0dz(s2+z2)3/2=μ0Kt2πs1−s2c2t2B_\varphi = \frac{\mu_0 K t\, s}{4\pi}\int_{-z_0}^{z_0}\frac{dz}{(s^2+z^2)^{3/2}} = \frac{\mu_0 K t}{2\pi s}\sqrt{1 - \frac{s^2}{c^2t^2}}, with z0=c2t2−s2z_0 = \sqrt{c^2t^2 - s^2}. Numbers at t=1 μt = 1\ \mus (I=0.1I = 0.1 mA, ct=300ct = 300 m): the square root is 1−1.4×10−61 - 1.4\times10^{-6}, so B=μ0(10−4)2π(0.5)=40.0 pTB = \frac{\mu_0(10^{-4})}{2\pi(0.5)} = \mathbf{40.0\ pT} — Biot–Savart with the present current, to a part per million, three hundred meters of causal horizon notwithstanding.

Step 4. The electric field has no such luck. From (1), only the J˙\dot{\mathbf{J}} term contributes (ρ=0\rho = 0 on a neutral wire):

Ez=−μ0K2π ln⁡ ⁣(ct+c2t2−s2s)=−μ0(100)2πln⁡(1200)=−0.142 mV/m,E_z = -\frac{\mu_0 K}{2\pi}\,\ln\!\left(\frac{ct + \sqrt{c^2t^2 - s^2}}{s}\right) = -\frac{\mu_0(100)}{2\pi}\ln(1200) = -\mathbf{0.142\ mV/m},

anti-parallel to the growing current (Lenz would approve) and growing logarithmically forever as the reporting window widens — a term with no static partner and no cancellation, the germ of the inductive electric field.

The lesson: the magnetic field forgets the dynamics (cancellation), the electric field remembers them (the lonely J˙\dot{\mathbf{J}} term). One source, both behaviors, straight from the anatomy of (1)–(2).

Going further

References. Griffiths, Introduction to Electrodynamics, 4th ed., §10.2.2 · Jackson, Classical Electrodynamics, 3rd ed., §6.5 · Zangwill, Modern Electrodynamics, §20.4 · Jefimenko, Electricity and Magnetism, 2nd ed. (1989) · Griffiths & Heald, Am. J. Phys. 59, 111 (1991)

Provenance. Equations follow Griffiths §10.2.2 (Eqs. 10.36–10.39) and Jackson §6.5. A verifier should: (1) re-derive E\mathbf{E} from −∇V−∂tA-\nabla V - \partial_t\mathbf{A} on the retarded potentials, confirming the three-term structure; (2) re-derive B=∇×A\mathbf{B} = \nabla\times\mathbf{A} confirming the two-term structure with the J˙×r^/cr\dot{\mathbf{J}}\times\hat{\boldsymbol{\mathscr{r}}}/c\mathscr{r} term; (3) check the static limits (Coulomb, Biot–Savart); (4) verify the quasistatic expansion showing first-order cancellation (Eq. 10.39's discussion). Constants: μ0ε0=1/c2\mu_0\varepsilon_0 = 1/c^2 exact. Addendum for this edition: the gradient of the retarded VV is recalled from Retarded Potentials (its verification derivation, Step 1) rather than recomputed — division of labor within the causality spine; this edition owns the field assembly (Derivation A), the quasistatic vindication with its exact cancellation (Derivation B, delivering the "why statics worked" payoff), the term-anatomy/range audit, and the Coulomb-gauge receipts connection promised by Coulomb Gauge's Derivation B. The point-charge specialization is deferred by name to Liénard–Wiechert Potentials and Fields of a Moving Point Charge (the potential route being calculationally superior — see the warning). Worked-example numbers (ωr/c\omega\mathscr{r}/c dial 10−610^{-6}/2.1/2.1; wire B=40.0B = 40.0 pT with 1−1.4×10−61-1.4\times10^{-6} horizon factor, E=−0.142E = -0.142 mV/m at ln⁡1200=7.09\ln 1200 = 7.09) verified symbolically and numerically for this edition.

Builds on: Retarded Potentials · Inhomogeneous Wave Equations for Potentials · Lorenz Gauge · Coulomb Gauge (causality receipts) Leads to: Liénard–Wiechert Potentials · Fields of a Moving Point Charge · Electric Dipole Radiation · Magnetic Dipole Radiation

Why are the exact field equations of electrodynamics named after a 20th-century textbook author?

Because for ninety years nobody needed to write them down. The potential route (Lorenz 1867 onward) was calculationally superior for every practical problem — multipole expansions, point charges, antennas all flow more easily from VV and A\mathbf{A} — and the quasistatic cancellation meant experimenters saw no anomalies that demanded exact field formulas. The integrals appeared in scattered forms (Heaviside's moving-charge fields, Schott's treatise, research literature), reached textbook print in Panofsky & Phillips (1962), and became a named teaching centerpiece with Oleg Jefimenko's 1966 Electricity and Magnetism, whose insistent causal reading — fields determined by past sources, full stop — struck a chord. Griffiths & Heald (1991) then fixed the modern canon. The naming is genuinely anachronistic, a fact worth knowing: priority in physics tracks use, and the community only reached for these equations when conceptual questions (causality, what radiates, what quasistatics really assumes) became teaching priorities.

Do Jefimenko's equations prove that electric fields are never produced by changing magnetic fields?

They prove something narrower and more useful: that the fields of a given source history can be written without fields appearing as sources — E\mathbf{E} and B\mathbf{B} each expressed directly in ρ\rho and J\mathbf{J}. Jefimenko himself read this ontologically ("Faraday induction" as correlation, not causation, both fields being co-effects of the sources), and the reading has real pedagogical value against loose "the changing B creates E, which creates B..." talk. But the strong metaphysical claim overreaches: the formulas presuppose retarded boundary data (a free incoming wave appears in neither integral yet is perfectly physical), and in the self-consistent problem — sources responding to fields responding to sources — the causal graph is genuinely entangled. The defensible summary: given the source history and quiescent past, sources suffice to determine the fields; whether that licenses "only sources cause fields" is philosophy the equations neither need nor settle.

Why is there a ρ̇ term in E but no ρ term at all in B?

Trace the parentage in Derivation A. B=∇×A\mathbf{B} = \nabla\times\mathbf{A}, and A\mathbf{A} is built from J\mathbf{J} alone — charge density never enters the vector potential, so no manipulation can plant a ρ\rho in B\mathbf{B}. (Deeper reading: static charge produces no magnetism, and where charge sits is exactly the information ρ\rho carries.) The electric field descends from both potentials: the ρ\rho and ρ˙\dot\rho terms from −∇V-\nabla V, the J˙\dot{\mathbf{J}} term from −∂tA-\partial_t\mathbf{A}. The ρ˙\dot\rho term is the interesting middle child — radial like Coulomb, 1/r1/\mathscr{r}-ranged like radiation. Its main structural role is revealed by the quasistatic expansion: it is the term that repairs the staleness of the retarded Coulomb piece (and via continuity, ρ˙=−∇′⋅J\dot\rho = -\nabla'\cdot\mathbf{J}, it is secretly current bookkeeping too — some authors rewrite (1) to eliminate it entirely in favor of J\mathbf{J} terms).

If these equations are exact, why does the chapter now go back to potentials for the moving point charge?

Because exactness and convenience are different virtues. For a point charge, ρ\rho and J\mathbf{J} are delta functions of a moving argument, and every term of (1)–(2) requires integrating deltas whose arguments contain the integration variable through trt_r — the collapse produces Jacobian factors, and doing this correctly for five terms, then combining, is an error-factory (published wrong answers exist). The potential route runs the treacherous delta collapse once, at the level of two scalar integrals (the Liénard–Wiechert derivation, next article), and then differentiates — heavy algebra but honest algebra. The two routes provably agree; the chapter's division is: Jefimenko for structure and doctrine, Liénard–Wiechert for point-charge computation. This is also the general lesson of the pair of articles: the fields-first formulation is conceptually supreme and calculationally third-best.

The J̇ term in E points along the current's acceleration, not along the line of sight. Why is that so important?

Because radial fields cannot radiate. In the far zone, energy transport needs E\mathbf{E} and B\mathbf{B} transverse to the propagation direction (the Poynting flux must point outward along r^\hat{\boldsymbol{\mathscr{r}}}), and the two radial terms of (1) — Coulomb and ρ˙\dot\rho — can never supply a transverse electric field. The J˙\dot{\mathbf{J}} term is the only one free of the line of sight: far away, its component perpendicular to r^\hat{\boldsymbol{\mathscr{r}}}, which is −J˙⊥/c2r-\dot{\mathbf{J}}_\perp/c^2\mathscr{r}, is the entire radiation field. This is also the Jefimenko-level origin of the sin⁡θ\sin\theta in dipole radiation: the transverse projection of a linear J˙\dot{\mathbf{J}} scales as sin⁡θ\sin\theta — the same projection the Coulomb-gauge article extracted from the transverse-current structure, and the pattern the dipole article will integrate into its power formula. One term, projected: that is what every antenna in the world is emitting.

How do Jefimenko's equations relate to Maxwell's equations — are they equivalent?

One-way traffic with a toll. Every Jefimenko field solves Maxwell's equations with the given ρ\rho, J\mathbf{J} (they descend from the retarded potentials, which were verified against the wave equations, which encode Maxwell in Lorenz gauge). But Maxwell's equations have more solutions: add any source-free field (an incoming wave from infinity, a cavity mode set ringing before your data began) and Maxwell is still satisfied while Jefimenko's integrals know nothing of it. The precise statement: Jefimenko = Maxwell + localized sources + quiescent past (no incoming radiation). That toll is exactly the retarded boundary condition prosecuted in the advanced-potentials article, and it is why Jefimenko's equations cannot be "the fundamental equations" — they are the fundamental solution of the fundamental equations, for the class of data nature usually supplies.

Can I see energy conservation directly in Jefimenko's equations?

Not easily, and the difficulty is instructive. Poynting's theorem is a local identity among fields, so it holds automatically for Jefimenko fields — but verifying it directly from (1)–(2) means multiplying retarded integrals against retarded integrals, producing double integrals over pairs of source elements with different retarded times: correlations between what one piece of the source did and what another piece did earlier. That structure is real physics (radiation is precisely the part of the field energy that escapes the mutual bookkeeping of the source's parts), but as an exercise it is brutal, and nobody audits energy that way. Practice: compute fields via potentials, form S=E×B/μ0\mathbf{S} = \mathbf{E}\times\mathbf{B}/\mu_0, integrate over a far sphere — the route the dipole-radiation article takes to its ω4\omega^4 power formula. The moral repeats: Jefimenko for structure, potentials for computation, Poynting for the ledger.

Before this

Where this leads