1.4 Brillouin zones

Concept of Brillouin zones

Brillouin zones define the fundamental region of reciprocal space you need to describe all the unique electronic and vibrational states in a crystal. They connect a crystal's real-space periodicity to the allowed wave vectors for electrons and phonons, making them essential for predicting everything from electrical conductivity to thermal behavior.

Definition and significance

A Brillouin zone is the Wigner-Seitz cell of the reciprocal lattice. That means it's the region of reciprocal space containing all points closer to one reciprocal lattice point than to any other. Every unique electron wave vector $\mathbf{k}$ in the crystal can be represented within this zone, so it acts as the "unit cell" of reciprocal space.

Why does this matter? Because the periodicity of a crystal means that wave vectors differing by a reciprocal lattice vector $\mathbf{G}$ are physically equivalent. The Brillouin zone captures exactly one copy of each distinct $\mathbf{k}$-state, which dramatically simplifies the analysis of electronic bands and phonon dispersions.

Relationship to reciprocal lattice

Brillouin zones are built directly from the reciprocal lattice. The zone boundaries are perpendicular bisector planes (Bragg planes) of the vectors connecting a central reciprocal lattice point to its neighbors. This means:

• The shape and symmetry of the Brillouin zone mirror the symmetry of the reciprocal lattice (and therefore the real-space crystal structure)
• Reciprocal lattice points sit at the vertices and edges of the zone boundaries
• A cubic real-space lattice produces a Brillouin zone with cubic-related symmetry, a hexagonal lattice produces a hexagonal zone, and so on

First vs. higher-order zones

The first Brillouin zone is the innermost region, enclosed by the nearest set of Bragg planes. It contains all the unique $\mathbf{k}$-vectors you need.

Higher-order zones are the successive shells you get by crossing additional Bragg planes outward from the first zone. The second Brillouin zone consists of all points reached by crossing exactly one Bragg plane from the first zone, the third zone by crossing two, and so on. Each higher zone has the same volume as the first and can be "folded back" into the first zone through translation by reciprocal lattice vectors $\mathbf{G}$. This folding is what produces multiple energy bands within the reduced zone scheme.

Construction of Brillouin zones

Building a Brillouin zone is a geometric procedure carried out in reciprocal space. Two equivalent approaches are commonly used.

Wigner-Seitz cell method

This is the standard construction:

1. Pick a reciprocal lattice point as the origin.
2. Draw vectors from the origin to all neighboring reciprocal lattice points.
3. Construct the perpendicular bisector plane of each vector at its midpoint.
4. The smallest enclosed region around the origin is the first Brillouin zone.

Every point inside this region is closer to the origin than to any other reciprocal lattice point. The procedure is identical to constructing a Wigner-Seitz cell in real space, just applied to the reciprocal lattice instead.

Bragg plane construction

This approach frames the same geometry in terms of diffraction. A Bragg plane is a plane in reciprocal space where the condition for constructive interference (Bragg reflection) is satisfied. For any reciprocal lattice vector $\mathbf{G}$, the Bragg plane is the set of wave vectors $\mathbf{k}$ satisfying:

$2\mathbf{k} \cdot \mathbf{G} = |\mathbf{G}|^2$

This is equivalent to the perpendicular bisector of $\mathbf{G}$. The first Brillouin zone is the region enclosed by the nearest Bragg planes. At these boundaries, electron waves undergo Bragg reflection, which is precisely where energy gaps open up in the band structure.

Symmetry considerations

The Brillouin zone inherits the full point-group symmetry of the crystal. This has a practical payoff: you only need to calculate electronic properties in the irreducible Brillouin zone, which is the smallest wedge of the full zone that can't be mapped onto a smaller piece by any symmetry operation (rotation, reflection, or inversion).

For example, in an FCC crystal the irreducible zone is roughly 1/48th of the full zone. High-symmetry points and lines within this wedge correspond to special $\mathbf{k}$-vectors where degeneracies or band crossings often occur.

Properties of Brillouin zones

Zone boundaries and symmetry points

Zone boundaries are where Bragg reflection occurs, and they're where energy band gaps typically appear. Within the zone, certain high-symmetry points have standard labels:

• $\Gamma$: the zone center ($\mathbf{k} = 0$), present in every Brillouin zone
• X, L, K, W, etc.: points on faces, edges, or corners of the zone, with specific labels depending on the lattice type

Band structure diagrams plot energy $E(\mathbf{k})$ along paths connecting these high-symmetry points. The choice of path (e.g., $\Gamma \to X \to W \to L \to \Gamma$) is dictated by the zone's symmetry and highlights the most physically important features like band gaps and band crossings.

