Skip to content

Landau Levels

A uniform magnetic field quantizes the transverse kinetic energy of a charged particle into equally spaced Landau levels. For a spinless particle of mass mm and signed charge qq in

B=Bz^,B>0,\mathbf B = B\hat{\mathbf z}, \qquad B>0,

the two-dimensional orbital energies are

En=ℏωc(n+12),n=0,1,2,…,E_n = \hbar\omega_c \left(n+\frac12\right), \qquad n=0,1,2,\ldots,

where

ωc=∣q∣Bm.\omega_c = \frac{\lvert q\rvert B}{m}.

Each level has a macroscopic guiding-center degeneracy. This card collects the energy, length, state-counting, gauge, and spin-extension formulas. The full oscillator reduction is at Landau Levels, and the careful finite-area counting argument is at Degeneracy of Landau Levels.

TaskFormula
Cyclotron frequencyωc=∣q∣B/m\omega_c=\lvert q\rvert B/m
Magnetic lengthℓB=ℏ/(∣q∣B)\ell_B=\sqrt{\hbar/(\lvert q\rvert B)}
Spinless two-dimensional energyEn=ℏωc(n+1/2)E_n=\hbar\omega_c(n+1/2)
Three-dimensional energyEn,kz=ℏωc(n+1/2)+ℏ2kz2/(2m)E_{n,k_z}=\hbar\omega_c(n+1/2)+\hbar^2k_z^2/(2m)
Degeneracy per areaNΦ/A=1/(2πℓB2)=∣q∣B/hN_\Phi/A=1/(2\pi\ell_B^2)=\lvert q\rvert B/h
Charge-qq flux quantumΦ0=h/∣q∣\Phi_0=h/\lvert q\rvert
Cyclotron radius squared⟨η2⟩n=(2n+1)ℓB2\langle\eta^2\rangle_n=(2n+1)\ell_B^2
Filling factorν=N/NΦ=2πℓB2n2D\nu=N/N_\Phi=2\pi\ell_B^2n_{2\mathrm D}

The cyclotron frequency and magnetic length are

ωc=∣q∣Bm,ℓB=ℏ∣q∣B.\omega_c = \frac{\lvert q\rvert B}{m}, \qquad \ell_B = \sqrt{ \frac{\hbar}{\lvert q\rvert B} }.

They obey the useful identity

mωcℓB2=ℏ.m\omega_c\ell_B^2 = \hbar.

The sign of qq controls the direction of cyclotron motion and the orientation of magnetic translation algebra, but not the positive energy spacing ℏωc\hbar\omega_c. The magnetic length is not the classical orbit radius. It is the basic quantum width from which the level-dependent orbit scale is built.

For charge magnitude ∣q∣\lvert q\rvert, define

Φ0=h∣q∣.\Phi_0 = \frac{h}{\lvert q\rvert}.

For an electron this is h/eh/e. It should not be confused with the superconducting flux quantum h/(2e)h/(2e), which reflects charge-2e2e pairs.

With minimal coupling,

H=π22m,π=p−qA,H = \frac{\boldsymbol\pi^2}{2m}, \qquad \boldsymbol\pi = \mathbf p-q\mathbf A,

where ∇×A=B\nabla\times\mathbf A=\mathbf B. In the transverse plane,

[πx,πy]=iℏqB.[\pi_x,\pi_y] = i\hbar qB.

Let

s=sgn⁡(qB)s = \operatorname{sgn}(qB)

and define

a=ℓB2ℏ(πx+isπy),a†=ℓB2ℏ(πx−isπy).a = \frac{\ell_B}{\sqrt2\hbar} \left( \pi_x+is\pi_y \right), \qquad a^\dagger = \frac{\ell_B}{\sqrt2\hbar} \left( \pi_x-is\pi_y \right).

Then

[a,a†]=1[a,a^\dagger]=1

and

H⊥=ℏωc(a†a+12).H_\perp = \hbar\omega_c \left( a^\dagger a+\frac12 \right).

Therefore

