Landau quantization

Landau quantization in quantum mechanics is the quantization of the cyclotron orbits of charged particles in magnetic fields. As a result, the charged particles can only occupy orbits with discrete energy values, called Landau levels. The Landau levels are degenerate, with the number of electrons per level directly proportional to the strength of the applied magnetic field. Landau quantization is directly responsible for oscillations in electronic properties of materials as a function of the applied magnetic field. It is named after the Jewish Soviet physicist Lev Landau.

Derivation

Consider a two-dimensional system of non-interacting particles with charge q and spin S confined to an area A = L_xL_y in the x-y plane.

Apply a uniform magnetic field $\mathbf {B} ={\begin{pmatrix}0\\0\\B\end{pmatrix}}$ along the z-axis. In CGS units, the Hamiltonian of this system is

${\hat {H}}={\frac {1}{2m}}({\hat {\mathbf {p} }}-q{\hat {\mathbf {A} }}/c)^{2}.$

Here, p̂ is the canonical momentum operator and  is the electromagnetic vector potential, which is related to the magnetic field by

$\mathbf {B} =\mathbf {\nabla } \times {\hat {\mathbf {A} }}.\,$

There is some gauge freedom in the choice of vector potential for a given magnetic field. The Hamiltonian is gauge invariant, which means that adding the gradient of a scalar field to  changes the overall phase of the wave function by an amount corresponding to the scalar field. But physical properties are not influenced by the specific choice of gauge. For simplicity in calculation, choose the Landau gauge, which is

${\hat {\mathbf {A} }}={\begin{pmatrix}0\\Bx\\0\end{pmatrix}}.$

where B=|B| and x̂ is the x component of the position operator.

In this gauge, the Hamiltonian is

${\hat {H}}={\frac {{\hat {p}}_{x}^{2}}{2m}}+{\frac {1}{2m}}\left({\hat {p}}_{y}-{\frac {qB{\hat {x}}}{c}}\right)^{2}.$

The operator ${\hat {p}}_{y}$ commutes with this Hamiltonian, since the operator ŷ is absent by the choice of gauge. Thus the operator ${\hat {p}}_{y}$ can be replaced by its eigenvalue ħk_y .

The Hamiltonian can also be written more simply by noting that the cyclotron frequency is ω_c = qB/mc, giving

${\hat {H}}={\frac {{\hat {p}}_{x}^{2}}{2m}}+{\frac {1}{2}}m\omega _{c}^{2}\left({\hat {x}}-{\frac {\hbar k_{y}}{m\omega _{c}}}\right)^{2}.$

This is exactly the Hamiltonian for the quantum harmonic oscillator, except with the minimum of the potential shifted in coordinate space by x₀ = ħk_y/mω_c .

To find the energies, note that translating the harmonic oscillator potential does not affect the energies. The energies of this system are thus identical to those of the standard quantum harmonic oscillator,

$E_{n}=\hbar \omega _{c}\left(n+{\frac {1}{2}}\right),\quad n\geq 0~.$

The energy does not depend on the quantum number k_y, so there will be degeneracies.

For the wave functions, recall that ${\hat {p}}_{y}$ commutes with the Hamiltonian. Then the wave function factors into a product of momentum eigenstates in the y direction and harmonic oscillator eigenstates $|\phi _{n}\rangle$ shifted by an amount x₀ in the x direction:

$\Psi (x,y)=e^{ik_{y}y}\phi _{n}(x-x_{0})~.$

In sum, the state of the electron is characterized by two quantum numbers, n and k_y .

Landau levels

Each set of wave functions with the same value of n is called a Landau level. Effects of Landau levels are only observed when the mean thermal energy is smaller than the energy level separation, kT ≪ ħω_c, meaning low temperatures and strong magnetic fields.

Each Landau level is degenerate because of the second quantum number k_y, which can take the values

$k_{y}={\frac {2\pi N}{L_{y}}}$,

where N is an integer. The allowed values of N are further restricted by the condition that the center of force of the oscillator, x₀, must physically lie within the system, 0 ≤ x₀ < L_x. This gives the following range for N,

$0\leq N<{\frac {m\omega _{c}L_{x}L_{y}}{2\pi \hbar }}.$

For particles with charge q = Ze, the upper bound on N can be simply written as a ratio of fluxes,

${\frac {ZBL_{x}L_{y}}{(hc/e)}}=Z{\frac {\Phi }{\Phi _{0}}},$

where Φ₀ = h/2e is the fundamental quantum of flux and Φ = BA is the flux through the system (with area A = L_xL_y).

Thus, for particles with spin S, the maximum number D of particles per Landau level is