Brillouin zone folding

When you map electronic states from higher Brillouin zones back into the first zone using reciprocal lattice translations, you get zone folding. A free-electron parabola $E = \hbar^2 k^2 / 2m$ in the extended zone scheme becomes multiple overlapping bands in the reduced scheme.

Where these folded bands would cross at zone boundaries, the crystal potential lifts the degeneracy and opens band gaps. This is the origin of the electronic band structure: the interplay between free-electron-like dispersion and Bragg reflection at zone boundaries.

Reduced zone scheme

The reduced zone scheme confines all $\mathbf{k}$-states to the first Brillouin zone. Any wave vector outside the first zone is translated back by an appropriate $\mathbf{G}$ vector. This is the standard representation used in most band structure plots and calculations because:

• It displays all bands compactly in one region
• It makes band gaps and avoided crossings immediately visible
• It facilitates direct comparison between different materials

The alternative extended zone scheme spreads bands across multiple zones, which is useful for seeing how the band structure evolves from the free-electron limit.

Applications in solid-state physics

Electronic band structure

Band structure calculations plot $E(\mathbf{k})$ along high-symmetry paths through the Brillouin zone. The position of band gaps relative to the Fermi level determines whether a material is a metal, semiconductor, or insulator. For instance, silicon has an indirect band gap of about 1.1 eV, visible as a gap between the valence band maximum near $\Gamma$ and the conduction band minimum near the $X$ point.

The density of states, integrated over the entire Brillouin zone, tells you how many electronic states are available at each energy and directly influences electrical conductivity, optical absorption, and thermoelectric properties.

Phonon dispersion curves

Lattice vibrations (phonons) are also described by wave vectors within the Brillouin zone. Phonon dispersion curves plot vibrational frequency $\omega(\mathbf{k})$ along high-symmetry directions. Two types of branches appear:

• Acoustic branches start at $\omega = 0$ at $\Gamma$ and correspond to atoms moving in phase (sound waves)
• Optical branches appear at higher frequencies when the unit cell contains more than one atom, corresponding to atoms moving out of phase

These curves are critical for understanding thermal conductivity, specific heat, and the speed of sound in a material.

X-ray diffraction analysis

Brillouin zone boundaries correspond to the conditions for Bragg diffraction. When X-rays scatter from crystal planes satisfying $2d\sin\theta = n\lambda$, the scattered wave vectors map onto Brillouin zone boundaries in reciprocal space. Systematic absences in diffraction patterns (missing reflections predicted by geometry) arise from the specific symmetry of the lattice and help determine the space group and lattice parameters of an unknown crystal.

Brillouin zones in different lattices

The shape of the first Brillouin zone depends entirely on the reciprocal lattice, which in turn depends on the real-space Bravais lattice.

Simple cubic lattice

The reciprocal lattice of a simple cubic (SC) lattice is also simple cubic, so the first Brillouin zone is a cube. Key high-symmetry points:

• $\Gamma$: zone center $(0, 0, 0)$
• $X$: center of a face
• $M$: center of an edge
• $R$: corner of the cube

This is the simplest Brillouin zone and serves as a useful reference case, though few real materials adopt the SC structure.

Body-centered cubic lattice

The reciprocal of a BCC lattice is an FCC lattice, so the BCC Brillouin zone is a truncated octahedron (a shape with 8 regular hexagonal faces and 6 square faces). Key points include:

• $\Gamma$: center
• $H$: center of a hexagonal face
• $P$: corner where three hexagons meet
• $N$: center of a square face

Many alkali metals (Na, K) and transition metals (Fe, Cr, W) have BCC structures, making this zone shape practically important.

Face-centered cubic lattice

The reciprocal of an FCC lattice is a BCC lattice, giving a Brillouin zone that is also a truncated octahedron, but with different proportions than the BCC case. Key points:

• $\Gamma$: center
• $X$: center of a square face
• $L$: center of a hexagonal face
• $W$: corner where a square and two hexagons meet
• $K$: midpoint of an edge between two hexagonal faces

This zone is relevant for many common metals (Cu, Al, Au, Ag) and semiconductors (Si, Ge, GaAs), making it one of the most frequently encountered Brillouin zones in practice.

Experimental techniques

Angle-resolved photoemission spectroscopy (ARPES)

ARPES directly maps the occupied electronic band structure. Photons (typically UV) eject electrons from a sample surface, and by measuring each electron's kinetic energy and emission angle, you can reconstruct its binding energy and crystal momentum $\mathbf{k}_\parallel$ within the Brillouin zone. The result is a direct image of $E(\mathbf{k})$ for occupied states below the Fermi level.

