emii2sol_13.pdf

7/30/2019 emii2sol_13.pdf

1/24


2/24

after noting that ii = = 1 since is an arbitrary real unit vector. We now observethat mjqqmk = 0 since mq is symmetric under the interchange of m and q whereas qmk isantisymmetric under the same interchange of indices. Thus, eq. (2) yields

(S)

3

jk = qqjk = q(Sq)jk = ( S)jk ,which establishes eq. (1). Note that we could have achieved the same result by writing outexplicitly

S=

0 0 0 00 0 3 20 3 0 10 2 1 0

,

and then computing ( S)3 via matrix multiplication. Indeed,

(S)

2

=

0 0 0 00

22

23 12 13

0 12 21 23 230 13 23 21 22

,and

( S)3 = ( S)2 S= (21 + 22 + 23)

0 0 0 00 0 3 20 3 0 10 2 1 0

= S,

after using the fact that is a real unit 3-vector.To establish eq. (2), we can follow either of the above two methods. Employing the first

method above, we make use of the block-matrix form for K

( K)jk =

0 k

j 0jk

, (3)

where 0jk stands for the matrix elements of the 3 3 zero matrix. In particular, j labels therow and k labels the column. Then,

( K)3jk =

0

j 0j

0 i

0i

0 k

i 0ik

=

0i

0Tj ji

0 k

i 0ik

= 1 0i

0Tj ji

0 ki 0ik

= 0 kj

0jk

= 0 kj 0jk

= ( K)jk ,after using the fact that is is real unit vector. Thus eq. (2) is established. Alternatively, wecan write out the explicit form for the matrix,

K=

0 1 2

3

1 0 0 02 0 0 03 0 0 0

,

2


3/24

and then compute ( K)3 via matrix multiplication. Indeed,

( K)2 = 21 +

22 +

23 0 0 0

0 21 12

13

0 12 22 230 1

3

23

23

,

and

( K)3 = ( K)2 K= ( 21 + 22 +

23 )

0 1 2

3

1 0 0 02 0 0 03 0 0 0

= K,

after using the fact that is a real unit 3-vector.

(b) Use the result of part (a) to show that:

exp K

= I K sinh + ( K)2 [cosh 1] ,

where I is the 4 4 identity matrix.

We employ the series expansion for the exponential (which defines the matrix exponential),

exp K =

n=0

()nn!

( K)n . (4)

In part (a), we established the following result: ( K)3 = K. Hence, it follows that

( K)2n = ( K)2 , ( K)2n+1 = K, for n = 1, 2, 3, . . . .

Thus, we can rewrite the series given in eq. (4) as

exp K

= I K

n oddn1

n

n!+ ( K)2

n evenn2

n

n!, (5)

after using the fact that ( K)

0

= I is the 4 4 identity matrix. Using,n=0

2n+1

(2n + 1)!= sinh ,

n=0

2n

(2n)!= cosh ,

and noting that the last summation in eq. (5) starts at n = 2, we end up with

exp K

= I K sinh + ( K)2 [cosh 1] , (6)

which is the desired result.

3


4/24

REMARKS:

To understand the significance of eq. (6), let us write it explicitly in matrix form. It isconvenient to use the block matrix form of eq. (3), where j labels the row and k labels the

column,

I =

1 0k

0Tj jk

, ( K)jk =

0 kj 0jk

, ( K)2jk =

1 0k

0Tj jk

.(7)

Then, eq. (6) yields

exp K

=

cosh k sinh

j sinh jk + jk(cosh

1)

.

In class, we identified = tanh1 as the rapidity, which satisfies

=1

1 2 = cosh , = sinh .

Hence, after writing = = (1 , 2 , 3), it follows that

exp K =

k

j jk + ( 1)jk

2

,

which we recognize as the boost matrix defined in eq. (11.98) of Jackson.

2. [Jackson, problem 11.13] An infinitely long straight wire of negligible cross-sectional areais at rest and has a uniform linear charge density q0 in the inertial frame K

. The frame K

(and the wire) move with velocity v parallel to the direction of the wire with respect to thelaboratory frame K.

(a) Write down the electric and magnetic fields in cylindrical coordinates in the rest frame

of the wire. Using the Lorentz transformation properties of the fields, find the components ofthe electric and magnetic fields in the laboratory.

