Skip to content

Problem Irodov 3.392

In a right-handed xyzxyz frame a uniform electric field of strength E\vect E points along the positive yy axis and a uniform magnetic field of induction B\vect B along the positive zz axis. A particle of specific charge q/mq/m is at rest at the origin OO at the moment t=0t = 0. Treating the motion as non-relativistic, find:

  1. the law of motion x(t)x(t), y(t)y(t) of the particle, and the shape of its trajectory;

  2. the length of the path between two successive points at which the particle’s velocity vanishes;

  3. the mean value of the xx component of the particle’s velocity (the drift velocity).

The frame with z toward the viewer: E along +y and B along +z, both uniform; the particle starts at rest at the origin O.

The frame with zz toward the viewer: E\vect E along +y+y and B\vect B along +z+z, both uniform; the particle starts at rest at the origin OO.

Given
  • E∥+y\vect E \parallel +y
  • B∥+z\vect B \parallel +z
  • q/mq/m
  • at rest at OO when t=0t = 0
Find
(a) x(t)x(t), y(t)y(t) and the trajectory’s shape · (b) the path between successive stops · (c) the drift velocity ⟨vx⟩\langle v_x\rangle

Solution

Step 1The equations of motion

The force is q(E+v×B)q(\vect E + \vect v\times\vect B). Neither field has a zz component of force on a particle moving in the xyxy plane, and the particle starts at rest, so it stays in that plane. With B=B k\vect B = B\,\vect k the magnetic part has the components

v×B=(vyB,  −vxB,  0)\vect v\times\vect B = \left(v_y B,\; -v_x B,\; 0\right)

Dividing by mm and writing ω\omega for the cyclotron frequency:

dvxdt=ωvy,dvydt=qEm−ωvx,ω=qBm(1) \begin{aligned} \frac{dv_x}{dt} &= \omega v_y, \\ \frac{dv_y}{dt} &= \frac{qE}{m} - \omega v_x, \\ \omega &= \frac{qB}{m} \end{aligned} \tag{1}
Step 2The velocity components

Differentiating the first equation of (1) and substituting the second eliminates vyv_y:

d2vxdt2=ω dvydt=ω2(EB−vx)\frac{d^2v_x}{dt^2} = \omega\,\frac{dv_y}{dt} = \omega^2\left(\frac{E}{B} - v_x\right)

So vx−E/Bv_x - E/B oscillates at the frequency ω\omega. The start from rest gives vx(0)=0v_x(0) = 0 and, through the first equation of (1), dvx/dt=0dv_x/dt = 0 at t=0t = 0; the same equation then returns vyv_y:

vx=EB (1−cos⁡ωt),vy=1ω dvxdt=EB sin⁡ωt(2) \begin{aligned} v_x &= \frac{E}{B}\,(1 - \cos\omega t), \\ v_y &= \frac{1}{\omega}\,\frac{dv_x}{dt} = \frac{E}{B}\,\sin\omega t \end{aligned} \tag{2}
Step 3(a) The law of motion and the trajectory

Integrating (2) from OO at t=0t = 0:

