emii4sol_13.pdf

7/30/2019 emii4sol_13.pdf

1/22

Physics 214 Solution Set 4 Winter 2013

1. [Jackson, problem 12.3] A particle with mass m and charge e moves in a uniform, static,electric field E0.

(a) Solve for the velocity and position of the particle as explicit functions of time,assuming that the initial velocity v0 was perpendicular to the electric field.

Using eqs. (12.1) and (12.2) of Jackson and setting B = 0, we have:

dp

dt= eE,

dW

dt= ev E,

where W is the total mechanical energy (usually called E, but we have renamed this W in

order to better distinguish it from the electric field) and v is the particle velocity (which isdenoted as u by Jackson).

Clearly, the motion takes place in a plane containing the E-field. Without loss ofgenerality, we assume that

E= Ex ,

and assume that the motion takes place in the xy plane. By assumption, v E = 0 att = 0, in which case px = 0 at t = 0. Solving the equations,

dpxdt

= eE ,dpydt

= 0 ,

in follows thatpx = eEt , py = p0 ,

where p0 is a constant.Using p = mv and E = mc2, it follows that

v =c2p

W=

c2p|p|2c2 + m2c4 .Hence,

vx =

c2eEt(p20 + e2E2t2)c2 + m2c4 , vy =

c2p0(p20 + e2E2t2)c2 + m2c4 .

Since v = dx/dt, it follows that

x = c2eE

tdt

W20 + (ceEt)2

, y = c2p0

dt

W20 + (ceEt)2

, (1)

where W20 = p20c

2 + m2c4.

1


2/22

We shall define the origin of the coordinate system to coincide with t = 0. Thencomputing the integrals in eq. (1) yields

x(t) =1

eEW20 + (ceEt)2 W0 , y(t) =p0c

eEsinh1

ceEt

W0 . (2)

(b) Eliminate the time to obtain the trajectory of the particle in space. Discuss theshape of the path for short and long times (define short and long times).

We can eliminate t from eq. (2),

t =W0ceE

sinh

eEy

p0c

.

Inserting this into the equation for x(t) and using the identity cosh2 z

sinh2 z = 1, it

follows thatx =

W0eE

cosh

eEy

p0c

1

,

which is the equation for a catenary curve.To describe the shape of the path for short and long times, we note that W0/(ceE)

has units of time. This we can define short and long times relative to this quantity. Fort W0/(ceE), we have

W20 + (ceEt)2 W0 + (ceEt)

2

2W0, sinh1

ceEt

W0

ceEt

W0.

Hence the approximate form of eq. (2) is

x(t) c2eEt2

2W0, y(t) p0c

2t

W0.

Solving for t and inserting the result back into the above equations yields

x eEW0y2

2p20c2

.

Since v0 = c2p0/W0, we can eliminate W0 from the above expression to obtain,

x

eEy2

2p0v0. (3)

That is, as short times, the motion is parabolic.1

1The result of eq. (3) also coincides with the non-relativistic limit (in which case p0 = mv0). To verifythis assertion, we can perform a formal expansion in powers of 1/c. In this limit, W0 mc2 and

t W0ceE

mceE

,

which is always true in the limit ofc (which is equivalent to taking the non-relativistic limit).

2


3/22

For t W0/(ceE), eq. (2) yields:

x(t) ct , y(t) p0ceE

ln

2ceEt

W0

.

In the latter case, we used:

sinh1 z = ln

z+

z2 + 1

ln 2z , for z 1 .

Hence, to a good approximation,

y p0ceE

ln

2eEx

W0

,

or equivalently,

x W02eE

exp

eEy

p0c .

That is, at long times the motion is exponential.

2. [Jackson, problem 12.9] The magnetic field of the earth can be represented approximatelyby a magnetic dipole of magnetic moment M = 8.1 1025 gauss-cm3. Consider the motionof energetic electrons in the neighborhood of the earth under the action of this dipole field(Van Allen electron belts). [Note that M points south.]

(a) Show that the equation for a line of magnetic force is r = r0 sin2 , where is the

usual polar angle (colatitude) measured from the axis of the dipole, and find an expression

