# Charged Particle in a Magnetic Field

Michael Fowler 1/16/08

### Introduction

Classically, the force on a charged particle in electric and magnetic fields is given by the Lorentz force law:

$\stackrel{\to }{F}=q\left(\stackrel{\to }{E}+\frac{\stackrel{\to }{v}×\stackrel{\to }{B}}{c}\right)$

This velocity-dependent force is quite different from the conservative forces from potentials that we have dealt with so far, and the recipe for going from classical to quantum mechanics$—$replacing momenta with the appropriate derivative operators$—$has to be carried out with more care.  We begin by demonstrating how the Lorentz force law arises classically in the Lagrangian and Hamiltonian formulations.

### Laws of Classical Mechanics

Recall first (or look it up in Shankar, Chapter 2) that the Principle of Least Action leads to the Euler-Lagrange equations for the Lagrangian L:

The canonical momentum pi is defined by the equation

${p}_{i}=\frac{\partial L}{\partial {\stackrel{˙}{q}}_{i}}$

and the Hamiltonian is defined by performing a Legendre transformation of the Lagrangian:

$H\left({q}_{i},{p}_{i}\right)=\sum {p}_{i}{\stackrel{˙}{q}}_{i}-L\left({q}_{i},{\stackrel{˙}{q}}_{i}\right)$

It is straightforward to check that the equations of motion can be written:

These are known as Hamilton’s Equations.  Note that if the Hamiltonian is independent of a particular coordinate qi, the corresponding momentum pi remains constant.  (Such a coordinate is termed cyclic, because the most common example is an angular coordinate in a spherically symmetric Hamiltonian, where angular momentum remains constant.)

For the conservative forces we have been considering so far, L = T  VH = T + V, with T the kinetic energy, V the potential energy.

### Poisson Brackets

Any dynamical variable f in the system is some function of the qi’s and pi’s and (assuming it does not depend explicitly on time) its development is given by:

$\frac{d}{dt}f\left({q}_{i},{p}_{i}\right)=\frac{\partial f}{\partial {q}_{i}}{\stackrel{˙}{q}}_{i}+\frac{\partial f}{\partial {p}_{i}}{\stackrel{˙}{p}}_{i}=\frac{\partial f}{\partial {q}_{i}}\frac{\partial H}{\partial {p}_{i}}-\frac{\partial f}{\partial {p}_{i}}\frac{\partial H}{\partial {q}_{i}}=\left\{f,H\right\}.$

The curly brackets are called Poisson Brackets, and are defined for any dynamical variables as:

$\left\{A,B\right\}=\frac{\partial A}{\partial {q}_{i}}\frac{\partial B}{\partial {p}_{i}}-\frac{\partial A}{\partial {p}_{i}}\frac{\partial B}{\partial {q}_{i}}.$

We have shown from Hamilton’s equations that for any variable $\stackrel{˙}{f}=\left\{f,H\right\}.$

It is easy to check that for the coordinates and canonical momenta,

This was the classical mathematical structure that led Dirac to link up classical and quantum mechanics: he realized that the Poisson brackets were the classical version of the commutators, so a classical canonical momentum must correspond to the quantum differential operator in the corresponding coordinate.

### Particle in a Magnetic Field

The Lorentz force is velocity dependent, so cannot be just the gradient of some potential. Nevertheless, the classical particle path is still given by the Principle of Least Action.  The electric and magnetic fields can be written in terms of a scalar and a vector potential:

The right Lagrangian turns out to be:

$L=\frac{1}{2}m{\stackrel{\to }{v}}^{2}-q\phi +\frac{q}{c}\stackrel{\to }{v}.\stackrel{\to }{A}.$

(Note: if you’re familiar with Relativity, the interaction term here looks less arbitrary: the relativistic version would have the relativistically invariant $\left(q/c\right)\int {A}^{\mu }d{x}_{\mu }$ added to the action integral, where the four-potential ${A}_{\mu }=\left(\stackrel{\to }{A},\phi \right)$ and $d{x}_{\mu }=\left(d{x}_{1},d{x}_{2},d{x}_{3},cdt\right)$.  This is the simplest possible invariant interaction between the electromagnetic field and the particle’s four-velocity.  Then in the nonrelativistic limit, $\left(q/c\right)\int {A}^{\mu }d{x}_{\mu }$ just becomes $\int q\left(\stackrel{\to }{v}.\stackrel{\to }{A}/c-\phi \right)\text{\hspace{0.17em}}dt.$ )

The derivation of the Lorentz force from this Lagrangian is given by Shankar on page 84.  We give the (equivalent) derivation from the Hamilton equations below.

Note that for zero vector potential, the Lagrangian has the usual T  V  form.

For this one-particle problem, the general coordinates qi are just the Cartesian co-ordinates ${x}_{i}=\left({x}_{1},{x}_{2},{x}_{3}\right)$, the position of the particle, and the ${\stackrel{˙}{q}}_{i}$ are the three components ${\stackrel{˙}{x}}_{i}={v}_{i}$ of the particle’s velocity.