In the rest frame of the wire (i.e. frame K), choose the z-axis to point along the wire. Then,to compute the electric field, we draw a cylinder of length L and radius r, whose symmetryaxis coincides with the z-axis. Applying Gauss law in gaussian units,

S

E n da = 4Q , (8)

4


5/24

where Q is the total charge enclosed inside the cylinder. In cylindrical coordinates (r, , z),2

the symmetry of the problem implies that E(r, , z) = E(r), where E(r) depends onlyon the radial distance from the symmetry axis. Choosing the surface S to be the surface of thecylinder, we have n = r, and so eq. (8) reduces to

2LE(r) = 4Q .

Defining the linear charge density (i.e. charge per unit length) by q0 = Q/L, we conclude that3

E(r) =2q0rr . (9)

Since there are no moving charges in the rest frame of the wire, it follows that B = 0 .The transformation laws for the electric and magnetic field between reference frames K and

K are given by4

E= E B

2

+ 1( E) ,

B = B + E

2

+ 1( B ) .

For this problem, = z. Using the results of part (a), and noting that r = r (since theradial direction is perpendicular to the direction of the velocity of frame K with respect to K),it follows that

E=2q0

rr , B =

2q0r

, (10)

where we have used z r = 0 and z r = .

(b) What are the charge and current densities associated with the wire in its rest frame? Inthe laboratory?

In reference frame K there are no moving charges, so that J = 0. The correspondingcharge density is

(r) =q0

2r(r) . (11)

To check this, let us integrate over a cylinder of length L and arbitrary nonzero radius, whose

symmetry axis coincides with the z-axis. Then,(r) dV =

(r) r dr ddz = q0

dr(r)dz = q0L = Q .

2We denote the radial coordinate of cylindrical coordinates in frame K to be r rather than the moretraditional , since we reserve the letter for charge density.

3The direction of the unit vectors r, and z are the same in frames Kand K, so no extra primed-superscriptis required on these quantities.

4Eq. (11.149) of Jackson provides the equations to transform the fields from reference frame K to reference

frame K. To transform the fields from K to K, simply change the sign of .

5


6/24

Since J = (c ; J) is a four-vector, the relevant transformation law between frames K and K

are:

c = (c + J) , (12)

J= J + 1

2( J) + c . (13)

Plugging in J = 0 and the result of eq. (11), and noting that = z and r = r, it followsthat5

(r) =q02r

(r) , J=cq0

2rz (r) = (r)vz = (r)v , (14)

after using v c.

(c) From the laboratory charge and current densities, calculate directly the electric and

magnetic fields in the laboratory. Compare with the results of part (a).

This is an electrostatics and magnetostatics problem, so we can use Gauss law to computeE and Amperes law to compute B. The computation of E is identical to the one given inpart (a) with q0 replaced by q0. Hence, it immediately follows from eq. (9) that

E(r) =2q0

rr ,

in agreement with eq. (10). Amperes law in gaussian units is

C B d =

4I

c ,

where I is the current enclosed in the loop C. With J given by eq. (14),

I =

A

J n da =

(r)v r dr d = q0v ,

after noting that n = z points along the direction of the current flow and da = rdrd isthe infinitesimal area element perpendicular to the current flow. Using the symmetry of theproblem, B = B(r). Thus, evaluating Amperes law with a contour C given by a circle

centered at r = 0 that lies in a plane perpendicular to the current flow, d = rd and

2rB(r) =4I

c=

4q0v

c,

which yields

B(r) =2q0v

r ,

after using v = c, in agreement with eq. (10).

5We can interpret q q0 as the linear charge density as observed in reference frame K. This is notunexpected due to the phenomenon of length contraction.

6


7/24

3. [Jackson, problem 11.15] In a certain reference frame, a static uniform electric field E0 isparallel to the x-axis, and a static uniform magnetic field B0 = 2E0 lies in the xy plane,making an angle with respect to the x-axis. Determine the relative velocity of a referenceframe in which the electric and magnetic fields are parallel. What are the fields in this frame

for 1 and 12?

In frame K, we haveE= E0 x , B = Bx x + By y , (15)

withE B = |E| |B| cos = E0B0 cos = 2E20 cos , (16)

after writing |E| = E0 and | B| = B0 = 2E0. It follows thatBx = 2E0 cos , By = 2E0 sin . (17)

The reference frame K

, in which the electric and magnetic fields are parallel, is moving ata velocity v c with respect to reference frame K. That is,

E B = 0 . (18)

