Case Study 2: Phonons — Quantized Sound

The Physical Situation

Sound is one of the most familiar phenomena in everyday life. We hear it, feel it, and manipulate it constantly. Yet the quantum mechanics of sound — the physics of phonons — is responsible for some of the most important and surprising properties of solid matter: why metals conduct heat, why crystals shatter along specific planes, why some materials become superconductors, and why your diamond engagement ring sparkles.

This case study traces the phonon concept from its origins in Einstein's 1907 theory of specific heat through the Debye model to modern applications in condensed matter physics and materials science.

The Specific Heat Puzzle

The Classical Prediction: Dulong and Petit (1819)

In 1819, Pierre Dulong and Alexis Petit made an empirical observation: at room temperature, the molar heat capacity of most solid elements is approximately $3R \approx 25$ J/(mol$\cdot$K), where $R$ is the gas constant. The explanation came from classical statistical mechanics: each atom in a solid vibrates in three dimensions with a harmonic potential. By the equipartition theorem, each degree of freedom contributes $\frac{1}{2}k_BT$ to the kinetic energy and $\frac{1}{2}k_BT$ to the potential energy. For $N$ atoms in 3D: $E = 3Nk_BT$, so $C_V = 3Nk_B = 3R$ per mole.

This works beautifully at room temperature. But it fails catastrophically at low temperatures: experimentally, $C_V \to 0$ as $T \to 0$, which the classical theory cannot explain (it predicts $C_V = 3R$ at all temperatures).

Einstein's Model (1907): Quantized Oscillators

Einstein proposed the first quantum theory of solids, nine years before quantum mechanics was formally developed. His key assumption: each atom vibrates as a quantum harmonic oscillator with the same frequency $\omega_E$ (the "Einstein frequency").

The average energy of a single quantum oscillator at temperature $T$ is:

$$\langle E \rangle = \hbar\omega_E\left(\frac{1}{e^{\hbar\omega_E/k_BT} - 1} + \frac{1}{2}\right)$$

The mean occupation number is the Bose-Einstein distribution:

$$\langle n \rangle = \frac{1}{e^{\hbar\omega_E/k_BT} - 1}$$

This is exactly the phonon occupation number from Section 34.6. At high temperature ($k_BT \gg \hbar\omega_E$), the occupation number becomes large ($\langle n \rangle \approx k_BT/\hbar\omega_E$), and the classical result $\langle E \rangle \approx k_BT$ is recovered (Dulong-Petit). At low temperature ($k_BT \ll \hbar\omega_E$), the occupation number vanishes exponentially ($\langle n \rangle \approx e^{-\hbar\omega_E/k_BT}$), and $C_V \to 0$.

Einstein's model correctly predicts that $C_V \to 0$ at low $T$, but it gives the wrong functional form: $C_V \propto e^{-\hbar\omega_E/k_BT}$ (exponential decay) instead of the observed $C_V \propto T^3$ (power law).

Debye's Model (1912): A Spectrum of Frequencies

Peter Debye realized that a real crystal does not vibrate at a single frequency — it has a spectrum of phonon modes with frequencies ranging from zero (long-wavelength sound waves) up to a maximum $\omega_D$ (the Debye frequency). By treating the solid as a continuum with a linear dispersion $\omega = v|k|$ (where $v$ is the sound velocity), Debye derived:

$$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$$

where $\Theta_D = \hbar\omega_D/k_B$ is the Debye temperature and $x = \hbar\omega/(k_BT)$.

In the low-temperature limit ($T \ll \Theta_D$), the upper limit of integration goes to infinity, and:

$$C_V = \frac{12\pi^4}{5}Nk_B\left(\frac{T}{\Theta_D}\right)^3$$

This is the famous Debye $T^3$ law, which matches experiment beautifully.

Debye Temperatures for Selected Materials