En=ℏωc(n+12).E_n = \hbar\omega_c \left(n+\frac12\right).

The oscillator number nn labels cyclotron excitation. It does not label the many distinct guiding-center states within the same Landau level.

Choose

A=Bx y^.\mathbf A = Bx\,\hat{\mathbf y}.

Then pyp_y commutes with the Hamiltonian, and a convenient basis is

ψn,ky(x,y)=eikyyLyφn(x−x0),\psi_{n,k_y}(x,y) = \frac{e^{ik_yy}}{\sqrt{L_y}} \varphi_n(x-x_0),

where

x0=ℏkyqB.x_0 = \frac{\hbar k_y}{qB}.

The sign of x0x_0 therefore depends on the signed charge. The normalized oscillator function is

φn(ξ)=Hn(ξ/ℓB)π1/42nn!ℓBexp⁡ ⁣(−ξ22ℓB2).\varphi_n(\xi) = \frac{ H_n(\xi/\ell_B) }{ \pi^{1/4}\sqrt{2^n n!\ell_B} } \exp\!\left( -\frac{\xi^2}{2\ell_B^2} \right).

Substitution reduces the Hamiltonian to

H⊥=px22m+12mωc2(x−x0)2.H_\perp = \frac{p_x^2}{2m} + \frac12m\omega_c^2(x-x_0)^2.

The energy is independent of kyk_y. Different kyk_y values shift the guiding center and generate the degeneracy. The canonical momentum label kyk_y and the center coordinate x0x_0 depend on gauge; the level energy and total state count do not.

The symmetric gauge

A=12B×r\mathbf A = \frac12\mathbf B\times\mathbf r

organizes the same degenerate subspace by angular-momentum-like labels. Use Landau Gauge and Symmetric Gauge when translating between those bases.

The position separates into guiding-center and cyclotron pieces:

r=R+η,\mathbf r = \mathbf R+\boldsymbol\eta,

with

Rx=x+πyqB,Ry=y−πxqB,ηx=−πyqB,ηy=πxqB.\begin{aligned} R_x &= x+\frac{\pi_y}{qB}, & R_y &= y-\frac{\pi_x}{qB}, \\ \eta_x &= -\frac{\pi_y}{qB}, & \eta_y &= \frac{\pi_x}{qB}. \end{aligned}

These sectors commute with one another:

[Ri,ηj]=0.[R_i,\eta_j]=0.

Within each sector,

[Rx,Ry]=−isℓB2,[ηx,ηy]=isℓB2.[R_x,R_y] = -is\ell_B^2, \qquad [\eta_x,\eta_y] = is\ell_B^2.

The Hamiltonian depends only on η\boldsymbol\eta; the independent guiding center R\mathbf R labels states at fixed energy. This is the gauge-independent content behind the Landau-gauge parameter kyk_y.

The root-mean-square cyclotron coordinate is

⟨η2⟩n=2n+1 ℓB.\sqrt{\langle\eta^2\rangle_n} = \sqrt{2n+1}\,\ell_B.

Thus the lowest Landau level has a finite transverse quantum extent even though n=0n=0.

Consider a rectangle Lx×LyL_x\times L_y and impose periodic boundary conditions along yy. Then

ky=2πjLy,j∈Z.k_y = \frac{2\pi j}{L_y}, \qquad j\in\mathbb Z.

Adjacent guiding centers are separated in magnitude by

∣Δx0∣=2πℓB2Ly.\lvert\Delta x_0\rvert = \frac{2\pi\ell_B^2}{L_y}.

The number whose centers fit across the bulk width LxL_x is

NΦ=LxLy2πℓB2=∣q∣BAh=BAΦ0.N_\Phi = \frac{L_xL_y}{2\pi\ell_B^2} = \frac{\lvert q\rvert BA}{h} = \frac{BA}{\Phi_0}.

Hence the orbital degeneracy per area is

NΦA=12πℓB2=∣q∣Bh.\frac{N_\Phi}{A} = \frac{1}{2\pi\ell_B^2} = \frac{\lvert q\rvert B}{h}.

