Landau levels explained

In quantum mechanics, the energies of cyclotron orbits of charged particles in a uniform magnetic field are quantized to discrete values, thus known as Landau levels. These levels are degenerate, with the number of electrons per level directly proportional to the strength of the applied magnetic field. It is named after the Soviet physicist Lev Landau.^[1]

Landau quantization contributes towards magnetic susceptibility of metals, known as Landau diamagnetism. Under strong magnetic fields, Landau quantization leads to oscillations in electronic properties of materials as a function of the applied magnetic field known as the De Haas–Van Alphen and Shubnikov–de Haas effects.

Landau quantization is a key ingredient in explanation of the integer quantum Hall effect.

Derivation

Consider a system of non-interacting particles with charge and spin confined to an area in the plane. Apply a uniform magnetic field

B = \begin{pmatrix} 0 \\ 0 \\B\end{pmatrix} along the -axis. In SI units, the Hamiltonian of this system (here, the effects of spin are neglected) is$$\backslash hat\; =\; \backslash frac\; \backslash left(\backslash hat\; -\; q\backslash hat\backslash right)^2.$$Here, $\backslash hat$ is the canonical momentum operator and $\backslash hat$ is the operator for the electromagnetic vector potential $\backslash mathbf$ (in position space $\backslash hat\; =\backslash mathbf$).

The vector potential is related to the magnetic field by

B = \nabla x 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.

In the Landau gauge

From the possible solutions for A, a gauge fixing introduced by Lev Landau is often used for charged particles in a constant magnetic field.^[2]

When

B = \begin{pmatrix} 0 \ 0 \ B \end{pmatrix} then A = \begin{pmatrix} 0 \ B ⋅ x \ 0 \end{pmatrix} is a possible solution^[3] in the Landau gauge (not to be mixed up with the Landau R _{\xi} gauge).

In this gauge, the Hamiltonian is$$\backslash hat\; =\; \backslash frac\; +\; \backslash frac\; \backslash left(\backslash hat\_y\; -\; qB\backslash hat\backslash right)^2\; +\; \backslash frac.$$The operator

\hat{p} _{y} commutes with this Hamiltonian, since the operator \hat y is absent for this choice of gauge. Thus the operator \hat{p} _{y} can be replaced by its eigenvalue \hbar k _{y} . Since \hat{z} does not appear in the Hamiltonian and only the z-momentum appears in the kinetic energy, this motion along the z-direction is a free motion.

The Hamiltonian can also be written more simply by noting that the cyclotron frequency is

