Physics · Book 3 · Bachelor Year 1

University Physics — Year 1

University Physics — Year 1 · Bachelor Year 1

17Charged Particles in E and B Fields

Inside the oscilloscope that drew the traces of Chapter 7, a beam of electrons is steered across a screen by the voltage between two plates; inside the magnet of a mass spectrometer, ions of two isotopes separate by millimeters because one is eight percent heavier than the other; inside a hospital’s cyclotron, protons spiral outward through two hundred turns before leaving at a third of the speed of light to burn a tumor. All three are the same dynamics: a charge in a uniform field, electric or magnetic, and Newton’s law. This chapter works it out — the parabola in an electric field, the circle and the helix in a magnetic field — and turns it into the instruments that sort, accelerate and measure charged particles.

Lawrence’s 60-inch cyclotron at Berkeley, 1939 (US Department of Energy): between the poles of the magnet, protons spiral outward as they are kicked twice per turn, exactly as in the weekend problem.
Lawrence’s 60-inch cyclotron at Berkeley, 1939 (US Department of Energy): between the poles of the magnet, protons spiral outward as they are kicked twice per turn, exactly as in the weekend problem.

17.1 The Lorentz force

Definition 17.1 (Lorentz force)

A particle of charge qq moving at velocity v\vect v in an electric field E\vect E and a magnetic field B\vect B receives the force

F=q(E+vB).\vect F = q\big(\vect E + \vect v\wedge\vect B\big) .

The electric part is along E\vect E (or opposite, for q<0q < 0) and independent of the motion; the magnetic part is perpendicular both to v\vect v and to B\vect B, proportional to the speed, and vanishes for a particle at rest. Units: E\vect E in V/m\mathrm{V}/\mathrm{m} (or N/C\mathrm{N}/\mathrm{C}), B\vect B in teslas (T\mathrm{T}): 1T=1N/(Am)1\,\mathrm{T} = 1\,\mathrm{N}/(\mathrm{A}\,\mathrm{m}).

Proposition 17.2 (Work of the Lorentz force)

The magnetic force does no work: its power q(vB)vq(\vect v\wedge\vect B)\cdot\vect v is zero, so a magnetic field alone changes the direction of the velocity but never the speed or the kinetic energy. The electric force has power qEvq\vect E\cdot\vect v, and between two points its work is q(VAVB)=qUABq(V_A - V_B) = qU_{AB} (Chapter 27): a particle of charge qq crossing a voltage UU gains the kinetic energy qU\abs{qU}.

Proof. vBv\vect v\wedge\vect B \perp \vect v. The electrostatic field derives from the potential, E=gradV\vect E = -\overrightarrow{\operatorname{grad}}V, so qEq\vect E is conservative with Ep=qVE_p = qV.

Definition 17.3 (The electron-volt)

The electron-volt is the kinetic energy gained by an elementary charge crossing one volt: 1eV=1.60×1019J1\,\mathrm{eV} = 1.60 \times 10^{-19}\,\mathrm{J}. Chemistry lives at a few eV, X-ray tubes at tens of keV, nuclear physics at MeV, particle physics at GeV and TeV.

Remark 17.4 (Why gravity is ignored)

For an electron in a field of 1kV/m1\,\mathrm{kV}/\mathrm{m} — a modest laboratory value — the electric force is eE=1.6×1016NeE = 1.6 \times 10^{-16}\,\mathrm{N}, against a weight mg=9×1030Nmg = 9 \times 10^{-30}\,\mathrm{N}: 101310^{13} times smaller. Weight is dropped from every calculation of this chapter without a second thought. (The relativistic limit is another matter: the formulas below hold for 12mv2mc2\tfrac12 mv^2 \ll mc^2, i.e. below about 50keV50\,\mathrm{keV} for electrons and 100MeV100\,\mathrm{MeV} for protons; see the Year 3 volume.)

17.2 Uniform electric field

Proposition 17.5 (Acceleration and deflection)

In a uniform field E\vect E the acceleration a=qE/m\vect a = q\vect E/m is constant: the motion is that of a projectile, with qE/mq\vect E/m in place of g\vect g.

  • Released from rest and accelerated through a voltage UU, the particle reaches v=2qU/mv = \sqrt{2\abs{qU}/m}.
  • Entering at speed v0v_0 along xx a region of length LL where E=Eey\vect E = E\vect e_y, it follows the parabola y=(qE/2mv02)x2y = (qE/2mv_0^2)x^2 and leaves deflected by the angle α\alpha with tanα=qEL/mv02\tan\alpha = qEL/mv_0^2, as if it had come in a straight line from the middle of the plates.