The electric and magnetic fields in frame K are related to the corresponding fields in frame Kby eq. (11.149) of Jackson. These relations can be rewritten in the following form,

E = E , B = B , (19)

E =

E + B

, B =

B E

. (20)

In our notation above, fields with a subscript are parallel to and fields with a subscriptare perpendicular to . For example,

E = 0 and E = 0 ,which implies that

E =( E)

2and E = E (

E)

2.

The form of eqs. (19) and (20) suggests that the relative velocity v should lie in the z-direction. That is,

=

z,

in which case E = Ezz and B = Bzz. Since Ez = Bz = 0, it follows from eq. (19) thatEz = B

z = 0. Using eq. (20), the transverse fields are given by

Ex = (Ex By) = E0(1 2sin ) , Ey = (Ey + Bx) = 2 E0 cos , (21)Bx = (Bx + Ey) = 2E0 cos , B

y = (Ey + Bx) = E0(2 sin ) , (22)

after using eqs. (15)(17). Note that eq. (18) implies that

ExBy EyBx = (E B )z = 0 .

7


8/24

Inserting the above results into eqs. (21) and (22) yields,

2E20(1 2sin )(2 sin ) 42E20 cos2 = 0 .

Multiplying out the factors above and writing cos2

= 1 sin2

, the above equation simplifiesto

22 sin 5+ 2 sin = 0 .This is a quadratic equation in which is easily solved. The larger of the two roots is greaterthan 1, which we reject since 0 1 (i.e., 0 v c). The smaller of the two roots isnon-negative and less than 1. Thus, we conclude that

=v

c=

5

25 16sin2 4sin

. (23)

The two limiting cases are easily analyzed. In the case of

1, we can work to first order

in . From eq. (23) we find that 25. Since 1 it follows that 1, in which case = (1 2)1/2 1 + O(2) .

Since we are working to first order in , we also must work to first order in . In particular wecan neglect terms such as . Hence, in this limiting case, eqs. (21) and (22) yield

E = 12B = E0(x + 2y) , for 25 1 ,

where we have neglected terms that are second order (or higher) in . Finally, in the limit of 1

2, eq. (23) yields = 1

2. Then = 2/

3, and eqs. (21) and (22) yield

E = 0 , B = 3E0y , for = 12 .

REMARK:

In fact, it is easy to show that eq. (23) implies that 0 12 . If we multiply the numeratorand denominator of eq. (23) by 5 +

25 16sin2 , we obtain,

=4sin

5 +25 16sin2 .

Since the polar angle lies in the range 0 or equivalently 0 sin 1, it followsimmediately that 0 (where = 0 corresponds to sin = 0 as expected). Finally, it is easyto verify that

4sin

5 +

25 16 sin2 1

2. (24)

Since the denominator on the left hand side above is positive, we can rewrite eq. (24) as

4sin 12

5 +

25 16 sin2

. (25)

8


9/24

This inequality is manifestly true for sin = 0. For sin > 0, eq. (25) can be rearranged intothe following form

8sin 5

25 16 sin2 . (26)

Squaring both sides and simplifying the resulting expression then yields

sin (sin 1) 0 . (27)Dividing both sides of the equation by sin (which is assumed positive) yields sin 1, whichis valid for all polar angles . Hence, eq. (24) is established. The inequality becomes an equalityif sin = 1, in which case = 1

2.

4. [Jackson, problem 11.18] The electric and magnetic fields of a particle of charge q movingin a straight line with speed v = c, given by eq. (11.52) of Jackson, become more and more

concentrated as 1, as indicated in Fig. 11.9 on p. 561 of Jackson. Choose axes so thatthe charge moves along the z axis in the positive direction, passing the origin at t = 0. Letthe spatial coordinates of the observation point be (x, y, z ) and define the transverse vector r,with components x and y. Consider the fields and the source in the limit of = 1.

(a) Show that the fields can be written as

E= 2qr

r2(ct z) , B = 2qv r

r2(ct z) , (28)

where v is a unit vector in the direction of the particles velocity.

We begin with eq. (11.154) on p. 560 of Jackson,

E=qR

R32(1 2 sin2 )3/2 , (29)

where is the angle between the vectors v and R. I have modified Jacksons notation byemploying the symbol R for the vector that points from the charge q to the observation pointr = (x, y, z ) in reference frame K.6 Eq. (29) was also derived in class along with the corre-sponding result for the magnetic field,

