10.1: Liénard-Wiechert Potentials
( \newcommand{\kernel}{\mathrm{null}\,}\)
A convenient starting point for the discussion of radiation by relativistic charges is provided by Eqs. (8.17) for the retarded potentials. In the free space, these formulas, with the integration variable changed from r′ to r′′ for the clarity of what follows, are reduced to
ϕ(r,t)=14πε0∫ρ(r′′,t−R/c)Rd3r′′,A(r,t)=μ04π∫j(r′′,t−R/c)Rd3r′′, with R≡r−r′′.
As a reminder, Eqs. (1a) were derived from the Maxwell equations without any restrictions, and are very natural for situations with continuous distributions of the electric charge and/or current. However, for a single point charge, with
ρ(r,t)=qδ(r−r′),j(r,t)=quδ(r−r′), with u=˙r′,
where r′ is the instantaneous position of the charge, it is more convenient to recast Eqs. (1a) into an explicit form that would not require integration in each particular case. Indeed, as Eqs. (1) show, the potentials at a given observation point {r,t} are contributed by only one specific point {r′(tret),tret} of the particle’s 4D trajectory (called its world line), which satisfies the following condition:
tret≡t−Rretc,
where tret is called the retarded time, and Rret is the length of the following distance vector
Rret≡r(t)−r′(tret)
– physically, the distance covered by the electromagnetic wave from its emission to observation.
The reduction of Eqs. (1a) to such a simpler form, however, requires some care. Indeed, their naïve integration over r′′ would yield the following apparent but wrong results:
ϕ(r,t)=14πε0qRret , i.e. ϕ(r,t)c=μ04πqcRret ;A(r,t)=μ04πquret Rret , (WRONG!)
where uret is the particle’s velocity at the retarded point r′(tret). Eqs. (4) is a good example of how the relativity theory (even the special one :-) cannot be taken too lightly. Indeed, the strings (9.84)-(9.85), formed from the apparent potentials (4), would not obey the Lorentz transform rule (9.91), because according to Eqs. (2)-(3), the distance Rret also depends on the reference frame it is measured in.
In order to correct the error, we need, first of all, discuss the conditions (2)-(3). Combining them (by eliminating Rret ), we get the following equation for tret:
c(t−tret)=|r(t)−r′(tret)|.Retarded time
Figure 1 depicts the graphical solution of this self-consistency equation as the only1 point of intersection of the light cone of the observation point (see Fig. 9.9 and its discussion) and the particle’s world line.

In Eq. (5), just as in Eqs. (1)-(3), all variables have to be measured in the same inertial (“lab”) reference frame, in which the observation point r rests. Now let us write Eqs. (1) for a point charge in another inertial reference frame 0’, whose velocity (as measured in the lab frame) coincides, at the moment t′=tret, with the velocity uret of the charge. In that frame, the charge rests, so that, as we know from the electro- and magnetostatics,
ϕ′=q4πε0R′,A′=0.
(Remember that this R′ may not be equal to Rret, because the latter distance is measured in the “lab” reference frame.) Let us use the identity 1/ε0≡μ0c2 again to rewrite Eqs. (6a) in the form of components of a 4-vector similar in structure to Eq. (4):
ϕ′c=μ04πqcR′,A′=0.
Now it is easy to guess the correct answer for the 4-potential for an arbitrary reference frame:
Aα=μ04πqcuαuβRβ,
where (mostly as a reminder), Aα≡{ϕ/c,A},uα≡γ{c,u}, and Rα is a 4-vector of the inter-event distance, formed similarly to that of a single event – cf. Eq. (9.48):
Rα≡{c(t−t′),R′}≡{c(t−t′),r−r′}.
Indeed, we needed the 4-vector Aα that would:
(i) obey the Lorentz transform,
(ii) have its spatial components Aj scaling, at low velocity, as uj, and
(iii) be reduced to the correct result (6) in the reference frame moving with the charge.
Eq. (7) evidently satisfies all these requirements, because the scalar product in its denominator is just
uβRβ=γ{c,−u}⋅{c(t−t′),R}≡γ[c2(t−t′)−u⋅R]≡γc(R−β⋅R)≡γcR(1−β⋅n),
where n≡R/R is a unit vector in the observer’s direction, β≡u/c is the normalized velocity of the particle, and γ≡1/(1−u2/c2)1/2. In the reference frame of the charge (in which β=0 and γ=1), the
expression (9) is reduced to cR, so that Eq. (7) is correctly reduced to Eq. (6b). Now let us spell out the components of Eq. (7) for the lab frame (in which t′=tret and R=Rret ):
ϕ(r,t)=14πε0q(R−β⋅R)ret=14πε0q[1R(1−β⋅n)]ret,
Liénard-Wiechert potentials
A(r,t)=μ04πq(uR−β⋅R)ret=μ04πqc[βR(1−β⋅n)]ret≡ϕ(r,t)uretc2.
These formulas are called the Liénard-Wiechert potentials.2 In the non-relativistic limit, they coincide with the naïve guess (4), but in the general case include an additional factor 1/(1−β⋅n)ret. Its physical origin may be illuminated by one more formal calculation – whose result we will need anyway. Let us differentiate the geometric relation (5), rewritten as
Rret=c(t−tret),
over tret and then, independently, over t, assuming that r is fixed. For that, let us first differentiate, over tret, both sides of the identity Rret2=Rret⋅Rret:
2Rret∂Rret∂tret=2Rret⋅∂Rret∂tret.
If r is fixed, then ∂Rret/∂tret≡∂(r−r′)/∂tret=−∂r′/∂tret≡−uret, and Eq. (12) yields
∂Rret∂tret=RretRret⋅∂Rret∂tret=−(n⋅u)ret.
Now let us differentiate the same Rret over t. On one hand, Eq. (11) yields
∂Rret∂t=c−c∂tret∂t.
On the other hand, according to Eq. (5), at the partial differentiation over time, i.e. if r is fixed, tret is a function of t alone, so that (using Eq. (13) at the second step), we may write
∂Rret∂tret=∂Rret∂tret∂tret∂t=−(n⋅u)ret∂tret∂t.
Now requiring Eqs. (14) and (15) to give the same result, we get:3
∂tret∂t=cc−(n⋅u)ret=(11−β⋅n)ret.
This relation may be readily re-derived (and more clearly understood) for the simple particular case when the charge’s velocity is directed straight toward the observation point. In this case, its vector u resides in the same space-time plane as the observation point’s world line r=const – say, the plane [x,t], shown in Fig. 2.

