Heat Capacity of Solids: Einstein and Debye Models

Introduction

Classical theory predicts a lattice heat capacity near 3R per mole of atoms, but many solids show a strong decrease on cooling. Einstein's model introduced quantised oscillators and explained why heat capacity tends toward zero. Debye's model added a spectrum of low-frequency collective vibrations and captured the observed low-temperature T³ law in many insulating crystalline solids.

Core explanation

In the Einstein model, N atoms supply 3N independent harmonic modes all at one frequency ν E. With x E = θ E/T and θ E = hν E/k, the lattice heat capacity is C V,E = 3Nk x E²e^x E/(e^x E − 1)². For one mole of atoms, replace Nk by R. At high T, x E is small and C V,E approaches 3Nk, the Dulong–Petit limit. At low T, it falls approximately with an exponential factor e^(−θ E/T). The model captures quantum freeze-out but assigns every mode the same frequency.

Real crystals have collective normal modes called phonons. Long-wavelength acoustic modes have frequencies approaching zero, so some modes remain thermally accessible even when T is far below a characteristic high-frequency scale. The Debye model approximates the acoustic mode density as proportional to frequency squared up to a cutoff ν D chosen to give a total of 3N modes. Its heat capacity can be written C V,D = 9Nk(T/θ D)³ ∫₀^(θ D/T) [x⁴e^x/(e^x − 1)²]dx. At high T it also tends to 3Nk; at low T the integral approaches a constant and C V,D ≈ (12π⁴/5)Nk(T/θ D)³.

The two models are alternatives for lattice vibrations, not two contributions to add for the same set of modes. A real solid can have acoustic and optical branches and may be fitted with a combination representing distinct mode subsets, but double-counting all 3N modes would be wrong. Metals also have an electronic heat-capacity contribution, often approximately proportional to T at sufficiently low temperature, in addition to a lattice term. Magnetic ordering or defects can add further features.

Heat capacity at constant pressure differs from C V because thermal expansion does work. The difference is often modest for many solids at ordinary conditions but need not be zero. The idealised Einstein and Debye formulas here describe lattice C V under harmonic assumptions.

Step-by-step reasoning

Identify whether the model assumes one oscillator frequency or a continuum of acoustic modes. Convert frequency or cutoff to θ. At high T, check the 3Nk limit. At low T, use exponential suppression for the Einstein model or T³ for the Debye acoustic model. If comparing with data, separate lattice, electronic and other contributions and check whether C p or C V was measured.

Visual explanation

Draw frequency on a horizontal axis. Einstein's model is one narrow spike at ν E; Debye's is a smooth density rising roughly as ν² to a cutoff. On a second graph, plot C V against T: both approach 3R at high T, but the Einstein curve drops exponentially at low T while the Debye curve follows a T³ guide line.

Real-world analogy

An orchestra with every instrument constrained to one pitch would become nearly silent below the energy needed to play that pitch. An orchestra with many very low notes can still produce sound at low energy. A crystal's low-frequency acoustic vibrations give it thermal modes that a single-frequency Einstein picture lacks. The analogy describes mode availability, not literal audible music.

Real-world example

Low-temperature calorimetry of an insulating crystal can reveal an approximately cubic lattice heat-capacity trend, supporting a Debye-style acoustic description. At warmer temperatures, heat capacity approaches the classical 3R scale per mole of atoms. Deviations may reveal optical phonons, electronic carriers or phase changes.

Why?

Quantum oscillator occupation is small when kT is below its energy gap. Einstein's model has one nonzero gap, giving exponential suppression. Debye's acoustic spectrum extends toward arbitrarily low frequency in a large crystal; counting the thermally accessible low-frequency modes produces a number proportional to T³ in three dimensions.

Common misconception

Einstein and Debye curves should not simply be added as full 3N-mode models. Another error is to say Debye's low-T heat capacity is zero because all vibrations are frozen; low-frequency acoustic modes remain available. The Dulong–Petit value is a high-temperature limit, not a universal constant for every solid at every temperature.

Worked example

Suppose an insulating solid has Debye temperature θ D = 300 K and is measured at T = 15 K, so T/θ D = 0.050. The low-T Debye estimate per mole of atoms is C V ≈ (12π⁴/5)R(0.050)³. Since 12π⁴/5 ≈ 234, C V ≈ 234(8.314)(1.25 × 10⁻⁴) ≈ 0.243 J mol⁻¹ K⁻¹. This is far below 3R ≈ 24.9 J mol⁻¹ K⁻¹, as expected at T ≪ θ D.

Quick check

1. Which model naturally predicts a low-temperature C V proportional to T³ in a three-dimensional insulating crystal? Answer: The Debye model, because it includes a continuum of low-frequency acoustic modes with an appropriate frequency-dependent density of states.

Exam focus

Identify whether the sample's heat capacity is per mole of atoms or per mole of formula units. Use the correct low-T form for the model and do not double-count modes. Distinguish C V from C p and note possible electronic or magnetic additions in real materials. Check the high-T 3R per mole of atoms limit.

Advanced insight

The Debye T³ coefficient depends on acoustic sound speeds through θ D, linking calorimetry to elastic properties. In low-dimensional materials, the density of states and phonon dispersions can change the low-T power law. The simple three-dimensional isotropic Debye result should therefore not be transplanted without checking dimensionality and mode structure.

Summary

Einstein assigns one frequency to 3N crystal modes, explaining quantum freeze-out but predicting exponential low-T decay. Debye uses an acoustic frequency spectrum and gives C V ∝ T³ at low T. Both approach 3Nk at high T. Real solids may also contain electronic, optical, magnetic and expansion-related contributions.

Practice questions

1. What is the approximate high-temperature lattice C V per mole of atoms in either model? Answer: 3R ≈ 24.9 J mol⁻¹ K⁻¹, the Dulong–Petit limit under harmonic classical behaviour. 2. If T doubles while still deep in the Debye T³ regime, how does lattice C V change? Answer: It increases by 2³ = 8, assuming θ D and the low-T model remain applicable. 3. Why does an Einstein model miss the Debye T³ law? Answer: It assigns all modes one finite frequency and therefore lacks the continuum of arbitrarily low-frequency acoustic modes that dominate a crystal's low-temperature heat capacity.