This counts one spinless orbital per charge-qq flux quantum. On a torus, consistent magnetic boundary conditions quantize the total flux and make NΦN_\Phi an integer. In a finite sample with physical edges, edge states and boundary details modify the literal center-counting picture. The bulk density of states remains the central result.

Do not multiply this count by two unless an unresolved twofold spin degeneracy is actually present.

For free motion parallel to the field,

H=H⊥+pz22m.H = H_\perp+\frac{p_z^2}{2m}.

The spectrum is

En,kz=ℏωc(n+12)+ℏ2kz22m.E_{n,k_z} = \hbar\omega_c \left(n+\frac12\right) + \frac{\hbar^2k_z^2}{2m}.

The magnetic field discretizes transverse kinetic energy but leaves a continuous longitudinal dispersion on infinite space. A finite length or an additional potential can quantize kzk_z separately.

The orbital formulas above are spinless. If the particle has magnetic moment

μs=gq2mS,\boldsymbol\mu_s = g\frac{q}{2m}\mathbf S,

the Zeeman Hamiltonian is

HZ=−μs⋅B.H_Z = -\boldsymbol\mu_s\cdot\mathbf B.

For Sz∣s,ms⟩=ℏms∣s,ms⟩S_z\lvert s,m_s\rangle=\hbar m_s\lvert s,m_s\rangle,

ΔEZ=−gqB2mℏms.\Delta E_Z = -g\frac{qB}{2m}\hbar m_s.

The combined ideal energy is therefore

En,ms=ℏωc(n+12)−gqB2mℏms.E_{n,m_s} = \hbar\omega_c \left(n+\frac12\right) - g\frac{qB}{2m}\hbar m_s.

For a free electron, q=−eq=-e and ΔEZ=gμBBms\Delta E_Z=g\mu_BBm_s, with μB=eℏ/(2me)\mu_B=e\hbar/(2m_e). In solids, the effective mass controlling ωc\omega_c and the effective gg factor can differ strongly from their vacuum values. State those parameters before comparing orbital and Zeeman gaps.

Relativistic particles, graphene-like Dirac bands, and particles with strong spin-orbit coupling have different Landau spectra; they are not obtained by blindly appending this Zeeman term.

For NN spinless particles in area AA, define the number density n2D=N/An_{2\mathrm D}=N/A. The orbital filling factor is

ν=NNΦ=2πℓB2n2D=n2Dh∣q∣B.\nu = \frac{N}{N_\Phi} = 2\pi\ell_B^2n_{2\mathrm D} = \frac{n_{2\mathrm D}h}{\lvert q\rvert B}.

In a clean noninteracting picture, integer ν\nu means an integer number of spin-resolved orbital Landau levels is filled. Quantized Hall plateaus require the response, disorder localization, and edge structure developed at the Integer Quantum Hall Effect canonical page; flux counting alone is not a derivation of the plateau physics.

  1. Record the signed charge qq but use ∣q∣\lvert q\rvert in ωc\omega_c and ℓB\ell_B.
  2. Decide whether the model is two- or three-dimensional and whether spin is included.
  3. Compute ωc\omega_c, ℓB\ell_B, and the spinless orbital energies.
  4. Use A/(2πℓB2)A/(2\pi\ell_B^2) only after stating the sample geometry and boundary approximation.
  5. Keep gauge-dependent labels such as kyk_y separate from gauge-invariant observables.
  6. Add Zeeman, confinement, disorder, lattice, or interaction terms only as explicit extensions of the ideal model.
  • Using qB/mqB/m as a positive cyclotron frequency for a negatively charged particle.
  • Losing the signed qq in the Landau-gauge center x0=ℏky/(qB)x_0=\hbar k_y/(qB).
  • Calling kyk_y an ordinary gauge-invariant mechanical momentum.
  • Confusing ℓB\ell_B with the level-dependent radius 2n+1 ℓB\sqrt{2n+1}\,\ell_B.
  • Counting A/(2πℓB2)A/(2\pi\ell_B^2) as an exact finite-edge degeneracy without specifying boundary conditions.
  • Adding an automatic factor of two for spin after using a spinless formula.
  • Forgetting the free kzk_z term in three dimensions.
  • Treating different-looking Landau- and symmetric-gauge wavefunctions as different spectra.
  • Assuming a uniform magnetic field confines the guiding center in the bulk.
  • Applying the quadratic-band result to relativistic or Dirac particles.
  • L. D. Landau, “Diamagnetismus der Metalle,” Zeitschrift fur Physik 64, 629-637 (1930).
  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Pergamon, 1977.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • R. E. Prange and S. M. Girvin, eds., The Quantum Hall Effect, 2nd ed., Springer, 1990.
  1. Show that the kinetic-momentum ladder operators satisfy [a,a†]=1[a,a^\dagger]=1.