B =

q(v R)

cR32(1 2 sin2 )3/2 . (30)The reference frame K is exhibited in Fig. 1. It is evident from this figure that

R= r vt . (31)The velocity vector is taken to lie along the z-direction. That is, v = vz.

6Jackson denotes the vector that points from the charge q to the observation point (x,y,z) by r. However,I prefer to employ r to represent the vector that points from the origin of reference frame K to the observationpoint, as shown in Fig. 1.

9


10/24

q

(x,y,z)

R

vt

r

x

z

Figure 1: A charge q moving at constant velocity v in the z-direction as seen from referenceframe K. The angle is defined so that v R= cos .

It is convenient to introduce the notation where

r = x

x+ y

y,

r = z

z,

so that r v = 0 and r v = 0. Likewise, we can resolve the vector R into componentsparallel and perpendicular to the velocity vector,

R= R + R ,

whereR Rz = (z vt)z , R = r . (32)

after making use of eq. (31). In particular, note that | R| R = R sin . It follows that

R3(1 2 sin2 )3/2 = (R2 2R2 sin2 )3/2 = (R2 + R2 2R2)3/2

= [R2 + R2(1 2)]3/2 = (R2 + R2/2)3/2 . (33)

Note that in obtaining eq. (33) we used R2 = R2 + R2 and (1 2)1/2. Moreover, since

R = r [cf. eq. (32)], we may replace R with r |r| = (x2 + y2)1/2 in the above formulae.Eqs. (29), (32) and (33) then yield

E=qr + (z vt)z

(2R2 + r

2)

3/2. (34)

Likewise, eqs.(30), (32) and (33) yield

B =q(v r)

c(2R2 + r2)

3/2. (35)

Consider the limit of 1. In this limit, , and we see that

lim

(2R2 + r2)

3/2=

0 , if R = 0 , , if R = 0 .

10


11/24

This implies that

lim

(2R2 + r2)

3/2= k(R) , (36)

for some constant k. Note that in light of eq. (32),

lim

R = z ct ,

since in the limit of v c. To determine k, we integrate eq. (36) from to ,since R can be any real number (either positive, negative or zero) depending on the value ofthe time t. Thus, employing the substitution u = R,

k =

dR(2R2 + r

2)

3/2=

du

(u2 + r2)3/2

=u

r2(u2 + r2)

1/2

=2

r2.

Hence, we conclude that

lim

(2R2 + r

2)

3/2= 2

r2(z ct) .

Inserting this result back into eqs. (34) and (35) yields

E= 2qr

r2(z ct) , B = 2qv r

r2(z ct) . (37)

In obtaining Eabove, we noted that in the limit ofv c, the z-component of the electric fieldis proportional to (z ct)(z ct) = 0 due to the properties of the delta function. In obtainingB above, we noted that limvc v/c = v. Finally, since the delta function is an even function

of its argument, we can write (z ct) = (ct z), and eq. (28) is verified.

(b) Show by substitution into the Maxwell equations that these fields are consistent withthe 4-vector source density

J = qcv(2)(r)(ct z) ,where the 4-vector v = (1 ; v).

The four-vector current is given by J = (c ; J). Hence, using the Maxwell equations ingaussian units,

E= 4 =

4J0

c .Hence, using eq. (37) and noting that Ez = 0, it follows that

J0 =c

4 E=

c

4

E+

Ezz

=

qc

2(z ct)

r

r2

, (38)

where x /x + y /y. For r x x + y y = 0, an elementary computation yields

r

r2

=

x

x

x2 + y2

+

y

y

x2 + y2

=

y2 x2(x2 + y2)2

+x2 y2

(x2 + y2)2= 0 . (39)

11


12/24

To determine the behavior at r = 0, we consider the two-dimensional analogue of the diver-gence theorem,

A dxdy

r

r2

=Crd

r

r2

r =20 d = 2 , (40)

where A is a circular disk and C is the circular boundary of the disk. Note that r = r/ris the outward normal to the circular boundary.

Eqs. (39) and (40) imply that

r

r2

= 2(2)(r) , (41)

where (2)(r) is a two-dimensional delta function. Inserting this result into eq. (38), we endup with

J0 = qc(2)(r)(z ct) . (42)Next, we employ the Maxwell equation,

B 1c

E

t=

4

cJ, (43)

to evaluate J. First, we compute

