Alan asked me an interesting question about a proton moving in a magnetic field. I hope he will forgive me for paraphrasing:

Consider a charged particle moving parallel to the ground in a uniform gravitational field. The gravitational force is normal to the ground. The charge immediately enters a region of uniform magnetic field orientated parallel to the ground, perpendicular to *both* the gravitational field and the charge’s initial velocity vector.

Provided the charge and the magnetic field have the correct sign, the Lorentz force will steer the charge *up* through the gravitational field. The charge will move in a semicircle, and leave the magnetic field, but at a higher point than that at which it entered.

In this way, the charge appears to have gained gravitational potential energy from nowhere. Has the conservation of energy been violated?

In the preamble above, I made an unproven assertion. I claimed that, simultaneously, there is a non-zero gravitational field *and* the charge moves in a semicircle. It is true that a charge moving in a plane normal to a uniform magnetic field will adopt a circular path. How do I know this is true in the presence of an additional gravitational field?

So rich and interesting is this simple physical problem that it warrants some discussion. In this relatively lengthy post, we will look in detail at the trajectories an electric charge may assume in uniform crossed and fields. We will first solve the equations of motion that describe the evolution of the charge’s spatial coordinates with time. Then we will explore the variety of trajectories available.

The quality of the animations has suffered a little in their previews, so please click on them to see them in their full glory in a new tab. All hail GeoGebra!

The particle has mass and charge .

Let the gravitational field be parallel to the axis. Let the magnetic field be parallel to the axis.

The particle starts at the origin:

The charge enters the magnetic field horizontally, so at time its velocity is

that is, parallel to the axis. In vector form, the equation of motion reads

which in column vector notation is

This gives the system of three differential equations

The first equation is the simplest to solve. The particle does not accelerate in the direction, because both the gravitational and Lorentz forces act in the plane.

This means the velocity in the direction will not vary with time. Making the velocity non-zero will not add any interesting new physics, so we will take it to be zero. This means

The second and third differential equations are *coupled* – each contains derivatives of both and . We will eliminate this dependence by integrating the second equation to give an expression for in terms of . We will then substitute this expression into the third equation. This will yield a differential equation involving and its own derivatives only.

The second equation reads

With inexplicable foresight we will introduce the quantity

so

Integrating up,

We now substitute this expression into the third equation:

The most general solution to this equation is

where and are constants. The constant term on the right is the so-called *particular integral*, an educated guess. The trigonometric terms form the *complementary function*, which can be added on to the solution ‘for free’ since they satisfy .

That sets

and means

Hence

We then substitute this into our equation for ,

and integrate, imposing initial conditions:

So we have obtained the charge’s trajectory. The position vector of the charge as a function of time is

Let’s now explore the family of trajectories by varying the featured parameters. We’ll first consider the most interesting effect, that of changing the initial velocity .

Suppose the charge is released from rest, so . In this case, the trajectory is described as a *cycloid*. A cycloid is the curve traced out by a point fixed on the circumference of a rolling circle. The animation below shows the motion of the charge for in the plane and the force acting on it.

There are several comments to be made.

This first and most obvious: this is a beautiful result.

The second: the charge *will not fall. *As the charge begins to gain velocity downward, it starts to experience a Lorentz force in the positive direction. As soon as the charge has a non-zero rightward velocity, there is a component of the Lorentz force acting *upwards*, decreasing the rate at which the charge is falling. Eventually, the charge has been steered round to such an extent that it is travelling horizontally. At this point, the Lorentz force acts exclusively in the positive direction, dragging the charge upwards. The charge ultimately comes to rest at the same height as where it started. The pattern repeats.

The form of this trajectory is identical for all values of and ; only its *scale* varies. Increasing either or increases its size. Increasing the magnitude of or decreases it. Changing the sign of causes the charge to move in the negative direction instead.

The rate at which the charge completes a ‘cycle’ is dictated by the angular frequency , which we defined as

This result is quite intuitive. The greater the magnitude of , the greater the strength of the coupling between the charge and the magnetic field. Hence the rate at which the particle’s path deviates from a straight vertical line is increased. The greater the value of , the more the charge resists the effect of the magnetic field and the slower is deviation from the vertical line.

Let us look at the more general case.

We’ll gradually increase from zero to a positive value:

The trajectory is a called a *trochoid*, of the the cycloid is a special case. Note the positions at which the charge meets the axis is invariant.

The trajectory passes through several interesting states as is increased.

For small , the cusps of the cycloid soften into smoothed hills. The name for this particular curve is a curtate trochoid.

As the initial velocity is increased further, you may notice that there exists a value of at which the trajectory is a straight line. From the expression for , this occurs when

When this equation is satisfied, the force on the charge vanishes.

Next, the trajectory forms another curtate trochoid, eventually becoming another cycloid *above* the axis. This cycloid has a phase shift relative to that for which .

This occurs when there exist solutions to the equation

which can be shown to occur when

Above this velocity, the charge starts to draw loops. Increasing the launch velocity distends the loops. This looping pattern is called a prolate trochoid.

Increasing the initial velocity further yet causes the charge to take a path which allows it to *escape* the magnetic field, assuming for . Some algebra shows that this happens if

This equation cannot, as far as I know, be solved analytically. A numerical solution is given below. It is quite interesting that this fraction

has appeared again. In fact, we can predict which trochoid the charge’s motion will assume based purely on the value of this fraction.

We’ll define the quantity as

Combing all of the previously established results, we know that for

- , the trajectory is a cycloid below the axis
- , the trajectory is a curtate trochoid below the axis
- , the trajectory is a straight line along the axis
- , the trajectory is a curtate trochoid above the axis
- , the trajectory is a cycloid above the axis
- , the trajectory is a prolate trochoid in the first quadrant
- , the particle can escape into the region

where the number is computed numerically. It is remarkable that the value of this simple linear equation in the initial velocity tells us which member of the trochoid family the trajectory will assume. Here is another animation showing the trajectory evolution along with the value of in the top right corner.

The conditions on can be easily inverted to give conditions for in relation to and instead.

As a sanity check, we can also watch how the particle’s trajectory changes as the gravitational field is switched off:

The trajectory tends to a circle, as we expect. For small , the charge moves in an almost circular path which drifts slowly to the right.

So, applying a gravitational field *downwards* causes the charge to migrate *sideways*!

On a final note, the dependence of the force on the charge with time goes like

for . Every trajectory above is produced by a force vector which rotates at a fixed frequency, like a phasor.

It can be shown from the expression for that the total energy of the charge,

is equal to

that is, mechanical energy is conserved. This shows that the particle cannot gain energy for free.

It’s likely that I’ve only scratched the surface of the physics encoded in these equations. In my next post, I’ll be looking more at the mechanics of a charge in a magnetic field.

In the line, \displaystyle \frac{d\dot{y}}{dt} = \omega \frac{dz}{dz}

should it actually be

\displaystyle \frac{d\dot{y}}{dt} = \omega \frac{dz}{dt}

?

Or am I missing something?

You’re absolutely right: thank you for drawing this to my attention! The error has been corrected.