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 and — 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 ( — statics with stale news, after all), a middle term driven by the rate of change of charge density (), and a term driven by the rate of change of current (). The magnetic field carries two: retarded Biot–Savart plus a correction. The 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 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.
- 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 terms are the static laws with retarded sources; the terms — sourced by and — are new, and they are the radiation.
- The equations are manifestly causal and manifestly gauge-free: fields at 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 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 is also good to second order — but the term has no static partner and survives as genuine radiation.
- Radiation requires changing currents: steady and , however arranged, produce only fields; only (accelerating charge) reaches infinity.
- Key idea: Jefimenko's equations are the honest derivatives of the retarded potentials — retarded Coulomb and Biot–Savart plus and terms of range that carry all the radiation and enforce causality with no reference to any gauge.
Jefimenko's equations — obtained by honestly differentiating the retarded potentials — express the exact fields as retarded Coulomb and Biot–Savart integrals plus and terms of 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 ) — 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
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 , 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 and "), 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.
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, , pointing along the stale line of sight — the only reporter statics ever employed. The second reports how fast the charge count was changing: the term, also radial but falling only as — 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 term, , pointing along the (negative) current acceleration rather than the line of sight — the only correspondent whose dispatches survive to infinity, because a field carries a Poynting flux through spheres of area : finite power forever outward. Radiation is the third correspondent, and nothing else. The magnetic desk runs leaner: retarded Biot–Savart () plus one stringer. No charge-density reporter at all — magnetic fields never cared where charge sits, only where it flows.
Two doctrinal points before the algebra. First, gauge-freedom: these formulas contain , , 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 , 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: , so any retarded source obeys .
where: , , dots are time derivatives evaluated at . 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 range divides the world: near a source (, with the source timescale) the static-shaped terms dominate; far away () only the term (and 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 , no memory integrals, no fields-sourcing-fields. Given the source history, two integrals, done.
The goal is equation (1): 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 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 — trivial, since — 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):
Why two pieces: the second is the static gradient acting on the 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 has at fixed positions, the derivative slides through the integral onto the source:
Step 3. Assemble , using to put both contributions over the common constant:
Why each parent matters: the Coulomb term and the term both descend from — they are radial, tied to the line of sight; the term descends from 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. requires the curl of ; the product rule gives a kernel piece and a retardation piece, — the second term from — 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 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 , and first-order errors are not usually invisible. The true answer is a conspiracy.
Go deeper · Derivation B · 2 minWhy Biot–Savart never failed you — the first-order cancellation
The goal is the quasistatic expansion of Jefimenko's , 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 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 :
Why the expansion is legitimate: it requires the source to change smoothly over the light-crossing time to the field point — the quasistatic condition , i.e. field points well inside a source wavelength.
Step 2. Insert into the integrand of (2) and collect orders:
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 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:
— Biot–Savart with the present current, with no first-order error. At 50 Hz and one meter, the neglected terms are of order : 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 , also good through first order — but the third term, , has no static partner to hide in. In the quasistatic zone it is small ( 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
| case | result | meaning |
|---|---|---|
| statics (all dots zero) | (1)–(2) → Coulomb + Biot–Savart | the parent laws, recovered exactly — with retarded sampling now irrelevant |
| quasistatic zone () | static laws with present sources, errors | the vindication theorem: first-order retardation cancels |
| far zone () | only terms survive; transverse part of dominates | radiation: finite Poynting flux to infinity, dipole formula three articles ahead |
| point charge on a trajectory | delta-function sources; integrals collapse with the LW Jacobian | next two articles — Jefimenko is correct but the collapse is subtle, so the potential route is used instead |
| no time-varying current anywhere | only | no radiation without : 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 for oscillation at . Evaluate the dial across three regimes.
Step 1. Power engineering: Hz at m: . The Jefimenko corrections are parts per million of the field — and by the vindication theorem, the surviving error is this number squared: . Circuit theory is safe by twelve orders of magnitude.
Step 2. RF on a bench: GHz at cm: . 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: MHz at m: again — same dial, same verdict. An AM receiver 100 m from the tower already lives in Jefimenko's terms; the Coulomb–Biot–Savart pieces of the tower's field died within a wavelength (300 m) of the mast.
The lesson: one dimensionless dial, , decides which terms of (1)–(2) run the show — the same 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 for (zero before), with A/s. Compute and at distance m exactly from Jefimenko's equations, and evaluate at s. (Griffiths Prob. 10.12 territory — done here as a Jefimenko showcase.)
Step 1. Set up the integral. For a line current, (2) becomes with , the factor coming from . Only elements with , i.e. , have reported in.
Step 2. Watch the integrand collapse. With and :
The retardation of the Biot–Savart term and the entire term cancel exactly — Derivation B's conspiracy, exhibited non-perturbatively (a linear ramp has , so the cancellation has no remainder from the expansion; only the finite reporting window survives as a correction).
Step 3. Integrate: , with . Numbers at s ( mA, m): the square root is , so — 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 term contributes ( on a neutral wire):
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 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 from on the retarded potentials, confirming the three-term structure; (2) re-derive confirming the two-term structure with the term; (3) check the static limits (Coulomb, Biot–Savart); (4) verify the quasistatic expansion showing first-order cancellation (Eq. 10.39's discussion). Constants: exact. Addendum for this edition: the gradient of the retarded 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 ( dial /2.1/2.1; wire pT with horizon factor, mV/m at ) 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 and — 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 — and each expressed directly in and . 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. , and is built from alone — charge density never enters the vector potential, so no manipulation can plant a in . (Deeper reading: static charge produces no magnetism, and where charge sits is exactly the information carries.) The electric field descends from both potentials: the and terms from , the term from . The term is the interesting middle child — radial like Coulomb, -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, , it is secretly current bookkeeping too — some authors rewrite (1) to eliminate it entirely in favor of 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, and are delta functions of a moving argument, and every term of (1)–(2) requires integrating deltas whose arguments contain the integration variable through — 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 and transverse to the propagation direction (the Poynting flux must point outward along ), and the two radial terms of (1) — Coulomb and — can never supply a transverse electric field. The term is the only one free of the line of sight: far away, its component perpendicular to , which is , is the entire radiation field. This is also the Jefimenko-level origin of the in dipole radiation: the transverse projection of a linear scales as — 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 , (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 , integrate over a far sphere — the route the dipole-radiation article takes to its power formula. The moral repeats: Jefimenko for structure, potentials for computation, Poynting for the ledger.