Chapter 2 Particle Motion in Electric and Magnetic Fields

Considering E and B to be given, we study the trajectory of particles under the inuence of Lorentz force F = q (E + v B) (2.1)

2.1

Electric Field Alone

m

dv = qE (2.2) dt Orbit depends only on ratio q/m. Uniform E uniform acceleration. In one-dimension z, Ez trivial. In multiple dimensions directly analogous to particle moving under inuence of q gravity. Acceleration gravity g m E. Orbits are parabolas. Energy is conserved taking

z x

qE __ m

Figure 2.1: In a uniform electric eld, orbits are parabolic, analogous to gravity. into account potential energy P.E. = q electric potential (2.3)

34

[Proof if needed, regardless of E spatial variation,

dv d .v = q.v = q dt dt d d 1 2 = (q) mv dt dt 2 m i.e. 1 mv 2 + q = const.] 2 A particle gains kinetic energy q when falling through a potential drop-. So consider the acceleration and subsequent analysis of particles electrostatically: How much deection
z

(2.4) (2.5)

=0 d

Detecting +a /2 Screen E z0 L

Particle Source

a /2 Electrostatic Analyser

Figure 2.2: Schematic of electrostatic acceleration and analysis.

2 will there be? After acceleration stage KE = 1 mvx = qs 2

vx =

2qs m

(2.6)

Supposing Ea , eld of analyser, to be purely z, this velocity is subsequently constant. Within the analyser dvz q q x m = qEa vz = Ea t = Ea . (2.7) dt m m vx So q t2 q 1 x2 z = vz dt = Ea = Ea 2 . (2.8) m 2 m 2 vx Hence height at output of analyser is zo = q 1 L2 q 1 m Ea 2 = Ea L2 m 2 vx m 2 (2qs ) 1 Ea 2 1 a L2 = L =+ 4 s d 4 s

(2.9)

35

using Ea = a /d. Notice this is independent of q and m! We could see this directly by eliminating the time from our fundamental equations noting d d = v (= v.) dt d

with

v=

2q ( s ) m

or

2q Es v= s +

q (2.10)

if there is initial energy Es . So equation of motion is m q

2q ( s ) d 2q ( s ) d d dx x = 2 s = Ea = , (2.11) s m d m d d d

which is independent of q and m. Trajectory of particle in purely electrostatic eld depends only on the eld (and initial particle kinetic energy/q). If initial energy is zero, cant deduce anything about q, m.

2.2

Electrostatic Acceleration and Focussing

Accelerated charged particle beams are widely used in science and in everyday applications. Examples: X-ray generation from e-beams (Medical, Industrial) Electron microscopes Welding. (e-beam) Surface ion implantation Nuclear activation (ion-beams) Neutron generation Television and (CRT) Monitors For applications requiring < few hundred keV energy electrostatic acceleration is easiest, widest used. Schematically

Source of Electrons or Ions s Beam t

Sample or Target

Figure 2.3: Obtaining dened energy from electrostatic acceleration is straightforward in principle. Beam focussing and transport to the target is crucial.

36

Clearly getting the required energy is simple. Ensure the potential dierence is right and particles are singly charged: Energy (eV ) Potential V. More interesting question: How to focus the beam? What do we mean by focussing?

Parallel Rays

(a) Optical Focussing

Object Focal Point

Image

(b) Particle beam Focussing

Lens

Figure 2.4: Analogy between optical and particle-beam focussing. What is required of the Lens? To focus at a single spot we require the ray (particle path) deviation from a thin lens to be systematic. Specically, all initially parallel rays converge to a point if the lens deviates their direction by such that

Deviation Angle Distance from axis r Parallel Rays f

Figure 2.5: Requirement for focussing is that the angular deviation of the path should be a linear function of the distance from the axis. r = f tan (2.12)

and for small angles, , r = f . This linear dependence ( = r/f )) of the deviation , on distance from the axis, r, is the key property. Electrostatic Lens would like to have (e.g.) Ea r (2.13) a but the lens cant have charged solids in its middle because the beams must pass through so (initially) = 0 .E = 0. Consequently pure Er is impossible (0 = .E = 1 (rEr )/r = r Er = 37

2Ea /a Ea = 0). For an axisymmetric lens (/ = 0) we must have both Er and Ez . Perhaps the simplest way to arrange appropriate Er is to have an aperture between two