ARPES is particularly powerful for studying surface states, topological insulators, and high-temperature superconductors.

Neutron scattering

Thermal neutrons have wavelengths comparable to interatomic spacings and energies comparable to phonon energies, making them ideal probes for both structure and dynamics.

• Inelastic neutron scattering measures phonon dispersion curves throughout the Brillouin zone by tracking the energy and momentum transferred to the neutron
• Elastic neutron scattering determines crystal structure and, because neutrons carry a magnetic moment, magnetic ordering as well

Neutron scattering is especially valuable for studying magnetic materials and light elements that scatter X-rays weakly.

Electron diffraction methods

• LEED (Low-Energy Electron Diffraction): Low-energy electrons (20–200 eV) scatter from a crystal surface, producing a diffraction pattern that maps the surface reciprocal lattice and surface Brillouin zone. Useful for determining surface reconstructions.
• RHEED (Reflection High-Energy Electron Diffraction): High-energy electrons at grazing incidence monitor thin-film growth in real time. Oscillations in RHEED intensity correspond to layer-by-layer deposition, making it a standard tool in molecular beam epitaxy (MBE).

Computational methods

Density functional theory (DFT)

DFT is the workhorse of modern electronic structure calculations. It replaces the many-electron Schrödinger equation with a set of single-particle Kohn-Sham equations that yield the ground-state electron density and total energy. From a DFT calculation, you can extract band structures, densities of states, and total energies across the Brillouin zone.

DFT calculations sample the Brillouin zone on a discrete $\mathbf{k}$-point grid (e.g., a Monkhorst-Pack grid). Convergence with respect to grid density is an important practical consideration: too few $\mathbf{k}$-points give inaccurate results, while too many waste computational resources.

Tight-binding approximation

The tight-binding model builds electronic states from linear combinations of localized atomic orbitals. You specify a small number of hopping parameters (overlap integrals between neighboring orbitals), and the resulting Hamiltonian matrix can be diagonalized at each $\mathbf{k}$-point to produce band structures.

Tight-binding is computationally cheap compared to DFT and captures the essential physics of band formation. It's widely used for building intuition about how orbital overlap and crystal geometry shape band structures, and for modeling large systems where full DFT is too expensive.

Pseudopotential calculations

Core electrons are tightly bound and don't participate much in bonding or conduction. Pseudopotentials replace the true electron-ion potential (which has deep, rapidly oscillating features near the nucleus) with a smoother effective potential that reproduces the correct behavior of valence electrons outside the core region.

This dramatically reduces the computational cost of plane-wave DFT calculations by eliminating the need for many high-frequency basis functions. Modern norm-conserving and ultrasoft pseudopotentials are standard in codes like Quantum ESPRESSO and VASP.

Advanced concepts

Extended vs. reduced zone schemes

These are two ways of displaying the same physics:

• Extended zone scheme: Each band occupies its own Brillouin zone. The first band lives in the first zone, the second in the second zone, and so on. This makes the connection to free-electron behavior transparent.
• Reduced zone scheme: All bands are folded into the first Brillouin zone. This is the standard representation because it shows all bands together, making gaps and interactions between bands easy to identify.

You can always convert between the two by adding or subtracting reciprocal lattice vectors $\mathbf{G}$.

Jones zone and Fermi surface

The Fermi surface is the constant-energy surface in $\mathbf{k}$-space at the Fermi energy $E_F$. It separates occupied from unoccupied states at zero temperature and governs electrical and thermal transport in metals.

When the Fermi surface approaches or touches a Brillouin zone boundary, the electrons near that boundary experience strong Bragg reflection, which can open gaps and distort the Fermi surface. The Jones zone is defined by the set of Bragg planes where the free-electron Fermi sphere would intersect zone boundaries. In Hume-Rothery alloys, the electron count is tuned so that the Fermi surface just touches certain zone boundaries, stabilizing particular crystal structures.

Umklapp processes

In a normal phonon-phonon scattering event, the total crystal momentum is conserved within the first Brillouin zone. In an Umklapp process, the resulting wave vector falls outside the first zone and must be folded back by a reciprocal lattice vector:

$\mathbf{k}_1 + \mathbf{k}_2 = \mathbf{k}_3 + \mathbf{G}$

This effectively reverses the direction of momentum flow, providing a mechanism for thermal resistance. Without Umklapp scattering, a perfect crystal would have infinite thermal conductivity. These processes become significant at higher temperatures where phonons with large wave vectors (near zone boundaries) are thermally populated, and they also contribute to electrical resistivity and certain nonlinear optical effects.