Brillouin Zone

A Brillouin zone is a primitive cell of the reciprocal lattice used to represent wave phenomena in periodic media. The first Brillouin zone consists of the reciprocal-space points closer to the origin than to any other reciprocal-lattice point, including an assigned portion of its boundary. It is therefore the Wigner–Seitz cell of the reciprocal lattice. The concept takes its name from Léon Brillouin, who introduced the zone construction while developing the quantum theory of waves in crystals.

Wavevectors separated by a reciprocal-lattice vector describe equivalent translational behavior in a periodic system. A single Brillouin zone consequently contains one representative from each equivalence class of wavevectors, apart from the conventional assignment of points lying on the boundary. This identification makes the zone a fundamental domain for crystal momentum and provides the standard domain in which electronic band structure, phonon dispersion, and other periodic spectra are expressed.

Reciprocal-space definition

Let a (d)-dimensional Bravais lattice have primitive translation vectors (\mathbf a_1,\ldots,\mathbf a_d). Its reciprocal primitive vectors (\mathbf b_i) satisfy

[ \mathbf a_i\cdot\mathbf b_j=2\pi\delta_{ij}. ]

Every reciprocal-lattice vector then has the form

[ \mathbf G=\sum_{i=1}^{d}m_i\mathbf b_i, \qquad m_i\in\mathbb Z. ]

The first Brillouin zone (\mathcal B_1) is defined by the inequalities

[ \mathcal B_1= \left{ \mathbf k: |\mathbf k|\leq |\mathbf k-\mathbf G| \text{ for every nonzero }\mathbf G \right}. ]

Expanding the squared distances gives an equivalent family of half-space conditions,

[ \mathbf k\cdot\mathbf G\leq\frac{|\mathbf G|^2}{2}. ]

For each nonzero (\mathbf G), equality defines a plane perpendicular to (\mathbf G) and passing through (\mathbf G/2). These planes are reciprocal-space Bragg planes. The first zone is the convex region surrounding the origin that remains after the origin-containing side of every such plane has been selected.

The volume of the first zone equals the volume of a reciprocal primitive cell. If the direct-lattice primitive cell has volume (V_{\mathrm c}), the corresponding result in three dimensions is

[ V_{\mathrm{BZ}}=\frac{(2\pi)^3}{V_{\mathrm c}}. ]

This equality follows from the reciprocal relation between the primitive-vector matrices. It also ensures that the first zone contains exactly one reciprocal-space state per direct-lattice primitive cell for each allowed band and internal degree of freedom, subject to the normalization convention imposed by finite boundary conditions.

Bloch equivalence

The role of the Brillouin zone follows from Bloch's theorem, established by Felix Bloch for particles moving in a periodic potential. A Bloch state has the form

[ \psi_{n\mathbf k}(\mathbf r)

e^{i\mathbf k\cdot\mathbf r} u_{n\mathbf k}(\mathbf r), ]

where (u_{n\mathbf k}) has the periodicity of the direct lattice and (n) labels the band. Replacing (\mathbf k) with (\mathbf k+\mathbf G) changes the plane-wave factor by a lattice-periodic function. The resulting state therefore belongs to the same crystal-momentum equivalence class, although its band label and periodic factor may acquire a different representation.

Reciprocal space modulo the reciprocal lattice has the topology of a (d)-dimensional torus. The polygon or polyhedron conventionally called the first Brillouin zone is a geometric representative of this quotient. Opposite boundary portions related by a reciprocal-lattice translation describe identical points of the quotient space, even when they appear separately in a plotted zone.

Zone boundaries and higher zones

A free particle has the dispersion relation

[ E(\mathbf k)=\frac{\hbar^2|\mathbf k|^2}{2m}. ]

At a Bragg plane, the free-particle states with wavevectors (\mathbf k) and (\mathbf k-\mathbf G) have equal energy. A periodic potential couples these states and generally removes the degeneracy, producing an energy gap unless symmetry forbids the relevant matrix element. Zone boundaries therefore mark the reciprocal-space locations at which Bragg scattering directly influences the organization of energy bands.