Proof. Energy: 12mv2=qU\tfrac12 mv^2 = \abs{qU}. Deflection: x=v0tx = v_0t, y=12(qE/m)t2y = \tfrac12(qE/m)t^2; at the exit,  ⁣dy/ ⁣dx=(qE/m)(x/v02)\dd y/\dd x = (qE/m)(x/v_0^2) evaluated at LL; the tangent at the exit, y=(qEL/mv02)(xL/2)y = (qEL/mv_0^2)(x - L/2), crosses the axis at x=L/2x = L/2.

Electrostatic deflection: between the plates the electron (here q < 0, deflected toward the positive plate) follows a parabola; beyond them it flies straight, as if coming from the midpoint of the plates — the geometry of the cathode-ray oscilloscope.
Electrostatic deflection: between the plates the electron (here q<0q < 0, deflected toward the positive plate) follows a parabola; beyond them it flies straight, as if coming from the midpoint of the plates — the geometry of the cathode-ray oscilloscope.

Example 17.6 (The oscilloscope tube)

Electrons accelerated through U0=2.0kVU_0 = 2.0\,\mathrm{kV} (v0=2.65×107m/sv_0 = 2.65 \times 10^{7}\,\mathrm{m}/\mathrm{s}, a tenth of cc — just non-relativistic) cross 4.0cm4.0\,\mathrm{cm} plates 1.0cm1.0\,\mathrm{cm} apart carrying 100V100\,\mathrm{V}: E=1.0×104V/mE = 1.0 \times 10^{4}\,\mathrm{V}/\mathrm{m}, tanα=eEL/mv02=EL/2U0=0.10\tan\alpha = eEL/mv_0^2 = EL/2U_0 = 0.10, and on a screen 30cm30\,\mathrm{cm} past the plates’ center the spot moves by 3.0cm3.0\,\mathrm{cm}0.3mm0.3\,\mathrm{mm} per volt, the tube’s sensitivity, set purely by geometry and the accelerating voltage.

17.3 Uniform magnetic field

Theorem 17.7 (Motion in a uniform magnetic field)

A particle of charge qq, mass mm, velocity v0\vect v_0 perpendicular to a uniform B\vect B, describes a circle at constant speed v0v_0, of radius and angular frequency

R=mv0qB,ωc=qBm,R = \frac{mv_0}{\abs q B}, \qquad \omega_c = \frac{\abs q B}{m},

the cyclotron frequency ωc\omega_c being independent of the speed; the sense of rotation depends on the sign of qq. If v0\vect v_0 also has a component vv_\parallel along B\vect B, that component is unchanged and the trajectory is a helix of radius mv/qBmv_\perp/\abs qB and pitch 2πmv/qB2\pi mv_\parallel/\abs qB wound around the field lines.

Proof. The force qvBq\vect v\wedge\vect B is perpendicular to v\vect v and to B\vect B; the speed is constant (no work) and the component along B\vect B sees no force. In the plane perpendicular to B\vect B the force has the constant norm qvB\abs qv_\perp B, always perpendicular to the velocity: uniform circular motion (Theorem 11.13) with mv2/R=qvBmv_\perp^2/R = \abs qv_\perp B. Then ωc=v/R\omega_c = v_\perp/R; the period 2πm/qB2\pi m/\abs qB times vv_\parallel is the pitch.

Left: a charge whose velocity is perpendicular to B circles at the cyclotron frequency qB/m, the magnetic force supplying the centripetal acceleration. Right: with a velocity component along B, a helix winds around the field line — the motion of electrons trapped in the Earth’s field.
Left: a charge whose velocity is perpendicular to B\vect B circles at the cyclotron frequency qB/m\abs qB/m, the magnetic force supplying the centripetal acceleration. Right: with a velocity component along B\vect B, a helix winds around the field line — the motion of electrons trapped in the Earth’s field.

Example 17.8 (Orders of magnitude)

A 1keV1\,\mathrm{keV} electron (v=1.9×107m/sv = 1.9 \times 10^{7}\,\mathrm{m}/\mathrm{s}) in 1mT1\,\mathrm{mT}: R=11cmR = 11\,\mathrm{cm}, ωc=1.8×108rad/s\omega_c = 1.8 \times 10^{8}\,\mathrm{rad}/\mathrm{s}, one turn in 36ns36\,\mathrm{ns} — the same period at any energy, which is what makes the cyclotron possible. A proton of the solar wind (1keV1\,\mathrm{keV}, 4.4×105m/s4.4 \times 10^{5}\,\mathrm{m}/\mathrm{s}) in the Earth’s 50µT50\,\text{µ}\mathrm{T}: R=92mR = 92\,\mathrm{m}; it cannot cross the field lines, only spiral along them toward the poles — the aurora.

17.4 Instruments

Proposition 17.9 (Velocity selector and mass spectrometer)