x=∫0tEB (1−cos⁡ωt′) dt′=EBω (ωt−sin⁡ωt),y=∫0tEB sin⁡ωt′ dt′=EBω (1−cos⁡ωt)\begin{aligned} x &= \int_0^t \frac{E}{B}\,(1 - \cos\omega t')\,dt' = \frac{E}{B\omega}\,(\omega t - \sin\omega t), \\ y &= \int_0^t \frac{E}{B}\,\sin\omega t'\,dt' = \frac{E}{B\omega}\,(1 - \cos\omega t) \end{aligned}

The common factor is a length; with ω\omega from (1),

a=EBω=mEqB2a = \frac{E}{B\omega} = \frac{mE}{qB^2}

which gives

x=a (ωt−sin⁡ωt),y=a (1−cos⁡ωt),a=mEqB2,ω=qBm\begin{gathered} x = a\,(\omega t - \sin\omega t), \quad y = a\,(1 - \cos\omega t), \\ a = \frac{mE}{qB^2}, \quad \omega = \frac{qB}{m} \end{gathered}
(3)

These are the equations of a cycloid: the particle moves as a point on the rim of a circle of radius aa that rolls without slipping along the xx axis, its centre advancing at the speed aω=E/Ba\omega = E/B.

The trajectory from rest at O: a cycloid of height 2a with cusps on the x axis. The dashed rolling circle of radius a turns clockwise as its centre moves along x.

The trajectory from rest at OO: a cycloid of height 2a2a with cusps on the xx axis. The dashed rolling circle of radius aa turns clockwise as its centre moves along xx.

Step 4(b) The path between two stops

The speed follows from (2), using 1−cos⁡ωt=2sin⁡2(ωt/2)1 - \cos\omega t = 2\sin^2(\omega t/2):

v2=vx2+vy2=E2B2(2−2cos⁡ωt)=4E2B2 sin⁡2ωt2v^2 = v_x^2 + v_y^2 = \frac{E^2}{B^2}\left(2 - 2\cos\omega t\right) = \frac{4E^2}{B^2}\,\sin^2\frac{\omega t}{2}

So v=(2E/B) ∣sin⁡(ωt/2)∣v = (2E/B)\,\lvert\sin(\omega t/2)\rvert, which vanishes at the cusps, ωt=0,2π,4π,…\omega t = 0, 2\pi, 4\pi, \ldots Integrating the speed over one arch, from t=0t = 0 to t=2π/ωt = 2\pi/\omega:

s=∫02π/ω2EB sin⁡ωt2 dt=2EB[−2ωcos⁡ωt2]02π/ω=2EB⋅4ωs = \int_0^{2\pi/\omega} \frac{2E}{B}\,\sin\frac{\omega t}{2}\,dt = \frac{2E}{B}\left[-\frac{2}{\omega}\cos\frac{\omega t}{2}\right]_0^{2\pi/\omega} = \frac{2E}{B}\cdot\frac{4}{\omega}
s=8mEqB2=8as = \frac{8mE}{qB^2} = 8a
(4)
Step 5(c) The drift velocity

Over one period T=2π/ωT = 2\pi/\omega the motion repeats itself shifted along xx by one arch, whose length along the axis is x(T)=2πax(T) = 2\pi a from the answer of (a). The mean of vxv_x is that shift over the period:

⟨vx⟩=x(T)T=2πa2π/ω=aω\langle v_x\rangle = \frac{x(T)}{T} = \frac{2\pi a}{2\pi/\omega} = a\omega
⟨vx⟩=EB\langle v_x\rangle = \frac{E}{B}
(5)
Check

Units: mE/(qB2)mE/(qB^2) is kg⋅(V/m)/(C⋅T2)\text{kg}\cdot(\text{V}/\text{m})/(\text{C}\cdot\text{T}^2); with V⋅C=J\text{V}\cdot\text{C} = \text{J} and T=kg/(C⋅s)\text{T} = \text{kg}/(\text{C}\cdot\text{s}) this is J⋅s2/(kg⋅m)=m\text{J}\cdot\text{s}^2/(\text{kg}\cdot\text{m}) = \text{m}, a length, and E/BE/B is (V/m)/(V⋅s/m2)=m/s(\text{V}/\text{m})/(\text{V}\cdot\text{s}/\text{m}^2) = \text{m}/\text{s}. The top speed, at ωt=π\omega t = \pi, is 2E/B2E/B, twice the drift, as for a point on the top of a rolling wheel. For a negative charge both aa and ω\omega change sign: y=a(1−cos⁡ωt)y = a(1 - \cos\omega t) becomes negative, so the arches hang below the xx axis, while x=a(ωt−sin⁡ωt)x = a(\omega t - \sin\omega t) stays positive, so the drift is still along +x+x, the direction of E×B\vect E\times\vect B.

Remark

In a frame moving along xx at u=E/Bu = E/B the electric field disappears, since E+u×B=E j−E j=0\vect E + \vect u\times\vect B = E\,\vect j - E\,\vect j = 0. There the particle starts with velocity −u-u in a pure magnetic field and simply circles with radius mu/(qB)=amu/(qB) = a; carried along at uu, that circle traces the cycloid. The drift E/BE/B does not depend on qq or mm, so in crossed fields electrons and ions of a plasma drift together and carry no net current along E×B\vect E\times\vect B.