Heat Capacity of Classical Gases

The heat capacity measures how much energy a substance absorbs per unit temperature rise: CV=(∂U/∂T)VC_V=(\partial U/\partial T)_V at fixed volume, CP=(∂H/∂T)PC_P=(\partial H/\partial T)_P at fixed pressure. For a classical ideal gas the internal energy depends only on temperature, so CVC_V follows directly from how energy is partitioned among a molecule's microscopic degrees of freedom. The tool is the equipartition theorem: every coordinate or momentum entering the Hamiltonian quadratically carries a mean energy 12kBT\tfrac12 k_BT, independent of its coupling constant. Counting these quadratic terms gives a strikingly simple result. Three translational momenta yield 32kBT\tfrac32 k_BT per molecule for any gas. A linear molecule adds two rotational degrees of freedom (rotation about the symmetry axis is inert); a nonlinear molecule adds three. Each vibrational mode contributes two quadratic terms — one kinetic, one potential — and so carries a full kBTk_BT. Summing and multiplying by Avogadro's number gives molar heat capacities that are simple rational multiples of the gas constant R=NAkBR=N_A k_B: 32R\tfrac32 R (monatomic), 52R\tfrac52 R (diatomic with rotation), 72R\tfrac72 R (with vibration active). Two universal relations round out the classical picture. First, Mayer's relation CP−CV=RC_P-C_V=R (per mole) for any ideal gas, a direct consequence of H=U+PV=U+NkBTH=U+PV=U+Nk_BT. Second, the adiabatic index γ=CP/CV\gamma=C_P/C_V, which equals 53\tfrac53, 75\tfrac75, and 97\tfrac97 for the three cases and sets the speed of sound. The triumph and the failure of this picture are equally instructive: it correctly predicts 32R\tfrac32 R for noble gases and 52R\tfrac52 R for N2\mathrm N_2 and O2\mathrm O_2 near room temperature, but it also predicts these values should never change with temperature — which experiment flatly contradicts. That failure, the freezing-out of degrees of freedom, is the doorway to quantum statistical mechanics and is treated in detail in the companion article on the failure of equipartition.

Key takeaways
  • Classical equipartition assigns 12kBT\tfrac12 k_BT per quadratic degree of freedom, making molar heat capacity a sum of temperature-independent integer/half-integer multiples of RR.
  • Monatomic gas: f=3f=3 translational ⇒CV=32R\Rightarrow C_V=\tfrac32 R; diatomic with rotation: f=5⇒CV=52Rf=5\Rightarrow C_V=\tfrac52 R; with vibration active: f=7⇒CV=72Rf=7\Rightarrow C_V=\tfrac72 R.
  • A vibrational mode counts as two quadratic terms (kinetic + potential), contributing kBTk_BT, not 12kBT\tfrac12 k_BT.
  • Mayer's relation: CP−CV=RC_P-C_V=R per mole for any ideal gas, because H=U+NkBTH=U+Nk_BT.
  • Adiabatic index γ=CP/CV=53,75,97\gamma=C_P/C_V=\tfrac53,\tfrac75,\tfrac97 for the three cases; it sets the speed of sound cs=γRT/Mc_s=\sqrt{\gamma RT/M}.
  • Key idea: count the active quadratic degrees of freedom ff; then CV=f2RC_V=\tfrac f2 R, CP=CV+RC_P=C_V+R, and γ=(f+2)/f\gamma=(f+2)/f — but "active" is temperature-dependent, and that is where classical physics breaks.
The core insight

Heat capacity is bookkeeping on degrees of freedom. Count the quadratic terms ff in the molecular Hamiltonian that are thermally active; each contributes 12R\tfrac12 R per mole. Then CV=f2RC_V=\tfrac f2 R, CP=CV+RC_P=C_V+R (Mayer), and γ=CP/CV=(f+2)/f\gamma=C_P/C_V=(f+2)/f. Monatomic f=3f=3, diatomic f=5f=5 (with rotation) or 77 (with vibration).

Where this lives today