Potential Contours

E1

E2

Figure 2.6: Potential variation near an aperture between two regions of dierent electric eld gives rise to focussing. regions of unequal electric eld. The potential contours bow out toward the lower eld region: giving Er . Calculating focal length of aperture Radial acceleration.

r z

Figure 2.7: Coordinates near an aperture.

dvr q = Er dt m 38

(2.14)

So

dvr 1 dvr q Er = = dz vz dt m vz .E = 0

(2.15)

But

1 (rEr ) Ez + =0 z r r Near the axis, only the linear part of Er is important i.e. Er Er (r, z) r r
r=0 So
1 Er rEr 2 r r r
r=0 and thus
Er Ez 2
+
=0 r
r=0 z and we may write Er 1 rEz /z. Then 2 dvr qr Ez
= , dz 2mvz z

(2.16)

(2.17)

(2.18)

(2.19)

(2.20)

which can be integrated approximately assuming that variations in r and vz can be neglected in lens to get qr vr = [vr ]f inal = [Ez ]2 (2.21) initial 1 2mvz The angular deviation is therefore = and the focal length is f = r/
2 2mvz 4E f= = q (Ez2 Ez1 ) q (Ez2 Ez1 )

vr +qr = [Ez2 Ez1 ] 2 vz 2mvz

(2.22)

(2.23)

When E1 is an accelerating region and E2 is zero or small the lens is diverging. This means that just depending on an extractor electrode to form an ion beam will give a diverging beam. Need to do more focussing down stream: more electrodes.

2.2.1

Immersion Lens

Two tubes at dierent potential separated by gap In this case the gap region can be thought of as an aperture but with the electric elds E1 , E2 the same (zero) on both sides. Previous eect is zero. However two other eects, neglected previously, give focussing: 1. vz is not constant. 2. r is not constant. Consider an accelerating gap: q(2 1 ) < 0. 39

Source

Diverging

Extraction Electrode

Figure 2.8: The extraction electrode alone always gives a diverging beam.

Region 1 Tube

Region 2

Figure 2.9: An Immersion Lens consists of adjacent sections of tube at dierent potentials. Eect (1) ions are converged in region 1, diverged in region 2. However because of zacceleration, vz is higher in region 2. The diverging action lasts a shorter time. Hence overall converging. Eect (2) The electric eld Er is weaker at smaller r. Because of deviation, r is smaller in diverging region. Hence overall converging. For a decelerating gap you can easily convince yourself that both eects are still converging. [Time reversal symmetry requires this.] One can estimate the focal length as 1 3
q
f 16
E 2
2

2 z

dz
r=0

(for weak focussing)

(2.24)

but numerical calculations give the values in gure 2.10 where 1 = E/q. Here E is the energy in region 1. Eect (2) above, that the focussing or defocussing deviation is weaker at points closer to the axis, means that it is a general principle that alternating lenses of equal converging and diverging power give a net converging eect. This principle can be considered to be the basis for

40

Image removed due to copyright restrictions.

Figure 2.10: Focal length of Electrostatic Immersion Lenses. Dependence on energy per unit charge () in the two regions, from S.Humphries 1986

2.2.2

Alternating Gradient Focussing

Idea is to abandon the cylindrically symmetric geometry so as to obtain stronger focussing. Consider an electrostatic conguration with Ez = 0 and Ex = Since .E = 0, we must have
dEx dEy dEy + =0 = const dx dx dy This situation arises from a potential = x y

2 2

dEx x with dx

dEx = const. dx

(2.25)

Ey =

dEx y dx

(2.26)

1 dEx 2 dx

(2.27)

so equipotentials are hyperbolas x2 y 2 = const. If qdEx /dx is negative, then this eld is converging in the x-direction, but dEy /dy = dEx /dx, so it is, at the same time, diverging in the y-direction. By using alternating sections of +ve and -ve dEx /dx a net converging focus can be obtained in both the x and y directions. This alternating gradient approach is very important for high energy particle accelerators, but generally magnetic, not electrostatic, elds are used. So well go into it more later.

