Skip to content

Problem Irodov 3.395

In a right-handed xyzxyz frame a constant uniform magnetic field of induction BB points along the positive zz axis, and a uniform electric field along the yy axis varies with time as E=Emcos⁡ωtE = E_m\cos\omega t, where ω=qB/m\omega = qB/m. A particle of specific charge q/mq/m starts from rest at the origin OO at the moment t=0t = 0. Treating the motion as non-relativistic, find the law of motion x(t)x(t), y(t)y(t) of the particle. What is the approximate shape of its trajectory?

The frame with z toward the viewer: the oscillating E along y (drawn for ω t > 0) and the constant B along +z; the particle starts at rest at the origin O.

The frame with zz toward the viewer: the oscillating E\vect E along yy (drawn for cos⁡ωt>0\cos\omega t > 0) and the constant B\vect B along +z+z; the particle starts at rest at the origin OO.

Given
  • B∥+z\vect B \parallel +z
  • Ey=Emcos⁡ωtE_y = E_m\cos\omega t
  • ω=qB/m\omega = qB/m
  • q/mq/m
  • at rest at OO, t=0t = 0
Find
x(t)x(t), y(t)y(t) and the approximate shape of the trajectory

Solution

Step 1The equations of motion

The force q(E+v×B)q(\vect E + \vect v\times\vect B) has no zz component for motion in the xyxy plane, so the particle stays in that plane. With B=B k\vect B = B\,\vect k, v×B=(vyB,−vxB,0)\vect v\times\vect B = (v_y B, -v_x B, 0), and writing aa for the amplitude of the electric acceleration:

dvxdt=ωvy,dvydt=acos⁡ωt−ωvx,a=qEmm(1) \begin{aligned} \frac{dv_x}{dt} &= \omega v_y, \\ \frac{dv_y}{dt} &= a\cos\omega t - \omega v_x, \\ a &= \frac{qE_m}{m} \end{aligned} \tag{1}
Step 2One complex equation

Combine the components into w=vx+ivyw = v_x + i v_y. Adding the first equation of (1) to ii times the second:

dwdt=ωvy+i acos⁡ωt−iωvx=−iω(vx+ivy)+i acos⁡ωt,dwdt+iωw=i acos⁡ωt\begin{align*} \frac{dw}{dt} &= \omega v_y + i\,a\cos\omega t - i\omega v_x = -i\omega\left(v_x + i v_y\right) + i\,a\cos\omega t, \\ \frac{dw}{dt} + i\omega w &= i\,a\cos\omega t \tag{2}\end{align*}
Step 3The velocity

Multiplying (2) by eiωte^{i\omega t} makes the left side a derivative, and the cosine splits into exponentials:

ddt(w eiωt)=i acos⁡ωt  eiωt=ia2(1+e2iωt)\frac{d}{dt}\left(w\,e^{i\omega t}\right) = i\,a\cos\omega t\;e^{i\omega t} = \frac{i a}{2}\left(1 + e^{2i\omega t}\right)

Integrating from rest, w=0w = 0 at t=0t = 0, and multiplying back by e−iωte^{-i\omega t}:

w eiωt=ia2 t+a4ω(e2iωt−1),w=iat2(cos⁡ωt−isin⁡ωt)+ia2ω sin⁡ωt\begin{aligned} w\,e^{i\omega t} &= \frac{i a}{2}\,t + \frac{a}{4\omega}\left(e^{2i\omega t} - 1\right), \\ w &= \frac{i a t}{2}\left(\cos\omega t - i\sin\omega t\right) + \frac{i a}{2\omega}\,\sin\omega t \end{aligned}

The real and imaginary parts are the two velocity components:

vx=at2 sin⁡ωt,vy=at2 cos⁡ωt+a2ω sin⁡ωt(3) \begin{aligned} v_x &= \frac{a t}{2}\,\sin\omega t, \\ v_y &= \frac{a t}{2}\,\cos\omega t + \frac{a}{2\omega}\,\sin\omega t \end{aligned} \tag{3}
Step 4The law of motion

Integrating (3) from OO, with the integral ∫tsin⁡ωt dt=−tcos⁡ωt/ω+sin⁡ωt/ω2\int t\sin\omega t\,dt = -t\cos\omega t/\omega + \sin\omega t/\omega^2 and its cosine companion, both found by parts:

x=a2[−tcos⁡ωtω+sin⁡ωtω2],y=a2[tsin⁡ωtω+cos⁡ωt−1ω2]+a2ω2(1−cos⁡ωt)\begin{aligned} x &= \frac{a}{2}\left[-\frac{t\cos\omega t}{\omega} + \frac{\sin\omega t}{\omega^2}\right], \\ y &= \frac{a}{2}\left[\frac{t\sin\omega t}{\omega} + \frac{\cos\omega t - 1}{\omega^2}\right] + \frac{a}{2\omega^2}\left(1 - \cos\omega t\right) \end{aligned}

In yy the two non-growing terms cancel, leaving

x=a2ω2(sin⁡ωt−ωtcos⁡ωt),y=a2ω tsin⁡ωt,a=qEmm\begin{gathered} x = \frac{a}{2\omega^2}\left(\sin\omega t - \omega t\cos\omega t\right), \\ y = \frac{a}{2\omega}\,t\sin\omega t, \\ a = \frac{qE_m}{m} \end{gathered}
(4)
Step 5The shape of the trajectory

Once ωt≫1\omega t \gg 1 the terms growing with tt dominate, and the motion is a circling whose radius rr grows in proportion to the time:

x≈−at2ωcos⁡ωt,y≈at2ωsin⁡ωt,r≈at2ωx \approx -\frac{a t}{2\omega}\cos\omega t, \qquad y \approx \frac{a t}{2\omega}\sin\omega t, \qquad r \approx \frac{a t}{2\omega}

Each turn adds πa/ω2\pi a/\omega^2 to the radius: the trajectory is an unwinding spiral, close to an Archimedean one, turning clockwise (for q>0q > 0) about the origin.

The law of motion plotted over two and a half turns: from rest at O the particle sets off along y and spirals outward clockwise, the gap between turns staying the same.

The law of motion plotted over two and a half turns: from rest at OO the particle sets off along yy and spirals outward clockwise, the gap between turns staying the same.

Check

Early times: for ωt≪1\omega t \ll 1, sin⁡ωt−ωtcos⁡ωt≈(ωt)3/3\sin\omega t - \omega t\cos\omega t \approx (\omega t)^3/3, so x≈aωt3/6x \approx a\omega t^3/6 and y≈at2/2y \approx a t^2/2; the particle starts off along yy as if the electric force qEmqE_m acted alone, which it does before any speed builds up for BB to act on. Units: aa is an acceleration, so at/ωa t/\omega and a/ω2a/\omega^2 are lengths. Speed: from (3), vx2+vy2≈(at/2)2v_x^2 + v_y^2 \approx (a t/2)^2 at late times, and rω≈at/2r\omega \approx a t/2 is the same speed, consistent with circling at the frequency ω\omega.

Remark

The field Emcos⁡ωtE_m\cos\omega t along a line is the sum of two fields of strength Em/2E_m/2 rotating in opposite senses at the frequency ω\omega. The one that rotates with the particle stays in step with its velocity and pushes it forward on every turn, so the speed grows steadily at the rate qEm/(2m)qE_m/(2m); the other rotates against it and its push averages out. This is cyclotron resonance, the principle of the cyclotron and of the resonant heating of plasmas; off resonance the speed would only beat, never grow.