In crossed uniform fields EB\vect E \perp \vect B, a particle moving perpendicular to both passes undeflected iff v=E/Bv = E/B, whatever its charge and mass (velocity selector). In a magnetic field alone a particle of known speed describes a semicircle of diameter 2mv/qB2mv/\abs qB: measuring it gives m/qm/q (mass spectrometer). Ions accelerated through a voltage UU then bent in BB land at

R=1B2mUq:R = \frac{1}{B}\sqrt{\frac{2mU}{\abs q}} :

two isotopes of masses mm and m+Δmm + \Delta m separate by ΔRRΔm/2m\Delta R \approx R\,\Delta m/2m.

Proof. qE+qvB=0q\vect E + q\vect v\wedge\vect B = \vect 0 requires vB=EvB = E with the right orientation. R=mv/qBR = mv/qB with v=2qU/mv = \sqrt{2qU/m}; differentiate lnR=12lnm+const\ln R = \tfrac12\ln m + \text{const}.

Example 17.10 (Carbon isotopes)

Singly charged 12^{12}C and 13^{13}C ions accelerated through 20kV20\,\mathrm{kV} in 0.50T0.50\,\mathrm{T}: R=14.1cmR = 14.1\,\mathrm{cm} and 14.7cm14.7\,\mathrm{cm}, 5.8mm5.8\,\mathrm{mm} apart — two separate spots on the detector, and the ratio of their intensities is the isotopic abundance, the raw material of radiocarbon dating and of every isotope analysis.

Proposition 17.11 (Cyclotron)

Two hollow half-cylinders (“dees”) in a uniform B\vect B, with an alternating voltage UU between them at the cyclotron frequency: the particle circles inside each dee (no field there), is accelerated by qU\abs qU at each of the two gap crossings per turn, and spirals outward with RvR \propto v; at the extraction radius RmaxR_{\max} its kinetic energy is

Ek=q2B2Rmax22m,E_k = \frac{q^2B^2R_{\max}^2}{2m},

independent of the voltage, which only sets the number of turns.

Proof. Because ωc\omega_c does not depend on the speed, a fixed-frequency voltage stays in step with the particle turn after turn; v=qBR/mv = \abs qBR/m at radius RR.

The cyclotron: two dees in a uniform magnetic field; the particle circles inside the dees and is kicked by the gap voltage twice per turn, always in phase because its period does not depend on its speed. The energy at the rim depends only on B and the radius.
The cyclotron: two dees in a uniform magnetic field; the particle circles inside the dees and is kicked by the gap voltage twice per turn, always in phase because its period does not depend on its speed. The energy at the rim depends only on BB and the radius.

Example 17.12 (A hospital cyclotron)

Protons, B=1.5TB = 1.5\,\mathrm{T}, Rmax=0.50mR_{\max} = 0.50\,\mathrm{m}: Ek=(1.6×1019×1.5×0.5)2/(2×1.67×1027)=4.3×1012J=27MeVE_k = (1.6\times10^{-19} \times 1.5 \times 0.5)^2/(2 \times 1.67\times10^{-27}) = 4.3 \times 10^{-12}\,\mathrm{J} = 27\,\mathrm{MeV}, at fc=ωc/2π=23MHzf_c = \omega_c/2\pi = 23\,\mathrm{MHz}. With 50kV50\,\mathrm{kV} on the dees, 100keV100\,\mathrm{keV} per turn: 270270 turns in 12µs12\,\text{µ}\mathrm{s}. Above a few tens of MeV the relativistic increase of the mass (the Year 3 volume) desynchronizes the protons; therapy machines at 230MeV230\,\mathrm{MeV} modulate the frequency or the field to compensate — the weekend problem’s subject.

Proposition 17.13 (Hall effect)

A ribbon of thickness dd carrying a current II of carriers of charge qq and density nn, in a field BB perpendicular to the ribbon, develops across its width the Hall voltage

UH=IBnqd:U_H = \frac{IB}{nqd} :

proportional to BB (a magnetometer) and inversely to nn (a probe of the carrier density and sign).

Proof. The magnetic force qvdBqv_dB on the carriers (drift speed vdv_d) pushes them to one edge until the transverse electric field EHE_H balances it: qEH=qvdBqE_H = qv_dB, and UH=EHwU_H = E_Hw across the width ww; with I=nqvdwdI = nqv_dwd, EHw=vdBw=IB/(nqd)E_Hw = v_dBw = IB/(nqd).

17.5 Exercises

Exercise 17.1

Convert 1eV1\,\mathrm{eV}, 13.6eV13.6\,\mathrm{eV}, 1MeV1\,\mathrm{MeV} to joules. Speed of an electron and of a proton accelerated from rest through 1.0kV1.0\,\mathrm{kV}.