41

2.3

Uniform Magnetic eld

dv = q (v B) (2.28) dt Take B in z-direction. Never any force in z -dir. vz = constant. Perpendicular dynamics are separate.

z x

v^B

Figure 2.11: Orbit of a particle in a uniform magnetic eld.

2.3.1

Brute force solution:

vx = vx q q vy B vy = vx B m m qB 2 qB 2 = vx vy = vy m m qB vy = v cos t m m qB y=v sin t + y0 qB m (2.29) (2.30)

Solution

qB vx = v sin t m m qB x = v cos t + x0 qB m

(2.31)

the equation of a circle. Center (x0 , y0 ) and radius (vm/qB) are determined by initial conditions.

2.3.2

Physics Solution

1. Magnetic eld force does no work on particle because F v. Consequently total |v| is constant. 2. Force is thus constant, to v. Gives rise to a circular orbit. 3. Centripetal acceleration gives vr = Larmor (or gyro) Radius.
2

F orce mass

= q vB i.e. r = m

mv . qB

This radius is called the

4. Frequency of rotation v = qB is called the Cyclotron frequency (angular fre r m quency, s1 , not cycles/sec, Hz).

42

When we add the constant vz we get a helical orbit. Cyclotron frequency = qB/m depends only on particle character q, m and B-strength not v (non relativistically, see aside). Larmor Radius r = mv/qB depends on particle momentum mv. All (non-relativistic) particles with same q/m have same . Dierent energy particles have dierent r. This variation can be used to make momentum spectrometers.

2.3.3

Relativistic Aside
d p = q ([E+] v B) dt

Relativistic dynamics can be written (2.32)

where relativistic momentum is m0 v p = mv = 2 1 v2 c Mass m is increased by factor v2 = 1 2 c

1
2

(2.33)

(2.34)

relative to rest mass m0 . Since for E = 0 the velocity |v| = const, is also constant, and so is m. Therefore dynamics of a particle in a purely magnetic eld can be calculated as if it were non-relativistic: m dv/dt = q(v B), except that the particle has mass greater by factor than its rest mass.

2.3.4

Momentum Spectrometers
Unform B v Slit x Detection Plane

Particles passing vertically through slit take dierent paths depending on mv/q. By meaDifferent Momenta

Figure 2.12: Dierent momentum particles strike the detection plane at dierent positions. suring where a particle hits the detection plane we measure its momentum/q : 2 mv =x qB : 43
mv Bx = q 2 . (2.35)

Why make the detection plane a diameter? Because detection position is least sensitive to velocity direction. This is a form of magnetic focussing. Of course we dont need to make the full 360 , so analyser can be reduced in size.
(a) v Different Angles Slit Focus Unform B (b) Unform B v Slit Entrance Detection Plane Detector Array Different Momenta

Figure 2.13: (a) Focussing is obtained for dierent input angles by using 180 degrees of orbit. (b) The other half of the orbit is redundant. Even so, it may be inconvenient to produce uniform B of sucient intensity over suciently large area if particle momentum is large.

2.3.5

Historical Day Dream (J.J. Thomson 1897)

Cathode rays: how to tell their charge and mass? Electrostatic Deection Tells only their energy/q = E/q and we have no independent way to measure E since the same quantity E/q just equals accelerating potential, which is the thing we measure. Magnetic Deection The radius of curvature is r= mv qB (2.36)

So combination of electrostatic and electromagnetic gives us

1 mv 2 2

q Hence

= M1

and

mv = M2 q

(2.37)

2M1 q = . 2 M2 m

(2.38)

We can measure the charge/mass ratio. In order to complete the job an independent measure of q (or m) was needed. Millikan (1911-13). [Actually Townsend in J.J. Thomsons lab had an experiment to measure q which was within factor 2 correct.]

44

2.3.6

Practical Spectrometer

In fusion research fast ion spectrum is often obtained by simultaneous electrostatic and electromagnetic analysis E parallel to B. This allows determination of E/q and q/m velocity of particle [E = 1 mv 2 ]. Thus e.g. deuterons and protons can be distinguished. 2

Particle E

y x x

