Phonons are quantized lattice vibrations obeying Bose-Einstein statistics. The number of phonons in mode k with energy ℏω_k is n_k = 1/[exp(ℏω_k/kT)−1]. Total thermal energy U = Σ_k ℏω_k n_k and heat capacity C_V = ∂U/∂T follow from the phonon distribution. Acoustic branches have linear dispersion at small k; optical branches are nearly flat.
From your study of Bose-Einstein statistics, you know that bosons — particles with integer spin — can pile into the same quantum state without restriction. Phonons are the quantum mechanical description of collective lattice vibrations in a solid: they behave exactly like bosons, and their thermal statistics determine how a solid stores and conducts heat. The key insight is that instead of thinking of a crystal as 3N coupled oscillators (where N is the number of atoms), you think of it as a gas of phonon "particles," each with a well-defined wavevector k and frequency ω(k).
Each phonon mode behaves like a quantum harmonic oscillator: the energy in mode k is E_k = ℏω_k (n_k + 1/2), where n_k is the number of phonons in that mode. Since phonon number is not conserved (phonons are created and destroyed as the lattice vibrates), the chemical potential μ = 0, and the Bose-Einstein distribution simplifies to n_k = 1 / [exp(ℏω_k / k_BT) − 1]. This is the Planck distribution — the same formula that describes photons in a cavity. The average thermal energy is U = Σ_k ℏω_k n_k (plus zero-point energy), and the heat capacity is C_V = ∂U/∂T.
The crucial ingredient is the dispersion relation ω(k), which determines how phonon frequency varies with wavevector. Acoustic branches have ω → 0 as k → 0 (long-wavelength sound waves), with linear dispersion ω = v_s k at small k (v_s is the speed of sound). At low temperatures, only the long-wavelength acoustic phonons are thermally excited, and the resulting heat capacity scales as T³ — the Debye T³ law. Optical branches appear in crystals with multiple atoms per unit cell; they have finite frequency at k = 0 (neighboring atoms vibrating against each other, like in ionic crystals) and are nearly dispersionless. Because optical phonons have higher energies, they "freeze out" below their characteristic temperature and contribute little to C_V at low temperatures.
The Debye model approximates all phonon branches as acoustic with linear dispersion, cutting off at the Debye frequency ω_D (or equivalently the Debye temperature θ_D = ℏω_D / k_B). This gives the correct low-temperature T³ dependence and converges to the classical Dulong-Petit value 3Nk_B at high temperatures (T >> θ_D), where all phonon modes are thermally activated. Real solids deviate from the Debye model at intermediate temperatures because the true dispersion relation is more complex, but the phonon picture — bosons following Bose-Einstein statistics with a mode-dependent frequency — is the correct quantum mechanical foundation for understanding thermal properties of solids.