Solution

Solution of Exercise 17.1.

1eV=1.60×1019J1\,\mathrm{eV} = 1.60 \times 10^{-19}\,\mathrm{J}, 13.6eV=2.18×1018J13.6\,\mathrm{eV} = 2.18 \times 10^{-18}\,\mathrm{J}, 1MeV=1.60×1013J1\,\mathrm{MeV} = 1.60 \times 10^{-13}\,\mathrm{J}. v=2eU/mv = \sqrt{2eU/m}: electron 1.9×107m/s1.9 \times 10^{7}\,\mathrm{m}/\mathrm{s}; proton 4.4×105m/s4.4 \times 10^{5}\,\mathrm{m}/\mathrm{s} (1836\sqrt{1836} times slower).

Exercise 17.2

A 1.0keV1.0\,\mathrm{keV} electron enters a 1.0mT1.0\,\mathrm{mT} field perpendicularly. Radius, cyclotron angular frequency, period.

Solution

Solution of Exercise 17.2.

v=1.88×107m/sv = 1.88 \times 10^{7}\,\mathrm{m}/\mathrm{s}; R=mv/eB=9.11×1031×1.88×107/(1.6×1019×103)=0.11mR = mv/eB = 9.11\times10^{-31} \times 1.88\times10^7/ (1.6\times10^{-19} \times 10^{-3}) = 0.11\,\mathrm{m}; ωc=eB/m=1.8×108rad/s\omega_c = eB/m = 1.8 \times 10^{8}\,\mathrm{rad}/\mathrm{s}; T=2π/ωc=36nsT = 2\pi/\omega_c = 36\,\mathrm{ns}.

Exercise 17.3

A proton at 1.0×105m/s1.0 \times 10^{5}\,\mathrm{m}/\mathrm{s} in the Earth’s field (50µT50\,\text{µ}\mathrm{T}): radius of its circle. What happens if its velocity makes 3030^\circ with the field?

Solution

Solution of Exercise 17.3.

R=mv/qB=1.67×1027×105/(1.6×1019×5×105)=21mR = mv/qB = 1.67\times10^{-27} \times 10^5/(1.6\times10^{-19} \times 5\times 10^{-5}) = 21\,\mathrm{m}. At 3030^\circ: a helix of radius Rsin30=10mR\sin 30^\circ = 10\,\mathrm{m} and pitch 2π(m/qB)vcos30=114m2\pi(m/qB)v\cos 30^\circ = 114\,\mathrm{m}, drifting along the field line.

Exercise 17.4

A velocity selector has E=10kV/mE = 10\,\mathrm{kV}/\mathrm{m} and B=20mTB = 20\,\mathrm{mT}. Which speed passes? What happens to faster ions, and to the same speed with the opposite charge?

Solution

Solution of Exercise 17.4.

v=E/B=104/0.020=5.0×105m/sv = E/B = 10^4/0.020 = 5.0 \times 10^{5}\,\mathrm{m}/\mathrm{s}. Faster ions: the magnetic force wins, they curve toward the magnetic side; opposite charge at the same speed: both forces reverse, still balanced — the selector is blind to the sign and to the mass.

Exercise 17.5 ★★

Oscilloscope: electrons at 2.0kV2.0\,\mathrm{kV}, plates 4.0cm4.0\,\mathrm{cm} long, 1.0cm1.0\,\mathrm{cm} apart, 100V100\,\mathrm{V} between them, screen 30cm30\,\mathrm{cm} past the plates’ center. Deflection on the screen and sensitivity in mm/V. Why does a higher accelerating voltage reduce the sensitivity?

Solution

Solution of Exercise 17.5.

tanα=EL/2U0=104×0.04/4000=0.10\tan\alpha = EL/2U_0 = 10^4 \times 0.04/4000 = 0.10; deflection 0.30×0.10=3.0cm0.30 \times 0.10 = 3.0\,\mathrm{cm}; sensitivity 0.30mm/V0.30\,\mathrm{mm}/\mathrm{V}. Faster electrons spend less time between the plates: tanα1/U0\tan\alpha \propto 1/U_0.

Exercise 17.6 ★★

Singly charged 12^{12}C and 13^{13}C ions, accelerated through 20kV20\,\mathrm{kV}, enter a 0.50T0.50\,\mathrm{T} field (1u=1.66×1027kg1\,\mathrm{u} = 1.66 \times 10^{-27}\,\mathrm{kg}). Radii and separation of the two spots after a half turn.

Solution

Solution of Exercise 17.6.