Figure 2.14: E parallel to B analyser produces parabolic output locus as a function of input velocity. The loci are dierent for dierent q/m. However, He4 and D2 have the same of ion orbits.
q m

so one cant distinguish their spectra on the basis

2.4

Dynamic Accelerators

In addition to the electrostatic accelerators, there are several dierent types of accelerators based on time-varying elds. With the exception of the Betatron, these are all based on the general principle of arranging for a resonance between the particle and the oscillating elds such that energy is continually given to the particle. Simple example

1+ 1

1-

2+ 2

2-

3+ 3

3-

Figure 2.15: Sequence of dynamically varying electrode potentials produces continuous ac celeration. Values at 3 times are indicated. Particle is accelerated through sequence of electrodes 3 at times (1) (2) (3). The potential of electrode is raised from negative to +ve while particle is inside electrode. So at each gap it sees an accelerating Ez . Can be thought of as a successive moving potential hill: Wave of potential propagates at same speed as particle so it is continuously accelerated. Historically earliest widespread accelerator based on this principle was the cyclotron.

45

z
Figure 2.16: Oscillating potentials give rise to a propagating wave.

+ Electric Poles Orbit B

Plan

Magnetic Poles

Elevation

Figure 2.17: Schematic of a Cyclotron accelerator.

2.4.1

Cyclotron

Take advantage of the orbit frequency in a uniform B-eld = qB . Apply oscillating m potential to electric poles, at this frequency. Each time particle crosses the gap (twice/turn) it sees an accelerating electric eld. Resonant frequency f= qB = = 1.52 107 B 2 m2 Hz (2.39)

15.2 MHz/T for protons. If magnet radius is R particle leaves accelerator when its Larmor radius is equal to R mv 1 1 q2 2 2 = R mv 2 = B R (2.40) qB 2 2m If iron is used for magnetic pole pieces then B < 2T (where it saturates). Hence larger accelerator is required for higher energy E R2 . [But stored energy in magnet R2 R3 ].

46

2.4.2 Limitations of Cyclotron Acceleration: Relativity

Mass increase (1 v 2 /c2 ) 2 breaks resonance, restricting maximum energy to 25MeV (protons). Improvement: sweep oscillator frequency (downward). Synchrocyclotron al lowed energy up to 500MeV but reduced ux. Alternatively: Increase B with radius. Leads to orbit divergence parallel to B. Compensate with azimuthally varying eld for focussing AVF-cyclotron. Advantage continuous beam.
1

2.4.3 Synchrotron
Vary both frequency and eld in time to keep beam in resonance at constant radius. High energy physics (to 800 GeV).

2.4.4 Linear Accelerators

Avoid limitations of electron synchrotron radiation. Come in 2 main types. (1) Induction (2) RF (linacs) with dierent pros and cons. (RF for highest energy electrons). Electron acceleration: v = c dierent problems from ion.

2.5 Magnetic Quadrupole Focussing (Alternating Gra dient)

Magnetic focussing is preferred at high particle energy. Why? Its force is stronger.
Magnetic force on a relativistic particle qcB.
Electric force on a relativistic particle qE.
E.g. B = 2T cB = 6 108 same force as an electric eld of magnitude 6 108 V /m = 0.6MV/mm! However magnetic force is perpendicular to B so an axisymmetric lens would like to have purely azimuthal B eld B = B . However this would require a current right
B

v v^B

Figure 2.18: Impossible ideal for magnetic focussing: purely azimuthal magnetic eld. where the beam is:

B.d = o I.

(2.41)

Axisymmetric magnetic lens is impossible. However we can focus in one cartesian direction (x, y) at a time. Then use the fact that successive combined focus-defocus has a net focus.

47

2.5.1

Preliminary Mathematics

Consider z = 0 purely transverse eld (approx) Bx , By . This can be represented by B = A with A = zA so A = (zA) = z A (since z = 0). In the vacuum region j = 0 (no current) so

0 = B = (z A) = z2 A + (z.) A

=0

(2.42)

i.e. 2 A = 0. A satises Laplaces equation. Notice then that solutions of electrostatic problems, 2 = 0 are also solutions of (2-d) vacuum magnetostatic problems. The same solution techniques work.

2.5.2