Let us consider an elementary time interval dtret≡dt′, during which the particle would travel the space interval dxret=uretdtret. In Fig. 2, the corresponding segment of its world line is shown with a solid vector. The dotted vectors in this figure show the world lines of the radiation emitted by the particle in the beginning and at the end of this interval, and propagating with the speed of light c. As it follows from the drawing, the time interval dt between the instants of the arrival of the radiation from these two points to any time-independent spatial point of observation is
dt=dtret−dxretc=dtret−uretcdtret, so that dtretdt=11−uret/c≡11−βret.
This expression coincides with Eq. (16), because in our particular case when the directions of the vectors β≡u/c and n≡R/R (both taken at time tret) coincide, and hence (β⋅n)ret=βret. Now the general Eq. (16) may be interpreted by saying that the particle’s velocity in the transverse directions (normal to the vector n) is not important for this kinematic effect4 – the fact almost evident from Fig. 1.
So, the additional factor in the Liénard-Wiechert potentials is just the derivative ∂tret/∂t. The reason for its appearance in Eqs. (10) is usually interpreted along the following lines. Let the charge q be spread along the direction of the vector Rret (in Fig. 2, along the x-axis) by an infinitesimal speed-independent interval δxret, so that the linear density λ of its charge is proportional to 1/δxret. Then the time rate of charge’s arrival at some spatial point is λuret=λdxret/dtret. However, the rate of radiation’s arrival at the observation point scales as 1/dt, so that due to the non-zero velocity uret of the particle, this rate differs from the charge arrival rate by the factor of dtret/dt, given by Eq. (16). (If the particle moves toward the observation point, (β⋅n)ret>0, as shown in Fig. 2, this factor is larger than 1.) This radiation compression effect leads to the field change (at (β⋅n)ret>0, its enhancement) by the same factor (16) – as described by Eqs. (10).
So, the 4-vector formalism was very instrumental for the calculation of field potentials. It may be also used to calculate the fields E and B – by plugging Eq. (7) into Eq. (9.124) to calculate the field strength tensor. This calculation yields
Fαβ=μ0q4π1uγRγddτ[Rαuβ−RβuαuδRδ].
Now using Eq. (9.125) to identify the elements of this tensor with the field components, we may bring the result to the following vector form:5
E=q4πε0[n−βγ2(1−β⋅n)3R2+n×{(n−β)×˙β}(1−β⋅n)3cR]ret,
Relativistic particle’s fields
B=nret×Ec, i.e. H=nret×EZ0.
Thus the magnetic and electric fields of a relativistic particle are always proportional and perpendicular to each other, and related just as in a plane wave – cf. Eq. (7.6), with the difference that now the vector nret may be a function of time. Superficially, this result contradicts the electro- and magnetostatics, because for a particle at rest, B should vanish while E stays finite. However, note that according to the Coulomb law for a point charge, in this case E=Enret , so that B∝nret×E∝nret×nret=0. (Actually, in these relations, the subscript “ret” is unnecessary.)
As a sanity check, let us use Eq. (19) as an alternative way to find the electric field of a charge moving without acceleration, i.e. uniformly, along a straight line – see Fig. 9.11a reproduced, with minor changes, in Fig. 3. (This calculation will also illustrate the technical challenges of practical applications of the Liénard-Wiechert formulas for even simple cases.) In this case, the vector β does not change in time, so that the second term in Eq. (19) vanishes, and all we need to do is to spell out the Cartesian components of the first term.
Let us select the coordinate axes and the time origin in the way shown in Fig. 3, and make a clear distinction between the actual position, r′(t)={ut,0,0} of the charged particle at the instant t we are considering, and its position r′(tret) at the retarded instant defined by Eq. (5), i.e. the moment when the particle’s field had to be radiated to reach the observation point r at the given time t, propagating with the speed of light. In these coordinates
β={β,0,0},r={0,b,0},r′(tret)={utret,0,0},nret={cosθ,sinθ,0},
with cosθ=−utret/Rret, so that [(n−β)x]ret=−utret/Rret−β, and Eq. (19) yields, in particular:
Ex=q4πε0−utret/Rret−βγ2[(1−β⋅n)3R2]ret≡q4πε0−utret−βRretγ2[(1−β⋅n)3R3]ret.