v r = z (xx + yy) = xy yx , (44)

where we have used the fact that v points in the z direction. It then follows that

v r

r2(z ct)

= det

x y z

x

y

zy

x2 + y2(z ct) x

x2 + y2(z ct) 0

= xx + yyx2 + y2

(z ct) +

x

x

x2 + y2

+

y

y

x2 + y2

(z ct)

= rr2

(z ct) + z rr2

(z ct)

= rr2

(z ct) + 2z (2)(r)(z ct) , (45)

after making use of eq. (41). The prime refers to differentiation with respect to z. Finally, wecompute

t

r

r2(z ct)

= c

z

r

r2(z ct)

= cr

r2(z ct) . (46)

12


13/24

Inserting eq. (37) into eq. (43) and using eqs. (45) and (46), we obtain

J=qc

2

v r

r2(z ct)

q

2

t r

r2(z ct)

=qc

2

r

r2(z ct) + 2z (2)(r)(z ct) + r

r2(z ct)

= qcv (2)(r)(z ct) , (47)

after using the fact that v = z. Combining eqs. (42) and (47), we can write

J = qcv(2)(r)(z ct) ,

where the four-vector v = (1 ; v).

(c) Show that the fields of part (a) are derivable from either of the following 4-vectorpotentials:

A0 = Az = 2q(ct z)ln(r) , A = 0 , (48)or

A0 = Az = 0 , A = 2q(ct z) ln(r) , (49)where is an irrelevant parameter setting the scale of the logarithm. Show that the twopotentials differ by a gauge transformation and find the corresponding gauge function .

The four-vector potential is A

= ( ;A). Given the four-vector potential, the electromag-netic fields are determined by

E= A0 1c

A

t, B = A .

Inserting the scalar and vector potentials given in eq. (48),

E= 2q

(ct z)ln(r)

+2q

cz ln(r)

t(ct z)

= 2q(z

ct)x x

+ y

y 1

2

ln(x2 + y2) + ln + 2qz ln(r) z

+1

c

t (ct z)

= 2qr

r2(z ct) ,

after using r = xx + yy and r2 = x

2 + y2. In particular, note that

z+

1

c

t

f(ct z) = 0 ,

13


14/24


15/24

Finally, we demonstrate that eqs. (48) and (49) differ by a gauge transformation. Under agauge transformation (using gaussian units),

A

A = A+ , A0

A 0 = A0

1

c

t.

Denoting A by eq. (48) and A by eq. (49), it follows that

t= 2qc(ct z)ln(r) ,

= 2q(ct z) ln(r) ,

z= 2q(ct z)ln(r) .

The solution to these equations can be determined by inspection,

(x, t) = 2q(ct z)ln(r) ,

up to an overall additive constant.

5. [Jackson, problem 11.27]

(a) A charge density of zero total charge, but with a dipole moment p , exists in reference

frame K. There is no current density in K. The frame K moves with velocity v = c in theframe K. Find the charge and current densities and J in the frame K and show that there

is a magnetic dipole moment m = 12(p ), correct to first order in . What is the electricdipole moment p in K to the same order in ?

We shall assume that in frame K, we have J = 0 and the charge density, (x ) is time-independent. Furthermore, we assume that the total charge in frame K is zero, i.e.

(x ) d3x = 0 , (50)

whereas the electric dipole moment in frame K, denoted by p , is assumed to be nonzero. Bydefinition, the electric dipole moment is given by

p =

x (x ) d3x . (51)

Frame K is assumed to move with velocity v = c with respect to frame K. We shall employeqs. (12) and (13), which we repeat here:

c = (c + J) , (52)

J= J + 1

2( J) + c . (53)

15


16/24

Setting J = 0, it follows that

(x) = (x ) , J(x) = c(x ) = c(x) . (54)

Hence, the electric dipole moment in frame K is given by

p =

x (x) d3x , (55)

and the magnetic dipole moment is defined (in gaussian units) in frame K as

m =1

2c

x J(x) d3x = 1

2

x (x) d3x = 1

2

x (x) d3x

,

after employing eq. (54) in the second step above. In light of eq. (55), it follows that

m = 12 p . (56)

In order to relate p to p , we need to compare (x) and (x ). Eq. (11.19) of Jackson gives

the four-vector x in terms of x. The inverse relation is obtained by changing the sign of ,which yields

x0 = (x0 + x

) , x = x + 1