Multipole Expansion

Potential can be expanded about some point in space in a kind of Taylor expansion. Choose origin at point of expansion and use coordinates (r, ), x = r cos , y = r sin . 2 A = 1 A 1 2A r + 2 2 =0 r r r r (2.43)

Look for solutions in the form A = u(r).w(). These require d2 w = const. w d2 d du r = const. u. dr dr Hence w solutions are sines and cosines r w = cos n or sin n and (2.44)

(2.45)

(2.46)

where n2 is the constant in the previous equation and n integral to satisfy periodicity. Correspondingly u = rn or ln r , rn (2.47) These solutions are called cylindrical harmonics or (cylindrical) multipoles: 1 rn cos n rn sin n ln r rn cos n rn sin n (2.48)

If our point of expansion has no source at it (no current) then the right-hand column is ruled out because no singularity at r = 0 is permitted. The remaining multipoles are 1 r cos (= x) r2 cos 2 = r2 (cos2 sin2 ) = x2 y 2 Higher orders 48 constant irrelevant to a potential uniform eld, A x non-uniform eld neglected.

The second order solution, x2 y 2 is called a quadrupole eld (although this is something of a misnomer). [Similarly r3 cos 3 hexapole, r4 cos octupole.] We already dealt with this potential in the electric case. A = x2 y 2 = 2x 2y x y So z A = 2xy 2yx Force on longitudinally moving charge: F = qv B = qv ( A) = qv ( A) = q (v.z) A qvz A z Magnetic quadrupole force is identical to electric quadrupole force replacing Avz (2.53) (2.51) (2.52) (2.50)

(2.49)

Consequently focussing in x-direction defocussing in y-direction but alternating gradients give net focussing. This is basis of all strong focussing.

2.6

Force on distributed current density

F = q (E + v B) (2.54)

We have regarded the Lorentz force law

as fundamental. However forces are generally measured in engineering systems via the in teraction of wires or conducting bars with B-elds. Historically, of course, electricity and magnetism were based on these measurements. A current (I) is a ow of charge: Coulombs/s Amp. A current density j is a ow of charge per unit area A/m2 . The charge is carried by particles: j= ni vi qi (2.55)
species i

Hence total force on current carriers per unit volume is F=

ni qi (vi B) = j B

(2.56)

Also, for a ne wire carrying current I, if its area is , the current density averaged across the section is I j = (2.57) Volume per unit length is . And the force/unit length = j B. = I B perpendicular to the wire.

49

2.6.1

Forces on dipoles

We saw that the eld of a localized current distribution, far from the currents, could be approximated as a dipole. Similarly the forces on a localized current by an external magnetic eld that varies slowly in the region of current can be expressed in terms of magnetic dipole. [Same is true in electrostatics with an electric dipole]. Total force F=

j Bd3 x

(2.58)

where B is an external eld that is slowly varying and so can be approximated as B (x ) = B0 + (x .) B where the tensor B (Bj /xi ) is simply a constant (matrix). Hence F = =

(2.59)

j Bo + j (x .B) d3 x jd3 x Bo +

j (x .B) d3 x

(2.60)

The rst term integral is zero and the second is transformed by our previous identity, which can be written as x

(x j) d x = 2x.

jx d x = 2x.
B) xi j

x jd3 x

(2.61)