The mode-counting rules are the backbone of chemical thermodynamics: statistical-mechanics software computes molar heat capacities, entropies, and equilibrium constants by summing translational, rotational, and vibrational partition functions for each molecule, with the classical f2R\tfrac f2 R values as the high-temperature limit and quantum corrections layered on top. In atmospheric science the temperature dependence of air's γ\gamma (as vibrations begin to thaw at high altitude and in hot exhaust plumes) affects sound propagation and shock structure. The speed-of-sound-to-degree-of-freedom link is used in acoustic gas analysis and in calibrating the thermodynamics of exotic gases. And the Dulong-Petit value 3R3R remains the standard high-temperature benchmark against which the Debye model's low-temperature T3T^3 behavior — the subject of the companion failure-of-equipartition article — is measured.

The physical picture

Think of a gas molecule as a small machine with several ways to store energy: it can translate, rotate, and (as a spring connecting its atoms) vibrate. Equipartition says the thermal bath fills each quadratic storage slot to exactly 12kBT\tfrac12 k_BT on average, like water finding the same level in connected vessels. Heat capacity is then just the count of open slots times 12R\tfrac12 R per mole — the more ways a molecule can hold energy, the more heat it takes to warm it by a degree. But the slots are not all open at once. Quantum mechanics quantizes each mode, and a mode stays frozen — stuck in its ground state, storing nothing — until kBTk_BT climbs above its level spacing. Translation has an effectively continuous spectrum and is always active; rotation activates at low temperatures (a few kelvin for most molecules, but ∼85 K\sim85\,\mathrm K for light H2\mathrm H_2); vibration needs thousands of kelvin. So the heat capacity of a real gas rises in steps as temperature climbs and new slots open, tracing a staircase from 32R\tfrac32 R to 52R\tfrac52 R to 72R\tfrac72 R.

Figure
A diatomic molecule stores energy in translation (3 quadratic terms), rotation (2, about the two perpendicular axes), and vibration (2: kinetic + potential of the bond spring). Rotation about the bond axis is inert. Each active quadratic slot holds ½kT.
diatomic molecule translation: 3 rotation: 2 (⊥ axes) vibration: 2 (KE + PE) translation: 3 × ½kT = 3/2 kT (always) + rotation: 2 × ½kT = kT → 5/2 kT total + vibration: 2 × ½kT = kT → 7/2 kT total C_V = (f/2) R, C_P = C_V + R, γ = (f+2)/f

The payoff is the staircase itself, plotted for hydrogen. At the lowest temperatures only translation is active and CV=32RC_V=\tfrac32 R. As TT rises past Θrot\Theta_{\rm rot}, rotation switches on and CVC_V climbs to 52R\tfrac52 R. Far higher, past Θvib\Theta_{\rm vib}, vibration activates and CVC_V approaches 72R\tfrac72 R. Each flat tread is a classical equipartition plateau; each riser is a quantum mode thawing. The classical theory correctly predicts the heights of the treads but is silent about the risers — which is precisely its boundary.

Figure
Punchline: the molar heat capacity C_V/R of a diatomic gas climbs in steps — 3/2 (translation) → 5/2 (rotation) → 7/2 (vibration). Classical equipartition predicts the plateau heights exactly; the quantum freeze-out sets the risers.
temperature (log scale) → C_V / R 3/2 5/2 7/2 translation + rotation + vibration frozen out Θ_rot Θ_vib

The core equations

For a classical ideal gas the internal energy is the equipartition sum over active quadratic degrees of freedom ff, and the heat capacity is its temperature derivative. Because UU is linear in TT, CVC_V is a constant.

U=f2NkBT,CV=(∂U∂T)V=f2NkB=f2R (per mole)(1)U=\frac f2 N k_BT,\qquad C_V=\left(\frac{\partial U}{\partial T}\right)_V=\frac f2 N k_B=\frac f2 R\ \text{(per mole)} \qquad (1)

where: ff — number of active quadratic degrees of freedom (3 translational always; +2 or +3 rotational; +2 per active vibrational mode); R=NAkB=8.314 J mol−1K−1R=N_A k_B=8.314\,\mathrm{J\,mol^{-1}K^{-1}}. Reading: heat capacity counts storage slots — each active quadratic slot is worth 12R\tfrac12 R per mole. The constant-pressure heat capacity exceeds CVC_V by exactly RR, because at constant pressure some of the added heat does expansion work. This is Mayer's relation.

CP−CV=R,γ=CPCV=f+2f(2)C_P-C_V=R,\qquad \gamma=\frac{C_P}{C_V}=\frac{f+2}{f} \qquad (2)