for the magnitude of B along any line of force as a function of .

Let the z-axis point from the origin in the direction of the north pole. Then, themagnetic dipole moment (which points south) is given by M = Mz, where M | M|.The vector potential is given in gaussian units by:

A(x) =M x

|x|3 =1

r3det

x y z0 0 M

r sin cos r sin sin r cos

=Msin

r

2(x sin

y cos ) =

Msin

r

2 ,

where r |x|. Then,

B = A =1

r2 sin det

r r r sin

r

0 0 r sin A

,

3


4/22

where A = Msin /r2. Evaluating the above determinant yields:

B = 2Mr3

cos r Msin r3

. (4)

Given the magnetic field at every point in space, B(x), one can consider a relatedvector field, F(x) = qm B(x), which gives the force on a magnetic test charge qm due tothe magnetic field at the point x.2 If we choose our test charge to have qm = 1, then thereis no distinction between the lines of magnetic force and the magnetic field lines. Wechoose to follow this convention in what follows.

The lines of force follow a curve x(x), where the arclength s parameterizes the location

along the curve. By definition B(x) is tangent to the lines of force. That is,

dx

ds=

Bx(s)

B

, (5)

where B | B|. To understand the normalization on the right hand side above, we notethat eq. (5) is equivalent to the three equations,

dx

ds=

BxB

,dy

ds=

ByB

,dz

ds=

BzB

.

Squaring each equation and summing the three resulting equations yields

(ds)2 = (dx)2 + (dy)2 + (dz)2 ,

which is the well-known formula for the differential arclength.

It is convenience to work in spherical coordinates. Consider an infinitesimal displace-ment dx, where

x = r sin cos x+ r sin sin y + r cos z .

By the chain rule,

dx =x

drdr +

x

dd +

x

dd

= (cos sin x + sin sin y + cos z) dr + r (cos cos x+ cos sin y sin z) d+r ( sin sin x + sin cos y) d

= r dr + r d + r sin d . (6)

2Of course, magnetic charges do not exist in classical electromagnetism. But the concept of lines offorce were developed before this fact was understood. In the case of the electric field, we do have F = qE,so the terminology lines of force makes sense. In the case of magnetic fields, it would be better to refer tothe lines of force as the magnetic field lines. Nevertheless, following Jackson, we retain the old terminologyin this problem.

4


5/22

The tangent to the curve x(s) then takes the form

dx

ds= r

dr

ds+ r

d

ds+ r sin

d

ds. (7)

Using eq. (4), it follows that the line of magnetic force is determined by the equation,

dx

ds=B (x(s))

B= 2M

Br3cos r Msin

Br3 , (8)

where r, and are functions of s. Equating eqs. (7) and (8) yields three differentialequations,

dr

ds= 2Mcos

Br3, r

d

ds= Msin

Br3,

d

ds= 0 . (9)

Dividing the first two equations above yields,

dr

d=

2r cos

sin ,

which is easily integrated, dr

r= 2

cos

sin d ,

Evaluating the integrals and imposing the condition r = r0 at =1

2, we obtain

ln

r

r0

= 2lnsin ,

or equivalentlyr = r0 sin

2 , (10)

which we identify as the equation for the line of magnetic force. Note that the third

equation in eq. (9) implies that is a constant along the line of magnetic force.Finally, we evaluate the magnitude of B along the line of force. Since

B | B| =

B2r + B2 + B

2 =

M

r3

4cos2 + sin2 ,

We simply plug in eq. (10) to obtain B as a function of along the line of magnetic force,

B() =M

r30

1 + 3 cos2

sin6 , (11)

after using sin2 = 1 cos2 in the numerator above.

(b) A positively charged particle circles around a line of force in the equatorial planewith a gyration radius a and a mean radius R (where a R). Show that the particlesazimuthal position (east longitude) changes approximately linearly in time according to:

(t) = 0 32

aR

2B(t t0) .

where B is the frequency of gyration at radius R.

5


6/22

Assuming that a R, we can use eq. (12.55) of Jackson to obtain an approximateformula for the gradient drift velocity,