Diamond's extraordinarily high Debye temperature explains why its specific heat is anomalously low at room temperature — a puzzle that baffled 19th-century physicists. In the phonon picture, the strong carbon-carbon bonds (large spring constant $\kappa$) and low atomic mass combine to give very high phonon frequencies. At room temperature, most phonon modes are "frozen out" (their occupation number is exponentially small), so diamond behaves as if it were at "low temperature" even at 300 K.

Phonons in Experiment

Inelastic Neutron Scattering

The most direct way to observe phonons is through inelastic neutron scattering. A beam of thermal neutrons (with wavelengths $\sim 1$–$10$ \AA, comparable to lattice spacings) is directed at a crystal sample. When a neutron scatters, it can create or absorb a phonon:

Phonon creation: The neutron loses energy $\hbar\omega_k$ and momentum $\hbar\mathbf{k}$: $$E_{\text{out}} = E_{\text{in}} - \hbar\omega_k, \qquad \mathbf{p}_{\text{out}} = \mathbf{p}_{\text{in}} - \hbar\mathbf{k}$$

Phonon absorption: The neutron gains energy $\hbar\omega_k$ and momentum $\hbar\mathbf{k}$.

By measuring the energy and momentum of scattered neutrons, experimentalists map out the phonon dispersion relation $\omega(\mathbf{k})$ — the fundamental property of the crystal's vibrational spectrum.

In second-quantized language, the neutron-phonon interaction involves the coupling of the neutron field to the displacement operator $\hat{u}_n$, which contains both $\hat{a}_k$ (phonon absorption) and $\hat{a}_k^\dagger$ (phonon creation) terms. The scattering cross section is proportional to $|\langle f|\hat{a}_k^\dagger + \hat{a}_{-k}|i\rangle|^2$, where $|i\rangle$ and $|f\rangle$ are the initial and final phonon states.

Thermal Conductivity

In insulating crystals, heat is carried entirely by phonons. The thermal conductivity $\kappa$ can be estimated from the kinetic theory expression:

$$\kappa = \frac{1}{3}C_V v_s \ell$$

where $v_s$ is the phonon velocity (speed of sound) and $\ell$ is the phonon mean free path. At high temperatures, phonon-phonon scattering (arising from the anharmonic terms in the interatomic potential) limits $\ell$, and $\kappa$ decreases with increasing $T$. At low temperatures, phonon scattering from crystal boundaries or impurities dominates, and $\kappa \propto T^3$ (following $C_V$).

Diamond has the highest thermal conductivity of any natural material ($\kappa \approx 2{,}200$ W/(m$\cdot$K) at room temperature) — roughly 5 times that of copper. This is because diamond's strong bonds and low atomic mass give both high phonon velocities and long phonon mean free paths.

Phonons and Superconductivity

The BCS Mechanism

One of the most dramatic consequences of phonon physics is superconductivity — the complete disappearance of electrical resistance below a critical temperature. The BCS theory (Bardeen, Cooper, Schrieffer, 1957) explains superconductivity through a phonon-mediated attraction between electrons.

The mechanism, in second-quantized language:

1. An electron with wavevector $\mathbf{k}$ polarizes the lattice, creating a local concentration of positive charge (a virtual phonon $\hat{a}_\mathbf{q}^\dagger$).
2. A second electron with wavevector $\mathbf{k}'$ is attracted to this positive charge concentration, absorbing the virtual phonon ($\hat{a}_\mathbf{q}$).
3. The net effect is an attractive interaction between the two electrons, mediated by phonon exchange.

The interaction Hamiltonian for this process is:

$$\hat{H}_{\text{e-ph}} = \sum_{\mathbf{k}, \mathbf{q}, \sigma} g_\mathbf{q}\, \hat{c}_{\mathbf{k}+\mathbf{q},\sigma}^\dagger \hat{c}_{\mathbf{k},\sigma}(\hat{a}_\mathbf{q} + \hat{a}_{-\mathbf{q}}^\dagger)$$