Solution

Using [πx,πy]=iℏqB=iℏs∣q∣B[\pi_x,\pi_y]=i\hbar qB=i\hbar s\lvert q\rvert B,

[a,a†]=ℓB22ℏ2[πx+isπy,πx−isπy]=ℓB22ℏ2(−2is[πx,πy])=ℓB2∣q∣Bℏ=1.\begin{aligned} [a,a^\dagger] &= \frac{\ell_B^2}{2\hbar^2} [\pi_x+is\pi_y,\pi_x-is\pi_y] \\ &= \frac{\ell_B^2}{2\hbar^2} \left(-2is[\pi_x,\pi_y]\right) \\ &= \frac{\ell_B^2\lvert q\rvert B}{\hbar} =1. \end{aligned}

The final equality uses ℓB2=ℏ/(∣q∣B)\ell_B^2=\hbar/(\lvert q\rvert B).

  1. Derive the bulk degeneracy of one spinless Landau level in a rectangle.
Solution

Periodic boundary conditions along yy give

ky=2πjLy.k_y=\frac{2\pi j}{L_y}.

Because x0=ℏky/(qB)x_0=\hbar k_y/(qB), adjacent centers have magnitude spacing

∣Δx0∣=2πℏ∣q∣BLy=2πℓB2Ly.\lvert\Delta x_0\rvert = \frac{2\pi\hbar}{\lvert q\rvert BL_y} = \frac{2\pi\ell_B^2}{L_y}.

The number of bulk centers in width LxL_x is therefore

NΦ=Lx∣Δx0∣=LxLy2πℓB2=∣q∣BAh.N_\Phi = \frac{L_x}{\lvert\Delta x_0\rvert} = \frac{L_xL_y}{2\pi\ell_B^2} = \frac{\lvert q\rvert BA}{h}.
  1. A three-dimensional spinless particle is in the state n=2n=2 with longitudinal wave number kzk_z. Find its energy and transverse root-mean-square cyclotron radius.
Solution

The energy is

E2,kz=52ℏωc+ℏ2kz22m.E_{2,k_z} = \frac{5}{2}\hbar\omega_c + \frac{\hbar^2k_z^2}{2m}.

The cyclotron-coordinate expectation value is

⟨η2⟩2=(2⋅2+1)ℓB2=5ℓB2,\langle\eta^2\rangle_2 = (2\cdot2+1)\ell_B^2 = 5\ell_B^2,

so

⟨η2⟩2=5 ℓB.\sqrt{\langle\eta^2\rangle_2} = \sqrt5\,\ell_B.
  1. A spinless two-dimensional gas has number density n2Dn_{2\mathrm D}. Find the field at which the orbital filling factor is ν=1\nu=1.
Solution

Use

ν=n2Dh∣q∣B.\nu = \frac{n_{2\mathrm D}h}{\lvert q\rvert B}.

Setting ν=1\nu=1 gives

B=n2Dh∣q∣.B = \frac{n_{2\mathrm D}h}{\lvert q\rvert}.

This is a flux-counting result. By itself it does not establish a quantized Hall plateau, which also depends on the spectrum, occupied states, response, and localization physics.