1.7 Lattice vibrations and phonons

Fundamentals of lattice vibrations

Atoms in a crystal aren't frozen in place. They vibrate around their equilibrium positions, and because they're connected through interatomic forces, those vibrations are collective: one atom's motion affects its neighbors. These collective vibrations are the starting point for understanding phonons, which in turn explain a huge range of material properties from heat capacity to superconductivity.

Crystal lattice structure

A crystal is a periodic arrangement of atoms or molecules. The smallest repeating unit is the unit cell, defined by lattice parameters (lengths and angles). In three dimensions, there are 14 distinct Bravais lattices that capture every possible periodic geometry. Symmetry operations like translations, rotations, and reflections describe how the structure maps onto itself. This periodicity is what makes collective vibrations so well-behaved and mathematically tractable.

Vibrational modes in solids

When atoms vibrate collectively, the resulting oscillations can be classified by how the atoms move relative to the direction the wave travels:

• Longitudinal modes: atoms displace parallel to the wave propagation direction (like a compression wave on a spring)
• Transverse modes: atoms displace perpendicular to the propagation direction (like a wave on a string)

Each mode has a polarization vector describing the displacement direction. When vibrations satisfy specific boundary conditions (e.g., fixed ends of a finite crystal), they form standing waves.

Normal modes vs. phonons

Normal modes are the classical picture: each mode is a collective oscillation with a well-defined frequency $\omega$ and wave vector $\mathbf{k}$. Every possible vibration of the crystal can be decomposed into a superposition of these normal modes.

Phonons are the quantum version. Just as photons are quanta of electromagnetic waves, phonons are quanta of lattice vibrations. Each phonon carries a discrete energy:

$E = \hbar\omega$

This quantization lets you apply quantum statistical mechanics to lattice vibrations, which turns out to be essential for correctly predicting thermal properties at low temperatures where classical physics fails.

Phonon dispersion relations

The dispersion relation connects a phonon's frequency $\omega$ to its wave vector $\mathbf{k}$. It tells you how fast phonons of different wavelengths propagate and what energies are available in the crystal. Dispersion relations are central to predicting thermal conductivity, sound speeds, and optical properties.

Acoustic vs. optical phonons

If the unit cell contains more than one atom, two types of phonon branches appear:

• Acoustic phonons: neighboring atoms oscillate roughly in phase. Near the Brillouin zone center ($k \to 0$), the dispersion is linear: $\omega \propto k$. The slope gives the speed of sound in the material. These are the phonons responsible for carrying most of the heat in insulators.
• Optical phonons: neighboring atoms oscillate out of phase. At $k = 0$, optical phonons have a non-zero frequency. They're called "optical" because in ionic crystals (like NaCl), these modes create oscillating electric dipoles that can couple to infrared light.

A crystal with $p$ atoms per unit cell in 3D has 3 acoustic branches and $3(p-1)$ optical branches.

Brillouin zones

The Brillouin zone is the primitive cell of the reciprocal lattice in $\mathbf{k}$-space. The first Brillouin zone contains all the unique wave vectors you need to describe every phonon mode. Wave vectors outside this zone are equivalent to ones inside it (they differ by a reciprocal lattice vector), so the physics is periodic in $\mathbf{k}$-space.

At the Brillouin zone boundaries, phonon waves satisfy the Bragg condition and form standing waves. In superlattices or crystals with large unit cells, zone folding maps branches from extended zones back into the first zone, creating additional apparent branches.

Dispersion curves

Dispersion curves are plots of $\omega$ vs. $\mathbf{k}$ along high-symmetry directions in the Brillouin zone. Several features are worth noting:

• The slope at any point gives the phonon group velocity $v_g = d\omega/dk$
• Flat regions mean low group velocity and a high density of states (many modes packed into a narrow energy range)
• Van Hove singularities are critical points where the density of states diverges or has kinks, often at zone boundaries or saddle points
• Experimentally, dispersion curves are most directly measured by inelastic neutron scattering

Quantum theory of phonons

Treating lattice vibrations quantum mechanically is what turns normal modes into phonons. This framework is necessary for correctly describing thermal properties, phonon statistics, and interactions.

Second quantization

Second quantization rewrites the lattice Hamiltonian in terms of creation ($a_q^\dagger$) and annihilation ($a_q$) operators for each phonon mode $q$. The Hamiltonian becomes:

$H = \sum_q \hbar\omega_q \left(a_q^\dagger a_q + \frac{1}{2}\right)$

The $\frac{1}{2}$ term represents the zero-point energy: even at absolute zero, each mode retains energy $\frac{1}{2}\hbar\omega_q$. This formalism treats phonons as bosonic quasiparticles and greatly simplifies calculations of phonon-phonon and electron-phonon interactions.

Phonon creation and annihilation

