1.3 Reciprocal lattice

Definition and concept

The reciprocal lattice transforms a real-space crystal structure into a frequency-domain representation. Instead of working with atoms arranged periodically in space, you work with a set of wavevectors that encode the same structural information. This makes it far easier to analyze diffraction experiments and calculate electronic properties of solids.

Fourier transform of real lattice

Think of the reciprocal lattice as the Fourier transform of the real-space lattice. Any periodic function (like the electron density in a crystal) can be decomposed into a sum of plane waves, and the reciprocal lattice tells you which wavevectors appear in that sum.

A reciprocal lattice vector $\mathbf{G}$ is defined by the condition:

$e^{i\mathbf{G}\cdot\mathbf{R}} = 1$

for all real-space lattice vectors $\mathbf{R}$. This means $\mathbf{G} \cdot \mathbf{R} = 2\pi n$ for some integer $n$. Every vector satisfying this condition is a point in the reciprocal lattice, and the full set of these points forms a lattice of its own with its own periodicity and symmetry.

The payoff: complicated spatial patterns in real space become discrete points in reciprocal space, which are much easier to work with mathematically.

Significance in crystallography

Reciprocal space is where diffraction data naturally lives. When X-rays scatter off a crystal, each diffraction peak corresponds directly to a reciprocal lattice point. This means you can read off structural information (lattice parameters, symmetry, atomic positions) straight from the diffraction pattern without first reconstructing the real-space arrangement.

This connection between diffraction peaks and reciprocal lattice points is what makes the reciprocal lattice indispensable for experimental crystallography.

Properties of reciprocal lattice

Periodicity and symmetry

The reciprocal lattice inherits the point group symmetry of the real lattice. If the real crystal has a four-fold rotation axis, so does its reciprocal lattice.

The key inverse relationship: dimensions in reciprocal space scale inversely with real space. A crystal with a large unit cell produces a finely spaced reciprocal lattice, and vice versa. Specifically, reciprocal lattice spacings have units of inverse length (e.g., $\text{Å}^{-1}$), reinforcing the idea that you're working in "frequency space" rather than position space.

Relationship to Brillouin zones

Brillouin zones are constructed directly from the reciprocal lattice. The first Brillouin zone is the Wigner-Seitz cell of the reciprocal lattice: the region of reciprocal space closer to the origin than to any other reciprocal lattice point.

Its boundaries are Bragg planes, which bisect the lines connecting the origin to neighboring reciprocal lattice points. Higher-order Brillouin zones are the next shells outward, separated by additional Bragg planes.

Why this matters: electronic band structure and phonon dispersion relations are conventionally plotted within the first Brillouin zone. It's the fundamental domain for describing wave-like excitations in a periodic solid.

Mathematical formulation

Reciprocal lattice vectors

Given real-space primitive vectors $\mathbf{a}_1$, $\mathbf{a}_2$, $\mathbf{a}_3$, the reciprocal lattice vectors are:

$\mathbf{b}_1 = 2\pi \frac{\mathbf{a}_2 \times \mathbf{a}_3}{\mathbf{a}_1 \cdot (\mathbf{a}_2 \times \mathbf{a}_3)}$

$\mathbf{b}_2 = 2\pi \frac{\mathbf{a}_3 \times \mathbf{a}_1}{\mathbf{a}_1 \cdot (\mathbf{a}_2 \times \mathbf{a}_3)}$

$\mathbf{b}_3 = 2\pi \frac{\mathbf{a}_1 \times \mathbf{a}_2}{\mathbf{a}_1 \cdot (\mathbf{a}_2 \times \mathbf{a}_3)}$

The denominator $\mathbf{a}_1 \cdot (\mathbf{a}_2 \times \mathbf{a}_3)$ is the volume of the real-space unit cell, $V_{\text{cell}}$.

These vectors satisfy the orthonormality condition:

$\mathbf{a}_i \cdot \mathbf{b}_j = 2\pi \delta_{ij}$

where $\delta_{ij}$ is the Kronecker delta (1 if $i = j$, 0 otherwise). This condition is what guarantees $e^{i\mathbf{G}\cdot\mathbf{R}} = 1$.

Any general reciprocal lattice vector is then written as:

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