where: CPC_P — molar heat capacity at constant pressure; γ\gamma — adiabatic index. Reading: γ\gamma decreases as more degrees of freedom open (53→75→97\tfrac53\to\tfrac75\to\tfrac97), because the same RR of expansion work is spread over a larger CVC_V. The adiabatic index controls the speed of sound, tying an abstract mode count to a directly measurable macroscopic quantity.

cs=γRTM(3)c_s=\sqrt{\frac{\gamma R T}{M}} \qquad (3)

where: csc_s — adiabatic sound speed; MM — molar mass. Reading: measuring the speed of sound in a gas is, in effect, counting its molecular degrees of freedom.

Derivation AMonatomic and diatomic heat capacities by counting quadratic terms

The goal is to derive the molar heat capacities 32R\tfrac32 R, 52R\tfrac52 R, 72R\tfrac72 R from the equipartition count of quadratic terms in the molecular Hamiltonian, and to fix the subtle point that a vibrational mode counts as two quadratic terms while the axial rotation of a linear molecule counts as none. This is the calculation that turns molecular structure into a measurable heat capacity and adiabatic index.

Step 1. Monatomic gas. The Hamiltonian per atom is purely translational, H=(px2+py2+pz2)/2mH=(p_x^2+p_y^2+p_z^2)/2m: three quadratic momentum terms, f=3f=3. By equipartition ⟨H⟩=32kBT\langle H\rangle=\tfrac32 k_BT, so U=32NkBTU=\tfrac32 Nk_BT and CV=32RC_V=\tfrac32 R per mole.

Step 2. Diatomic gas, rotation active. Add rotation about the two axes perpendicular to the bond: Hrot=L12/2I1+L22/2I2H_{\rm rot}=L_1^2/2I_1+L_2^2/2I_2, two quadratic terms. Why only two, not three: rotation about the bond axis has a negligible moment of inertia (the mass is concentrated on the axis), so its rotational quantum is enormous and the mode is always frozen — it contributes nothing. Total f=3+2=5f=3+2=5, so CV=52RC_V=\tfrac52 R.

Step 3. Diatomic gas, vibration active. Add the bond vibration, modeled as a harmonic oscillator Hvib=p2/2μ+12μω2x2H_{\rm vib}=p^2/2\mu+\tfrac12\mu\omega^2 x^2: two quadratic terms (kinetic and potential). Why a vibration is worth kBTk_BT, not 12kBT\tfrac12 k_BT: both the kinetic and potential energies are quadratic, and each earns 12kBT\tfrac12 k_BT, so a fully active vibrational mode contributes kBTk_BT. Total f=3+2+2=7f=3+2+2=7, so CV=72RC_V=\tfrac72 R.

Step 4. Adiabatic indices: γ=(f+2)/f\gamma=(f+2)/f gives 53\tfrac53 (f=3f=3), 75=1.4\tfrac75=1.4 (f=5f=5), 97≈1.29\tfrac97\approx1.29 (f=7f=7). Why γ\gamma falls: γ=1+R/CV=1+2/f\gamma=1+R/C_V=1+2/f, so more degrees of freedom mean a smaller γ\gamma. Verified: the mode counts f=3,5,7f=3,5,7, the heat capacities 32R,52R,72R\tfrac32 R,\tfrac52 R,\tfrac72 R, and the indices γ=53,75,97\gamma=\tfrac53,\tfrac75,\tfrac97 are confirmed symbolically, along with the numeric 52R=20.79 J mol−1K−1\tfrac52 R=20.79\,\mathrm{J\,mol^{-1}K^{-1}}, in verify_aq-sm-heat-capacity-classical-gases.py.

Go deeper · Derivation B · 1 minMayer's relation C_P − C_V = R from the enthalpy
Derivation BMayer's relation C_P − C_V = R from the enthalpy

The goal is to prove that for any ideal gas the constant-pressure heat capacity exceeds the constant-volume one by exactly the gas constant RR per mole, regardless of the number of internal degrees of freedom, and to interpret the difference as the expansion work done when heating at constant pressure. This universal relation lets one measure γ\gamma and infer ff from either heat capacity.

Step 1. Definitions: CV=(∂U/∂T)VC_V=(\partial U/\partial T)_V and CP=(∂H/∂T)PC_P=(\partial H/\partial T)_P, where H=U+PVH=U+PV is the enthalpy. For an ideal gas U=U(T)U=U(T) depends only on temperature (no interactions), and PV=NkBTPV=Nk_BT.