where $g_\mathbf{q}$ is the electron-phonon coupling constant. The term $\hat{c}_{\mathbf{k}+\mathbf{q},\sigma}^\dagger\hat{c}_{\mathbf{k},\sigma}\hat{a}_\mathbf{q}$ describes an electron absorbing a phonon (scattering from $\mathbf{k}$ to $\mathbf{k}+\mathbf{q}$), and $\hat{c}_{\mathbf{k}+\mathbf{q},\sigma}^\dagger\hat{c}_{\mathbf{k},\sigma}\hat{a}_{-\mathbf{q}}^\dagger$ describes an electron emitting a phonon.

When two electrons with opposite momenta and spins ($\mathbf{k}\uparrow$ and $-\mathbf{k}\downarrow$) exchange virtual phonons, they form a Cooper pair. The Cooper pair is a boson (two fermions bound together have integer total spin), and below the critical temperature, Cooper pairs undergo Bose-Einstein condensation into a macroscopic quantum state — the superconducting ground state.

This remarkable chain — phonons (bosonic quasiparticles) mediate an attraction between electrons (fermions) that bind into Cooper pairs (composite bosons) that condense — beautifully illustrates the power of second quantization. Every step involves creation and annihilation operators acting on Fock space.

Phonons as a Template for Quasiparticles

Phonons are the prototype for a vast family of quasiparticles — collective excitations of many-body systems that behave like particles:

Every one of these quasiparticles is described by creation and annihilation operators acting on Fock space. The formalism of second quantization, developed for phonons and photons, is the universal language of collective excitations in condensed matter physics.

Quantitative Analysis: Phonon Mean Free Path in Silicon

Let us work through a specific calculation to see second quantization in action.

Problem: Estimate the thermal conductivity of silicon at $T = 300$ K.

Given: $\Theta_D = 645$ K, $v_s \approx 8{,}400$ m/s (average sound velocity), lattice constant $a = 5.43$ \AA.

Step 1: Specific heat at 300 K. Since $T/\Theta_D = 300/645 = 0.47$, we are in the intermediate regime. Numerical evaluation of the Debye integral gives $C_V/V \approx 1.6 \times 10^6$ J/(m$^3\cdot$K).

Step 2: Mean free path. At room temperature, the dominant scattering mechanism is three-phonon (Umklapp) scattering. For silicon, $\ell \approx 40$ nm at 300 K (determined from neutron scattering data).

Step 3: Thermal conductivity: $$\kappa = \frac{1}{3}C_V v_s \ell = \frac{1}{3}(1.6 \times 10^6)(8400)(40 \times 10^{-9}) \approx 180 \text{ W/(m$\cdot$K)}$$

The experimental value is $\kappa = 148$ W/(m$\cdot$K). Our estimate, which used only the simplest kinetic theory formula and Debye model phonon properties, is within 20% — a testament to the quantitative power of the phonon picture.

Questions for Reflection

1. Why does the Debye model work so well despite its crude approximation (linear dispersion with a sharp cutoff)? In particular, why does it give the correct $T^3$ law at low temperatures regardless of the detailed dispersion?

2. The electron-phonon coupling that produces superconductivity also causes electrical resistance at high temperatures (phonon scattering of electrons). How can the same interaction produce both resistance (destructive) and superconductivity (constructive)?

3. Phonons carry crystal momentum $\hbar\mathbf{k}$ but no real momentum (the center of mass of the crystal does not move). What is the physical meaning of phonon momentum, and why is it conserved (up to reciprocal lattice vectors)?

4. At very low temperatures, the phonon mean free path can exceed the size of the sample ($\ell > L$). In this regime, phonon transport is ballistic rather than diffusive — phonons travel across the sample without scattering. This has been observed experimentally and has practical implications for heat management in nanoscale devices. How would you modify the kinetic theory expression for $\kappa$ in this regime?