vGBa

=a

2B2 B B , (12)where a is the gyration radius and B is the field at the equator ( = 1

2). Using eq. (4),

this means that

B = Mr3

r=R

= MR3

, B = rB

r

r=R

= 3MR4

r , (13)

where R is the mean radius. In computing B, we used the fact that B | B| = M/r3and

= n = r

r+

1

r sin

,

where n B = 0. Inserting the results of eq. (13) into eq. (12), we end up with

vG =Ba

2

2

R6

M2

M

R3

3M

R4

r = 3Ba

2

2R . (14)

Finally, we can express vG in terms of the angular velocity d/dt by

vG = Rd

dt .

Comparing this equation with eq. (14), we conclude that

d

dt=

3a2

2R2B

.

Solving this differential equation, and imposing the initial condition (t0) = 0 , we end upwith

(t) = 0 3a2

2R2B(t t0) . (15)

(c) If, in addition to its circular motion of part (b), the particle has a small componentof velocity parallel to the lines of force, show that it undergoes small oscillations in around = 1

2 with frequency = (3/

2)(a/R)B. Find the change in longitude per cycle of

oscillation in latitude.

As discussed in Chapter 12, section 4 of Jackson, the transverse velocity of gyrationis v = Ba [cf. discussion below eq. (12.61) of Jackson]. If in addition, we now includethe small component of the velocity parallel to the lines of magnetic force, we may useeq. (12.72) of Jackson to write:

v2 = v20 v2 0

B(z)

B0.

6


7/22

Here, the subscript 0 refers to the equator z = 0 (or equivalently to = 12

). In particular,we can write v20 = v

2

0 + v2 0 so that

v2 = v2

0 + v2

01 B(z)

B0 . (16)

In part (a), we found that along the lines of magnetic force,

B() =M

r30

1 + 3 cos2

sin6 , (17)

where r0 r( = 12). In this problem, we are interested in the behavior of the particle atthe mean radius R, so we take r0 = R. To compute B(z), we expand about z = 0. Sincez = R cos , we expand about z = 0 by writing = 1

2 + . Then,

z = R cos = R cos 1

2 + = R sin R .

Hence, z/R and 12 z/R. It follows that

cos cos

2 z

R

= sin

z

R, sin sin

2

zR

= cos

z

R.

Using eq. (17),

B(z) Mr30

1 + 3 sin2(z/R)

cos6(z/R) M

r30

1 + 3z2/R2

[1 z2/(2R2)]6 M

r30

1 +

9z2

2R2

.

Plugging this result into eq. (16) yields

v2(z) = v2 0 9

2

Ba

R2

z2 . (18)

As discussed below eq. (12.72) of Jackson, this equation is equivalent to the conservationof energy of a one-dimensional non-relativistic mechanics problem with total mechanicalenergy,

E(z) = 12

mv2 + V(z) ,

where

V(z) = 12

m

92

Ba2

2R2

z2 , (19)

is the potential energy of a one-dimensional harmonic oscillator. Indeed, eq. (18) is equiv-

alent to the statement that E(z) = E(0), i.e. conservation of energy. If we write theharmonic oscillator potential in the standard form,

V(z) = 12

m2z2 ,

the eq. (19) implies that the effective oscillator frequency is given by

=3

2

Ba

R.

7


8/22

That is, the charged particle undergoes small oscillations in around = 12

with fre-quency .

One period T of oscillation is given by

T = 2 = 2

2R3Ba

. (20)

Using the results of part (b) [cf. eq. (15)], the change of longitude is

= 3a2

2R2Bt . (21)

Choosing t = T then yields the change of longitude per cycle of oscillation in latitude,

=

2a

R.

(d) For an electron of 10 MeV kinetic energy at a mean radius of R = 3 107 m, find and a, and so determine how long it takes to drift once around the earth and how longit takes to execute one cycle of oscillation in latitude. Calculate the same quantities for anelectron of 10 keV at the same radius.

Given M = 8.11025 gauss-cm3 and R = 3109 cm, the magnetic field at the equator is

B =M

R3= 3 103 gauss .