• $a_q^\dagger$ creates a phonon in mode $q$ (adds one quantum of vibrational energy)
• $a_q$ annihilates a phonon in mode $q$ (removes one quantum)
• They obey bosonic commutation relations: $[a_q, a_{q'}^\dagger] = \delta_{qq'}$
• The number operator $N_q = a_q^\dagger a_q$ gives the phonon occupation of mode $q$
• A general state is specified by occupation numbers: $|n_1, n_2, \ldots, n_q\rangle$

Phonon statistics

Phonons are bosons, so they follow Bose-Einstein statistics. The average number of phonons in mode $q$ at temperature $T$ is:

$\langle n_q \rangle = \frac{1}{e^{\hbar\omega_q / k_B T} - 1}$

At low temperatures, high-frequency modes are "frozen out" ($\langle n_q \rangle \approx 0$). At high temperatures ($k_B T \gg \hbar\omega_q$), each mode approaches $\langle n_q \rangle \approx k_B T / \hbar\omega_q$, recovering the classical equipartition result.

The phonon density of states $g(\omega)$ counts the number of modes per unit frequency interval. It's the key ingredient for computing thermodynamic quantities like heat capacity.

Thermal properties of phonons

Phonons dominate the thermal behavior of solids, especially insulators and semiconductors. The two classic models for phonon heat capacity are the Debye and Einstein models, and understanding where each works (and fails) is a common exam topic.

Specific heat capacity

The phonon contribution to specific heat $C_V$ shows two limiting behaviors:

• Low temperatures: $C_V \propto T^3$ (the Debye $T^3$ law), because only long-wavelength acoustic phonons are thermally excited
• High temperatures: $C_V \to 3Nk_B$ per mole, the classical Dulong-Petit value (about 25 J/mol·K for monatomic solids)

Anomalies in $C_V$ (sudden jumps or kinks) can signal phase transitions.

Debye model

The Debye model treats the phonon spectrum as a continuous distribution with a linear dispersion $\omega = v|\mathbf{k}|$, cut off at a maximum frequency $\omega_D$. The Debye temperature is:

$\Theta_D = \frac{\hbar\omega_D}{k_B}$

The resulting specific heat is:

$C_V = 9Nk_B\left(\frac{T}{\Theta_D}\right)^3 \int_0^{\Theta_D/T} \frac{x^4 e^x}{(e^x - 1)^2}\,dx$

This correctly gives the $T^3$ behavior at low $T$ and the Dulong-Petit limit at high $T$. However, it assumes only acoustic-like modes with a single speed of sound, so it misses optical phonon contributions in crystals with a multi-atom basis.

Einstein model

The Einstein model takes the opposite extreme: all phonons vibrate at a single frequency $\omega_E$, with Einstein temperature $\Theta_E = \hbar\omega_E / k_B$. The specific heat is:

$C_V = 3Nk_B\left(\frac{\Theta_E}{T}\right)^2 \frac{e^{\Theta_E/T}}{(e^{\Theta_E/T} - 1)^2}$

This captures the exponential freeze-out of modes at low $T$ but predicts $C_V \propto e^{-\Theta_E/T}$, which drops too fast compared to the actual $T^3$ behavior. The Einstein model works well for optical phonon branches (which are relatively flat in dispersion). In practice, you can combine both models: Debye for acoustic branches, Einstein for optical branches.

Phonon interactions

In a perfectly harmonic crystal, phonons wouldn't interact at all and thermal conductivity would be infinite. Real crystals have anharmonic terms in the interatomic potential, and these give rise to phonon scattering, finite thermal conductivity, and thermal expansion.

Phonon-phonon scattering

When phonons scatter off each other, two types of processes matter:

• Normal (N) processes: total crystal momentum is conserved. These redistribute phonons among modes but don't directly resist heat flow.
• Umklapp (U) processes: the total wave vector changes by a reciprocal lattice vector $\mathbf{G}$. These do resist heat flow and are the dominant source of intrinsic thermal resistance.

At low temperatures, Umklapp processes freeze out exponentially (they require phonons with large enough $\mathbf{k}$ to reach the zone boundary), so thermal conductivity rises sharply as $T$ drops. Three-phonon processes (two phonons combining into one, or one splitting into two) dominate at moderate temperatures. Higher-order processes become relevant at elevated temperatures.

Electron-phonon coupling

Electrons in a metal or semiconductor interact with lattice vibrations. This coupling is responsible for:

• Electrical resistivity in metals (phonon scattering limits electron mean free path)
• Conventional superconductivity (phonon exchange binds electrons into Cooper pairs)
• Polaron formation (an electron "dresses" itself with a cloud of virtual phonons, increasing its effective mass)

The interaction is described by the Fröhlich Hamiltonian, and the coupling strength is characterized by a dimensionless parameter $\lambda$. Larger $\lambda$ means stronger coupling and, in BCS theory, a higher superconducting transition temperature.

Anharmonic effects

The harmonic approximation assumes the interatomic potential is perfectly parabolic. Real potentials are anharmonic (asymmetric), which leads to:

• Thermal expansion: the average atomic spacing increases with temperature because the potential well is asymmetric
• Finite phonon lifetimes: phonons decay through scattering, giving spectral lines a finite width
• Phonon frequency shifts with temperature and pressure, quantified by the Grüneisen parameter $\gamma = -\frac{V}{\omega}\frac{\partial \omega}{\partial V}$

Anharmonicity is also essential for understanding phase transitions and nonlinear optical responses.

Experimental techniques

Three major techniques probe phonon properties, each with different strengths.

Neutron scattering

Inelastic neutron scattering is the most direct way to measure phonon dispersion relations across the entire Brillouin zone. Neutrons have wavelengths comparable to interatomic spacings and energies comparable to phonon energies, making them ideal probes.

• Energy and momentum conservation in the scattering event determine the phonon's $\omega$ and $\mathbf{k}$
• Common instruments: triple-axis spectrometers (scan specific points in $(\mathbf{k}, \omega)$ space) and time-of-flight spectrometers (measure broad regions simultaneously)
• Can probe both acoustic and optical branches

The main limitation is that neutron sources (reactors or spallation sources) are large-scale facilities.

Raman spectroscopy

Raman spectroscopy measures optical phonons near the zone center ($k \approx 0$) by detecting inelastically scattered laser light.

• Stokes peaks: the photon loses energy by creating a phonon
• Anti-Stokes peaks: the photon gains energy by absorbing a phonon
• Selection rules based on crystal symmetry determine which modes are Raman-active
• Temperature-dependent measurements reveal anharmonic effects (peak shifts and broadening)

Raman is a lab-scale technique, making it far more accessible than neutron scattering, though it only probes a small region of $\mathbf{k}$-space.

X-ray diffraction

While primarily used for crystal structure determination, X-ray techniques also provide phonon information:

• Thermal diffuse scattering (the diffuse background around Bragg peaks) reflects phonon populations
• Inelastic X-ray scattering (IXS) at synchrotron sources can measure dispersion relations with very high momentum resolution
• Useful for small samples where neutron scattering isn't practical

Applications of phonons

Thermal conductivity

In insulators and semiconductors, phonons are the primary heat carriers. Thermal conductivity $\kappa$ depends on phonon specific heat, group velocity, and mean free path:

$\kappa \approx \frac{1}{3} C_V v \ell$

where $v$ is an average phonon velocity and $\ell$ is the mean free path. Phonon engineering tailors $\kappa$ for specific applications:

• Phononic crystals create band gaps that block certain phonon frequencies, enabling precise thermal management
• Thermoelectric materials need low $\kappa$ (to maintain a temperature gradient) but high electrical conductivity, so the goal is to scatter phonons without scattering electrons

Superconductivity

In conventional (BCS) superconductors, phonons mediate the attractive interaction between electrons that forms Cooper pairs. Evidence for this phonon mechanism includes the isotope effect: replacing atoms with heavier isotopes lowers the superconducting transition temperature $T_c$, because heavier atoms vibrate at lower frequencies.

Strong electron-phonon coupling produces higher $T_c$ values. $\text{MgB}_2$ ($T_c \approx 39$ K) is a well-known example of a phonon-mediated superconductor with unusually strong coupling. In unconventional superconductors (like cuprates), other pairing mechanisms may compete with or replace phonon-mediated pairing.

Thermoelectric materials

Thermoelectric devices convert temperature differences into voltage (or vice versa). Their efficiency is governed by the figure of merit:

$ZT = \frac{S^2 \sigma T}{\kappa}$

where $S$ is the Seebeck coefficient, $\sigma$ is electrical conductivity, and $\kappa$ is thermal conductivity. The "phonon glass, electron crystal" concept captures the ideal: scatter phonons aggressively (low $\kappa$) while letting electrons move freely (high $\sigma$). Strategies include alloying, nanostructuring, and incorporating "rattler" atoms in cage-like crystal structures.

Phonons in low-dimensional systems

When one or more dimensions of a material shrink to the nanoscale, phonon behavior changes qualitatively due to confinement and surface effects.

Nanostructures

• Quantum confinement modifies phonon dispersion and density of states, potentially creating discrete phonon energy levels in quantum dots
• Surface and interface phonons become significant as the surface-to-volume ratio grows
• When the structure is smaller than the phonon mean free path, transport becomes ballistic rather than diffusive
• The phonon bottleneck effect in quantum dots slows carrier relaxation because the discrete energy spacing limits which phonons can participate

Two-dimensional materials

Graphene and similar 2D materials have distinctive phonon physics:

• Flexural phonons (out-of-plane vibrations with $\omega \propto k^2$ dispersion) are unique to 2D and contribute anomalously to thermal conductivity
• Phonon-limited carrier mobility in 2D differs from bulk because of the changed dimensionality and screening
• In van der Waals heterostructures, weak interlayer coupling introduces low-frequency interlayer breathing and shear modes

Phononic crystals

Phononic crystals are artificially structured materials designed to control phonon propagation, analogous to how photonic crystals control light.

• Periodic variations in density or elastic properties create phononic band gaps where certain frequencies cannot propagate
• Applications include acoustic cloaking, thermal diodes, phonon waveguides, and emerging quantum phonon circuits
• Acoustic metamaterials extend these ideas to achieve properties not found in natural materials (e.g., negative effective mass density)