Les vibrations thermiques du réseau cristallin jouent un rôle essentiel dans la compréhension de nombreux phénomènes physiques, notamment la capacité calorifique à volume constant du cristal et les propriétés optiques. Ces vibrations sont modélisées à l'aide de phonons, des quanta d'oscillations du réseau. Le comportement de la chaleur spécifique, qui dépend de ces vibrations, peut être analysé en fonction des interactions des phonons et des propriétés du cristal.

L'analyse de l'énergie vibratoire du cristal à partir de la densité d'états des phonons D(ω)D(\omega) permet d'obtenir une expression pour la chaleur spécifique CVC_V. Lorsque l'on insère la densité d'états dans l'expression de l'énergie, la chaleur spécifique du réseau cristallin à température constante est donnée par l'intégrale :

CV=09NkBz4ez(ez1)2dz,C_V = \int_0^\infty \frac{9Nk_B z^4 e^z}{(e^z - 1)^2} \, dz,

z=ωkBTz = \frac{\hbar \omega}{k_B T} et θ=ωDkB\theta = \frac{\hbar \omega_D}{k_B} est la température de Debye. À des températures faibles ( TθT \ll \theta ), on observe que la capacité calorifique suit une relation CVT3C_V \sim T^3, ce qui correspond aux résultats expérimentaux. À haute température ( TθT \gg \theta ), la capacité calorifique devient indépendante de la température et atteint une valeur constante, conformément à la loi classique de Dulong et Petit.

Le comportement des phonons dans le cristal peut être décrit à l'aide d'un modèle simple, où le réseau cristallin est représenté par des masses ponctuelles reliées par des ressorts élastiques. Dans ce modèle, la relation de dispersion des phonons est donnée par :

ω=4fmsin(Ka2),\omega = \sqrt{\frac{4f}{m}} \sin\left( \frac{K a}{2} \right),