\omega _{c=qB/m} , giving$$\backslash hat\; =\; \backslash frac\; +\; \backslash frac\; m\; \backslash omega\_^2\; \backslash left(\backslash hat\; -\; \backslash frac\; \backslash right)^2\; +\; \backslash frac.$$This is exactly the Hamiltonian for the quantum harmonic oscillator, except with the minimum of the potential shifted in coordinate space by x _{0=\hbar} k _{y/m\omega} _{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,^[4] $$E\_n=\backslash hbar\backslash omega\_\backslash left(n+\backslash frac\backslash right)\; +\; \backslash frac,\backslash quad\; n\backslash geq\; 0.$$The energy does not depend on the quantum number

k _{y} , so there will be a finite number of degeneracies (If the particle is placed in an unconfined space, this degeneracy will correspond to a continuous sequence of p _{y} ). The value of p _{z} is continuous if the particle is unconfined in the z-direction and discrete if the particle is bounded in the z-direction also. Each set of wave functions with the same value of n is called a Landau level.

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 _{0} in the x direction:$$\backslash Psi(x,y,z)\; =\; e^\; \backslash phi\_n(x-x\_0)$$where k _{z} = p _{z} / \hbar . In sum, the state of the electron is characterized by the quantum numbers, n , k _{y} and k _{z} .

In the symmetric gauge

The derivation treated and y as 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 and .

A more adequate choice of gauge, is the symmetric gauge, which refers to the choice$$\backslash hat\; =\backslash frac\; \backslash mathbf\backslash times\; \backslash hat\; =\; \backslash frac\backslash begin\; -By\backslash \backslash \; Bx\; \backslash \backslash 0\; \backslash end.$$

In terms of dimensionless lengths and energies, the Hamiltonian can be expressed as$$\backslash hat\; =\; \backslash frac\; \backslash left[\backslash left(-i\backslash frac\{\backslash partial\}\{\backslash partial\; x\}\; +\; \backslash frac\{y\}\{2\}\backslash right)^2\; +\; \backslash left(-i\; \backslash frac\{\backslash partial\}\{\backslash partial\; y\}\; -\; \backslash frac\{x\}\{2\}\backslash right)^2\; \backslash right]$$

The correct units can be restored by introducing factors of

q, \hbar, B and m .

Consider operators

These operators follow certain commutation relations$$[\backslash hat\{a\},\; \backslash hat\{a\}^\{\backslash dagger\}]\; =\; [\backslash hat\{b\},\backslash hat\{b\}^\{\backslash dagger\}]\; =\; 1.$$

In terms of above operators the Hamiltonian can be written as$$\backslash hat\; =\; \backslash hbar\backslash omega\_\backslash left(\backslash hat^\backslash hat\; +\; \backslash frac\backslash right),$$where we reintroduced the units back.

The Landau level index

n is the eigenvalue of the operator \hat{N} = \hat{a} ^{\dagger} \hat{a} .

The application of

\hat{b} ^{\dagger} increases m _{z} by one unit while preserving n , whereas \hat{a} ^{\dagger} application simultaneously increase n and decreases m _{z} by one unit. The analogy to quantum harmonic oscillator provides solutions$$\backslash hat\; |n,m\_z\backslash rangle\; =\; E\_n\; |n,m\_z\backslash rangle,$$where$$E\_n\; =\; \backslash hbar\backslash omega\_\backslash left(n\; +\; \backslash frac\backslash right)$$and$$|n,m\_z\backslash rangle\; =\; \backslash frac\; \backslash frac|0,0\backslash rangle.$$

One may verify that the above states correspond to choosing wavefunctions proportional to$$\backslash psi\_(x,\; y)\; =\; \backslash left(\backslash frac\; -\; \backslash frac\; \backslash right)^n\; w^\; e^$$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}} .

Degeneracy of the Landau levels

In the Landau gauge

The effects of Landau levels may only be observed when the mean thermal energy is smaller than the energy level separation,

kT\ll\hbar\omega _{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\; =\; \backslash frac,$$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 _{0} , must physically lie within the system, 0 \le x _{0<L} _{x} . This gives the following range for N ,

For particles with charge

q = Ze , the upper bound on N can be simply written as a ratio of fluxes,$$\backslash frac\; =\; Z\backslash frac,$$where \Phi _{0=h/e} is the fundamental magnetic flux quantum and \Phi = 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)\; \backslash frac,$$which for electrons (where Z = 1 and S =1 / 2 ) gives D = 2 \Phi/\Phi _{0} , two available states for each flux quantum that penetrates the system.

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^[5] an active area of research.

In general, Landau levels are observed in electronic systems. 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 = \Phi/\Phi _{0} . Zeeman splitting has a significant effect on the Landau levels because their energy scales are the same, 2 \mu _{BB=\hbar\omega} _{c} . 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.

Moreover, the above derivation in the Landau gauge 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, (\hbar k ^{2 /( 2 m)} _{z)} , 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.

In the symmetric gauge

Each Landau level has degenerate orbitals labeled by the quantum numbers

m _{z} in symmetric gauge. The degeneracy per unit area is the same in each Landau level.

The z component of angular momentum is$$\backslash hat\_z\; =\; -i\; \backslash hbar\; \backslash frac\; =\; -\; \backslash hbar\; (\backslash hat^\backslash hat\; -\; \backslash hat^\backslash hat)$$

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 _{z} \hbar , where it is clear that m _{z} \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.

Relativistic case

See also: Dirac matter. An electron following Dirac equation under a constant magnetic field, can be analytically solved.^{[6] [7]} The energies are given by$$E\_=\backslash pm\; \backslash sqrt$$

where c is the speed of light, the sign depends on the particle-antiparticle component and ν is a non-negative integer. Due to spin, all levels are degenerate except for the ground state at .

The massless 2D case can be simulated in single-layer materials like graphene near the Dirac cones, where the eigenergies are given by^[8] $$E\_=\backslash pm\; \backslash sqrt$$where the speed of light has to be replaced with the Fermi speed v_F of the material and the minus sign corresponds to electron holes.

Magnetic susceptibility of a Fermi gas