Using eq. (12.39) of Jackson,

B =eB

mc=

ecB

mc2. (22)

Although the last step above is rather trivial, it is convenient to write B in this form. Thenumerical value of the quantity ec is given by

ec = (4.8 1010 statcoulombs)(3 1010 cm s1) = 14.4 statcoulombs cm s1 . (23)

It is convenient to eliminate statcoulombs in favor of gauss. That is,

1 gauss = 1 dyne statcoulomb1 = 1 erg cm1 statcoulomb1 .

Using 1 eV = 1.6 1012 ergs, we can write:1 gauss = (1.6 1012)1 eV cm1 statcoulomb1 = 6.25 1011 eV cm1 statcoulomb1 .

Hence, it follows that

1 statcoulomb = 6.25 1011 eV cm1 gauss1 .

8


9/22

Inserting this result into eq. (23) yields

ec = 9 1012 eV gauss1 s1 .Therefore, the gyration frequency can be written as

B = 9 1012 s1B (gauss)

mc2 (eV). (24)

For the electron, we have mc2 = 511 keV. If the electron has a kinetic energy of K =10 MeV, then E = mc2 = mc2 + K, which yields K = ( 1)mc2. Hence,

= 1 +K

mc2= 1 +

10 MeV

0.511 MeV= 20.57 .

It follows from eq. (24) that

B = 9 1012

s1

3

103

(20.57)(5.11 105) = 2.57 103

s1

.

Next we use v v = Ba to determine a. Since 1, it follows that v c, so that

a =c

B=

3 1010 cm s12.57 103 s1 = 117 km .

To drift once around the earth requires the longitude (or azimuthal angle ) to change by2. Inserting = 2 in eq. (21) [the overall sign is not significant here], we obtain

t =4R2

3a2B

=4(3 109 cm)2

3(1.17

107 cm)2(2.57

103 s1)= 107 s .

Finally, the time it takes to execute one cycle of oscillation in latitude was obtained inpart (c) [cf. eq. (20)]:

T =2

2R

3Ba=

2

2(3 109 cm)3(1.17 107 cm)2(2.57 103 s1) = 0.3 s .

For an electron with kinetic energy of 10 keV,

= 1 +K

mc2= 1 +

10 keV

511 keV= 1.02 . (25)

It follows from eq. (24) that

B =(9 1012 s1)(3 103)

(1.02)(5.11 105) = 5.18 104 s1 .

To determine a, we first compute v using eq. (25):

11 v2/c2 = 1.02 =

v

c= 0.195 .

9


10/22

Hence,

a =v

B=

(0.195)(3 1010 cm s1)5.18 104 s1 = 1.13 km .

Finally, following the previous computation,

t =4R2

3a2B=

4(3 109 cm)23(1.13 105 cm)2(5.18 104 s1) = 5.7 10

4 s ,

and

T =2

2R

3Ba=

2

2(3 109 cm)3(1.13 105 cm)2(5.18 104 s1) = 1.52 s .

Note that in both computations above, we have a R, which implies that the gradientof the magnetic field is small over the orbit of the electrons. Hence, the approximationsintroduced in Chapter 12, sections 4 and 5 of Jackson are valid for the charged particlemotions examined in this problem.

3. [Jackson, problem 12.11] Consider the precession of the spin of a muon, initially lon-gitudinally polarized, as the muon moves in a circular orbit in a plane perpendicular to auniform magnetic field B.

(a) Show that the difference of the spin precession frequency and the orbital gyrationfrequency is

=eBa

mc,

independent of the muons energy, where a = 12

(g 2) is the magnetic moment anomaly.Find the equations of motion for the components of the spin along the mutually perpen-dicular directions defined by the particles velocity, the radius vector from the center of thecircle to the particle, and the magnetic field.

Our starting point is the Thomas equation, which Jackson writes in the following form[cf. eq. (11.170) of Jackson]:

ds

dt=

e

mcs

g

2 1 + 1

B

g2

1

+ 1( B)

g

2

+ 1

E

,

(26)where the time derivative of the velocity vector is given by [cf. eq. (11.168) of Jackson]:

ddt

= emc

E+ B ( E) . (27)

For a particle moving in a circular orbit in a plane perpendicular to a uniform magneticfield B, we have B = 0, where v c is the particle velocity. Hence, eqs. (26) and(27) reduce to

ds

dt=

e

mc

g

2 1 + 1

s B ,

dv

dt=

e

mcv B , (28)

10


11/22

since by assumption there is no electric field present ( E = 0). That is, eq. (28) can bewritten in the form of precession equations,

ds

dt= s ,

dv

dt= v B ,

where the spin precession frequency and the orbital gyration frequency B are given by:

emc

1 +

g 2

2

B , B

e

mcB .

The difference of these two frequencies is

B =e

mc

g 2

2

B ,

and the magnitude of this frequency difference is given by

=eBa

mc, where a = 1

2(g 2) .

To find the equations of motion for the components of the spin vector, we first decomposethis vector into longitudinal and transverse components with respect to the direction of thevelocity, /. That is, s = s + s, where

s = ( s) , s = s s .By construction,

s = 0 . (29)

We first work out ds/dt.

dsdt

=d

dt

( s)

=

d

dt

s

+ s

d

dt. (30)

Jackson gives the following result in his eq. (11.171),

d

dt

s

= e

mcs

g2

1 B +

g

2 1

E

.

Setting E= 0, we obtain

d

dt

s

= eB

mc

g 2

2

s ( B) . (31)

We also need to work out d/dt.

d

dt=

d

dt

=

1

d

dt

2d

dt. (32)

11


12/22

Using

d

dt=

d

dt

1/2

=1

2

1/2 d

dt

=

1

22

d

dt=

d

dt,

in eq. (32), we conclude that

d

dt=

1

d

dt

d

dt

.

From eq. (28), we obtain

d

dt=

e

mc B .

Hence d/dt = 0, and we end up with

ddt

= eBmc

B . (33)

Inserting eqs. (31) and (33) into eq. (30), we obtain

dsdt

= eBmc

g 2

2

[s ( B)] +

eB

mcs ( B) .

Since s (s ), it immediately follows that

s ( B) = sB .

We can further simplify the quantity [s ( B)] by using s = 0 [cf. eq. (29)] and

B = 0. First, consider the triple cross product

s ( B)

= [s ( B)] ( B)s = [s ( B)] .

However, ( B) = ( B) B = B. Hence,[s ( B)] = s B .

Inserting eqs. (35) and (36) into eq. (34) then yields

dsdt

=eB

mc

g 2

2

s +

1

s

B

Using this result, we can evaluate ds/dt.

dsdt

=d

dt

s s

=

eB

mc

g

2 1 + 1

(s + s)

B =eB

mc

g 2

2

s +

1

s

B ,

12


13/22

which simplifies to

dsdt

=eB

mc

g 2

2

s +

1

s

B

Finally, we need to further decompose s

into components along the direction of themagnetic field and along the radius vector r which points to the center of the circular path ofthe moving spin. In light of eq. (27) [with E= 0], dv/dt B. But for circular motion,r = 0 and the acceleration dv/dt points radially into the origin, i.e. dv/dt r. Itfollows that r = B , and we conclude that the unit vectors {B , , r} form a mutuallyorthogonal right-handed triad of vectors. Thus, we can write:

s sB + sr , where sB (s B)B and sr (s r)r . (34)Note that

dsBdt

=

B

ds

dt

B = 0 , (35)

since B is time-independent by assumption and

B ds

dt B (s B) = 0 ,

in light of eq. (28). Thus, sB is a constant in time, from which it follows that

dsrdt

=d

dt(s + sB) =

dsdt

. (36)

Hence, the equations of motion for the components of the spin vector are:

dsB

dt = 0 ,

dsrdt

=eB

mc

g 2

2

s +

1

sr

B ,

dsdt

=eB

mc

g 2

2

sr +

1

s

B ,

after using sB B = (s B)B B = 0.

(b) For the CERN Muon Storage Ring, the orbit radius is R = 2.5 meters and B =

