6.1: Charged Particle in a Magnetic Field
( \newcommand{\kernel}{\mathrm{null}\,}\)
Classically, the force on a charged particle in electric and magnetic fields is given by the Lorentz force law:
→F=q(→E+→v×→Bc)
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 that the Principle of Least Action leads to the Euler-Lagrange equations for the Lagrangian L:
ddt(∂L(qi,˙qi)∂˙qi)−∂L(qi,˙qi)∂qi=0
with qi and ˙qi being coordinates and velocities. The canonical momentum pi is defined by the equation
pi=∂L∂˙qi
and the Hamiltonian is defined by performing a Legendre transformation of the Lagrangian:
H(qi,pi)=∑i(pi˙qi−L(qi,˙qi))
It is straightforward to check that the equations of motion can be written:
˙qi=∂H∂pi,˙pi=−∂H∂qi
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−V
and
H=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:
ddtf(qi,pi)=∂f∂qi˙qi+∂f∂pi˙pi=∂f∂qi∂H∂pi−∂f∂pi∂H∂qi={f,H}.
The curly brackets are called Poisson Brackets, and are defined for any dynamical variables as:
{A,B}=∂A∂qi∂B∂pi−∂A∂pi∂B∂qi.
We have shown from Hamilton’s equations that for any variable ˙f={f,H}.
It is easy to check that for the coordinates and canonical momenta,
qi,qj=0=pi,pj,qi,pj=δij.
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.
Poisson brackets are the classical version of the commutators
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:
→B=→∇×→A
→E=−→∇φ−1c∂→A∂t.
The right Lagrangian turns out to be:
L=12m→v2−qφ+qc→v⋅→A.
Relativity Effects
If you’re familiar with Relativity, the interaction term here looks less arbitrary: the relativistic version would have the relativistically invariant (q/c)∫Aμdxμ added to the action integral, where the four-potential Aμ=(→A,φ) and dxμ=(dx1,dx2,dx3,cdt). This is the simplest possible invariant interaction between the electromagnetic field and the particle’s four-velocity. Then in the nonrelativistic limit, (q/c)∫Aμdxμ just becomes ∫q(→v⋅→A/c−φ) dt.
The derivation of the Lorentz force from the Hamilton equations is straightforward.
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 xi=(x1,x2,x3), the position of the particle, and the ˙qi are the three components ˙xi=vi of the particle’s velocity.
The important new point is that the canonical momentum pi=∂L∂˙qi=∂L∂˙xi=mvi+qcAi
is no longer mass × velocity—there is an extra term!
The Hamiltonian is
H(qi,pi)=∑pi˙qi−L(qi,˙qi)=∑(mvi+qcAi)vi−12m→v2+qφ−qc→v⋅→A=12m→v2+qφ
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=(→p−q→A(→x,t)/c)22m+qφ (→x,t)
where we have noted explicitly that the potentials mean those at the position →x of the particle at time t.
Let us now consider Hamilton’s equations
˙xi=∂H∂pi,˙pi=−∂H∂xi
It is easy to see how the first equation comes out, bearing in mind that
pi=mvi+qcAi=m˙xi+qcAi.
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,
˙pi=m¨xi+qc˙Ai=m¨xi+qc(∂Ai∂t+vj∇jAi).
The right-hand side of the second Hamilton equation ˙pi=−∂H∂xi is −∂H∂xi=(→p−q→A(→x,t)/c)m⋅qc⋅∂→A∂xi−q∂φ(→x,t)∂xi=qcvj∇iAj−q∇iφ.
Putting the two sides together, the Hamilton equation reads:
m¨xi=−qc(∂Ai∂t+vj∇jAi)+qcvj∇iAj−q∇iφ .
Using →v×(→∇×→A)=→∇(→v⋅→A)−(→v⋅→∇)→A, →B=→∇×→A, and the expressions for the electric and magnetic fields in terms of the potentials, the Lorentz force law emerges: m¨→x=q(→E+→v×→Bc)
Quantum Mechanics of a Particle in a Magnetic Field
We make the standard substitution:
→p=−iℏ→∇,sothat[xi,pj]=iℏδijasusual:butnowpi≠mvi.
This leads to the novel situation that the velocities in different directions do not commute. From mvi=−iℏ∇i−qAi/c
it is easy to check that [vx,vy]=iqℏm2cB
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, →A(x,y,z)=(−By,0,0)
giving a constant field B in the z direction. The equation is Hψ(x,y)=[12m(px+qBy/c)2+p2y2m]ψ(x,y)=Eψ(x,y).
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 eipxx/ℏ, so the common eigenstates must have the form: ψ(x,y)=eipxx/ℏχ(y).
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 py=−iℏd/dy, the y-component χ(y) of the wavefunction satisfies: −ℏ22md2dy2χ(y)+12m(qBmc)2(y−y0)2χ(y)=Eχ(y)
where y0=−cpx/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=(n+12)ℏω=(n+12)ℏ|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 mv2/r=qvB/c,ω=v/r=qB/mc ) .
Let us confine our attention to states corresponding to the lowest oscillator state, E=12ℏω. How many such states are there? Consider a square of conductor, area A=Lx×Ly, and, for simplicity, take periodic boundary conditions. The center of the oscillator wave function y0 must lie between 0 and Ly. But remember that y0=−cpx/qB, and with periodic boundary conditions eipxLx/ℏ=1, so px=2nπℏ/Lx=nh/Lx. This means that y0 takes a series of evenly-spaced discrete values, separated by Δy0=ch/qBLx.
So the total number of states N=Ly/Δy0, N=LxLy(hcqB)=A⋅BΦ0,
where Φ0 is called the “flux quantum”. So the total number of states in the lowest energy level E=12ℏω (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.
Writing m˙→v=qc→v×→B in components,
m¨x=qBc˙y,m¨y=−qBc˙x.
These equations integrate trivially to give:
m˙x=qBc(y−y0),m˙y=−qBc(x−x0).
Here (x0,y0) are the coordinates of the center of the classical circular motion (the velocity vector ˙→r=(˙x,˙y) is always perpendicular to (→r−→r0) ) , and →r0 is given by
y0=y−cmvx/qB=−cpx/qBx0=x+cmvy/qB=x+cpy/qB.
(Recall that we are using the gauge →A(x,y,z)=(−By,0,0), and px=∂L∂˙x=mvx+qcAx, etc.)
Just as y0 is a conserved quantity, so is x0: it commutes with the Hamiltonian since
[x+cpy/qB,px+qBy/c]=0.
However, x0 and y0 do not commute with each other: [x0,y0]=−iℏ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=√ℏcqB.
References: the classical mechanics at the beginning is similar to Shankar’s presentation, the quantum mechanics is closer to that in Landau.