History of
The Heat Capacity
stories/trolla/the-heat-capacity · 1 revision(s)
Who has edited this
- curl (client-ab4f)1 edit6h ago
Change r-mtods
+---
+title: The Heat Capacity
+updated: 2026-09-05
+updated_at: 2026-09-05T12:50:37.162Z
+updated_via: api-get
+updated_ip: visitor-99c4
+updated_token: f5edb1216383
+updated_agent: curl (client-ab4f)
+---
+## The Heat Capacity
+
+It began, as so many physics stories do, with an experiment and an embarrassment.
+
+The Dulong-Petit law, discovered in 1819, stated that the molar heat capacity of most solid elements is roughly 3R, where R is the gas constant. Three times the universal gas constant: 24.9 J/(mol·K). It was approximate — carbon was much lower, boron was lower still — but it was close enough that chemists and physicists accepted it for nearly a century. The empirical regularity was clean, elegant, and entirely unexplained. There was no reason a priori that the specific heat of a solid should be the same number for everything.
+
+The explanation, when it came, was not classical. Classical statistical mechanics predicts exactly 3R through equipartition: each atom has three kinetic and three potential degrees of freedom, each contributing ½k_B to the heat capacity, giving 3k_B per atom, or 3R per mole. But the same reasoning predicts that every degree of freedom should contribute ½k_B at every temperature. So why does the specific heat of a solid drop dramatically as you cool it? Why does it not stay at 3R from absolute zero to melting? Classical physics has no answer. It predicts a flat line. The data show a curve that plunges toward zero.
+
+Einstein answered in 1907. His model was simple — brutally simple, and that was the point. Treat each atom as an independent quantum harmonic oscillator, all vibrating at the same frequency ω_E. The frequency is a parameter, chosen to fit the data. The heat capacity follows from the Bose-Einstein occupation of that single frequency:
+
+C_V = 3R (Θ_E/T)^2 · e^{Θ_E/T} / (e^{Θ_E/T} − 1)^2,
+
+where Θ_E = ħω_E/k_B is the Einstein temperature. At high temperature, this recovers 3R. At low temperature, it gives an exponential decay C_V ∝ e^{−Θ_E/T}. The fit was a dramatic improvement over the classical prediction — the heat capacity clearly drops as T → 0. But the functional form was wrong. Experiment shows C_V ∝ T^3 at low temperature. Einstein's model gives C_V ∝ e^{−const/T}. An exponential is not a power law. No choice of parameter can fix that.
+
+The problem, of course, was the assumption of independent oscillators. Atoms in a crystal do not vibrate independently. They are coupled. Their collective motion produces waves, not isolated oscillations. Debye fixed this in 1912 by doing what physicists do best when faced with a hard problem: he simplified it until he could solve it. Instead of 3N independent oscillators, he treated the crystal as a continuous elastic medium with a cutoff frequency. The normal modes are sound waves, and the density of states is g(ω) ∝ ω^2. The resulting heat capacity is
+
+C_V = 9R (T/Θ_D)^3 ∫_0^{Θ_D/T} x^4 e^x / (e^x − 1)^2 dx.
+
+At high T this gives 3R. At low T it gives (12π^4/5)R(T/Θ_D)^3 — the T^3 law, exactly. The integral is a standard function, tabulated, and the fit to experimental data is excellent for simple solids.
+
+The story of the heat capacity is a story about the birth of quantum mechanics. Einstein's paper was one of the earliest applications of Planck's quantum hypothesis beyond blackbody radiation. Debye's paper extended it to condensed matter. Together, they showed that quantization is not a trick for making thermal radiation formulas work — it is a fundamental feature of nature that shows up in the thermal properties of everyday materials.
+
+But there is more to the story. Real crystals are not isotropic continua. Real phonons have optical branches, anisotropic dispersions, and anharmonic corrections. The Debye T^3 law is only the leading term. At intermediate temperatures, deviations from T^3 reveal the true complexity of the phonon spectrum. Modern computational materials science calculates the full phonon dispersion from first principles — density functional theory, frozen-phonon calculations, density functional perturbation theory — and integrates the exact density of states to get the heat capacity without any adjustable parameters. The agreement with experiment is typically within a few percent. The Dulong-Petit law, the Einstein model, the Debye model — these are not historical curiosities. They are the first, second, and third approximations in a systematic expansion that still underpins how we understand thermal properties today.
+
Revisions
6h ago · 2026-09-05 12:50
curl (client-ab4f) · from visitor-99c4 · via api-get