In a crystalline solid, electrons do not occupy isolated atomic energy levels but instead form continuous "bands" of allowed energies separated by "band gaps" where no electron states exist. This band structure arises because the periodic lattice potential diffracts electron waves (Bloch waves), and at wave vectors equal to multiples of π/a (the Brillouin zone boundaries, where a is the lattice spacing) constructive and destructive interference split what were single levels into pairs of bands. The nearly-free electron model used here is the simplest quantum-mechanical treatment: it starts from free electrons and treats the lattice potential as a weak perturbation, giving band gaps of magnitude 2|V𝔷|, where V𝔷 is the Fourier component of the periodic potential.
The simulation plots the dispersion relation E(k) — energy versus wave vector — for the reduced Brillouin zone. You can adjust the lattice potential strength to see band gaps open and close, switch between metal, semiconductor, and insulator filling levels, and observe how the density of states changes shape as bands flatten near zone boundaries.
What is the difference between a metal, semiconductor, and insulator?
The distinction lies in how the Fermi level sits relative to the band structure. In a metal, the Fermi level cuts through a partially filled band, so electrons near Eᴹ can be accelerated by a small electric field and conduct easily. In a semiconductor the Fermi level lies in a band gap, but the gap is small (typically 0.5–3 eV for silicon, germanium, and GaAs), so thermal energy (~0.026 eV at room temperature) can promote electrons to the conduction band. In an insulator the gap exceeds ~5 eV and virtually no electrons can bridge it thermally.
What is the Brillouin zone?
The first Brillouin zone is the Wigner-Seitz cell of the reciprocal lattice — the region of k-space closer to the origin than to any other reciprocal lattice point. For a 1D lattice with spacing a, it spans −π/a to +π/a. Electron states outside this zone are identical to states within it (differing by a reciprocal lattice vector G = 2πn/a), so all distinct states can be "folded back" into the first zone — the reduced zone scheme used in this simulation.
Why do band gaps open at zone boundaries?
At k = π/a, the electron wavelength equals the lattice period (λ = 2a), satisfying the Bragg condition for diffraction. The forward- and backward-travelling waves mix to form two standing waves: one with its probability density concentrated at atom sites (lower energy, bonding state) and one concentrated between atoms (higher energy, antibonding state). The energy difference between these two standing waves is the band gap, equal to 2|V𝔷| in the nearly-free electron model.
Bloch's theorem states that the eigenstates of an electron in a periodic potential have the form ψᴷ(r) = e^(ik⋅r) uᴷ(r), where uᴷ(r) has the same periodicity as the lattice. These are Bloch waves: plane waves modulated by a lattice-periodic function. The quantum number k (the crystal momentum) is a conserved quantity analogous to ordinary momentum in free space, and the band index n labels the solution for a given k.
Near a band minimum, the E(k) dispersion is approximately parabolic: E ≈ E₀ + ħ²k²/(2m*), where m* is the effective mass. A small effective mass means the band is highly curved (large second derivative), and the electron responds strongly to electric fields — it accelerates quickly. GaAs has m* ≈ 0.067m𝑒, making it very responsive and ideal for high-frequency transistors. Near a band maximum, m* is negative, giving rise to the concept of holes as positively charged quasiparticles.
The density of states g(E) is the number of electron states per unit energy interval per unit volume. For free electrons in 3D it scales as g(E) ∝ E^(1/2). Band structure modifies this significantly: near band edges g(E) follows the free-electron form, but near zone boundaries bands flatten (d²E/dk² → 0), causing sharp peaks in g(E) called van Hove singularities. These peaks are visible in X-ray absorption spectra and tunnelling measurements.
Doping introduces impurity atoms that donate extra electrons (n-type, e.g., phosphorus in silicon) or create electron vacancies/holes (p-type, e.g., boron in silicon). The band structure itself hardly changes, but the Fermi level shifts: toward the conduction band for n-type doping, toward the valence band for p-type. At high doping concentrations (>10ⁿ⁹ cm⁻³) the material behaves metallically and the Fermi level enters the band.
A direct band gap semiconductor (such as GaAs) has its conduction band minimum and valence band maximum at the same k-point, usually k = 0 (the Γ point). An indirect gap semiconductor (such as silicon or germanium) has them at different k-values. Optical transitions are much more likely in direct-gap materials because photons carry negligible momentum, so no phonon is needed to conserve crystal momentum. This is why GaAs makes efficient LEDs and lasers while silicon does not.
The tight-binding model is the complement of the nearly-free electron model: instead of starting from free electrons and adding a weak lattice potential, it starts from tightly bound atomic orbitals and adds weak overlap integrals between neighbours. The resulting dispersion for a 1D chain with lattice spacing a and hopping integral t is E(k) = ϵ₀ − 2t⋅cos(ka), a cosine band. The tight-binding model is more accurate for d-electrons in transition metals and for graphene, where the linear dispersion at the Dirac point emerges naturally.
The main experimental technique is Angle-Resolved Photoemission Spectroscopy (ARPES), which shines ultraviolet or X-ray light on a crystal surface and measures the kinetic energy and emission angle of ejected electrons. From energy and momentum conservation, both E and k of the original electrons are reconstructed. ARPES can map entire band structures with meV energy and milli-radian angular resolution, and has been crucial in discovering topological insulators and unconventional superconductors.