 Solution

The Debye model approximates a three-dimensional crystal as an elastic continuum supporting quantized normal modes, or phonons. At long wavelengths, the acoustic phonon branches have approximately linear dispersion relations ${\omega }_s(k)=v_{s}k$: one longitudinal branch and two transverse branches. Replace their complicated actual spectra by these linear forms, then introduce a cutoff chosen to preserve the total $3N$ vibrational degrees of freedom for $N$ atoms. This is an approximation to the spectrum, not a claim that a real crystal has no non-acoustic branches.

For periodic boundary conditions in volume $V$, each allowed wavevector occupies $(2\pi {)}^3/V$ in the space of wavevectors. Counting a spherical shell for each branch gives the density of states per unit angular frequencyAn effective speed of sound $v_D$ can be defined by $3/v_D^3=1/v_l^3+2/v_t^3$. In the single-speed Debye model, choose ${\omega }_D=v_D{k}_D$ so thatThus the Debye wavevector, Debye frequency and Debye temperature satisfyHere $g$ has units of seconds when angular frequency is measured in ${\mathrm{s}}^{-1}$; its integral is the dimensionless number of normal modes. The spherical cutoff replaces the real Brillouin zone by an equal-count continuum construction.

Each phonon mode is a quantum harmonic oscillator. Since phonon number is not conserved at thermal equilibrium, the chemical potential is zero and its Bose-Einstein distribution gives mean thermal energy $\hslash \omega /[e^{\hslash \omega /(k_{B}T)}-1]$. The zero-point energy is independent of temperature in the harmonic, fixed-volume approximation and does not contribute to the heat capacity at constant volume. Excluding that constant term,Differentiating the first integral at fixed volume avoids losing the temperature-dependent upper limit in the second one. The resulting Debye heat capacity isFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{r\ge 1}{e}^{-rx}$ and integrating term by term gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\zeta (4)={\pi }^4/15$. ThereforeThis $T^3$ law follows from the three-dimensional ${\omega }^2$ density of states and the availability only of low-frequency modes at low temperature. It is a central success of the Debye model.

For $T\gg {\Theta }_D$, the Bose-Einstein distribution reduces to the classical equipartition theorem: every mode contributes thermal energy $k_{B}T$. Consequentlythe Dulong-Petit law, or $3R$ per mole of atoms. The simple model omits detailed phonon dispersion relation, distinct optical phonon branches, and anharmonic effects such as thermal expansion. In a metal, the low-temperature total heat capacity additionally contains a term approximately proportional to $T$ from Electrons, so the lattice $T^3$ term should not be mistaken for the whole measurement.