But according to Eq. (5), the product βRret may be represented as βc(t−tret)≡u(t−tret). Plugging this expression into Eq. (22), we may eliminate the explicit dependence of Ex on time tret:
Ex=q4πε0−utγ2[(1−β⋅n)R]3ret.
The only non-zero transverse component of the field also has a similar form:
Ey=q4πε0[sinθγ2(1−β⋅n)3R2]ret=q4πε0bγ2[(1−β⋅n)R]3ret,
while Ez=0. From Fig. 3, β−nret=βcosθ=−βutret/Rret, so that (1−β⋅n)Rret≡Rret+βutret, and we may again use Eq. (5) to get (1−β⋅n)Rret=c(t−tret)+βutret≡ct−ctret/γ2. What remains is to calculate tret from the self-consistency equation (5), whose square in our current case (Fig. 3) takes the form
R2ret≡b2+(utret)2=c2(t−tret)2.
This is a simple quadratic equation for tret, which (with the appropriate negative sign before the square root, to get tret<t) yields:
tret=γ2t−[(γ2t)2−γ2(t2−b2/c2)]1/2≡γ2t−γc(u2γ2t2+b2)1/2,
so that the only retarded-function combination that participates in Eqs. (23)-(24) is
[(1−β⋅n)R]ret=cγ2(u2γ2t2+b2)1/2,
and, finally, the electric field components are
Ex=−q4πε0γut(b2+γ2u2t2)3/2,Ey=q4πε0γb(b2+γ2u2t2)3/2,Ez=0.
These are exactly Eqs. (9.139),6 which had been obtained in Sec. 9.5 by much simpler means, without the necessity to solve the self-consistency equation (5). However, that alternative approach was essentially based on the inertial motion of the particle, and cannot be used in problems in which particle moves with acceleration. In such problems, the second term in Eq. (19), dropping with distance more slowly, as 1/Rret, and hence describing wave radiation, is essential and most important.
Reference
1 As Fig. 1 shows, there is always another, “advanced” point {r′(tadv),tadv} of the particle’s world line, with tadv>t, which is also a solution of Eq. (5), but it does not fit Eqs. (1), because the observation, at the point {r,t<tadv}, of the field induced at the advanced point, would violate the causality principle.
2 They were derived in 1898 by Alfred-Marie Liénard and (independently) in 1900 by Emil Wiechert.
3 This relation may be used for an alternative derivation of Eqs. (10) directly from Eqs (1) – the exercise highly recommended to the reader.
4 Note that this effect (linear in β) has nothing to do with the Lorentz time dilation (9.21), which is quadratic in β. (Indeed, all our arguments above referred to the same, lab frame.) Rather, it is close in nature to the Doppler effect.
5 An alternative way of deriving these formulas (highly recommended to the reader as an exercise) is to plug Eqs. (10) into the general relations (9.121), and carry out the required temporal and spatial differentiations directly, using Eq. (16) and its spatial counterpart (which may be derived absolutely similarly):
∇tret=−[nc(1−β⋅n)]ret.
6 A similar calculation of magnetic field components from Eq. (20) gives the results identical to Eqs. (9.140).