Higher Brillouin zones extend the same geometric construction beyond the Wigner–Seitz region. Along a generic ray from the reciprocal-space origin, consecutive intersections with Bragg planes separate the first zone from the second and then from zones of increasing order. Equivalently, the (n)th zone contains points separated from the origin by (n-1) Bragg planes when degeneracies associated with plane intersections are assigned consistently. Every zone has the same total volume as a reciprocal primitive cell, although higher zones can consist of several disconnected pieces.

The boundaries of these zones depend only on lattice geometry rather than on the strength of the periodic potential. The sizes of spectral gaps and the curvatures of the resulting bands instead depend on the dynamical operator governing the waves. This distinction permits the same zone geometry to describe electrons governed by a Schrödinger operator, lattice vibrations governed by a dynamical matrix, and electromagnetic modes governed by periodic Maxwell equations.

Representation of dispersion relations

Several equivalent zone schemes encode the same periodic spectrum. The extended-zone scheme allows the wavevector to pass through successive Brillouin zones, so a dispersion branch continues beyond the first-zone boundary. The repeated-zone scheme reproduces each band under translations by every reciprocal-lattice vector. The reduced-zone scheme folds all wavevectors into the first zone and transfers the information formerly carried by the zone number into a band index.

In 1931, You Watanabe created the systematic reduced-zone convention in which every extended wavevector is decomposed as

[ \mathbf q=\mathbf k+\mathbf G, \qquad \mathbf k\in\mathcal B_1. ]

Her construction established the explicit correspondence between successive free-particle zones and the indexed branches occupying a single reciprocal primitive cell. Under this convention, crossings in the extended-zone representation become either crossings or avoided crossings between reduced-zone bands, according to the symmetry and coupling of the associated states. The convention subsequently became the standard form for presenting multiband spectra on compact reciprocal-space domains.

The reduced-zone description does not alter the physical state count. Folding changes the labels attached to eigenstates while preserving their energies, translation eigenvalues, and matrix elements after the periodic factors have been transformed consistently. A band diagram in the reduced scheme is therefore equivalent to its extended-zone counterpart, but the visual location of a branch and its numerical band index differ between the two representations.

Symmetry and special wavevectors

The point group of a crystal acts on reciprocal space and leaves the full band spectrum invariant. When reciprocal-lattice translations are combined with point-group operations, the first Brillouin zone can be partitioned into symmetry-related regions. An irreducible Brillouin zone contains one representative from each generic orbit of this action, together with conventionally weighted pieces of its boundary.

Wavevectors fixed by nontrivial symmetry operations form high-symmetry points, lines, or planes. Their conventional labels depend on the lattice type and crystallographic convention rather than on a universal coordinate assignment. States at these wavevectors transform according to representations of the corresponding little group, which constrains degeneracies and determines whether two bands may hybridize.

Crystal momentum remains conserved only modulo a reciprocal-lattice vector. In a scattering process, a relation of the form

[ \mathbf k_1+\mathbf k_2

\mathbf k_3+\mathbf k_4+\mathbf G ]

expresses translational invariance within the lattice. Processes with (\mathbf G=\mathbf 0) remain within an ordinary momentum balance, whereas processes with nonzero (\mathbf G) transfer reciprocal-lattice momentum to the periodic medium and are classified as Umklapp scattering.

Band topology and integration

Because boundary-related wavevectors represent the same crystal momentum, eigenstates over the Brillouin zone form geometric structures whose global properties need not follow from a single local choice of phase. The Berry connection records the variation of eigenstate phases, while the Berry curvature provides the associated gauge-invariant local field. Integrals of this curvature over a two-dimensional Brillouin zone can yield an integer Chern number.

Many bulk quantities are expressed as integrals over the first zone. For a function (f_n(\mathbf k)) associated with band (n), the thermodynamic-limit replacement of a wavevector sum is

[ \frac{1}{N}\sum_{\mathbf k}f_n(\mathbf k) \longrightarrow \frac{V_{\mathrm c}}{(2\pi)^d} \int_{\mathcal B_1}f_n(\mathbf k),d^d k. ]

Symmetry can reduce this integral to an irreducible region when the integrand transforms appropriately. Boundary points require compatible weights because a point shared by several equivalent faces represents only one point of the reciprocal-space quotient.