2( x ) + x0 . (57)

Working to first order in , we can set = (1 2)1/2 1, and

x0 x0 + x , x x + x0 . (58)In the same approximation, d3x d3x and (x) (x ). Hence, using eq. (50),

p =

x (x) d3x

(x + x0)

(x ) d3x =

x (x ) d3x = p . (59)

to first order in . Using eqs. (56) and (59), we conclude that to first order in ,

m = 12 p .

(b) Instead of the charge density, but no current density, in K, consider no charge density,

but a current density J that has a magnetic dipole moment m. Find the charge and currentdensities in K and show that to first order in there is an electric dipole moment p = m

in addition to the magnetic dipole moment.

Suppose that (x ) = 0 in frame K. Then, using the continuity equation

J +

t= 0 , (60)

16


17/24

and it follows that J is a steady current; that is,

J(x ) = 0 . (61)

The magnetic dipole moment in frame K

is nonzero and is denoted by m

. By definition, themagnetic dipole moment is given by

m =1

2c

x J(x ) d3x . (62)

Then, using eqs. (52) and (53),

(x) =

c J(x ) , (63)

J(x) = J(x ) + 1

2[ J(x )] . (64)

Taking the dot product of eq. (64) with ,

J(x) = J(x ) . (65)

Inserting this result into eq. (63) yields

(x) =

c J(x) . (66)

Using eq. (66), the electric dipole moment in frame K is given by

p =x (x) d3x = 1

cx J(x)

d3x . (67)

Working to first order in , we can set = (1 2)1/2 1, in which case d3x d3x and J(x) = J(x ). Hence,

p 1c

x J(x )

d3x .

We now apply using eq. (58) which yields,

p

1

cx J

(x

)

d

3

x

+

x0

c J

(x

) d

3

x

. (68)

This result can be simplified by using the following trick. One can express the ith componentof the current as follows,

J i = k(J kx i) x i J = k(J kx i) , (69)

after applying eq. (61) in the final step. Thus,J i d3x =

k(J kx i) d3x = 0 , (70)

17


18/24

after converting to a surface integral at infinity, which vanishes under the assumption of alocalized current. Hence, the last term in eq. (68) vanishes, and we are left with

p

1

cx J(x ) d3x . (71)

We can rewrite the integrand of eq. (71) by using the vector identity

x J(x )

=

x J(x )

+ J(x )( x ) .

Inserting this result into eq. (71) and using eq. (62), it follows that

p = 2 m +1

c

J(x )( x ) d3x . (72)

The remaining integral can evaluated by employing eq. (69) for J. In particular, the ith

component of the integral in eq. (72) is

1

c

J ijxj d3x =

j

c

xj k(J kx i) d3x =

j

c

J kx i k(xj) d3x =

j

c

Jjx i d3x ,

(73)after an integration by parts, where the surface integral at infinity can be dropped as it vanishesfor a localized current. In the final step of eq. (73), we used k(xj) = jk and summed overthe repeated index k. Thus, eq. (73) yields

J(x )( x ) d3x = 1

c x