$D=Z(2S+1){\frac {\Phi }{\Phi _{0}}}~.$

The above gives only a rough idea of the effects of finite-size geometry. Strictly speaking, using the standard solution of the harmonic oscillator is only valid for systems unbounded in the x-direction (infinite strips). If the size L_x is finite, boundary conditions in that direction give rise to non-standard quantization conditions on the magnetic field, involving (in principle) both solutions to the Hermite equation. The filling of these levels with many electrons is still^[1] an active area of research.

In general, Landau levels are observed in electronic systems, where Z=1 and S=1/2. As the magnetic field is increased, more and more electrons can fit into a given Landau level. The occupation of the highest Landau level ranges from completely full to entirely empty, leading to oscillations in various electronic properties (see de Haas–van Alphen effect and Shubnikov–de Haas effect).

If Zeeman splitting is included, each Landau level splits into a pair, one for spin up electrons and the other for spin down electrons. Then the occupation of each spin Landau level is just the ratio of fluxes D = Φ/Φ₀. Zeeman splitting has a significant effect on the Landau levels because their energy scales are the same, 2μ_BB = ħω . However, the Fermi energy and ground state energy stay roughly the same in a system with many filled levels, since pairs of split energy levels cancel each other out when summed.

Discussion

This derivation treats x and y as slightly asymmetric. However, by the symmetry of the system, there is no physical quantity which distinguishes these coordinates. The same result could have been obtained with an appropriate interchange of x and y.

Moreover, the above derivation assumed an electron confined in the z-direction, which is a relevant experimental situation — found in two-dimensional electron gases, for instance. Still, this assumption is not essential for the results. If electrons are free to move along the z direction, the wave function acquires an additional multiplicative term exp(ik_zz); the energy corresponding to this free motion, (ħ k_z)²/(2m), is added to the E discussed. This term then fills in the separation in energy of the different Landau levels, blurring the effect of the quantization. Nevertheless, the motion in the x-y-plane, perpendicular to the magnetic field, is still quantized.

Landau Levels in Symmetric Gauge

The symmetric gauge refers to the choice

${\hat {{\mathbf {A}}}}={\frac {1}{2}}{\begin{pmatrix}-By\\Bx\\0\end{pmatrix}}$

In terms of dimensionless lengths and energies, the Hamiltonian can be expressed as

${\hat {H}}={\frac {1}{2}}\left[\left(-i{\frac {\partial }{\partial x}}-{\frac {y}{2}}\right)^{2}+\left(-i{\frac {\partial }{\partial y}}+{\frac {x}{2}}\right)^{2}\right]$

The correct units can be restored by introducing factors of $q,c,\hbar ,{\mathbf {B}}$ and $m$

Consider operators

${\hat {a}}={\frac {1}{{\sqrt {2}}}}\left[\left({\frac {x}{2}}+{\frac {\partial }{\partial x}}\right)-i\left({\frac {y}{2}}+{\frac {\partial }{\partial y}}\right)\right]$

${\hat {a}}^{{\dagger }}={\frac {1}{{\sqrt {2}}}}\left[\left({\frac {x}{2}}-{\frac {\partial }{\partial x}}\right)+i\left({\frac {y}{2}}-{\frac {\partial }{\partial y}}\right)\right]$

${\hat {b}}={\frac {1}{{\sqrt {2}}}}\left[\left({\frac {x}{2}}+{\frac {\partial }{\partial x}}\right)+i\left({\frac {y}{2}}+{\frac {\partial }{\partial y}}\right)\right]$

${\hat {b}}^{{\dagger }}={\frac {1}{{\sqrt {2}}}}\left[\left({\frac {x}{2}}-{\frac {\partial }{\partial x}}\right)-i\left({\frac {y}{2}}-{\frac {\partial }{\partial y}}\right)\right]$

These operators follow certain commutation relations

$[{\hat {a}},{\hat {a}}^{{\dagger }}]=[{\hat {b}},{\hat {b}}^{{\dagger }}]=1$.

In terms of above operators the Hamiltonian can be written as

${\hat {H}}={\hat {a}}^{{\dagger }}{\hat {a}}+{\frac {1}{2}}$

Landau Level index $n$ is the eigenvalue of ${\hat {a}}^{{\dagger }}{\hat {a}}$

The z component of angular momentum is

${\hat {L}}_{z}=-i\hbar {\frac {\partial }{\partial \theta }}=-\hbar ({\hat {b}}^{{\dagger }}{\hat {b}}-{\hat {a}}^{{\dagger }}{\hat {a}})$