The important new point is that the canonical momentum

${p}_{i}=\frac{\partial L}{\partial {\stackrel{˙}{q}}_{i}}=\frac{\partial L}{\partial {\stackrel{˙}{x}}_{i}}=m{v}_{i}+\frac{q}{c}{A}_{i}$

is no longer mass $×$ velocity$—$there is an extra term!

The Hamiltonian is

$\begin{array}{c}H\left({q}_{i},{p}_{i}\right)=\sum {p}_{i}{\stackrel{˙}{q}}_{i}-L\left({q}_{i},{\stackrel{˙}{q}}_{i}\right)\\ =\sum \left(m{v}_{i}+\frac{q}{c}{A}_{i}\right){v}_{i}-\frac{1}{2}m{\stackrel{\to }{v}}^{2}+q\phi -\frac{q}{c}\stackrel{\to }{v}.\stackrel{\to }{A}\\ =\frac{1}{2}m{\stackrel{\to }{v}}^{2}+q\phi \end{array}$

Reassuringly, the Hamiltonian just has the familiar form of kinetic energy plus potential energy.  However, to get Hamilton’s equations of motion, the Hamiltonian has to be expressed solely in terms of the coordinates and canonical momenta.  That is,

$H=\frac{{\left(\stackrel{\to }{p}-q\stackrel{\to }{A}\left(\stackrel{\to }{x},t\right)/c\right)}^{2}}{2m}+q\phi \text{\hspace{0.17em}}\left(\stackrel{\to }{x},t\right)$

where we have noted explicitly that the potentials mean those at the position $\stackrel{\to }{x}$ of the particle at time t.

Let us now consider Hamilton’s equations

It is easy to see how the first equation comes out, bearing in mind that

${p}_{i}=m{v}_{i}+\frac{q}{c}{A}_{i}=m{\stackrel{˙}{x}}_{i}+\frac{q}{c}{A}_{i}.$

The second equation yields the Lorentz force law, but is a little more tricky.  The first point to bear in mind is that dp/dt is not the acceleration, the A term also varies in time, and in a quite complicated way, since it is the field at a point moving with the particle.  That is,

${\stackrel{˙}{p}}_{i}=m{\stackrel{¨}{x}}_{i}+\frac{q}{c}{\stackrel{˙}{A}}_{\text{ }i}=m{\stackrel{¨}{x}}_{i}+\frac{q}{c}\left(\frac{\partial {A}_{\text{ }i}}{\partial t}+{v}_{j}{\nabla }_{j}A{\text{ }}_{i}\right).$

The right-hand side of the second Hamilton equation ${\stackrel{˙}{p}}_{i}=-\frac{\partial H}{\partial {x}_{i}}$ is

$\begin{array}{c}-\frac{\partial H}{\partial {x}_{i}}=\frac{\left(\stackrel{\to }{p}-q\stackrel{\to }{A}\left(\stackrel{\to }{x},t\right)/c\right)}{m}.\frac{q}{c}.\frac{\partial \stackrel{\to }{A}}{\partial {x}_{i}}-q\frac{\partial \phi \left(\stackrel{\to }{x},t\right)}{\partial {x}_{i}}\\ =\frac{q}{c}{v}_{j}{\nabla }_{i}{A}_{j}-q{\nabla }_{i}\phi .\end{array}$

Putting the two sides together, the Hamilton equation reads:

$m{\stackrel{¨}{x}}_{i}=-\frac{q}{c}\left(\frac{\partial {A}_{\text{\hspace{0.17em}}i}}{\partial t}+{v}_{j}{\nabla }_{j}{A}_{i}\right)+\frac{q}{c}{v}_{j}{\nabla }_{i}{A}_{j}-q{\nabla }_{i}\phi \text{\hspace{0.17em}}.$

and the expressions for the electric and magnetic fields in terms of the potentials, the Lorentz force law emerges:

$m\stackrel{¨}{\stackrel{\to }{x}}=q\left(\stackrel{\to }{E}+\frac{\stackrel{\to }{v}×\stackrel{\to }{B}}{c}\right)$

### Quantum Mechanics of a Particle in a Magnetic Field

We make the standard substitution:

This leads to the novel situation that the velocities in different directions do not commute.  From

$m{v}_{i}=-i\hslash {\nabla }_{i}-q{A}_{i}/c$

it is easy to check that

$\left[{v}_{x},{v}_{y}\right]=\frac{iq\hslash }{{m}^{2}c}B$

To actually solve Schrödinger’s equation for an electron confined to a plane in a uniform perpendicular magnetic field, it is convenient to use the Landau gauge,

$\stackrel{\to }{A}\left(x,y,z\right)=\left(-By,0,0\right)$

giving a constant field B in the z direction.  The equation is

$H\psi \left(x,y\right)=\left[\frac{1}{2m}{\left({p}_{x}+qBy/c\right)}^{2}+\frac{{p}_{y}^{2}}{2m}\right]\psi \left(x,y\right)=E\psi \left(x,y\right).$