17 103 gauss. What is the momentum of the muon? What is the time dilation factor ?How many periods of precession T = 2/ occur per observed laboratory mean lifetimeof the muons? [Relevant data: m = 105.66 MeV, 0 = 2.2 106 s, a /(2) where 1/137.]

For circular motion,

a =dv

dt= v

2

Rr . (37)

13


14/22

Since the circular motion is in a plane that is perpendicular to the magnetic field B, itfollows that B, v and r are mutually orthogonal vectors. Moreover, eqs. (12.38) and (12.39)of Jackson yield

dv

dt

=e

mc

v B . (38)

Thus, if B points in the z-direction, then v = v and the circular motion is clockwise inthe xy plane. Combining eqs. (37) and (38), it follows that

mv =eBR

c, (39)

which we recognize as the relativistic momentum of the muon, p. Using eq. (12.42) ofJackson, we can rewrite eq. (39) as3

p (MeV/c) = 3 104 BR (gauss-cm) .

Hence,p = (3 104)(1.7 104)(250) MeV/c = 1.275 103 MeV/c .

The -factor is

=E

mc2=

(p2c2 + m2c4)1/2

mc2=

p2

m2c2+ 1

1/2.

The muon rest energy is mc2 = 105.66 MeV. Hence,

=

1 +

(1.275 103)2(105.66)2

1/2= 12.11 .

The number of periods of precession, T = 2/, occurring per observed mean muonlifetime, 0 = (2.2 106 s), is given by4

0T

=0

2=

0eBa

2mc=

20va

2R,

where eq. (39) was used to arrive at the final result above. Since 1, we can approximatev c. In addition, we take

a = 12

(g 2) 2

, where 1137

,

as predicted at lowest non-trivial order in quantum electrodynamics. Hence,

0T

20c

42R=

(12.11)2(2.2 106 s)(3 1010 cm s1)42(250 cm)(137)

= 7.156 .

3The factor of 3 104 arises as follows. In gaussian units, e = 4.8 1010 esu and 1 MeV= 1.6 106 ergs. Hence, the conversion factor between ergs and MeV is 4.8 1010/1.6 106 = 3 104.

4Note that in the laboratory frame, the observed muon lifetime is given by 0, where 0 is the muonlifetime in the muon rest frame.

14


15/22

(c) Express the difference frequency in units of orbital rotation frequency and com-pute how many precessional periods (at the difference frequency) occur per rotation for a300 MeV muon, a 300 MeV electron, a 5 GeV electron (this last typical of the e+e storagering at Cornell).

NOTE: The energy values above correspond to the total relativistic energies.

For a 300 MeV muon,

=E

mc2=

300

105.66= 2.839 ,

and

=eBa

mc= Ba

B

2= 3.3 103B .

One revolution occurs in time t = 2R/v. In this time, the number of periods of precession,T = 2/, is given by

t

T =2R

v

2

=R

v .

We can rewrite the above result using eq. (39), which yields

R

v=

mc

eB=

1

B.

Hence, for a 300 MeV muon, we have

t

T=

B

2= 3.3 103 .

For a 300 MeV electron, we use mec2 = 511 keV to obtain

=300

0.511= 587 .

Hence,t

T=

B

2= 0.682 .

Finally, for a 5 GeV electron, we have

=5000

0.511= 9.785

103 .

It follows thatt

T=

B

2= 11.37 .

15


16/22

4. [Jackson, problem 14.4] Using the Lienard-Wiechert fields, discuss the time-averagedpower radiated per unit solid angle in nonrelativistic motion of a particle with charge e,moving:

(a) along the z axis with instantaneous position z(t) = a cos 0(t) ,

(b) in a circle of radius R in the xy plane with constant angular frequency 0.

Sketch the angular distribution of the radiation of the radiation and determine the totalpower radiated in each case.

(a) Case 1: Non-relativistic motion of a particle with charge e moving along the z axiswith instantaneous position z(t) = a cos 0(t) .

We make use of eq. (14.20) of Jackson, which is relevant for non-relativistic motion,

dP

d =

e2

4cnn d

dt

2

, (40)

where =

v

c=

1