Exploiting the property $[{\hat {H}},{\hat {L}}_{z}]=0$ we chose eigenfunctions which diagonalize ${\hat {H}}$ and ${\hat {L}}_{z}$, The eigenvalue of ${\hat {L}}_{z}$ is denoted by $-m\hbar$, where it is clear that $m \ge -n$ in the $n$th Landau level. However, it may be arbitrarily large, which is necessary to obtain the infinite degeneracy (or finite degeneracy per unit area) exhibited by the system.

The application of $\hat{b}^{\dagger}$ increases $m$ by one unit while preserving $n$, whereas $\hat{a}^{\dagger}$ application simultaneously increase $n$ and decreases $m$ by one unit. The analogy to quantum harmonic oscillator provides solutions

${\hat {H}}|n,m\rangle =E_{n}|n,m\rangle$

$E_{n}=\left(n+{\frac {1}{2}}\right)$

$|n,m\rangle ={\frac {({\hat {b}}^{{\dagger }})^{{m+n}}}{{\sqrt {(m+n)!}}}}{\frac {({\hat {a}}^{{\dagger }})^{{n}}}{{\sqrt {n!}}}}|0,0\rangle$

Each Landau level has degenerate orbitals labeled by the quantum numbers k_y and $m$ in the Landau and symmetric gauges respectively. The degeneracy per unit area is the same in each Landau level.

One may verify that the above states correspond to choosing wavefunctions proportional to

$\psi_{n,m}(x, y) = \left( \frac{\partial}{\partial w} - \frac{\bar{w}}{4} \right)^n w^{n + m} e^{-|w|^2 / 4}$

where $w = x + i y$.

In particular, the lowest Landau level $n=0$ consists of arbitrary analytic functions multiplying a Gaussian, $\psi(x,y) = f(w) e^{-|w|^2/4}$.

Effects of Gauge Transformation

${\vec {A}}\to {\vec {A}}'={\vec {A}}+{\vec {\nabla }}\lambda ({\vec {x}})$

The definition for kinematical momenta is

${\hat {\pi }}={\hat {{\mathbf {p}}}}-q{\hat {{\mathbf {A}}}}/c$

where ${\hat {\mathbf {p} }}$ is the canonical momenta. The Hamiltonian is a gauge invariant so $\langle {\hat {\pi }}\rangle$ and $\langle {\hat {x}}\rangle$ will remain invariant under gauge transformations but $\langle {\hat {{\mathbf {p}}}}\rangle$ will depend upon gauge. For observing the effect of gauge transformation on the quantum state of the particle, consider the state with A and A' as Vector Potential, with states $|\alpha \rangle$ and $|\alpha '\rangle$.

As $\langle {\hat {x}}\rangle$ and $\langle {\hat {\pi }}\rangle$ is invariant under the gauge transformation we get

$\langle \alpha |{\hat {x}}|\alpha \rangle =\langle \alpha '|{\hat {x}}|\alpha '\rangle$

$\langle \alpha |{\hat {\pi }}|\alpha \rangle =\langle \alpha '|{\hat {\pi '}}|\alpha '\rangle$

$\langle \alpha |\alpha \rangle =\langle \alpha '|\alpha '\rangle$

Consider an operator ${\mathcal {G}}$ such that $|\alpha '\rangle ={\mathcal {G}}|\alpha \rangle$

from above relation we deduce that

${\mathcal {G}}^{{\dagger }}{\hat {x}}{\mathcal {G}}={\hat {x}}$

${\mathcal {G}}^{{\dagger }}\left({\hat {p}}-{\frac {e{\hat {A}}}{c}}-{\frac {e{\vec {\nabla }}\lambda (x)}{c}}\right){\mathcal {G}}={\hat {p}}-{\frac {e{\hat {A}}}{c}}$

${\mathcal {G}}^{{\dagger }}{\mathcal {G}}=1$

from this we conclude

$\mathcal{G} = \exp\left(\frac{ie\lambda(\vec{x})}{\hbar c}\right)$

See also

• Barkhausen effect
• de Haas–van Alphen effect
• Shubnikov–de Haas effect
• Quantum Hall effect
• Laughlin wavefunction
• Coulomb potential between two current loops embedded in a magnetic field

Further reading

• Landau, L. D.; and Lifschitz, E. M.; (1977). Quantum Mechanics: Non-relativistic Theory. Course of Theoretical Physics. Vol. 3 (3rd ed. London: Pergamon Press). ISBN 0750635398 .

1. Mikhailov, S. A. (2001). "A new approach to the ground state of quantum Hall systems. Basic principles". Physica B: Condensed Matter. 299: 6. arXiv:cond-mat/0008227. Bibcode:2001PhyB..299....6M. doi:10.1016/S0921-4526(00)00769-9.