R=2mU/e/BR = \sqrt{2mU/e}/B: m12=1.99×1026kgm_{12} = 1.99 \times 10^{-26}\,\mathrm{kg}, R12=2×1.99×1026×2×104/1.6×1019/0.5=14.1cmR_{12} = \sqrt{2 \times 1.99\times10^{-26} \times 2\times10^4/1.6\times10^{-19}}/0.5 = 14.1\,\mathrm{cm}; R13=R1213/12=14.7cmR_{13} = R_{12}\sqrt{13/12} = 14.7\,\mathrm{cm}; separation 2ΔR=1.2cm2\Delta R = 1.2\,\mathrm{cm} between the landing points (radii differ by 5.8mm5.8\,\mathrm{mm}, diameters by twice that).

Exercise 17.7 ★★

A proton cyclotron has B=1.5TB = 1.5\,\mathrm{T} and Rmax=0.50mR_{\max} = 0.50\,\mathrm{m}. Frequency of the dee voltage, maximal kinetic energy, number of turns for 50kV50\,\mathrm{kV} on the dees, time spent inside. By what fraction has the proton’s mass increased at the exit (γ1Ek/mc2\gamma - 1 \approx E_k/mc^2, mc2=938MeVmc^2 = 938\,\mathrm{MeV})?

Solution

Solution of Exercise 17.7.

fc=eB/2πm=23MHzf_c = eB/2\pi m = 23\,\mathrm{MHz}; Ek=e2B2R2/2m=4.3×1012J=27MeVE_k = e^2B^2R^2/2m = 4.3 \times 10^{-12}\,\mathrm{J} = 27\,\mathrm{MeV}; 100keV100\,\mathrm{keV} per turn: 270270 turns; time 270/fc=12µs270/f_c = 12\,\text{µ}\mathrm{s}; γ1=27/938=2.9%\gamma - 1 = 27/938 = 2.9\%.

Exercise 17.8 ★★

An electron at 1.0×106m/s1.0 \times 10^{6}\,\mathrm{m}/\mathrm{s} enters a 2.0mT2.0\,\mathrm{mT} field at 3030^\circ to the field lines. Radius and pitch of its helix; number of turns per meter along the field.

Solution

Solution of Exercise 17.8.

v=5.0×105m/sv_\perp = 5.0 \times 10^{5}\,\mathrm{m}/\mathrm{s}, v=8.7×105m/sv_\parallel = 8.7 \times 10^{5}\,\mathrm{m}/\mathrm{s}: R=mv/eB=9.11×1031×5×105/(3.2×1022)=1.4mmR = mv_\perp/eB = 9.11\times10^{-31} \times 5\times10^5/(3.2\times10^{-22}) = 1.4\,\mathrm{mm}; period 2πm/eB=18ns2\pi m/eB = 18\,\mathrm{ns}, pitch vT=1.5cmv_\parallel T = 1.5\,\mathrm{cm}: about 6464 turns per meter.

Exercise 17.9 ★★

A copper strip 1.0mm1.0\,\mathrm{mm} thick carries 5.0A5.0\,\mathrm{A} in 1.0T1.0\,\mathrm{T} (n=8.5×1028m3n = 8.5 \times 10^{28}\,\mathrm{m}^{-3}). Hall voltage. Why are Hall sensors made of semiconductors (n1×1022m3n \sim 1 \times 10^{22}\,\mathrm{m}^{-3})?

Solution

Solution of Exercise 17.9.

UH=IB/nqd=5.0/(8.5×1028×1.6×1019×103)=3.7×107VU_H = IB/nqd = 5.0/(8.5\times10^{28} \times 1.6\times10^{-19} \times 10^{-3}) = 3.7 \times 10^{-7}\,\mathrm{V}: sub-microvolt. With nn ten million times smaller, a semiconductor gives millivolts — a usable sensor.

Exercise 17.10 ★★★

Prove from Newton’s law that a charge in any magnetic field (uniform or not), with no electric field, keeps a constant speed. A 10keV10\,\mathrm{keV} electron enters a region of strong, non-uniform field: what can and cannot change?

Solution

Solution of Exercise 17.10.

m ⁣dv/ ⁣dt=qvBm\,\dd\vect v/\dd t = q\vect v\wedge\vect B; dot with v\vect v:  ⁣d(12mv2)/ ⁣dt=0\dd(\tfrac12 mv^2)/ \dd t = 0 whatever B(r,t)\vect B(\vect r, t), provided there is no E\vect E. The direction of the velocity, the radius of curvature, the pitch can all change; the speed, the kinetic energy (10keV10\,\mathrm{keV}) cannot.

Exercise 17.11 ★★★