for any x. Use the quantity B for x (i.e. xi

giving (2.62)

1 (x j) d3 x B = m B = j (x .) Bd3 x 2

This tensor identity is then contracted by an internal cross-product [ijk Tjk ] to give the vector identity (m ) B = j [(x .) B] d3 x (2.63) Thus F = (m ) B = (m.B) m (.B) (remember the operates only on B not m). This is the force on a dipole: F = (m.B) Total Torque (Moment of force) is M = =

(2.64)

(2.65)

x (j B) d3 x j (x .B) B (x .j) d3 x 50

(2.66) (2.67)

B here is (to lowest order) independent of x : B0 so second term is zero since

x .jd3 x =

1 2 . |x | j |x |2 .j d3 x = 0. 2

(2.68)

The rst term is of the standard form of our identity. M = B.

1 x jd x = B (x j) d3 x 2
3

(2.69) (2.70)

M=mB

Moment on a dipole.

y x z dy dx I

Figure 2.19: Elementary circuit for calculating magnetic force.

2.6.2

Force on an Elementary Magnetic Moment Circuit

Consider a plane rectangular circuit carrying current I having elementary area dxdy = dA. Regard this as a vector pointing in the z direction dA. The force on this current in a eld B(r) is F such that Fx = Idy [Bz (x + dx) Bz (x)] = Idydx Fy Fz Bz x Bz = Idx [Bz (y + dy) Bz (y)] = Idydx y = Idx [By (y + dy) By (y)] Idy [Bx (x + dx) Bx (x)] Bx By Bz = Idxdy + = Idydx x y z (2.71) (2.72)

(2.73)

(using .B = 0). Hence, summarizing: F = IdydxBz . Now dene m = IdA = Idydxz and take it constant. Then clearly the force can be written F = (B.m) or strictly (B).m. (2.74)

51

N m S M

Figure 2.20: Moment on a bar magnet in a uniform eld.

|B| Net Force: S m N N m S B

Figure 2.21: A magnetic moment in the form of a bar magnet is attracted or repelled toward the stronger eld region, depending on its orientation.

2.6.3

Example

Small bar magnet: archetype of dipole. In uniform B feels just a torque aligning it with B. In a uniform eld, no net force. Non-uniform eld: If magnet takes its natural resting direction, m parallel to B, force is F = m|B| (2.75)

A bar magnet is attracted to high eld. Alternatively if m parallel to minus B the magnet points other way F = m|B| repelled from high |B|. (2.76) Same would be true for an elementary circuit dipole. It is attracted/repelled according to whether it acts to increase or decrease B locally. A charged particle moving in its Larmor orbit is always diamagnetic: repelled from high |B|.

2.6.4

Intuition

There is something slightly non-intuitive about the natural behavior of an elementary wire circuit and a particle orbit considered as similar to this elementary circuit. Their currents ow in opposite directions when the wire is in its stable orientation. The reason is that the strength of the wire sustains it against the outward magnetic expansion force, while the particle needs an inward force to cause the centripetal acceleration.

52

 |B| Net Force: Paramagnetic Attracted Diamagnetic Repelled B

Figure 2.22: Elementary circuit acting as a dipole experiences a force in a non-uniform magnetic eld.

Stable Orientation j^B outwards. Wire Loop

To balance acceleration j^B inwards. vq Orbiting Particle

Figure 2.23: Dierence between a wire loop and a particle orbit in their natural orientation.

2.6.5 Angular Momentum

If the local current is made up of particles having a constant ratio of charge to mass: q/M say (Notational accident m is magnetic moment). Then the angular momentum is L = 1 qi xi vi . So i Mi xi v and magnetic moment is m = 2 q m = L. Classical (2.77) 2M This would also be true for a continuous body with constant (charge density)/(mass density) (/m ). Elementary particles, e.g. electrons etc., have spin with moments m, L. However they do not obey the above equation. Instead q m=g L (2.78) 2M with the Land g-factor ( 2 for electrons). This is attributed to quantum and relativistic e eects. However the classical value might not occur if /m were not constant. So we should not be surprised that g is not exactly 1 for particles spin.

2.6.6 Precession of a Magnetic Dipole (formed from charged par ticle)

The result of a torque m B is a change in angular momentum. Since m = gLq/2M we have dL q =mB=g (L B)) (2.79) dt 2M 53

m^B

m^B L m

Figure 2.24: Precession of an angular momentum L and aligned magnetic moment m about the magnetic eld.
q This is the equation of a circle around B. [Compare with orbit equation dv = m v B]. The dt direction of L precesses like a tilted top around direction of B with a frequency

=g

qB 2M

(2.80)

For an electron (g = 2) this is equal to the cyclotron frequency. For protons g = 2 2.79 [Written like this because spin is 1 ]. For neutrons g = 2 (1.93). 2 Precession frequency is thus f= electron = (28GHz) (B/Tesla) 2 proton = (43MHz) (B/Tesla) 2 (2.81) (2.82)

This is the (classical) basis of Nuclear Magnetic Resonance but of course that really needs QM.

54