Step 2. Write the enthalpy of one mole: H=U+PV=U(T)+RTH=U+PV=U(T)+RT. Why PV=RTPV=RT contributes: at constant pressure, warming the gas by dT\mathrm dT expands it, and the term RTRT is the work-plus-flow-energy bookkeeping that enthalpy carries.

Step 3. Differentiate at constant pressure: CP=dHdT=dUdT+R=CV+RC_P=\dfrac{\mathrm dH}{\mathrm dT}=\dfrac{\mathrm dU}{\mathrm dT}+R=C_V+R. Since UU depends only on TT, dU/dT=CV\mathrm dU/\mathrm dT=C_V regardless of the process. Hence CP−CV=RC_P-C_V=R, equation (2), independent of ff.

Step 4. Physical reading: at constant volume all added heat raises UU; at constant pressure the gas also does expansion work P dV=R dTP\,\mathrm dV=R\,\mathrm dT per mole, so more heat is needed for the same dT\mathrm dT — exactly RR more. Verified: CP=CV+RC_P=C_V+R from H=U+RTH=U+RT (symbolic), and the resulting γ\gamma values, are confirmed in the verify script.

Limiting cases

CaseCVC_Vγ\gammaExample
Monatomic (f=3f=3)32R=12.47\tfrac32 R=12.475/3≈1.6675/3\approx1.667He, Ar, Ne
Diatomic, rotation (f=5f=5)52R=20.79\tfrac52 R=20.797/5=1.4007/5=1.400N2,O2\mathrm N_2,\mathrm O_2 at 300 K
Diatomic, +vibration (f=7f=7)72R=29.10\tfrac72 R=29.109/7≈1.2869/7\approx1.286N2\mathrm N_2 above ∼3000 K\sim3000\,\mathrm K
Nonlinear polyatomic (f=6f=6, rot only)3R=24.943R=24.944/3≈1.3334/3\approx1.333H2O,CO2\mathrm{H_2O,CO_2} (bent) rot part
Solid (Dulong-Petit, f=6f=6)3R=24.943R=24.94—Cu, Ag, Pb above ΘD\Theta_D

Worked examples

Example 1 — Predicting the heat capacity and adiabatic index of nitrogen at room temperature (count, apply, check). N2\mathrm N_2 at T=300 KT=300\,\mathrm K: at this temperature rotation is fully active but vibration (Θvib≈3340 K\Theta_{\rm vib}\approx3340\,\mathrm K) is frozen.

Step 1. Count active quadratic terms: 3 translational + 2 rotational (linear molecule) + 0 vibrational = f=5f=5.

Step 2. Apply equation (1): CV=52R=52(8.314)=20.79 J mol−1K−1C_V=\tfrac52 R=\tfrac52(8.314)=\mathbf{20.79\,J\,mol^{-1}K^{-1}}; CP=CV+R=72R=29.10 J mol−1K−1C_P=C_V+R=\tfrac72 R=29.10\,\mathrm{J\,mol^{-1}K^{-1}}; γ=CP/CV=7/5=1.400\gamma=C_P/C_V=7/5=\mathbf{1.400}.

Step 3. Boundary check: the measured CVC_V of N2\mathrm N_2 at 300 K300\,\mathrm K is 20.8 J mol−1K−120.8\,\mathrm{J\,mol^{-1}K^{-1}} and γ=1.40\gamma=1.40 — an essentially exact match, confirming that vibration is frozen out and only 5 degrees of freedom are active. If vibration were active, CVC_V would be 72R=29.1\tfrac72 R=29.1 and γ\gamma would drop to 1.291.29.

Example 2 — The speed of sound in air from the adiabatic index. Air is mostly diatomic (N2,O2\mathrm N_2,\mathrm O_2), so γ=1.40\gamma=1.40 and mean molar mass M=29 g/mol=0.029 kg/molM=29\,\mathrm{g/mol}=0.029\,\mathrm{kg/mol}.

Step 1. From equation (3), cs=γRT/Mc_s=\sqrt{\gamma RT/M}. At T=293 KT=293\,\mathrm K (20 °C), γRT=1.40(8.314)(293)=3410 J/mol\gamma RT=1.40(8.314)(293)=3410\,\mathrm{J/mol}.

Step 2. Divide by MM: cs=3410/0.029=1.176×105=343 m/sc_s=\sqrt{3410/0.029}=\sqrt{1.176\times10^5}=\mathbf{343\,m/s}.