A charge qq starts from rest at the origin in E=Eey\vect E = E\vect e_y and B=Bez\vect B = B\vect e_z. Write the equations of motion, show that the particle drifts on average at the velocity E/BE/B along ex\vect e_x (hint: look for a solution of the form v=u+w\vect v = \vect u + \vect w with u\vect u constant), and describe the trajectory (a cycloid). Numbers for an electron with E=1.0×104V/mE = 1.0 \times 10^{4}\,\mathrm{V}/\mathrm{m}, B=0.10TB = 0.10\,\mathrm{T}.

Solution

Solution of Exercise 17.11.

mv˙x=qvyBm\dot v_x = qv_yB, mv˙y=qEqvxBm\dot v_y = qE - qv_xB, vz=0v_z = 0. With v=(E/B)ex+w\vect v = (E/B)\vect e_x + \vect w: mw˙x=qwyBm\dot w_x = qw_yB, mw˙y=qwxBm\dot w_y = -qw_xB — a pure cyclotron rotation of w\vect w at ωc\omega_c. Starting from rest, w(0)=(E/B)ex\vect w(0) = -(E/B)\vect e_x: w\vect w turns at constant norm E/BE/B, so the velocity is the drift E/BE/B plus a rotating vector of the same norm — a cycloid, with cusps where v=0\vect v = \vect 0. Electron: drift E/B=1.0×105m/sE/B = 1.0 \times 10^{5}\,\mathrm{m}/\mathrm{s}, ωc=1.8×1010rad/s\omega_c = 1.8 \times 10^{10}\,\mathrm{rad}/\mathrm{s}, arch height 2mE/qB2=11µm2mE/qB^2 = 11\,\text{µ}\mathrm{m}.

Exercise 17.12 ★★★

In the oscilloscope of Exercise 17.5, replace the plates by a magnetic deflection coil producing 1.0mT1.0\,\mathrm{mT} over the same 4.0cm4.0\,\mathrm{cm}. Deflection angle (small-angle: the arc of radius RR over a length LL), and deflection on the screen. Why do television tubes, with their large angles, use magnetic rather than electric deflection?

Solution

Solution of Exercise 17.12.

R=mv0/eB=9.11×1031×2.65×107/(1.6×1022)=0.15mR = mv_0/eB = 9.11\times10^{-31} \times 2.65\times10^7/(1.6\times10^{-22}) = 0.15\,\mathrm{m}; αL/R=0.04/0.15=0.27\alpha \approx L/R = 0.04/0.15 = 0.27 (about 1515^\circ): deflection 0.30×0.27=8cm\approx 0.30 \times 0.27 = 8\,\mathrm{cm} for a modest millitesla. Large angles with plates would need kilovolts and long plates; a coil bends the beam through tens of degrees with no voltage at all.

17.6 Problem: Protons against a tumor

Problem 17.1

Weekend problem — the cyclotron behind the wall of a hospital: how protons are spun up to a third of the speed of light, why the simple machine stops working before the energy a therapist needs, and how many protons it takes to treat a tumor

Protons: m=1.67×1027kgm = 1.67 \times 10^{-27}\,\mathrm{kg}, q=e=1.60×1019Cq = e = 1.60 \times 10^{-19}\,\mathrm{C}, mc2=938MeVmc^2 = 938\,\mathrm{MeV}; c=3.00×108m/sc = 3.00 \times 10^{8}\,\mathrm{m}/\mathrm{s}. Cyclotron: B=1.5TB = 1.5\,\mathrm{T}, extraction radius Rmax=0.50mR_{\max} = 0.50\,\mathrm{m}, dee voltage amplitude U=50kVU = 50\,\mathrm{kV}.

Part I — One turn.

  1. Write the Lorentz force on a proton moving at v\vect v perpendicular to B\vect B, and show that its speed cannot change inside a dee.
  2. Derive the radius of its circle and its period; show that the period does not depend on the speed.
  3. Compute the cyclotron frequency fcf_c for these protons.
  4. Energy gained at each gap crossing; per turn.
  5. Why must the dee voltage alternate at exactly fcf_c? What happens if it is slightly off?
  6. The protons are injected at the center with negligible speed: radius after the first turn? After nn turns?

Part II — To the rim.

  1. Kinetic energy at extraction (in J and MeV), and the speed, as a fraction of cc.
  2. Number of turns needed and time spent in the machine.
  3. Total length of the spiral path (sum of the circumferences — approximate the sum by an integral).
  4. Does the result of question 7 depend on UU? On what does UU act?
  5. The vacuum inside must be good: estimate the distance a proton travels between collisions with residual gas molecules if the mean free path must exceed the spiral length; the mean free path scales as 1/p1/p and is about 70nm70\,\mathrm{nm} at atmospheric pressure. What pressure is needed?
  6. Compute γ1Ek/mc2\gamma - 1 \approx E_k/mc^2 at extraction and the resulting relative change of the period. Over the number of turns found above, by what fraction of a period has the proton slipped relative to the voltage? Comment.

