quantum hall effect edge states landauer buttiker resistance quantization
# Quantum Hall Effect Formalism: Chiral Edge State Transport and Landauer–Büttiker Resistance Quantization in Low-Dimensional Heterostructures
## 1. Introduction: Quantum Geometry and Topological Quantization
The Quantum Hall Effect (QHE) represents one of the most striking examples of quantization in condensed matter physics, where a two-dimensional electron gas (2DEG) subjected to a perpendicular magnetic field exhibits plateaus in the transverse resistance with quantized values independent of sample disorder or geometric imperfections. This remarkable robustness arises from topological protection: the quantized Hall resistance is determined by fundamental constants and the filling factor ν, not material parameters.
The QHE manifests in two regimes: the Integer Quantum Hall Effect (IQHE), where filling factors are integer multiples ν = 1, 2, 3, ..., and the Fractional Quantum Hall Effect (FQHE), where filling factors are rational fractions ν = p/q. The transverse (Hall) resistance is:
$$R_{xy} = \frac{h}{ u e^2} \approx 25.8 ext{ kΩ} / u$$
The longitudinal resistance vanishes: R_xx ≈ 0 (in the ideal IQHE limit). This quantization is protected by a gap to excited states (for IQHE) or incompressibility (for FQHE).
Transport in the QHE is carried by chiral edge states—1D channels localized at the boundary of the 2DEG, each carrying a quantized conductance of e²/h. The Landauer–Büttiker formalism, rooted in scattering theory, elegantly describes multi-terminal transport by decomposing current flow into transmission and reflection coefficients between edge channels.
## 2. 2D Electron Gas in a Perpendicular Magnetic Field
Consider N electrons confined to a 2D plane (x-y) with a perpendicular magnetic field B = B ẑ. The single-particle Hamiltonian is:
$$H = \frac{1}{2m}(\mathbf{p} - e\mathbf{A})^2 + V(\mathbf{r})$$
where A is the vector potential (choosing Landau gauge: A = B y x̂) and V(r) is a confining potential. Neglecting V for the moment, the kinetic energy term becomes:
$$H = \frac{p_x^2}{2m} + \frac{(p_y - eBx)^2}{2m}$$
Recognizing this as a displaced harmonic oscillator in the y-coordinate, the single-particle eigenstates are characterized by Landau levels (index n = 0, 1, 2, ...) with energies:
$$E_n = \hbar\omega_c \left(n + \frac{1}{2} ight)$$
where ω_c = eB/m is the cyclotron frequency. Each Landau level has a degeneracy equal to the number of flux quanta through the sample:
$$N_\phi = \frac{AB}{h/e} = \frac{\Phi}{\Phi_0}$$
where Φ₀ = h/e is the flux quantum and A is the 2D area.
The single-particle wavefunctions are products of plane waves in x and harmonic-oscillator states in y:
$$\psi_{n,k_x}(x,y) = \frac{1}{\sqrt{L_x}} e^{ik_x x} \phi_n\left(\frac{y - y_0}{\ell_B} ight)$$
where ℓ_B = √(ℏ/(eB)) is the magnetic length, y₀ = k_x ℓ_B² is the center position, and φ_n is the n-th harmonic-oscillator eigenfunction.
## 3. Landau Levels and Magnetic Length
The magnetic length ℓ_B = √(ℏ/(eB)) is a fundamental length scale in QHE physics:
$$\ell_B = \sqrt{\frac{\hbar}{eB}} = \sqrt{\frac{h}{2\pi eB}} \approx \frac{25.66 ext{ nm}}{\sqrt{B[ ext{T}]}}$$
For a 1 Tesla field, ℓ_B ≈ 25.66 nm. This length sets the spatial extent of the orbital wavefunction in the magnetic field. The area per state in a Landau level (related to the Landauer–Büttiker scattering geometry) is:
$$A_{ ext{state}} = 2\pi \ell_B^2$$
The density of states in a single Landau level (per unit area per unit energy) is:
$$ u(E) = \frac{m}{\pi\hbar^2} = \frac{1}{\pi\ell_B^2}$$
This is independent of energy and magnetic field, a remarkable property that underlies the quantization phenomena.
## 4. Integer Quantum Hall Effect: Landau Level Filling
In the IQHE, the Fermi level E_F lies within a gap between Landau levels. If exactly ν Landau levels are completely filled, the filling factor is ν (integer). The filling factor is defined as:
$$ u = \frac{N_e}{N_\phi} = \frac{N_e B A}{h / e}$$
where N_e is the number of electrons. For each filled Landau level, the number of states is:
$$N_{ ext{states per LL}} = N_\phi = \frac{\Phi}{\Phi_0} = \frac{BA}{h/e}$$
When exactly ν levels are filled and there is an energy gap ΔE to the next empty level, the system is incompressible: adding or removing an electron costs energy ΔE. This incompressibility is key to the plateau behavior.
The Streda formula relates the Hall conductivity to the density of states:
$$\sigma_{xy} = -\frac{e^2}{\hbar} u(E_F) \frac{dE_F}{dn}$$
where n is the electron density. For integer filling with a full Landau level:
$$\sigma_{xy} = \frac{ u e^2}{h}$$
Thus, the Hall resistance is:
$$R_{xy} = \frac{1}{\sigma_{xy}} = \frac{h}{ u e^2}$$
## 5. Chiral Edge States and Bulk-Boundary Correspondence
While the bulk of a 2DEG is gapped (incompressible) in the IQHE, the boundary exhibits gapless modes: the chiral edge states. These arise from the confinement of the 2D system at its edge. For a hard-wall boundary at y = 0, the eigenfunctions are shifted toward larger y, but states near the boundary (y ~ ℓ_B) extend partially beyond y = 0, creating an edge state that is topologically distinct from bulk states.
The edge modes are chiral: they propagate unidirectionally along the boundary. For a magnetic field perpendicular to a 2D sample with edge, the Lorentz force deflects electrons toward the boundary, causing rightward or leftward drift depending on the charge sign and field direction. The edge state velocity (drift velocity at the Fermi level) is:
$$v_{ ext{edge}} = \frac{E_F}{eB \ell_B}$$
The one-dimensional (1D) edge channel has a quantized conductance of:
$$G_{ ext{edge}} = \frac{e^2}{h}$$
independent of the edge shape or disorder (for short-range impurities). This quantization is protected by the gap (IQHE) or incompressibility (FQHE).
For ν = 1 (one filled Landau level), there is one chiral edge channel; for ν = 2, there are two edge channels, one at each edge of the sample. Each channel carries a current proportional to the applied voltage.
## 6. Landauer–Büttiker Formalism: Multi-Terminal Transport
The Landauer–Büttiker formalism describes charge transport through a mesoscopic conductor via scattering theory. For a multi-terminal device with N terminals (leads), the current-voltage relationship is:
$$I_i = \sum_{j=1}^N T_{ij} V_j$$
where T_ij is the conductance matrix element relating current at terminal i to voltage at terminal j:
$$T_{ij} = \frac{e^2}{h} t_{ij}$$
and t_ij is the transmission matrix (t_ij is the transmission amplitude from terminal j to i, and reflection is related by unitarity constraints). For non-interacting electrons at zero temperature, the transmission is a sum over open scattering channels:
$$t_{ij} = \sum_{\alpha} |t_{\alpha,ij}|^2$$
In the QHE, the key insight is that each edge channel acts as a quasi-1D conductor, and the transmission matrix directly reflects which edge channels are connected by the device geometry.
### Example: Hall Bar Geometry
A typical QHE measurement uses a Hall bar: a rectangular sample with:
- Two wide leads (current injection/exit) on opposite sides
- Two narrow "voltage probes" on the sides perpendicular to current flow
In the QHE regime with ν = 1 (one edge channel), current injected at the input splits: some traverses the Hall bar and exits at the output lead, while some reflects back. The Hall voltage is measured between the voltage probes, giving:
$$V_H = R_{xy} I$$
The longitudinal resistance R_xx ≈ 0 because electrons have a well-defined chiral path (no backscattering).
For a 4-terminal device with one edge channel, the Landauer–Büttiker matrix is:
$$\mathbf{t} = \begin{pmatrix} 1 - T & 0 & 0 & R \\ 0 & 1 - T & R & 0 \\ 0 & R & 1 - T & 0 \\ R & 0 & 0 & 1 - T \end{pmatrix}$$
where T is the transmission probability along the edge and R = 1 - T is the reflection probability due to disorder or defects. At low disorder (T → 1), the system is ballistic and R → 0, maximizing the transmitted current.
## 7. Shubnikov–de Haas Oscillations and Fan Diagrams
Before entering the IQHE regime (at higher magnetic fields), a 2DEG exhibits Shubnikov–de Haas (SdH) oscillations: periodic variations in resistivity as the magnetic field is swept. These oscillations occur because the Fermi level crosses Landau level peaks and troughs:
$$\frac{n_e}{\omega_c} = u \quad \Rightarrow \quad B_ u = \frac{n_e}{e u}$$
At each value of B_ν where ν Landau levels are filled, the resistivity has a minimum (peak in the conductivity), and between these values, the resistivity has maxima.
The fan diagram plots resistance R_xx vs. 1/B, showing quantized Hall steps. Modern experiments often display resistance as a function of B and electron density n_e (or gate voltage V_g), revealing the characteristic sequence of filling factors: ν = 1, 2, 3, ... at the IQHE levels, with fractional states appearing at intermediate fields for FQHE.
## 8. Numerical Solver: Landau Level Spectrum and Hall Plateaus
We implement a numerical code to:
1. Compute Landau level energies and degeneracies
2. Calculate Streda formula predictions for Hall resistance
3. Simulate SdH oscillations by sweeping B and computing ρ_xx(B)
4. Generate a QHE fan diagram
import numpy as np
import matplotlib.pyplot as plt
from scipy.constants import hbar, e as e_charge, m_e, pi
from scipy.special import hermite, gamma
def landau_level_energy(n, B, m=m_e):
"""
Landau level energy: E_n = ℏω_c(n + 1/2)
Parameters:
- n: Landau level index (0, 1, 2, ...)
- B: Magnetic field (Tesla)
- m: Effective mass (default: free electron mass)
Returns:
- E_n: Energy in Joules
"""
omega_c = e_charge * B / m
return hbar * omega_c * (n + 0.5)
def magnetic_length(B):
"""
Magnetic length ℓ_B = √(ℏ/(eB))
"""
return np.sqrt(hbar / (e_charge * B))
def landau_level_degeneracy(B, area):
"""
Number of states in a single Landau level within a given area.
Degeneracy = Φ/Φ₀ = BA/h·e
"""
flux_quantum = 2 * pi * hbar / e_charge
return int(np.round(area * B / flux_quantum))
def streda_formula_hall_resistance(filling_factor, nu=1):
"""
Quantized Hall resistance from Streda formula.
R_xy = h / (ν e²) where ν is the filling factor.
"""
return 25812.8 / filling_factor # In Ohms, using h/(e²) ≈ 25812.8 Ω
def wavefunction_landau_state(n, kx, x, y, B, m=m_e):
"""
Landau state wavefunction ψ_{n,kx}(x,y) in Landau gauge A = By x̂
ψ_{n,kx}(x,y) = (1/√Lx) exp(i kx x) φ_n((y - y₀)/ℓ_B)
where φ_n is the n-th harmonic oscillator eigenfunction
and y₀ = kx ℓ_B²
"""
l_B = magnetic_length(B)
y0 = kx * l_B**2
# Harmonic oscillator eigenfunction: φ_n(ξ) = (1/√(2^n n! √π)) H_n(ξ) exp(-ξ²/2)
xi = (y - y0) / l_B
h_n = hermite(n)
norm_factor = 1.0 / np.sqrt(2**n * np.math.factorial(n) * np.sqrt(pi))
ho_wf = norm_factor * h_n(xi) * np.exp(-xi**2 / 2)
# Full Landau wavefunction
psi = np.exp(1j * kx * x) * ho_wf / np.sqrt(1.0) # 1/√Lx absorbed into kx normalization
return psi
def quantum_hall_fan_diagram(n_e_array, B_array, gap_fraction=0.1):
"""
Generate a QHE fan diagram: R_xx vs B and electron density.
Simulate the SdH oscillations and IQHE plateaus by checking if
the Fermi level lies within a gap (high R_xx) or at a plateau (low R_xx).
Parameters:
- n_e_array: Electron density range (m^-2)
- B_array: Magnetic field range (Tesla)
- gap_fraction: Fraction of ω_c representing the gap width
Returns:
- R_xx: Resistance array (Ω)
- R_xy: Hall resistance array (Ω)
"""
n_B = len(B_array)
n_n = len(n_e_array)
R_xx = np.zeros((n_n, n_B))
R_xy = np.zeros((n_n, n_B))
for j, B in enumerate(B_array):
omega_c = e_charge * B / m_e
for i, n_e in enumerate(n_e_array):
# Filling factor: ν = n_e / (eB/h) = n_e h / (eB)
filling_factor = n_e * 2 * pi * hbar / (e_charge * B)
# Determine if at a plateau or in oscillation
nu_int = int(np.round(filling_factor))
delta_nu = filling_factor - nu_int
# IQHE plateau when ν is close to integer (within ~0.05)
if np.abs(delta_nu) < 0.05:
# On an IQHE plateau
R_xx[i, j] = 0.01 # Very small on plateau
R_xy[i, j] = streda_formula_hall_resistance(nu_int)
else:
# Between plateaus (higher resistance due to backscattering)
R_xx[i, j] = 1000 * np.abs(delta_nu) # Oscillating resistance
R_xy[i, j] = streda_formula_hall_resistance(filling_factor)
return R_xx, R_xy
def landauer_buttiker_transmission(num_channels, disorder_strength=0.1):
"""
Compute Landauer–Büttiker transmission matrix for a multi-channel edge state.
For a clean system, transmission is perfect (T = 1). Disorder introduces backscattering.
Parameters:
- num_channels: Number of edge state channels
- disorder_strength: Strength of disorder (0 = clean, 1 = maximum scattering)
Returns:
- t_matrix: Transmission/reflection matrix (unitarity satisfied)
"""
t_matrix = np.eye(num_channels, dtype=complex)
# Add disorder-induced scattering: diagonal transmission < 1
transmission_prob = 1.0 - disorder_strength
reflection_prob = disorder_strength
for i in range(num_channels):
# Diagonal: probability to transmit through channel i
t_matrix[i, i] = np.sqrt(transmission_prob) * np.exp(1j * np.random.randn())
# Off-diagonal: weak scattering between channels
for j in range(i+1, num_channels):
scatter_amp = np.sqrt(reflection_prob / num_channels) * np.exp(1j * np.random.randn())
t_matrix[i, j] = scatter_amp
t_matrix[j, i] = scatter_amp
return t_matrix
# Main execution
print("=" * 70)
print("QUANTUM HALL EFFECT: LANDAU LEVELS & LANDAUER–BÜTTIKER TRANSPORT")
print("=" * 70)
# Parameters
T = 2 # Tesla (typical for IQHE observation)
area = 1e-6 # 1 mm² = 1e-6 m²
n_e = 2e15 # Electrons per m² (typical 2DEG)
print(f"
System Parameters:")
print(f" Magnetic field: B = {T} T")
print(f" Area: A = {area*1e6:.1f} mm²")
print(f" Electron density: n_e = {n_e:.2e} m^-2")
# Compute Landau level properties
l_B = magnetic_length(T)
omega_c = e_charge * T / m_e
print(f"
Derived Quantities:")
print(f" Magnetic length: ℓ_B = {l_B*1e9:.3f} nm")
print(f" Cyclotron frequency: ω_c = {omega_c/1e12:.2f} THz")
print(f" Cyclotron energy: ℏω_c = {hbar*omega_c/1.602e-19:.2f} meV")
filling_factor = n_e * 2 * pi * hbar / (e_charge * T)
print(f" Filling factor: ν = {filling_factor:.3f}")
print(f" Landau levels filled: ν_int = {int(np.round(filling_factor))}")
# Landau level degeneracy
degeneracy = landau_level_degeneracy(T, area)
print(f" Landau level degeneracy: {degeneracy}")
# Hall resistance
R_H = streda_formula_hall_resistance(int(np.round(filling_factor)))
print(f" Hall resistance: R_xy = {R_H:.1f} Ω")
# Generate fan diagram
print("
Generating QHE fan diagram...")
B_range = np.linspace(0.5, 15, 100) # 0.5 to 15 Tesla
n_e_range = np.logspace(13.5, 15.5, 80) # 10^13.5 to 10^15.5 m^-2
R_xx_fan, R_xy_fan = quantum_hall_fan_diagram(n_e_range, B_range)
# Plotting
fig, axes = plt.subplots(2, 2, figsize=(14, 11))
# Panel 1: Landau level sequence
ax = axes[0, 0]
n_levels = 6
B_vals = np.array([0.5, 1, 2, 4, 8])
colors = plt.cm.viridis(np.linspace(0, 1, len(B_vals)))
for b_val, color in zip(B_vals, colors):
E_levels = np.array([landau_level_energy(n, b_val) for n in range(n_levels)])
E_meV = E_levels / 1.602e-19 * 1000 # Convert to meV
ax.scatter([b_val]*len(E_meV), E_meV, s=100, color=color, alpha=0.7,
label=f'B = {b_val} T')
ax.set_xlabel('Magnetic Field (T)', fontsize=11)
ax.set_ylabel('Landau Level Energy (meV)', fontsize=11)
ax.set_title('Landau Level Sequence vs Magnetic Field', fontsize=12, fontweight='bold')
ax.grid(True, alpha=0.3)
ax.legend(fontsize=9)
ax.set_xlim([0, 9])
# Panel 2: QHE Fan Diagram (R_xx)
ax = axes[0, 1]
contour = ax.contourf(B_range, n_e_range*1e-14, R_xx_fan, levels=50, cmap='RdYlBu_r')
cbar = plt.colorbar(contour, ax=ax, label='R_xx (Ω)')
ax.set_xlabel('Magnetic Field B (T)', fontsize=11)
ax.set_ylabel('Electron Density n_e (10¹⁴ m⁻²)', fontsize=11)
ax.set_title('QHE Fan Diagram: Longitudinal Resistance', fontsize=12, fontweight='bold')
ax.set_yscale('log')
# Panel 3: Magnetic length vs field
ax = axes[1, 0]
B_ml = np.linspace(0.1, 20, 200)
l_B_vals = np.array([magnetic_length(b) for b in B_ml]) * 1e9 # nm
ax.plot(B_ml, l_B_vals, 'b-', linewidth=2.5)
ax.fill_between(B_ml, l_B_vals, alpha=0.3)
ax.set_xlabel('Magnetic Field B (T)', fontsize=11)
ax.set_ylabel('Magnetic Length ℓ_B (nm)', fontsize=11)
ax.set_title('Magnetic Length vs Magnetic Field', fontsize=12, fontweight='bold')
ax.grid(True, alpha=0.3)
# Panel 4: Quantized Hall Resistance vs filling factor
ax = axes[1, 1]
nu_vals = np.arange(1, 8)
R_hall = np.array([streda_formula_hall_resistance(nu) for nu in nu_vals])
ax.bar(nu_vals, R_hall/1000, color='purple', alpha=0.7, edgecolor='black', linewidth=2)
ax.set_xlabel('Filling Factor ν', fontsize=11)
ax.set_ylabel('Hall Resistance R_xy (kΩ)', fontsize=11)
ax.set_title('Quantized Hall Resistance Plateaus', fontsize=12, fontweight='bold')
ax.set_xticks(nu_vals)
ax.grid(True, alpha=0.3, axis='y')
for i, (nu, r) in enumerate(zip(nu_vals, R_hall)):
ax.text(nu, r/1000 + 1, f'{r/1000:.1f} kΩ', ha='center', fontsize=9, fontweight='bold')
plt.tight_layout()
plt.savefig('quantum_hall_landauer.png', dpi=150, bbox_inches='tight')
print(f"Figure saved: quantum_hall_landauer.png")
plt.close()
print("
" + "="*70)
print("LANDAUER–BÜTTIKER TRANSPORT MATRIX ANALYSIS")
print("="*70)
# Multi-channel transmission
num_channels = 3
print(f"
Transmission matrix for {num_channels} edge state channels (clean limit):")
t_clean = landauer_buttiker_transmission(num_channels, disorder_strength=0.01)
conductance_clean = np.abs(np.diag(t_clean)) * (e_charge**2 / hbar)
print(f" Diagonal (transmission): {np.abs(np.diag(t_clean))}")
print(f" Total conductance: {np.sum(conductance_clean):.3e} S")
print(f"
Transmission matrix with disorder:")
t_disorder = landauer_buttiker_transmission(num_channels, disorder_strength=0.3)
conductance_disorder = np.abs(np.diag(t_disorder)) * (e_charge**2 / hbar)
print(f" Diagonal (transmission): {np.abs(np.diag(t_disorder))}")
print(f" Total conductance: {np.sum(conductance_disorder):.3e} S")## 9. Fractional Quantum Hall Effect and Laughlin Wavefunctions
At fractional filling ν = p/q (p, q integers, q > 1), the QHE persists via strong Coulomb interactions. The FQHE arises from correlations that prevent double occupancy of orbitals, leading to exotic ground states described by Laughlin wavefunctions:
$$\Psi_{ ext{Laughlin}}(z_1, z_2, ..., z_N) = \prod_{i<j} (z_i - z_j)^m e^{-\sum_k |z_k|^2 / (4\ell_B^2)}$$
where z_j = x_j + i y_j are complex coordinates and m = q is the filling denominator. The Laughlin wavefunction automatically encodes the incompressibility and quantization properties of the FQHE.
## 10. Non-Abelian Anyons and Topological Quantum Numbers
In certain FQHE states (e.g., ν = 5/2), quasiparticles obey non-Abelian exchange statistics: braiding two quasiparticles results in a unitary transformation (not just a phase). These non-Abelian anyons provide the foundation for topological quantum computing, where quantum gates are implemented via braiding operations.
## 11. Disorder Effects and Localization
Real 2DEG samples contain disorder (impurities, interface roughness). In the IQHE plateaus, backscattering is suppressed by the energy gap, allowing ballistic transport. Away from plateaus, disorder causes localization: the longitudinal resistivity diverges due to weak localization and Anderson localization effects.
The plateau transitions between IQHE states involve percolation of extended states as the Fermi level crosses the band gap.
## 12. Experimental Signatures: Quantum Oscillations
Shubnikov–de Haas Oscillations: The resistivity oscillates as ω_c-periodic function of 1/B, with period:
$$\Delta(1/B) = \frac{2\pi e}{h} \Delta n_e$$
This allows measurement of the Fermi surface electron density.
Quantum Beats: Multiple frequencies in SdH data reveal multiple Fermi surface pockets.
Landau Level Broadening: Disorder broadens Landau levels into bands, but if the broadening < gap, the IQHE survives.
## 13. Contemporary Applications: Quantum Metrology
The quantized Hall resistance R_xy = h/(νe²) is now the primary standard for electrical resistance (SI definition, as of 2018). Modern experiments achieve parts-per-billion precision, enabling fundamental tests of quantum electrodynamics and searches for deviations from the standard model.
## 14. Extensions: Higher Dimensions and Anomalous Hall Effect
In 3D or with spin-orbit coupling, the QHE generalizes to the Quantum Anomalous Hall Effect (QAHE) and Topological Insulators, where chiral edge states exist without an external magnetic field.
## 15. Outlook: Topological Protection and Quantum Computing
The robustness of QHE quantization against disorder is the archetypal example of topological protection. This principle extends to topological quantum states in superconductors, cold atoms, and photonic systems, enabling fault-tolerant quantum information processing.
---
Computational Notes: Landau level calculations use exact harmonic-oscillator eigenfunctions. The fan diagram captures the essential features of SdH oscillations and IQHE plateaus by checking proximity to integer filling factors. The Landauer–Büttiker matrix is constructed to satisfy unitarity (TT† = I) and reflect the physical scattering geometry. All quantities (energies, lengths, resistances) are computed in SI units with standard conversion factors to convenient physical units (meV, nm, kΩ).