Phonons

A phonon is one quantum of a stable vibrational normal mode of a material. In a crystal, the modes are collective displacement patterns labeled by wave vector $\mathbf q$ , branch $\nu$ , and polarization. Their frequencies $\omega_{\nu}(\mathbf q)$ form a phonon dispersion, and their eigenvectors determine which atoms move, in which directions, and with what relative phases.

The word phonon does three jobs in materials physics:

it names the harmonic quanta used to count vibrational energy;
it labels sharp, renormalized vibrational peaks when interactions remain weak enough;
it provides the kinematic language for heat capacity, thermal transport, structural instabilities, and scattering experiments.

Those uses are related but not identical. A calculated harmonic mode can be invisible in a chosen probe. A measured broad feature need not define a long-lived phonon. An imaginary harmonic frequency is not a negative-energy particle.

Phonons as Many-Body Excitations is the canonical home for the universal derivation from coupled coordinates to bosonic operators, one-phonon states, coherent lattice waves, and phonon propagators. This page owns the material-facing lattice-dynamics problem: force constants, real-crystal dispersions and polarizations, phonon density of states, thermodynamics, probe-dependent spectra, computational workflows, and stability interpretation.

Required background. Crystals and Lattices supplies primitive cells, basis atoms, and equilibrium positions. Reciprocal Lattice supplies reciprocal vectors and the Fourier-phase language used in the dynamical matrix. Brillouin Zones supplies reduced-wave-vector zones and reciprocal-space paths used to read a phonon dispersion.

Helpful background. Phonons as Many-Body Excitations supplies normal-mode quantization and the bosonic creation and annihilation operators used here without rederivation.

From a Crystal Structure to a Spectrum

Let $\mathbf R_l$ label primitive cells and $\boldsymbol\tau_\kappa$ the equilibrium position of basis atom $\kappa$ within a cell. Its position is

\mathbf r_{l\kappa} = \mathbf R_l + \boldsymbol\tau_\kappa + \mathbf u_{l\kappa},

where $\mathbf u_{l\kappa}$ is a displacement from the chosen reference structure.

Within a clamped-electron or Born–Oppenheimer description, expand the potential-energy surface:

\begin{aligned} V ={}& V_0 + \frac{1}{2} \sum_{\substack{l\kappa\alpha\\l'\kappa'\beta}} u_{l\kappa\alpha} \Phi_{l\kappa\alpha,l'\kappa'\beta} u_{l'\kappa'\beta} \\ &+ V^{(3)} + V^{(4)} + \cdots. \end{aligned}

The linear term vanishes at a stationary reference structure. The harmonic interatomic force constants are

\Phi_{l\kappa\alpha,l'\kappa'\beta} = \frac{\partial^2V} {\partial u_{l\kappa\alpha} \partial u_{l'\kappa'\beta}} \biggr\rvert_{\mathbf u=0}.

They are a real-space Hessian, not merely pairwise spring constants. Electronic screening, chemical bonding, long-range electrostatics, and the chosen electronic-structure approximation are all encoded in them.

Dynamical matrix

This page uses a Fourier phase containing the full equilibrium atomic position. Translation symmetry lets one fix the first cell at $l=0$ and define

\begin{aligned} D_{\kappa\alpha,\kappa'\beta}(\mathbf q) ={}& \frac{1} {\sqrt{M_\kappa M_{\kappa'}}} \sum_l \Phi_{0\kappa\alpha,l\kappa'\beta} \\ &\times \exp \left[ i\mathbf q\cdot \left( \mathbf R_l + \boldsymbol\tau_{\kappa'} - \boldsymbol\tau_\kappa \right) \right]. \end{aligned}

Some codes omit the basis positions from the Fourier phase and move the corresponding phases into the eigenvectors and structure factors. Frequencies are unchanged, but raw polarization vectors from the two conventions cannot be compared component by component without conversion.

The mass-weighted eigenproblem is