Part III — 230 MeV. Proton therapy needs about 230MeV230\,\mathrm{MeV} to reach tumors 30cm30\,\mathrm{cm} deep.

  1. With the non-relativistic formula, what product BRBR would give 230MeV230\,\mathrm{MeV}? With B=2.0TB = 2.0\,\mathrm{T}, what radius?
  2. At that energy γ1=0.245\gamma - 1 = 0.245: what happens to the cyclotron period as the proton gains energy, and why does a fixed-frequency cyclotron fail?
  3. Two cures exist: vary the voltage frequency during the acceleration (synchrocyclotron), or make BB grow with radius so that B/γB/\gamma stays constant (isochronous cyclotron). Explain each in one sentence and name a drawback of the first.
  4. The relativistic momentum is p=γmvp = \gamma mv and the radius R=p/qBR = p/qB still holds: with γ=1.245\gamma = 1.245 and v=0.60cv = 0.60c, compute RR for B=2.0TB = 2.0\,\mathrm{T} and compare with question 13.
  5. The extracted beam passes between two plates 1.0m1.0\,\mathrm{m} long with E=1.0×105V/mE = 1.0 \times 10^{5}\,\mathrm{V}/\mathrm{m} to be steered by a small angle: estimate the deflection angle (non-relativistically, with pp from the previous question as mvmv replacement). Is electric steering practical here?
  6. The same steering with a magnet of 1.0T1.0\,\mathrm{T} over 1.0m1.0\,\mathrm{m}: angle? Conclude on why beam lines use magnets.

Part IV — The dose.

  1. A 230MeV230\,\mathrm{MeV} proton deposits most of its energy near the end of its range (the Bragg peak). If all of it is absorbed in a tumor of mass 0.50kg0.50\,\mathrm{kg}, how many protons deliver a dose of 2.0Gy2.0\,\mathrm{Gy} (1Gy=1J/kg1\,\mathrm{Gy} = 1\,\mathrm{J}/\mathrm{kg})?
  2. The beam current is 1.0nA1.0\,\mathrm{nA}: protons per second, and duration of the session.
  3. Power carried by the beam; compare with a light bulb.
  4. The protons are slowed to a stop inside the patient: what is the energy deposited per proton in the last centimeter, roughly, if a quarter of the energy is lost there? In electron-volts per micrometer?
  5. Each ionization in tissue costs about 30eV30\,\mathrm{eV}: how many ionizations does one proton cause along its whole path?
  6. Why is a proton beam preferable to X-rays for a deep tumor, in one sentence about where the energy goes?
  7. Summarize the chain from the Lorentz force to the dose: which equations fixed the energy, the radius, the number of turns and the treatment time.
Solution

Solution of Problem 17.1.

1. F=evBv\vect F = e\vect v\wedge\vect B \perp \vect v: no power, constant speed.

2. mv2/R=evBmv^2/R = evB: R=mv/eBR = mv/eB; T=2πR/v=2πm/eBT = 2\pi R/v = 2\pi m/eB, free of vv.

3. fc=eB/2πm=1.6×1019×1.5/(2π×1.67×1027)=23MHzf_c = eB/2\pi m = 1.6\times10^{-19} \times 1.5/(2\pi \times 1.67 \times10^{-27}) = 23\,\mathrm{MHz}.

4. eU=50keVeU = 50\,\mathrm{keV} per crossing, 100keV100\,\mathrm{keV} per turn.

5. The gap must be at the accelerating polarity each time the proton arrives, half a period apart; off frequency, the phase slips and the kicks eventually decelerate.

6. After one turn Ek=100keVE_k = 100\,\mathrm{keV}, v=2Ek/m=4.4×106m/sv = \sqrt{2E_k/m} = 4.4 \times 10^{6}\,\mathrm{m}/\mathrm{s}, R=mv/eB=3.0cmR = mv/eB = 3.0\,\mathrm{cm}; after nn turns Ek=n×100keVE_k = n \times 100\,\mathrm{keV} and R=3.0cmnR = 3.0\,\mathrm{cm}\sqrt n.

7. Ek=e2B2R2/2m=4.3×1012J=27MeVE_k = e^2B^2R^2/2m = 4.3 \times 10^{-12}\,\mathrm{J} = 27\,\mathrm{MeV}; v=2Ek/m=7.2×107m/s=0.24cv = \sqrt{2E_k/m} = 7.2 \times 10^{7}\,\mathrm{m}/\mathrm{s} = 0.24c.

8. 270270 turns; 270×43.5ns=12µs270 \times 43.5\,\mathrm{ns} = 12\,\text{µ}\mathrm{s}.