Step 3. Boundary check: the measured speed of sound in dry air at 20 ∘C20\,^\circ\mathrm C is 343 m/s343\,\mathrm{m/s} — exact. Had we (wrongly) used the isothermal relation γ=1\gamma=1, we would get 290 m/s290\,\mathrm{m/s}, the value Newton originally computed and which puzzled physicists for a century until Laplace supplied the adiabatic γ\gamma. The audible speed of sound is a direct measurement of the molecular degree-of-freedom count.

Origins · 1 mindulong-petit law and how this was discovered
DL
dulong-petit law
statistical-mechanics concept · Dulong & Petit 1819; Mayer 1842

The Dulong-Petit law (Pierre Louis Dulong and Alexis Thérèse Petit, 1819) states that the molar heat capacity of a crystalline solid is 3R≈24.9 J mol−1K−13R\approx24.9\,\mathrm{J\,mol^{-1}K^{-1}}, independent of the material — the equipartition result for 3N3N oscillators, each with two quadratic terms. It was one of the first quantitative successes of the kinetic-molecular picture and let chemists estimate atomic weights. Its low-temperature failure (heat capacities falling toward zero as T3T^3) was, like the diatomic-gas freeze-out, a decisive clue that classical equipartition was incomplete, resolved only by the Einstein and Debye quantum models. For gases the analogous universal relation is Mayer's (CP−CV=RC_P-C_V=R, Julius Robert von Mayer 1842), which follows from the ideal-gas enthalpy and holds for any number of internal degrees of freedom.

History
1819
— Dulong and Petit find the molar heat capacity of solids is nearly 3R3R, independent of substance.
1842
— Mayer states CP−CV=RC_P-C_V=R and estimates the mechanical equivalent of heat.
1859
— Maxwell derives equipartition of translational and rotational energy.
1876
— Boltzmann formalizes equipartition for arbitrary quadratic degrees of freedom.
1877
— the anomalously low room-temperature heat capacity of diatomic gases is recognized as an open puzzle.

Going further

References. Pathria & Beale, Statistical Mechanics, 3rd ed., §3.7, §6.5 · Kardar, Statistical Physics of Particles, Ch. 4 · Reif, Fundamentals of Statistical and Thermal Physics, §7.5–7.7 · Sethna, Statistical Mechanics: Entropy, Order Parameters, and Complexity, Ch. 6 · Landau & Lifshitz, Statistical Physics Part 1, §44–§47.

Provenance. All numbers are standard textbook values: R=8.314 462 J mol−1K−1R=8.314\,462\,\mathrm{J\,mol^{-1}K^{-1}}, kB=1.380649×10−23 J K−1k_B=1.380649\times10^{-23}\,\mathrm{J\,K^{-1}}, NA=6.022 140 76×1023 mol−1N_A=6.022\,140\,76\times10^{23}\,\mathrm{mol^{-1}} (CODATA/SI). Heat-capacity values 32R,52R,72R\tfrac32 R,\tfrac52 R,\tfrac72 R and indices γ=5/3,7/5,9/7\gamma=5/3,7/5,9/7 are derivable by mode counting (Reif §7.5, Pathria §3.7). Characteristic temperatures (Θrot≈85 K\Theta_{\rm rot}\approx85\,\mathrm K for H2\mathrm H_2, ≈2.9 K\approx2.9\,\mathrm K for N2\mathrm N_2; Θvib≈3340 K\Theta_{\rm vib}\approx3340\,\mathrm K for N2\mathrm N_2) are standard molecular-spectroscopy values (Pathria §6.5, Kardar Ch. 4). Sound speed in air ≈343 m s−1\approx343\,\mathrm{m\,s^{-1}} at 20 ∘C20\,^\circ\mathrm C follows from cs=γRT/Mc_s=\sqrt{\gamma RT/M} with γ=1.40\gamma=1.40, M=29 g mol−1M=29\,\mathrm{g\,mol^{-1}}. A verifier should recompute the Gaussian integral behind equipartition, confirm Mayer's relation from H=U+NkBTH=U+Nk_BT, and check the numeric heat capacities. All checks are in verify_aq-sm-heat-capacity-classical-gases.py.