c

dx

dt.

In this case, we havex(t) = z a cos 0t ,

which yieldsd

dt= z a

20

ccos 0t .

Working out the absolute square of the triple product in eq. (40),nn

d

dt

2

=

nn

d

dt

d

dt

2

=

d

dt

2

n

d

dt

2

(41)

=a240

c2cos2 0t

1 (n z)2 = a240

c2cos2 0t sin

2 .

In obtaining the final result above, we chose to work in a coordinate system in which theorigin corresponds to the instantaneous position of the charged particle, and the unit vectorn has polar angle and azimuthal angle with respect to the z-axis,

n = x sin cos + y sin sin + z cos . (42)

The time-averaged power is easily obtained by noting that5

cos2 0t = 12 .5To compute the time-average of cos2 0t, note that the time averages satisfy cos2 0t = sin2 0t,

and cos2 0t + sin2 0t = 1.

16


17/22

Hence, it follows that dP

d

=

e2a2040

8c3sin2 . (43)

In Figure 1, the angular distribution of the radiated power is exhibited as a polar plot.

Figure 1: A polar plot of the angular distribution of the power radiated by a charged particlemoving non-relativistically along the z axis with instantaneous position z(t) = a cos0(t). Theangular distribution is given by eq. (43) and is proportional to sin2 . This plot was created withMaple 15 software.

Integrating over the solid angle yields the total radiated power,

P = e2

a2

40

3c3.

(b) Case 2: Non-relativistic motion of a particle with charge e moving in a circle ofradius R in the xy plane with constant angular frequency 0.

For circular motion in the xy plane, the trajectory of the particle is given by

x(t) = R(x cos 0t + y sin 0t) .

Then, we easily compute

ddt

=1c

d2xdt2

= 20cx(t) .

We again choose to work in a coordinate system in which the origin corresponds to theinstantaneous position of the charged particle, and the unit vector n given by eq. (42) haspolar angle and azimuthal angle with respect to the z-axis. Consequently,

n d

dt=

20R

c(cos 0t sin cos + sin 0t sin sin ) .

17


18/22

Evaluating the absolute square of the triple cross product as in part (a) [cf. eq. (41)], weobtain:

nn d

dt

2

=40R

2

c2 1 sin

2 (cos cos 0t + sin sin 0t)2

=40R

2

c2

1 sin2 cos2(0t )

.

Using eq. (40), it follows that

dP

d=

e240R2

4c3

1 sin2 cos2(0t )

.

The time-averaged power is easily obtained by noting that cos2(0t ) = 12 . Employingthe trigonometric identity, 1 1

2sin2 = 1

2(1 + cos2 ) , it follows that

dP

d

=

e240R2

8c3

1 + cos2

. (44)

In Figure 2, the angular distribution of the radiated power is exhibited as a polar plot.

Figure 2: A polar plot of the angular distribution of the power radiated by a charged particlemoving non-relativistically in a circle of radius R in the xy plane with constant angular fre-

quency 0. The angular distribution is given by eq. (44) and is proportional to 1 + cos2 . Thisplot was created with Maple 15 software.

Integrating over solid angles yields the total radiated power,

P = 2e240R

2

3c3.

18


19/22

5. [Jackson, problem 14.5] A nonrelativistic particle of charge ze, mass m and kineticenergy E makes a head-on collision with a fixed central force field of finite range. Theinteraction is repulsive and described by a potential V(r), which becomes greater than Eat close distances.

(a) Show that the total energy radiated is given by

W =4

3

z2e2

m2c3

m

2

rmin

dVdr2

drV(rmin) V(r)

,

where rmin is the closest distance of approach in the collision.

Consider a particle with kinetic energy E at t = that is initially an infinite distanceaway and is headed in a radial direction toward the origin. Because the potential V(r) isrepulsive and becomes greater than E at close distances, there is a distance of closestapproach, rmin, where the particles radial velocity drops to zero. At this point, the particlecollides head on with the central force field and is turned around. It now travels backalong its original radial path until it reaches its original starting point (an infinite distanceaway) at t = .