Note that x does not appear in this Hamiltonian, so it is a cyclic coordinate, and px is conserved. In other words, this H commutes with px, so H and  px have a common set of eigenstates.  We know the eigenstates of px are just the plane waves ${e}^{i{p}_{x}x/\hslash }$, so the common eigenstates must have the form:

$\psi \left(x,y\right)={e}^{i{p}_{x}x/\hslash }\chi \left(y\right).$

Operating on this wavefunction with the Hamiltonian, the operator px appearing in H simply gives its eigenvalue. That is, the  px  in H just becomes a number!  Therefore, writing ${p}_{y}=-i\hslash d/dy,$ the y-component $\chi \left(y\right)$ of the wavefunction satisfies:

$-\frac{{\hslash }^{2}}{2m}\frac{{d}^{2}}{d{y}^{2}}\chi \left(y\right)+\frac{1}{2}m{\left(\frac{qB}{mc}\right)}^{2}{\left(y-{y}_{0}\right)}^{2}\chi \left(y\right)=E\chi \left(y\right)$

where

${y}_{0}=-c{p}_{x}/qB.$

We now see that the conserved canonical momentum px  in the x-direction is actually the coordinate of the center of a simple harmonic oscillator potential in the y-direction!  This simple harmonic oscillator has frequency ω = |q|B/mc, so the allowed values of energy for a particle in a plane in a perpendicular magnetic field are:

$E=\left(n+\frac{1}{2}\right)\hslash \omega =\left(n+\frac{1}{2}\right)\hslash |q|B/mc.$

The frequency is of course the cyclotron frequency$—$that of the classical electron in a circular orbit in the field (given by $m{v}^{2}/r=qvB/c,\text{ }\omega =v/r=qB/mc$ ) .

Let us confine our attention to states corresponding to the lowest oscillator state, $E=\frac{1}{2}\hslash \omega$.  How many such states are there?  Consider a square of conductor, area $A={L}_{x}×{L}_{y}$, and, for simplicity, take periodic boundary conditions.  The center of the oscillator wave function y0 must lie between 0 and Ly.  But remember that ${y}_{0}=-c{p}_{x}/qB$, and with periodic boundary conditions   This means that y0 takes a series of evenly-spaced discrete values, separated by

$\Delta {y}_{0}=ch/qB{L}_{x}.$

So the total number of states $N={L}_{y}/\Delta {y}_{0},$

$N=\frac{{L}_{x}{L}_{y}}{\left(\frac{hc}{qB}\right)}=A.\frac{B}{{\Phi }_{0}},$

where ${\Phi }_{0}$ is called the “flux quantum”.  So the total number of states in the lowest energy level $E=\frac{1}{2}\hslash \omega$ (usually referred to as the lowest Landau level, or LLL) is exactly equal to the total number of flux quanta making up the field B penetrating the area A.

It is instructive to find y0 from a purely classical analysis.

These equations integrate trivially to give:

$\begin{array}{l}m\stackrel{˙}{x}=\frac{qB}{c}\left(y-{y}_{0}\right),\\ m\stackrel{˙}{y}=-\frac{qB}{c}\left(x-{x}_{0}\right).\end{array}$

Here (x0, y0) are the coordinates of the center of the classical circular motion (the velocity vector $\stackrel{˙}{\stackrel{\to }{r}}=\left(\stackrel{˙}{x},\stackrel{˙}{y}\right)$ is always perpendicular to $\left(\stackrel{\to }{r}-{\stackrel{\to }{r}}_{0}\right)$ ) , and ${\stackrel{\to }{r}}_{0}$ is given by

$\begin{array}{l}{y}_{0}=y-cm{v}_{x}/qB=-c{p}_{x}/qB\\ {x}_{0}=x+cm{v}_{y}/qB=x+c{p}_{y}/qB.\end{array}$

(Recall that we are using the gauge $\stackrel{\to }{A}\left(x,y,z\right)=\left(-By,0,0\right)$, and ${p}_{x}=\frac{\partial L}{\partial \stackrel{˙}{x}}=m{v}_{x}+\frac{q}{c}{A}_{x}$, etc.)

Just as y0 is a conserved quantity, so is x0: it commutes with the Hamiltonian since

$\left[x+c{p}_{y}/qB,{p}_{x}+qBy/c\right]=0.$

However, x0 and y0 do not commute with each other:

$\left[{x}_{0},{y}_{0}\right]=-i\hslash c/qB.$

This is why, when we chose a gauge in which y0 was sharply defined, x0 was spread over the sample. If we attempt to localize the point (x0, y0) as well as possible, it is fuzzed out over an area essentially that occupied by one flux quantum.  The natural length scale of the problem is therefore the magnetic length defined by

$l=\sqrt{\frac{\hslash c}{qB}}.$

References: the classical mechanics at the beginning is similar to Shankar’s presentation, the quantum mechanics is closer to that in Landau.