9. 2πRn=2π×0.030n2π×0.030×23×2703/2=560m\sum 2\pi R_n = 2\pi \times 0.030\sum\sqrt n \approx 2\pi \times 0.030 \times \tfrac23 \times 270^{3/2} = 560\,\mathrm{m}.

10. No: EkE_k depends on BB and RmaxR_{\max} only; UU sets the number of turns and the time.

11. Path 560m560\,\mathrm{m}: need λ1km\lambda \gtrsim 1\,\mathrm{km}, i.e. p105 Pa×7×108/103=7×106Pap \lesssim 10^5\ \mathrm{Pa} \times 7\times10^{-8}/10^3 = 7 \times 10^{-6}\,\mathrm{Pa} — a high vacuum.

12. γ1=27/938=0.029\gamma - 1 = 27/938 = 0.029; the period grows by 2.9%2.9\%. Over 270270 turns the slip accumulates to several periods — the last turns would be out of phase; in practice BB and the injection phase are trimmed, and 27MeV27\,\mathrm{MeV} is near the limit of the simple machine.

13. Ek=(eBR)2/2mE_k = (eBR)^2/2m: BR=2mEk/e=2×1.67×1027×3.68×1011/1.6×1019=2.2TmBR = \sqrt{2mE_k}/e = \sqrt{2 \times 1.67\times 10^{-27} \times 3.68\times10^{-11}}/1.6\times10^{-19} = 2.2\,\mathrm{T}\,\mathrm{m}; R=1.1mR = 1.1\,\mathrm{m} at 2T2\,\mathrm{T}.

14. T=2πγm/eBT = 2\pi\gamma m/eB grows by 24%24\% from center to rim: a fixed frequency falls hopelessly out of step.

15. Synchrocyclotron: lower the frequency as the bunch gains energy — only one bunch at a time, low average current. Isochronous: shape B(R)γB(R) \propto \gamma so the period stays constant — continuous beam, the modern choice.

16. p=γmv=1.245×1.67×1027×1.8×108=3.7×1019kgm/sp = \gamma mv = 1.245 \times 1.67\times10^{-27} \times 1.8\times10^8 = 3.7 \times 10^{-19}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}; R=p/eB=1.2mR = p/eB = 1.2\,\mathrm{m}, a little more than the non-relativistic 1.1m1.1\,\mathrm{m}.

17. Transverse impulse eEL/v=1.6×1019×105×1.0/1.8×108=8.9×1023kgm/seEL/v = 1.6\times10^{-19} \times 10^5 \times 1.0/1.8\times10^8 = 8.9 \times 10^{-23}\,\mathrm{kg}\,\mathrm{m}/\mathrm{s}; angle 8.9×1023/3.7×1019=2.4×104rad\approx 8.9\times10^{-23}/ 3.7\times10^{-19} = 2.4 \times 10^{-4}\,\mathrm{rad}: hopeless for steering.

18. Magnetic: R=p/eB=2.3mR = p/eB = 2.3\,\mathrm{m}, angle L/R=0.43rad25L/R = 0.43\,\mathrm{rad} \approx 25^\circ: a thousand times more — magnets it is.

19. Energy per proton 3.68×1011J3.68 \times 10^{-11}\,\mathrm{J}; dose 2.0×0.50=1.0J2.0 \times 0.50 = 1.0\,\mathrm{J}: 2.7×10102.7\times10^{10} protons.

20. 109/1.6×1019=6.2×109protons/s10^{-9}/1.6\times10^{-19} = 6.2 \times 10^{9}\,\mathrm{protons}/\mathrm{s}: about 4.4s4.4\,\mathrm{s}.

21. P=230MeV×6.2×109s1=0.23WP = 230\,\mathrm{MeV} \times 6.2 \times 10^{9}\,\mathrm{s}^{-1} = 0.23\,\mathrm{W}: a night light, delivered exactly where wanted.

22. 58MeV58\,\mathrm{MeV} in 1cm1\,\mathrm{cm}: about 6keV/µm6\,\mathrm{keV}/\text{µ}\mathrm{m} — dense ionization along the last track.

23. 230×106/308×106230\times10^6/30 \approx 8 \times 10^{6} ionizations per proton — the molecular damage that kills the cell.

24. X-rays deposit energy all along their path, most near the entrance; protons deposit most of theirs at the end of their range, sparing the tissue in front and leaving nothing behind.

25. R=mv/eBR = mv/eB and T=2πm/eBT = 2\pi m/eB (the circle), Ek=(eBR)2/2mE_k = (eBR)^2/2m (the energy), Ek/2eUE_k/2eU (the turns), dose == (protons ×\times energy)/mass and current =e×= e \times rate (the time).

Terms defined in this chapter

See all 393 terms in the glossary