If the particle does not radiate, then we can use energy conservation to compute theinstantaneous velocity of the particle at all points along its trajectory. In particular, theconservation of the sum of the kinetic and potential energy yields

E = 12

m[v(r)]2 + V(r) , for rmin r < , (45)since the potential is assumed to be of finite range which means that limr V(r) = 0. Atthe point of closes approach to the origin, v(rmin) = 0. Hence, it follows from eq. (45) that

E = V(rmin) . (46)

Writing v(r) = dr/dt, we can solve for the velocity using eq. (45),

v(t) =dr

dt=

2

m

E V(r) , (47)

where we employ the minus sign as the particle moves toward the origin and the plus sign asthe particle moves away from the origin. The acceleration can be obtained by differentiatingeq. (47) with respect to t. However, a more direct computation uses Newtons second law,

F = ma = V(r) = rdVdr

,

which yields

a = r 1m

dV

dr. (48)

During the period of acceleration, the particle radiates and hence loses energy, W.Thus, it is not justified to ignore this energy loss in the energy balance equation given in

19


20/22

eq. (45). Nevertheless, if W E, then it is justified in first approximation to ignore theradiated energy loss in deriving the acceleration given in eq. (48).6 Thus, we shall assumethat one can neglect the energy loss due to radiation and check for consistency at the endof the computation. In this case, we can use the Larmor formula [cf. eq. (14.22) of Jackson]

for the instantaneous power emitted by a nonrelativistic, accelerated charge,

P =2

3

z2e2

c3|a|2 = 2

3

z2e2

m2c3

dVdr2

. (49)

where eq. (48) has been used for the acceleration. To compute W, we first consider theenergy emitted by the radiation from the initial position of the particle at r = until thepoint of closest approach to the origin, rmin. Then,

7

W =

P

r(t)

dt =

rmin

P(r)dt

drdr =

m

2

rmin

P(r)dr

E V(r) . (50)

Since we are neglecting the energy loss in computing the instantaneous acceleration of theparticle, the energy loss of the particle as it moves from the distance of closest approach backout to infinity again yields eq. (50). Hence the total energy radiated during < t < is just twice that of eq. (50),

W = 2

m

2

rmin

P(r)dr

V(rmin) V(r) ,

after employing eq. (46). Finally, we substitute for P using eq. (49), which yields

W =4

3

z2e2

m2c3m

2

rmindV

dr2

drV(rmin) V(r)

, (51)

To justify this computation, we would have to show that

W 12

mv20 , (52)

where E 12

mv20 and v0 is the initial velocity at time t = . Since the motion is non-relativistic, we have v0 c, and one can check that for for reasonable potentials, eq. (52)is satisfied [cf. eq. (56)].

(b) If the interaction is a Coulomb potential V(r) = zZ e2/r, show that the total energy

radiated isW =

8

45

zmv50Zc3

,

6To properly take the radiation loss into account, we must address the question of radiation reaction,which is treated in Chapter 16 of Jackson. It turns out that this topic involves numerous subtleties, notall of which are completely understood. For further details, check out the first few sections of Chapter 16.

7Note that we use the minus sign in eq. (47) when the particle moves toward the origin. This minussign is then used to reverse the limits of integration in eq. (50).

20


21/22


22/22

It follows that rmin

dr

r7/21

r

rmin 1

=16

15r5/2min

.

Plugging this result back into eq. (53), we end up with

W =4

3

z2e5

m2c3(zZ)3/2

m

2

16

15r5/2min

. (55)

We can simplify this expression by writing E = 12

mv20 , where v0 is the initial velocity ofthe particle at t = . Using eq. (46),

1

2mv20 = V(rmin) =

zZ e2

rmin.

Solving for rmin and inserting this result back into eq. (55) yields our final result,

W =64

45

z2e5

m2c3(zZ)3/2

m

2

mv20

2zZ e2

5/2=

8

45

zmv50Zc3

.

The condition of eq. (52) then implies that

16z

45Z

v0c

3 1 , (56)

which is always satisfied for non-relativistic motion [assuming that z/Z O(1)].

22

emii4sol_13.pdf

Documents