Einstein solid ($C_V/(3N{k}_B)=x^2{e}^x/(e^x-1{)}^2$)

The Einstein model of a solid represents the vibrations of $N$ chemical atoms by $3N$ independent one-dimensional quantum harmonic oscillators, all with the same positive angular frequency $\omega$. Equivalently, each chemical atom has three independent oscillator directions. The common angular frequency is a simplifying model assumption, not a claim about every vibration of a real solid. Hold the volume, $N$ and $\omega$ fixed. Define the characteristic temperature ${\theta }_E=\hslash \omega /k_B$ and the dimensionless ratio $x={\theta }_E/T$. The vibrational internal energy isThe first term is the zero-point energy $3N\hslash \omega /2$. It is constant at fixed $\omega$, so it contributes no heat capacity at constant volume. The temperature-dependent term givesThis is the vibrational heat capacity at constant volume of the model, with no expansion or other contributions included. ^[ref]

For a worked comparison, choose the model parameter ${\theta }_E=300\mathrm{K}$. The normalized heat capacity at constant volume depends only on $T/{\theta }_E$:

At $T={\theta }_E$, one mole of chemical atoms therefore has vibrational heat capacity at constant volume $3R(0.92067)\simeq 22.965\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, using the unrounded model value. Even at twice the characteristic temperature, the model remains about $2.06\%$ below the Dulong-Petit law limit. These are calculations for the ideal Einstein solid, not material measurements.

A short derivation uses the thermal partition function of a quantum harmonic oscillator. The energy eigenvalues $E_n=\hslash \omega (n+1/2)$ and their Boltzmann factors form a geometric series, giving the one-oscillator canonical partition function $z=e^{-x/2}/(1-e^{-x})$. The mean occupation is $1/(e^x-1)$, hence the displayed $U$ after summing over $3N$ independent quantum harmonic oscillators. Taking its derivative with $dx/dT=-x/T$ gives the displayed heat capacity at constant volume. ^[ref]

The two limits explain the shape. At high temperature, a Taylor series in small $x$ givesAt low temperature, $1-e^{-x}\to 1$, soThe zero-point energy remains even when the thermal response is tiny: freeze-out means that changing temperature scarcely changes the occupations, not that the chemical atoms have zero energy. This is the useful lesson of the Einstein model of a solid: the classical counting rule reappears when many levels are accessible, while discrete excitation spacings explain its failure upon cooling.

The exponential low-temperature decay is also the model's limitation. Real three-dimensional crystalline solids have acoustic phonon branches with arbitrarily small angular frequencies in the bulk limit. The Debye model includes these collective normal modes and predicts a low-temperature vibrational contribution $C_V\propto T^3$, which the single-frequency Einstein solid cannot reproduce. This qualification concerns the vibrational contribution; additional excitations can change a material's total heat capacity. ^[ref]

Einstein introduced the quantum treatment of solid vibrations in his 1907 paper Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Its historical importance is that quantum excitation energies explained why the classical Dulong-Petit law could fail; the modern quantum harmonic oscillator formulation above includes the zero-point energy explicitly. ^[ref]