\sum_{\kappa'\beta} D_{\kappa\alpha,\kappa'\beta}(\mathbf q) e_{\kappa'\beta,\nu}(\mathbf q) = \omega_{\nu}^2(\mathbf q) e_{\kappa\alpha,\nu}(\mathbf q).

For a Hermitian dynamical matrix, the eigenvectors can be normalized as

\sum_{\kappa\alpha} e_{\kappa\alpha,\nu}^*(\mathbf q) e_{\kappa\alpha,\nu'}(\mathbf q) = \delta_{\nu\nu'}.

The eigenvalue is $\omega^2$ , not $\omega$ . A negative eigenvalue therefore gives exponential growth in the harmonic equations of motion. Plotting such a mode at a negative ordinate is a visualization convention for an imaginary frequency.

Translation invariance

A rigid translation cannot change the energy of a free bulk crystal. The force constants must obey the acoustic sum rule

\sum_{l\kappa'} \Phi_{0\kappa\alpha,l\kappa'\beta} = 0.

Consequently, an ordinary stable three-dimensional crystal has three zero-frequency translation modes at $\mathbf q=0$ . Small nonzero values in a calculation can come from finite convergence, interpolation, or imperfect enforcement of invariance. Large violations are evidence of a defective force-constant model, not a harmless plotting issue.

Rotational invariance imposes additional constraints. They are especially important for low-dimensional materials, where a freestanding membrane has a flexural acoustic branch that is quadratic at sufficiently long wavelength in the ideal harmonic continuum.

Quantization Without Repeating the Derivation

For every stable nonzero normal mode, canonical quantization gives

H_{\mathrm{ph}} = \sum_{\mathbf q\nu} \hbar\omega_{\nu}(\mathbf q) \left( b_{\mathbf q\nu}^{\dagger} b_{\mathbf q\nu} + \frac{1}{2} \right).

The displacement operator is

\begin{aligned} \hat u_{l\kappa\alpha} ={}& \frac{1} {\sqrt{N_cM_\kappa}} \sum_{\mathbf q\nu} \sqrt{ \frac{\hbar} {2\omega_\nu(\mathbf q)} } e_{\kappa\alpha,\nu}(\mathbf q) \\ &\times e^{i\mathbf q\cdot (\mathbf R_l+\boldsymbol\tau_\kappa)} \left( b_{\mathbf q\nu} + b_{-\mathbf q\nu}^{\dagger} \right), \end{aligned}

up to the paired phase convention for the polarization vectors.

A one-phonon state is obtained by applying $b_{\mathbf q\nu}^{\dagger}$ . It is not one atom oscillating, and it need not have nonzero mean displacement. A classical-looking lattice wave is generally a coherent state involving many occupation numbers.

Phonons carry energy $\hbar\omega_\nu(\mathbf q)$ and crystal momentum $\hbar\mathbf q$ modulo a reciprocal vector. Their number is not generally conserved: cubic and higher anharmonic interactions can create and annihilate phonons, and equilibrium phonons ordinarily have zero chemical potential.

Reading a Phonon Dispersion

A dispersion plot samples $\omega_\nu(\mathbf q)$ along a chosen path through the Brillouin zone. It is analogous to an electronic band plot, but the eigenvalues are squared frequencies and the eigenvectors are mass-weighted displacement patterns.

Branch counting

With $r$ basis atoms in a three-dimensional bulk crystal, there are

3r

branches at each $\mathbf q$ . Near the zone center, a stable free crystal ordinarily has:

\begin{aligned} 3 &\quad \text{acoustic branches}, \\ 3r-3 &\quad \text{optical branches}. \end{aligned}

The names refer to long-wavelength displacement character, not directly to audibility or optical brightness. A monatomic primitive cell has no optical branches in the harmonic branch-counting sense. A two-dimensional monolayer embedded in three-dimensional space still has three displacement directions per atom, even though $\mathbf q$ is two-dimensional.

Acoustic branches and sound

For a regular acoustic branch at small wave vector,

\omega_\nu(\mathbf q) \approx v_\nu(\hat{\mathbf q}) |\mathbf q|.

The slope is a direction- and polarization-dependent sound velocity. In continuum elasticity, the Christoffel matrix

\Gamma_{ij}(\hat{\mathbf n}) = \frac{1}{\rho} C_{ikjl} n_k n_l

has eigenvalues $v_\nu^2(\hat{\mathbf n})$ . Here $C_{ikjl}$ is the elastic tensor, $\rho$ is mass density, and $\hat{\mathbf n}$ is the propagation direction.

Longitudinal and transverse labels are exact in isotropic media and along suitable symmetry directions. At a generic wave vector in a low-symmetry crystal, a mode can have mixed polarization. The eigenvector, not the branch color chosen by plotting software, decides.

Optical branches

At a zone-center optical mode, basis atoms move against one another while the cell’s center-of-mass translation cancels. The frequency can remain finite as $\mathbf q\to0$ .

An optical branch is:

infrared active if its mode effective charge is nonzero;
Raman active if the mode changes the polarizability tensor in an allowed symmetry channel;
silent in a given first-order optical probe if neither selection rule is met.

These are symmetry and matrix-element statements. “Optical phonon” does not mean “visible in every optical spectrum.”

Polar crystals and directional limits

In a polar insulator, a longitudinal optical displacement creates a macroscopic electric field. The long-range Coulomb interaction adds a nonanalytic term as $\mathbf q\to0$ . If the Born effective charges $\mathsf Z_\kappa^*$ are expressed in units of $e$ , a representative SI form is

\begin{aligned} D_{\kappa\alpha,\kappa'\beta}^{\mathrm{NA}} (\mathbf q) ={}& \frac{e^2} {\epsilon_0\Omega_c \sqrt{M_\kappa M_{\kappa'}}} \\ &\times \frac{ \left( \mathbf q\cdot\mathsf Z_\kappa^* \right)_\alpha \left( \mathbf q\cdot\mathsf Z_{\kappa'}^* \right)_\beta }{ \mathbf q\cdot \boldsymbol\epsilon_\infty \mathbf q }. \end{aligned}

The limit depends on the direction of approach. This produces longitudinal–transverse optical splitting and means that “the frequency at $\Gamma$ ” can require a boundary condition and direction, not just one number.

Crossings, avoided crossings, and branch identity

Eigenvalues can cross when symmetry forbids mixing. Modes with compatible symmetry generally repel and exchange character through an avoided crossing. Sorting branches only by frequency can then create discontinuous colors or labels.

Robust mode tracking uses some combination of:

eigenvector overlap between neighboring $\mathbf q$ points;
irreducible representations at high-symmetry points and lines;
projections onto atoms, directions, or molecular units;
continuity of experimentally relevant structure factors.

At a degeneracy, individual eigenvectors may rotate arbitrarily inside the degenerate subspace. The subspace is meaningful even when a particular polarization vector is not.

Phonon Density of States

For $N_q$ wave vectors sampling a primitive Brillouin zone, define the density of states per primitive cell by

g(\omega) = \frac{1}{N_q} \sum_{\mathbf q\nu} \delta \left[ \omega-\omega_\nu(\mathbf q) \right].

For a three-dimensional crystal with $r$ basis atoms,

\int_0^\infty d\omega\, g(\omega) = 3r.

This normalization counts modes per cell per angular-frequency interval. A density per physical volume is $g_V(\omega)=g(\omega)/\Omega_c$ . If energy $E=\hbar\omega$ is used instead, the Jacobian must be included:

g_E(E) = \frac{1}{\hbar} g_\omega \left( \frac{E}{\hbar} \right).

Constant-frequency surfaces

In the continuum limit,

g(\omega) = \frac{\Omega_c}{(2\pi)^3} \sum_\nu \int_{\mathrm{BZ}} d^3q\, \delta \left[ \omega-\omega_\nu(\mathbf q) \right].

At a regular frequency,

g(\omega) = \frac{\Omega_c}{(2\pi)^3} \sum_\nu \int_{\mathcal S_{\nu,\omega}} \frac{dS} { \left| \nabla_{\mathbf q} \omega_\nu \right| }.

Large DOS appears where constant-frequency surfaces have large area or where group velocity is small. Critical points of $\omega_\nu(\mathbf q)$ produce van Hove features. A sharp DOS peak does not by itself imply a localized vibration or a long lifetime.

Projected and weighted densities

For mass-weighted eigenvectors, an atom-projected density can be defined as

g_\kappa(\omega) = \frac{1}{N_q} \sum_{\mathbf q\nu} \sum_\alpha \left| e_{\kappa\alpha,\nu}(\mathbf q) \right|^2 \delta \left[ \omega-\omega_\nu(\mathbf q) \right].

With the stated normalization,

\sum_\kappa g_\kappa(\omega) = g(\omega).

Other projections can resolve Cartesian direction, chemical species, layer, or molecular unit. Their meaning depends on the eigenvector convention and projector.

An experimentally extracted “phonon DOS” is often a generalized or probe-weighted density. In incoherent neutron measurements, species weights can scale roughly as scattering cross section divided by mass. Multiphonon contributions, Debye–Waller factors, detector response, and instrumental resolution must be removed or modeled before comparing with the unweighted theoretical $g(\omega)$ .

Low-frequency laws

For three linearly dispersing acoustic branches in a three-dimensional continuum,

g_V(\omega) = \frac{\omega^2}{2\pi^2} \sum_{s=1}^{3} \frac{1}{v_s^3}.

In an isotropic solid this becomes

g_V(\omega) = \frac{\omega^2}{2\pi^2} \left( \frac{1}{v_L^3} + \frac{2}{v_T^3} \right).

The $\omega^2$ law is not universal across dimension or dispersion. A two-dimensional linear branch gives $g(\omega)\propto\omega$ , while an ideal quadratic flexural branch gives a constant low-frequency contribution.

Vibrational Thermodynamics

For harmonic phonons,

n_{\mathrm B}(\omega) = \frac{1} {e^{\beta\hbar\omega}-1}.

The phonon internal energy is

U_{\mathrm{ph}} = \sum_{\mathbf q\nu} \hbar\omega_{\nu}(\mathbf q) \left[ n_{\mathrm B} \left( \omega_\nu \right) + \frac{1}{2} \right].

The free energy is

\begin{aligned} F_{\mathrm{ph}} = \sum_{\mathbf q\nu} \biggl\{ &\frac{\hbar\omega_\nu}{2} \\ &+ k_{\mathrm B}T \ln \left[ 1-e^{-\beta\hbar\omega_\nu} \right] \biggr\}. \end{aligned}

The mode heat capacity is

C_{\mathbf q\nu} = k_{\mathrm B} \frac{x_{\mathbf q\nu}^2e^{x_{\mathbf q\nu}}} {\left(e^{x_{\mathbf q\nu}}-1\right)^2}, \qquad x_{\mathbf q\nu} = \frac{\hbar\omega_{\mathbf q\nu}} {k_{\mathrm B}T}.

Summing over all modes gives $C_V$ in the fixed-frequency harmonic approximation.

Low and high temperature

For a three-dimensional crystal with linear acoustic modes,

\frac{C_V}{V} \approx \frac{2\pi^2k_{\mathrm B}^4T^3} {15\hbar^3} \sum_{s=1}^{3} \frac{1}{v_s^3}.

This is the Debye $T^3$ law. Optical modes and high-frequency acoustic states are exponentially suppressed at sufficiently low temperature.

At temperatures large compared with all relevant $\hbar\omega/k_{\mathrm B}$ ,

C_V \longrightarrow 3r k_{\mathrm B}

per primitive cell in a three-dimensional harmonic crystal. This is the Dulong–Petit limit. Electronic, magnetic, defect, and anharmonic contributions can prevent measured heat capacity from being purely phononic.

Zero-point energy affects total energies, isotope shifts, equilibrium structures, and phase competition, but it does not contribute to harmonic heat capacity when the frequencies are temperature independent.

Volume dependence and thermal expansion

The quasiharmonic approximation lets frequencies depend on volume while retaining harmonic occupation at each volume. Define the mode Grüneisen parameter

\gamma_{\mathbf q\nu} = - \frac{\partial \ln\omega_{\mathbf q\nu}} {\partial\ln V}.

The vibrational contribution to volumetric expansion is schematically

\alpha_V = \frac{1} {B_TV} \sum_{\mathbf q\nu} \gamma_{\mathbf q\nu} C_{\mathbf q\nu},

where $B_T$ is the isothermal bulk modulus. Negative-Grüneisen modes can drive negative thermal expansion over a temperature range.

Quasiharmonicity includes thermal expansion from volume-dependent harmonic spectra. It does not include every explicit phonon–phonon frequency shift or linewidth at fixed volume.

Thermal transport is not fixed by the DOS

In a relaxation-time description,

\kappa_{ij} = \frac{1}{V} \sum_{\mathbf q\nu} C_{\mathbf q\nu} v_{\mathbf q\nu,i} v_{\mathbf q\nu,j} \tau_{\mathbf q\nu}.

The DOS supplies heat-capacity weighting, but conductivity also requires group velocities and lifetimes. Two materials can have similar phonon DOS and very different thermal conductivity.

How Scattering Sees Phonons

In a one-phonon event, the probe transfers momentum $\mathbf Q$ and angular frequency $\Omega$ . Crystal momentum conservation gives

\mathbf Q = \mathbf G+\mathbf q,

while energy conservation gives

\Omega = \pm\omega_\nu(\mathbf q).

Here positive $\Omega$ means energy lost by the probe and gained by the sample. The plus sign creates a phonon; the minus sign annihilates a thermally occupied phonon.

Neutron and nonresonant x-ray scattering

For coherent one-phonon neutron scattering, a representative mode structure factor is

\begin{aligned} F_\nu(\mathbf Q) = \sum_\kappa &\frac{ b_\kappa e^{-W_\kappa(\mathbf Q)} }{ \sqrt{M_\kappa} } e^{i\mathbf Q\cdot\boldsymbol\tau_\kappa} \\ &\times \mathbf Q\cdot \mathbf e_{\kappa\nu}(\mathbf q). \end{aligned}

$b_\kappa$ is the nuclear scattering length and $W_\kappa$ the Debye–Waller exponent. Omitting overall kinematic prefactors, the harmonic one-phonon contribution has the form

\begin{aligned} S^{(1)}(\mathbf Q,\Omega) \propto \sum_\nu \frac{|F_\nu(\mathbf Q)|^2} {2\omega_\nu(\mathbf q)} \biggl\{ & \left[ n_\nu+1 \right] \delta \left( \Omega-\omega_\nu \right) \\ &+ n_\nu \delta \left( \Omega+\omega_\nu \right) \biggr\}. \end{aligned}

This equation explains several experimental facts:

the creation side remains at zero temperature through $n+1$ ;
the annihilation side vanishes as $n\to0$ ;
modes polarized perpendicular to $\mathbf Q$ can be dark through $\mathbf Q\cdot\mathbf e$ ;
basis atoms can interfere constructively or destructively;
heavier atoms have smaller displacement amplitudes for the same mass-weighted eigenvector;
intensity changes between reciprocal zones even though $\omega_\nu(\mathbf q)$ is periodic.

Nonresonant inelastic x-ray scattering has closely related displacement and interference factors, with electronic form factors replacing neutron scattering lengths and different instrumental kinematics. Neutrons are often favorable for light atoms and isotope contrast; hard x-rays can access very small samples and extreme sample environments. Neither probe displays “the dispersion” without matrix-element and resolution effects.

Neutron Scattering owns the neutron beam, nuclear and magnetic channel separation, instrument coverage, resolution, background, absolute normalization, and reciprocal-zone measurement strategy.

X-Ray Scattering owns hard-X-ray momentum and energy kinematics, form factors, attenuation, analyzer resolution, and the nonresonant IXS forward model.

Detailed balance requires

S(\mathbf Q,-\Omega) = e^{-\beta\hbar\Omega} S(-\mathbf Q,\Omega).

When symmetry makes the $\mathbf Q$ and $-\mathbf Q$ responses equivalent, the ratio of annihilation to creation intensity is $e^{-\beta\hbar\omega}$ .

Raman, infrared, and Brillouin probes

Visible and infrared photons carry little crystal momentum compared with a Brillouin-zone dimension. First-order Raman and infrared spectra therefore probe modes near $\mathbf q=0$ .

Infrared absorption couples to a changing dipole and requires a nonzero mode effective charge.
Raman scattering couples to a changing polarizability tensor and obeys polarization-dependent symmetry selection rules.
Brillouin light scattering resolves long-wavelength acoustic modes and elastic properties.
Ultrafast pump–probe measurements can launch coherent zone-center or finite-wave-vector lattice motion, depending on the excitation mechanism and geometry.

A Raman peak is not automatically an infrared peak, and a zone-boundary phonon is not ordinarily visible in first-order optical spectroscopy unless disorder, a superstructure, resonance, or a multiphonon process relaxes momentum selection.

Frequencies, linewidths, and spectral weight

An ideal harmonic mode gives a delta-function line. Real spectra are broadened by:

anharmonic phonon–phonon decay;
isotope and defect disorder;
boundaries and finite domains;
electron–phonon or spin–phonon coupling;
hybridization with other collective modes;
instrumental energy and momentum resolution.

The observed peak position can differ from the harmonic frequency, and its width can depend on temperature and wave vector. A lifetime inferred from a width must state whether the width is a half width or full width and whether the lifetime refers to amplitude or population decay. When the width is comparable to the mode frequency or neighboring separation, a sharp-phonon description can fail.

Structure Factors owns exact correlation-function and detailed-balance conventions. Spectral Functions owns pole, continuum, linewidth, residue, and sum-rule interpretation.

Computing Phonons in Materials

Two common first-principles routes are:

finite displacement: displace atoms in a supercell, calculate forces, infer real-space force constants, and Fourier interpolate;
density-functional perturbation theory: calculate the linear response to periodic displacements directly at selected $\mathbf q$ points.

Both routes approximate the electronic ground-state energy surface and require convergence tests. Agreement between software packages is not meaningful unless structure, pseudopotentials, exchange-correlation approximation, magnetic order, spin–orbit treatment, supercell or $\mathbf q$ mesh, and nonanalytic corrections are aligned.

A reproducible workflow

Relax the intended structure and state pressure, strain, magnetic order, charge state, and symmetry constraints.
Converge electronic basis, Brillouin-zone sampling, smearing, and self-consistency more tightly than needed for a rough total energy.
Calculate force constants by finite displacement or linear response.
Check permutation symmetry, acoustic sum rules, and, where important, rotational invariance.
Include long-range dipole corrections for polar materials using Born charges and $\boldsymbol\epsilon_\infty$ .
Interpolate to a dense mesh for the DOS and to a declared path for visualization.
Test supercell size, displacement amplitude, $\mathbf q$ mesh, and integration broadening.
Compare frequencies and intensities with the geometry and resolution of the relevant experiment.
Add quasiharmonic, anharmonic, electronic, magnetic, or disorder effects only when the observable requires them.

Finite differences need a window

Displacements that are too small amplify force noise. Displacements that are too large mix in anharmonicity. A reliable finite-displacement calculation repeats the extraction over a reasonable amplitude range and checks that the force constants and frequencies are stable.

A path is not the Brillouin zone

A standard high-symmetry path can miss:

an instability at a generic wave vector;
a narrow avoided crossing;
anisotropic sound velocity away from displayed directions;
a DOS feature generated by a broad off-path surface.

Stability and thermodynamics require a converged full-zone mesh, not only an attractive line plot.

Harmonic stability is conditional

If every harmonic $\omega^2$ is nonnegative, the reference structure is a local minimum within the modeled degrees of freedom. This is dynamical stability, not proof of global thermodynamic stability. Another structure may have lower free energy.

Conversely, an imaginary zero-temperature harmonic mode can be stabilized at finite temperature by anharmonic fluctuations, disorder, or entropy. The correct response is to examine the eigenvector, convergence, reference symmetry, and temperature-dependent effective potential rather than deleting the point or replacing $\omega$ by $|\omega|$ .

Interpreting Characteristic Features

Soft mode

A branch whose renormalized frequency decreases toward zero can signal proximity to a structural instability. The eigenvector identifies the distortion pattern and the soft wave vector identifies the translation symmetry of the emerging structure. Not every structural transition is a simple displacive soft-mode transition.

Kohn anomaly

Electron screening can produce a cusp or rapid softening at wave vectors connecting important electronic states. Such a Kohn anomaly is evidence of momentum-dependent electron–phonon response, but nesting alone does not determine a transition. Electronic matrix elements, dimensionality, interactions, and phonon self-energy all matter.

Avoided crossing

When two compatible modes hybridize, their frequencies repel and their atomic or polarization character exchanges. Fitting each side as an unrelated branch can create false discontinuities in group velocity or coupling.

Acoustic–optical gap

A frequency gap between acoustic and optical sectors can restrict three-phonon decay phase space and alter thermal transport. It does not by itself guarantee high thermal conductivity; velocities, anharmonic vertices, disorder, and boundaries remain decisive.

Flat branch

A flat branch has small group velocity and often contributes a DOS peak. It can arise from weakly coupled local units, heavy atoms, geometric constraints, or hybridization. Flatness alone does not establish localization.

A Monatomic-Chain Reading Example

The Harmonic Chain Model Dossier derives the periodic nearest-neighbor chain. Its dispersion is

\omega(q) = 2\sqrt{\frac{K}{M}} \left| \sin\frac{qa}{2} \right|.

This one formula illustrates how to read a material dispersion:

there is one branch because the primitive cell has one scalar displacement;
$\omega(0)=0$ follows from translation invariance;
the sound speed is

v_s = a\sqrt{\frac{K}{M}};

the group velocity vanishes at the zone boundary;
the one-dimensional DOS per cell is

g(\omega) = \frac{2} {\pi \sqrt{\omega_{\max}^2-\omega^2}}, \qquad 0<\omega<\omega_{\max},

where

\omega_{\max} = 2\sqrt{\frac{K}{M}};

the square-root divergence at $\omega_{\max}$ is a one-dimensional van Hove singularity caused by the flat zone-boundary dispersion.

The model is a benchmark, not a generic material fit. Real crystals have vector displacements, basis atoms, longer-range force constants, anharmonicity, and probe-dependent intensities.

Common Mistakes

Treating a phonon as one atom moving independently.
Calling every normal mode a phonon without checking harmonic or quasiparticle stability.
Confusing $\omega$ , ordinary frequency $\nu$ , energy $\hbar\omega$ , and spectroscopic wavenumber.
Comparing polarization vectors from different Fourier or mass-normalization conventions directly.
Assuming every optical branch is infrared or Raman active.
Reading longitudinal or transverse labels from a plot without checking the eigenvector.
Treating a high-symmetry path as a full stability test.
Replacing imaginary frequencies by absolute values and continuing as if the structure were stable.
Comparing an unweighted theoretical DOS with a neutron- or x-ray-weighted spectrum.
Inferring the absence of a branch from zero intensity in one scattering geometry.
Calling numerical broadening a physical linewidth.
Using the Dulong–Petit limit or Debye $T^3$ law outside its temperature regime.
Inferring thermal conductivity from the phonon DOS alone.

Exercises

1. Branch and state-counting audit

A three-dimensional crystal has five atoms in its primitive cell. How many phonon branches occur at each $\mathbf q$ ? How many are acoustic near $\Gamma$ , how many are optical, and what is the integral of the DOS per primitive cell?

Solution

There are three displacement coordinates per atom, so

3r = 3\times5 = 15

branches. An ordinary stable bulk crystal has three acoustic branches associated with rigid translations. The remaining

15-3 = 12

are optical branches near $\Gamma$ . With the DOS normalization used here,

\int_0^\infty d\omega\, g(\omega) = 15.

2. Normalize the chain density of states

Show that

g(\omega) = \frac{2} {\pi\sqrt{\omega_{\max}^2-\omega^2}}

contains one mode per cell over $0<\omega<\omega_{\max}$ .

Solution

Integrate using $\omega=\omega_{\max}\sin\theta$ :

\begin{aligned} \int_0^{\omega_{\max}} d\omega\, g(\omega) &= \frac{2}{\pi} \int_0^{\pi/2} d\theta \\ &= 1. \end{aligned}

The divergence at the upper edge is integrable. It produces a sharp feature without creating extra states.

3. Which sound branch dominates low-temperature heat capacity?

An isotropic solid has $v_L=2v_T$ . Compare the longitudinal contribution to the combined contribution of the two transverse branches in the Debye $T^3$ coefficient.

Solution

Each branch contributes in proportion to $v^{-3}$ . The longitudinal weight is

\frac{1}{v_L^3} = \frac{1}{8v_T^3}.

The two transverse branches contribute

\frac{2}{v_T^3}.

Their ratio is

\frac{ 2/v_T^3 }{ 1/(8v_T^3) } = 16.

The slower transverse branches dominate because low-frequency state counting scales as $v^{-3}$ .

4. A dark optical mode

Two identical atoms in a primitive cell have an optical eigenvector with equal magnitudes and opposite directions. At a scattering vector where their basis phases are equal, assume equal $b_\kappa/\sqrt{M_\kappa}$ . What happens to the one-phonon structure factor?

Solution

The two terms are equal in magnitude. Their opposite polarizations give opposite signs in $\mathbf Q\cdot\mathbf e_{\kappa\nu}$ , while the basis phases and species weights are equal. Therefore they cancel:

F_\nu(\mathbf Q) = 0.

The mode is dark in that geometry even though it exists. Moving to a reciprocal zone with a different basis phase, changing polarization geometry, or using another probe can restore intensity.

5. Creation and annihilation intensity

For a sharp mode of frequency $\omega$ , find the ratio of the creation-side intensity to the annihilation-side intensity at temperature $T$ , assuming all other structure factors are equal.

Solution

The creation side is proportional to $n_{\mathrm B}+1$ , while the annihilation side is proportional to $n_{\mathrm B}$ . Since

\frac{n_{\mathrm B}+1} {n_{\mathrm B}} = e^{\beta\hbar\omega},

the ratio is

\frac{I_{\mathrm{create}}} {I_{\mathrm{annihilate}}} = e^{\beta\hbar\omega}.

At low temperature, annihilation is strongly suppressed because few thermal phonons are available.

6. Negative thermal expansion from one mode

A mode frequency increases as the crystal volume increases. What is the sign of its Grüneisen parameter, and what sign does its weighted contribution make to $\alpha_V$ ?

Solution

Because

\frac{\partial\ln\omega} {\partial\ln V} > 0,

the definition gives

\gamma = - \frac{\partial\ln\omega} {\partial\ln V} < 0.

Since the mode heat capacity is positive, its term $\gamma C$ contributes negatively to the quasiharmonic thermal-expansion sum. Whether the total expansion is negative depends on all modes and the elastic response.

7. Audit an imaginary branch

A first-principles calculation shows a small imaginary acoustic frequency near $\Gamma$ and a large imaginary optical branch at a generic wave vector. Give a defensible interpretation workflow.

Solution

First separate likely numerical drift from a physical instability:

check relaxation forces, stress, electronic convergence, and the intended magnetic or charge state;
enforce or diagnose the acoustic sum rule and rotational invariance;
converge supercell size or $\mathbf q$ mesh, displacement amplitude, basis, and Brillouin-zone sampling;
inspect whether the small near- $\Gamma$ acoustic error shrinks under those tests;
inspect the eigenvector of the large generic- $\mathbf q$ mode and build a commensurate supercell if possible;
displace along that eigenvector and map the energy surface rather than replacing $\omega$ by $|\omega|$ ;
compare lower-symmetry relaxed structures and, when relevant, finite-temperature anharmonic stabilization.

A tiny sum-rule-sensitive acoustic artifact and a robust large optical instability should not be assigned the same physical meaning.

Connections

Phonons as Many-Body Excitations derives quantization, one-phonon states, coherent waves, propagators, and anharmonic number-changing processes.
Harmonic Chain Model Dossier supplies the exact periodic-chain benchmark and finite-size audits.
Crystals and Lattices defines cells and basis atoms; Reciprocal Lattice and Brillouin Zones define the $\mathbf q$ geometry.
Quantum Matter Conventions fixes angular-frequency, energy, Fourier-phase, cell, and reciprocal-space conventions.
Quasiparticles Overview explains lifetime, residue, and breakdown of a particle-like mode.
Polarons Preview treats a mobile carrier dressed by lattice motion.
BCS Theory explains how retarded electron–phonon attraction can produce Cooper pairing and where the reduced interaction gives way to Eliashberg theory.
Boltzmann Transport uses phonon absorption, emission, detailed balance, and momentum transfer in electronic collision integrals.
Heat Capacity and Thermodynamics owns calorimeter models, addenda subtraction, experimental Debye-coefficient extraction, and entropy accounting.
Math Needed for Quantum Matter collects the Fourier, eigenvalue, and numerical tools used here.
Condensed Matter Roadmap places phonons between band theory and material response.

References

M. Born and K. Huang, Dynamical Theory of Crystal Lattices, Oxford University Press, 1954.
N. W. Ashcroft and N. D. Mermin, Solid State Physics, Holt, Rinehart and Winston, 1976.
M. T. Dove, Introduction to Lattice Dynamics, Cambridge University Press, 1993, doi:10.1017/CBO9780511619885.
G. L. Squires, Introduction to the Theory of Thermal Neutron Scattering, 3rd ed., Cambridge University Press, 2012.
S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Giannozzi, “Phonons and Related Crystal Properties from Density-Functional Perturbation Theory,” Reviews of Modern Physics 73, 515–562 (2001), doi:10.1103/RevModPhys.73.515.
X. Gonze and C. Lee, “Dynamical Matrices, Born Effective Charges, Dielectric Permittivity Tensors, and Interatomic Force Constants from Density-Functional Perturbation Theory,” Physical Review B 55, 10355–10368 (1997), doi:10.1103/PhysRevB.55.10355.
A. Togo and I. Tanaka, “First Principles Phonon Calculations in Materials Science,” Scripta Materialia 108, 1–5 (2015), doi:10.1016/j.scriptamat.2015.07.021.
B. Fultz, “Vibrational Thermodynamics of Materials,” Progress in Materials Science 55, 247–352 (2010), doi:10.1016/j.pmatsci.2009.05.002.
A. Q. R. Baron, “Introduction to High-Resolution Inelastic X-Ray Scattering,” in Synchrotron Light Sources and Free-Electron Lasers, Springer, 2016, doi:10.1007/978-3-319-14394-1_52.
W. Cochran, “Crystal Stability and the Theory of Ferroelectricity,” Advances in Physics 9, 387–423 (1960), doi:10.1080/00018736000101229.
F. Giustino, “Electron–Phonon Interactions from First Principles,” Reviews of Modern Physics 89, 015003 (2017), doi:10.1103/RevModPhys.89.015003.
National Institute of Standards and Technology, “Neutron Spectroscopy”, overview of energy- and momentum-resolved neutron probes for phonons and other material excitations.