Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

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,tR/c)Rd3r,A(r,t)=μ04πj(r,tR/c)Rd3r, with Rrr.

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δ(rr),j(r,t)=quδ(rr), 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:

 trettRretc,

where  tret is called the retarded time, and  Rret  is the length of the following distance vector

 Rretr(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(ttret)=|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.

Screen Shot 2022-03-02 at 11.31.02 PM.pngFig. 10.1. Graphical solution of Eq. (5).

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(tt),R}{c(tt),rr}.

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(tt),R}γ[c2(tt)uR]γc(RβR)γcR(1βn),

where  nR/R is a unit vector in the observer’s direction,  βu/c is the normalized velocity of the particle, and  γ1/(1u2/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(ttret),

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=RretRret:

 2RretRrettret=2RretRrettret.

If r is fixed, then  Rret/tret(rr)/tret=r/treturet, and Eq. (12) yields

 Rrettret=RretRretRrettret=(nu)ret.

Now let us differentiate the same  Rret  over  t. On one hand, Eq. (11) yields

 Rrett=cctrett.

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

 Rrettret=Rrettrettrett=(nu)rettrett.

Now requiring Eqs. (14) and (15) to give the same result, we get:3

 trett=cc(nu)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.

Screen Shot 2022-03-03 at 12.48.51 AM.png
Fig. 10.2. Deriving Eq. (16) for the case  βn=β.

Let us consider an elementary time interval  dtretdt, 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=dtretdxretc=dtreturetcdtret, so that dtretdt=11uret/c11βret.

This expression coincides with Eq. (16), because in our particular case when the directions of the vectors  βu/c and  nR/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  Bnret×Enret×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πε0utret/Rretβγ2[(1βn)3R2]retq4πε0utretβRretγ2[(1βn)3R3]ret.

Screen Shot 2022-03-03 at 1.28.33 AM.pngFig. 10.3. The linearly moving charge problem.

But according to Eq. (5), the product  βRret may be represented as  βc(ttret)u(ttret). Plugging this expression into Eq. (22), we may eliminate the explicit dependence of  Ex on time  tret:

 Ex=q4πε0utγ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)RretRret+βutret, and we may again use Eq. (5) to get  (1βn)Rret=c(ttret)+βutretctctret/γ2. What remains is to calculate  tret from the self-consistency equation (5), whose square in our current case (Fig. 3) takes the form

 R2retb2+(utret)2=c2(ttret)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(t2b2/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).


This page titled 10.1: Liénard-Wiechert Potentials is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Konstantin K. Likharev via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?