The Fermi gas (an ensemble of non-interacting fermions) is part of the basis for understanding of the thermodynamic properties of metals. In 1930 Landau derived an estimate for the magnetic susceptibility of a Fermi gas, known as Landau susceptibility, which is constant for small magnetic fields. Landau also noticed that the susceptibility oscillates with high frequency for large magnetic fields,^[9] this physical phenomenon is known as the De Haas–Van Alphen effect.

Two-dimensional lattice

See also: Hofstadter's butterfly. The tight binding energy spectrum of charged particles in a two dimensional infinite lattice is known to be self-similar and fractal, as demonstrated in Hofstadter's butterfly. For an integer ratio of the magnetic flux quantum and the magnetic flux through a lattice cell, one recovers the Landau levels for large integers.^[10]

Integer quantum Hall effect

See main article: Quantum Hall effect. The energy spectrum of the semiconductor in a strong magnetic field forms Landau levels that can be labeled by integer indices. In addition, the Hall resistivity also exhibits discrete levels labeled by an integer . The fact that these two quantities are related can be shown in different ways, but most easily can be seen from Drude model: the Hall conductivity depends on the electron density as

\rho _{xy} = (B)/(n e) .

Since the resistivity plateau is given by

\rho _{xy} = (2 \pi\hbar)/(e ^{2}) (1)/(\nu) ,

the required density is

n = (B)/(\Phi _{0}) \nu,

which is exactly the density required to fill the Landau level. The gap between different Landau levels along with large degeneracy of each level renders the resistivity quantized.

See also

• Laughlin wavefunction

External links

• Web site: Lev Landau (1930). de. Diamagnetismus der Metalle.

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). .

Notes and References

1. Landau . L. . Lev Landau . Diamagnetismus der Metalle . Diamagnetism of Metals. Zeitschrift für Physik . Springer Science and Business Media LLC . 64 . 9–10 . 1930 . 1434-6001 . 10.1007/bf01397213 . 629–637 . 1930ZPhy...64..629L . 123206025 . de.
2. Web site: Charge in Magnetic Field . 11 March 2023 . courses.physics.illinois.edu.
3. An equally correct solution in the Landau gauge would be: A = \begin{pmatrix} - B y & 0 & 0 \end{pmatrix} ^{T} .
4. Book: Landau . L. D. . Lev Landau. Lifshitz. E. M.. Evgeny Lifshitz. Quantum mechanics : non-relativistic theory . 3rd. Butterworth Heinemann . Amsterdam . 1977 . 978-0-7506-3539-4 . 846962062 . 424–426.
5. 10.1016/S0921-4526(00)00769-9. A new approach to the ground state of quantum Hall systems. Basic principles. Physica B: Condensed Matter. 299. 6–31. 2001. Mikhailov . S. A. . 1–2. cond-mat/0008227 . 2001PhyB..299....6M . 118500817.
6. Rabi. I. I.. 1928. Das freie Elektron im homogenen Magnetfeld nach der Diracschen Theorie. Zeitschrift für Physik. de. 49. 7–8. 507–511. 10.1007/BF01333634. 1928ZPhy...49..507R . 121121095. 1434-6001.
7. Book: Berestetskii. V. B.. Quantum Electrodynamics: Volume 4. Pitaevskii. L. P.. Lifshitz. E. M.. 2012-12-02. Elsevier. 978-0-08-050346-2. en.
8. Yin. Long-Jing. Bai. Ke-Ke. Wang. Wen-Xiao. Li. Si-Yu. Zhang. Yu. He. Lin. 2017. Landau quantization of Dirac fermions in graphene and its multilayers. Frontiers of Physics. en. 12. 4 . 127208. 10.1007/s11467-017-0655-0. 1703.04241 . 2017FrPhy..12l7208Y . 2095-0462. free.
9. Book: Landau, L. D.. Lifshitz, E. M.. Statistical Physics: Volume 5. 22 October 2013. Elsevier. 978-0-08-057046-4. 177.
10. Analytis. James G.. Blundell. Stephen J.. Ardavan. Arzhang. May 2004 . Landau levels, molecular orbitals, and the Hofstadter butterfly in finite systems. American Journal of Physics. en. 72. 5. 613–618. 10.1119/1.1615568 . 2004AmJPh..72..613A . 0002-9505.