J(x d

3x p ,

where we have used eq. (71) in the final step. Inserting this result back into eq. (72) yields

p = 2 m p .Solving for p, we end up with

p = m . (74)

For completeness, we shall evaluate m in terms of m, keeping only terms to first orderin . In this approximation, we have d3x d3x and J(x) = J(x ). Using these results alongwith eq. (58), it follows that

m = 12c

x J(x) d3x = 1

2c

x J(x ) d3x = 1

2c

(x + x0) J

(x ) d3x

= m +x02c

J(x ) d3x ,

after using eq. (62). However, in light of eq. (70), the remaining integral above vanishes. Hence,to first order in ,

m = m . (75)

18


19/24

ADDENDUM to Problem 5

Although Jackson only asks for the the results to first order in , it is not too hard to derivethe exact results. This derivation will be given below.

(a) By assumption, (x ) is a time-independent charge density in frame K which satisfieseq. (50) [i.e. the total charge vanishes]. We are asked to compute p and m in frame K. Byassumption, this means that we should evaluate these quantities at a fixed time x0 in frame K.Using eq. (57), it follows that

x0 =x0

x .

Inserting this result back into eq. (57) then yields

x = x + 1

2 ( x

) + x0 .

We can simplify this expression by noting that = (1 2)1/2 implies that 22 = 2 1.Substituting for 2 above,

12

= + 1

. (76)

Hence,

x = x + 1

( x ) + x0 . (77)

We can evaluate the Jacobian of the transformation, eq. (77),7

xixj= ij + 1 ij , (78)

which we evaluate at fixed x0. To compute the determinant of the Jacobian, we use the followinggeneral result which is easily proved (see the Appendix following this addendum):

det(ij aiaj) = 1 |a|2 .

Hence,

det

xixj

= 1

2

+ 1=

1

+ 1

+ 1

2 1

=

1

.

It follows thatd3x = det

xixj

d3x =

d3x

, (79)

which is just the well-known length contraction in special relativity. Using eqs. (54) and (79),it follows that the charge dq located inside an infinitesimal volume element is

dq = (x) d3x = (x ) d3x , (80)

7In eq. (78), we do not distinguish between lowered and raised indices, as all involved quantities are three-dimensional.

19


20/24

since the factors of cancel out. Eq. (80) is expected since the electric charge is a Lorentzscalar, which must be independent of the reference frame used to evaluate it.

We now can compute p in frame K. Using eqs. (77) and (80),

p =x (x) d3x =

x (x ) d3x =

x

+ 1( x ) + x0

(x ) d3x . (81)

Using eq. (51), it follows that

p = p + 1

( p ) (82)

where we have used eq. (50) [i.e., the total charge in frame K vanishes] to eliminate the lastterm in eq. (81). Indeed, if we work to first order in , then the second term on the right handside of eq. (82) can be dropped and we recover the result of eq. (59).

To obtain m, we start with the exact result obtained in eq. (56). Inserting eq. (82) for p

and using = 0, it follows that8

m = 12p , (83)

which is exact to all orders in .

REMARK:

One can derive eq. (82) more directly as follows. Using eqs. (54) and (57),

p =

x (x) d3x =

x +

12

( x ) + x0

(x ) d3x

=

x + 12 ( x ) + x0

(x ) (x0 ct)d4x .

In the last step, I inserted the integral

(x0 ct)dx0 = 1 which does not change the result.This is useful, since d4x = d4x = d3xdx0. Using eq. (57) to express x0 in terms of x

0,

p =

x +

12

( x ) + x0

(x )

x0 + x

ct

d3xdx0 ,

after writing (x0 ct) = ([x0 ct]/) = (x0 + x ct/). Integrating over x0 yields

p = x + 12

( x ) (x ) d3x + ct (x ) d3x .Using eqs. (50) and (76) along with p =

x (x ) d3x, we end up with

p = p + 1

( p ) ,

which reproduces eq. (82). Of course, the two derivations are equivalent.

8In the literature, one often finds the result of eq. (83) quoted without the factor of 12

. For a discussionof this discrepancy, see V. Hnizdo, Magnetic dipole moment of a moving electric dipole, American Journal ofPhysics 80, 645647 (2012).

20


21/24

(b) By assumption, J(x ) is a time-independent, steady current in frame K which satisfieseq. (61). We are asked to compute p and m in frame K. In this computation, it will be quiteuseful to derive a formula for the following integral,

x iJj d3x .

This can be accomplished using a identity similar to that of eq. (69). First we write

x iJj = 12(x iJj + xjJ i) + 12(x

iJj xjJ i) .We then observe that

x iJj + xjJ i = k(x ixjJ k) + x ixj J = k(x ixjJ k) ,after applying eq. (61), and

x iJj xjJ i = ijk(x J)k .Consequently,

x iJj d3x = 12

(x iJj + xjJ i) d3x + 12

(x iJj xjJ i) d3x

= 12

k(x ixjJ k) d3x + 12ijk

(x J)k d3x . (84)

The first term in the second line of eq. (84) can be converted into a surface integral at infinity,

which vanishes under the assumption that the current is localized. The second term is relatedto the magnetic moment via eq. (62). Hence,

1

c

x iJj d3x = ijkm k . (85)

We are now ready to compute m. We again perform the computation at fixed x0, in whichcase we can use eqs. (77) and (79). Employing eq. (64),

m =1

2c

x J(x) d3x

=1

2c

x

+ 1 ( x

) + x0J

(x

) +

1

2 [ J

(x

)]

d3

x

=m

+

1

2c

1

2[ J(x )]x

+ 1( x ) J(x ) + x0 J

(x )

d3x ,

where some of the cross-terms vanish due to = 0. It is convenient to rewrite the aboveexpression in component form,

mk =m k

+

1

2cijk

12

x iJ j + 1

x Jji + x0Jji

d3x .

21


22/24

Using eqs. (70) and (85), the integrals above are easily evaluated,

mk =m k

+

1

2ijk

12

inmnj + 1

jnmni

.

Summing over the product of -tensors,

mk =m k

+

1

2mn

12

(jkn jnk)j + 1

(ikn ink)i

=m k

+

1

2mn(2kn kn)

1

2

+ 1

=m k

+

12(+ 1)

2m k ( m )k

,

after using 2 = (2

1)/2. We can also use the latter to combine the first two terms above,

1

+

2( 1)2(+ 1)

=1

2

1 +

1

2

.

Thus, we end up with

m =1

2

1 +

1

2

m 1

2(+ 1)( m ) .

To first order in , we can approximate 1, and it follows that m m , which confirms theresult of eq. (75).

Next, we evaluate p at fixed x0 starting from eq. (67),

p =1

c

x J(x)

d3x =

1

c

x J(x )

d3x ,

after using eqs. (65) and (79). We now employ eq. (77) and obtain

p =1

c

x

+ 1 ( x ) + x0

J(x )

d3x ,

which we rewrite in component form,

pi =1

c x i

+ 1ikx k + x0

i

jJj d3x .

Using eqs. (70) and (85), the integrals above are easily evaluated,

pi = j

ijm + 1

ikkjm

= ( m )k ,

after noting that kjjk = 0. Hence,

p = m ,

which is exact to all orders in .

22


23/24

REMARK:

The results of parts (a) and (b) can be combined by linear superposition into the followingresult. Suppose a time-independent localized charge density (x ) of zero total charge and a

time-independent, steady localized current density J

(x

) exist in a reference frame K

. In thisreference frame, the system of charges and currents possesses an electric dipole moment p anda magnetic dipole moment m . The frame K moves with velocity v = c in the frame K.Then, in frame K there is a charge density and current density given by:

c(x) =

c(x ) + J(x )

, J(x) = J(x ) + 1

2 J(x )

+ c(x ) ,

where = (1 2)1/2. Moreover, the electric and magnetic dipole moments in frame K are:

p = p + 1

( p ) + m ,

m =1

2

1 +

1

2

m 1

2(+ 1)( m ) + 12 p

.

The above results are exact to all orders in . If one sets m = 0 [or p = 0], then one recoversthe results of part (a) [part (b)] above, respectively.

An equivalent result can be found by defining the unit vector /. Then, using2 = (2 1)/2, we end up with9

p = p 1 1

( p) + m , (86)

m =1

2

1 +

1

2

m

1 1

2( m ) + p

. (87)

9Eqs. (86) and (87) were obtained in George P. Fisher, The electric dipole moment of a moving magneticdipole, American Journal of Physics 39, 15281533 (1971). See also Marijan Ribaric and Luka Sustersic, Movingpointlike charges and electric and magnetic dipoles, American Journal of Physics 60, 513519 (1992).

23


24/24

APPENDIX: Proof of a determinantal formula

Theorem: Let A be an n n matrix, whose matrix elements are given by:

Aij = ij aiaj . (88)Then,

det A = 1 |a|2 , where |a|2 =ni=1

a2i . (89)

Proof: The determinant of A is defined as

det A = i1i2in A1i1A2i2 Anin ,

where there is an implicit sum over repeated indices. Plugging in eq. (88),

det A = i1i2iN(1i1 a1ai1)(2i2 a2ai2) (nin anain) .

Expanding the product above, and using the Kronecker deltas to perform the sums, we obtain

det A = 123n a1ai1i123n a2ai21i23n anain123in= 123n(1 a21 a22 a2n) , (90)

after performing the final set of summations and using the fact that i1i2in vanishes unless allof its indices are distinct. Note the absence of any terms in eq. (90) that are quartic (or higher

order) in the ai, since all such terms will be symmetric under the interchange of two indicesthat are summed against two corresponding indices of the Levi-Civita tensor. For example,

i1i2iNa1ai1a2ai2 = 0 ,

since ai1ai2 is symmetric under the interchange of i1 and i2 whereas i1i2iN is antisymmetricunder this interchange of indices.

Since 123n = 1, eq. (90) yields

det A = 1 n

i=1 a2i ,

which completes the proof.

emii2sol_13.pdf

Documents