where $h, k, l$ are integers. These are the same Miller indices used to label crystal planes, which is not a coincidence: the reciprocal lattice vector $\mathbf{G}_{hkl}$ is perpendicular to the $(hkl)$ family of planes, and its magnitude is $2\pi / d_{hkl}$, where $d_{hkl}$ is the interplanar spacing.

Calculation methods

For orthogonal lattices (cubic, tetragonal, orthorhombic), the calculation is straightforward since the cross products simplify. For non-orthogonal lattices (monoclinic, triclinic), matrix methods are more practical:

1. Construct the matrix $A$ whose columns are the real-space vectors $\mathbf{a}_1, \mathbf{a}_2, \mathbf{a}_3$
2. Compute $B = 2\pi (A^{-1})^T$
3. The columns of $B$ are the reciprocal lattice vectors

Fourier transform algorithms (particularly FFT) handle the conversion between real-space functions (like electron density) and their reciprocal-space representations computationally.

Applications in solid state physics

X-ray diffraction analysis

Diffraction occurs when the scattering vector $\Delta \mathbf{k}$ equals a reciprocal lattice vector $\mathbf{G}$. This is the Laue condition, and it's equivalent to Bragg's law:

$2d\sin\theta = n\lambda$

Here $d$ is the spacing between crystal planes, $\theta$ is the scattering angle, and $\lambda$ is the X-ray wavelength. The reciprocal lattice makes this relationship geometric: a diffraction peak appears whenever the scattering geometry connects two reciprocal lattice points.

From the positions of diffraction peaks, you extract the lattice parameters and symmetry. From the intensities, you get information about atomic positions within the unit cell (through the structure factor, covered below).

Electronic band structure

Electron wavevectors in a periodic potential are only unique within one Brillouin zone. States outside the first zone can always be "folded back" by subtracting a reciprocal lattice vector, thanks to Bloch's theorem.

This is why band diagrams are plotted along high-symmetry paths in the first Brillouin zone (e.g., $\Gamma \to X \to M \to \Gamma$ for a square lattice). The reciprocal lattice framework also underpins calculations of conductivity, optical absorption, and electron-phonon coupling.

Reciprocal space vs real space

These two spaces contain the same information, just represented differently. Real space shows you where atoms sit. Reciprocal space shows you the spatial frequencies present in the structure. Switching between them (via Fourier transforms) is routine in both experiment and computation.

Visualization techniques

Reciprocal lattice points are typically plotted in 2D cross-sections or 3D projections. Point size or color often encodes diffraction intensity. Stereographic projections are useful for representing high-symmetry directions and zone axes on a 2D diagram, which helps when orienting crystals in the lab.

Interpretation of data

Reading a diffraction pattern in reciprocal space:

• Peak positions give you lattice parameters and crystal system
• Peak intensities encode atomic positions within the unit cell
• Systematic absences (missing peaks where you'd expect them) reveal the space group, since glide planes and screw axes cause specific reflections to vanish
• Diffuse scattering (broad features between Bragg peaks) indicates disorder, thermal motion, or short-range correlations

Experimental techniques

Neutron scattering

Thermal neutrons have wavelengths comparable to interatomic spacings (around 1–2 Å), making them suitable for probing crystal structures in reciprocal space. Two advantages over X-rays:

• Neutrons scatter from nuclei rather than electron clouds, so they can distinguish light elements (like hydrogen) and different isotopes
• Neutrons carry a magnetic moment, so they can probe magnetic ordering and spin-wave (magnon) dispersions

Inelastic neutron scattering maps out phonon dispersion relations across the Brillouin zone, giving direct information about lattice dynamics.

Electron diffraction

High-energy electrons (typically 100–300 keV in a TEM) produce diffraction patterns from very small sample volumes, down to nanometer-scale regions. This makes electron diffraction ideal for studying nanostructures, thin films, and local variations in crystal structure.

• Kikuchi patterns arise from inelastically scattered electrons and reveal precise crystal orientation
• RHEED (reflection high-energy electron diffraction) monitors surface structure during thin-film growth in real time

Reciprocal lattice in 2D systems

Graphene and other 2D materials

Graphene's real-space honeycomb lattice produces a hexagonal reciprocal lattice. The first Brillouin zone is also hexagonal, with high-symmetry points labeled $\Gamma$ (center), $K$ and $K'$ (corners), and $M$ (edge midpoints).

The $K$ and $K'$ points are where graphene's famous Dirac cones appear: the valence and conduction bands meet at these points with a linear dispersion, giving electrons an effective mass of zero. Stacking multiple layers (as in bilayer graphene or twisted bilayer graphene) modifies the reciprocal lattice and can produce flat bands and correlated electron phenomena.

Surface science applications

Surfaces often reconstruct, meaning the surface atoms rearrange into a periodicity different from the bulk. These reconstructions show up as additional spots in the surface reciprocal lattice.

LEED (low-energy electron diffraction) is the standard technique for probing surface reciprocal lattices. Electrons with energies of 20–200 eV have short mean free paths, so they scatter primarily from the top few atomic layers. The resulting diffraction pattern directly maps the 2D surface reciprocal lattice.

Advanced concepts

Ewald construction

The Ewald construction is a geometric tool for determining which reciprocal lattice points satisfy the diffraction condition for a given experimental setup.

1. Place the origin of reciprocal space at the point where the incident beam would end
2. Draw a sphere (the Ewald sphere) of radius $1/\lambda$ centered on the crystal, passing through the reciprocal lattice origin
3. Any reciprocal lattice point that lies on the surface of this sphere satisfies the Laue condition and produces a diffraction peak

For a perfect monochromatic beam and an ideal crystal, the sphere must pass exactly through a reciprocal lattice point. In practice, peak broadening (from finite crystal size, beam divergence, and wavelength spread) relaxes this condition. Systematic absences appear as reciprocal lattice points that the sphere intersects but that have zero structure factor.

Structure factor calculations

The structure factor determines the intensity of each diffraction peak. It accounts for how the atoms within a single unit cell contribute to scattering:

$F(\mathbf{G}) = \sum_j f_j \, e^{i\mathbf{G} \cdot \mathbf{r}_j}$

Here $f_j$ is the atomic form factor of atom $j$ (which depends on the element and the scattering angle), and $\mathbf{r}_j$ is its position within the unit cell. The measured diffraction intensity is proportional to $|F(\mathbf{G})|^2$.

A classic example: in BCC structures, the structure factor vanishes when $h + k + l$ is odd, producing systematic absences. In FCC structures, reflections vanish unless $h, k, l$ are all even or all odd. These rules follow directly from evaluating the sum above for the specific basis atom positions.

Computational methods

Simulation of reciprocal lattice

Modern computational approaches to reciprocal-space analysis include:

• Density functional theory (DFT) calculates electronic band structures and charge densities in reciprocal space, using plane-wave basis sets defined by reciprocal lattice vectors
• FFT algorithms convert between real-space grids and reciprocal-space representations efficiently, scaling as $O(N \log N)$ rather than $O(N^2)$
• Molecular dynamics simulations can generate diffuse scattering patterns by Fourier-transforming atomic displacement snapshots

Software tools for analysis

Common tools in the field:

• VESTA and CrystalMaker for visualizing crystal and reciprocal lattice structures
• Python (with NumPy/SciPy) and MATLAB for custom reciprocal-space calculations
• ICDD PDF-4 database for matching experimental diffraction patterns to known phases
• FullProf and GSAS for Rietveld refinement of diffraction data

Importance in material characterization

Crystal structure determination

Solving a crystal structure from diffraction data is essentially a problem in reciprocal space. The workflow:

1. Collect diffraction data (peak positions and intensities)
2. Index the peaks to determine the unit cell and space group
3. Use direct methods or Patterson synthesis to get initial atomic positions from the structure factors
4. Refine the model using least-squares fitting against the measured $|F(\mathbf{G})|^2$ values

Charge flipping is a more recent algorithm that can solve structures ab initio without prior assumptions about atom types, working iteratively between real and reciprocal space.

Defect analysis

Reciprocal space is also sensitive to imperfections in crystals:

• Peak broadening indicates small crystallite size (via the Scherrer equation) or microstrain
• Diffuse streaks along specific reciprocal-space directions signal stacking faults or planar defects
• Satellite peaks near Bragg reflections can indicate modulated structures or compositional ordering
• Asymmetric peak shapes may reveal dislocation distributions or strain gradients

These features, which would be difficult to interpret in real space, become systematic and quantifiable in reciprocal space.