Builds on: Equipartition Theorem · Classical Ideal Gas · Generalized Equipartition Theorem · Maxwell-Boltzmann Velocity Distribution Leads to: Sackur-Tetrode Equation · Virial Equation of State · Thermal de Broglie Wavelength

Why is a diatomic gas 5/2 R and not 7/2 R at room temperature?

Because vibration is frozen out. A diatomic molecule has 3 translational, 2 rotational, and 2 vibrational quadratic degrees of freedom, which would give f=7f=7 and CV=72RC_V=\tfrac72 R if all were active. But the vibrational quantum ℏω\hbar\omega corresponds to a characteristic temperature of thousands of kelvin (e.g. Θvib≈3340 K\Theta_{\rm vib}\approx3340\,\mathrm K for N2\mathrm N_2), far above room temperature, so at 300 K almost no molecules are vibrationally excited and the mode contributes nothing. Only translation and rotation are active, giving f=5f=5 and CV=52RC_V=\tfrac52 R. Heat a diatomic gas to a few thousand kelvin and CVC_V does climb toward 72R\tfrac72 R.

Why does a linear molecule have only two rotational degrees of freedom?

Because rotation about the molecular (bond) axis has a vanishingly small moment of inertia — nearly all the mass sits on that axis, so spinning about it barely stores any energy. Quantum-mechanically, a tiny moment of inertia means a huge rotational level spacing (ℏ2/2I\hbar^2/2I is enormous), so the mode's characteristic temperature is astronomically high and it is always frozen. Only the two rotations about axes perpendicular to the bond are active. A nonlinear molecule, by contrast, has three genuine rotational axes with comparable moments of inertia, so it gets three rotational degrees of freedom.

Where does Mayer's relation C_P − C_V = R come from physically?

At constant volume, all the heat you add goes into raising the internal energy, so CV=dU/dTC_V=\mathrm dU/\mathrm dT. At constant pressure, the gas also expands as it warms, doing work P dVP\,\mathrm dV against its surroundings; you must supply that work in addition to raising UU. For an ideal gas PV=RTPV=RT per mole, so P dV=R dTP\,\mathrm dV=R\,\mathrm dT, meaning you need exactly RR more joules per kelvin at constant pressure than at constant volume. Hence CP=CV+RC_P=C_V+R. The relation is universal for ideal gases because it depends only on PV=RTPV=RT, not on the internal structure of the molecule.

How does the speed of sound reveal the number of degrees of freedom?

Sound is an adiabatic (fast, no heat exchange) compression wave, so its speed depends on the adiabatic index γ=CP/CV\gamma=C_P/C_V, via cs=γRT/Mc_s=\sqrt{\gamma RT/M}. Since γ=(f+2)/f\gamma=(f+2)/f depends only on the active degree-of-freedom count ff, measuring the sound speed at known TT and MM pins down ff. This is exactly how the puzzle of Newton's too-low sound-speed prediction was resolved: Newton assumed isothermal compression (γ=1\gamma=1) and got 290 m/s290\,\mathrm{m/s}; Laplace realized the compression is adiabatic and used γ=1.4\gamma=1.4, obtaining the correct 343 m/s343\,\mathrm{m/s}.

Does the heat capacity of a monatomic gas ever change with temperature?

Not for its translational motion — a point atom has only 3 translational degrees of freedom, and translation has an essentially continuous energy spectrum, so it is always active and CV=32RC_V=\tfrac32 R from near absolute zero up to where the atom ionizes or dissociates. This is why noble gases are the cleanest test of equipartition: CV=32RC_V=\tfrac32 R and γ=5/3\gamma=5/3 hold over an enormous temperature range. Electronic excitation can add a contribution at very high temperatures, but for ordinary conditions a monatomic gas has a genuinely temperature-independent heat capacity — the one case where the classical prediction of constancy is actually correct.

Why does classical equipartition predict constant heat capacity when experiments show it rising with temperature?

Because equipartition is purely classical: it treats every quadratic mode as always available, with no notion of a minimum excitation energy. Real molecular modes are quantized, and a mode contributes to the heat capacity only once kBTk_BT exceeds its level spacing. Below that threshold the mode is "frozen" and contributes nothing; as temperature rises it thaws and its 12R\tfrac12 R or RR switches on. So the true heat capacity is a staircase rising with temperature, while classical equipartition sees only the flat treads and predicts each is permanent. The risers — the temperature dependence — are a quantum effect, treated fully in the failure-of-equipartition article.

Before this