Vibrational Heat Capacity and Einstein Modes

Introduction

A bond vibration can store energy only in quantised steps. At low temperature, thermal energy may be too small to populate its first excited state, so the vibration contributes almost nothing to heat capacity even though it retains zero-point motion. At high temperature, many levels populate and classical equipartition returns. A single Einstein mode provides an exact model for this crossover.

Core explanation

For one harmonic vibrational mode of frequency ν, level energies are ε v = hν(v + 1/2), with v = 0, 1, 2, … . Define θ vib = hν/k and x = θ vib/T. The partition sum is q vib = Σ v exp[−x(v + 1/2)] = exp(−x/2)/(1 − exp(−x)). If energy is measured relative to the zero-point level, the thermal part is q thermal = 1/(1 − e^(−x)); both conventions give the same heat capacity.

Differentiating ln q gives mean mode energy U vib = hν/2 + hν/(e^x − 1). The first term is zero-point energy and does not change with T if ν is fixed. Differentiation gives C V,vib = kx²e^x/(e^x − 1)². At x ≫ 1, excitation is suppressed and C V,vib tends toward zero approximately as kx²e^(−x). At x ≪ 1, the expression tends to k, the classical contribution of one quadratic kinetic plus one quadratic potential term.

For a molecule with several approximately independent normal modes, vibrational partition factors multiply and the corresponding heat capacities add under the harmonic approximation. Different modes can have very different frequencies and therefore become thermally active at different temperatures. A stiff bond with a high ν has a large θ vib and remains frozen over a wider temperature range than a soft low-frequency bend. Real anharmonicity and mode coupling become important at high excitation and near dissociation.

The zero-point term matters to absolute internal energy and reaction energy if vibrational frequencies differ between reactants and products. It does not contribute to heat capacity at fixed ν because heat capacity is the temperature derivative. A temperature-dependent effective frequency or structural change would require additional care.

Step-by-step reasoning

Calculate θ vib from ν or wavenumber, keeping units consistent. Form x = θ vib/T. Decide whether energy includes zero-point, then use the corresponding q consistently. Compute U and C V from the formulas or derivatives. Check low- and high-T limits, and sum contributions only for approximately separable modes.

Visual explanation

Draw equally spaced vibrational levels, starting at hν/2 above the potential minimum. At low T show nearly all dots on v = 0; at high T spread dots across several levels. Plot C V/k against T/θ vib, rising smoothly from zero to one. Mark zero-point energy as a constant vertical offset that does not change the plotted heat-capacity curve.

Real-world analogy

A device that accepts coins only in fixed denominations may stay at its lowest setting when available money is far below the first price jump. When the budget grows, more settings become affordable. A vibrational mode's Boltzmann populations behave similarly, although molecules occupy quantum states probabilistically rather than making purchases.

Real-world example

High-frequency stretching modes in simple diatomic gases often contribute little to room-temperature heat capacity, while lower-frequency vibrations in larger flexible molecules can contribute more. Infrared spectral frequencies provide the ν values needed to estimate characteristic temperatures and calculate mode contributions.

Why?

Heat capacity measures the change in mean energy when T changes. A large energy gap relative to kT prevents populations from shifting, so C V is small. Once kT is large relative to the gap, occupation changes readily and the discrete sum approaches a continuous classical oscillator with energy kT above its fixed zero-point reference.

Common misconception

Zero heat-capacity contribution at low T does not mean zero vibrational energy; the oscillator still has hν/2. Another error is to add only kT/2 at high T because there is one vibrational coordinate; its kinetic and potential quadratic terms each contribute kT/2. The energy-zero convention changes q but not C V.

Worked example

Take a mode with θ vib = 1000 K at T = 300 K. Then x = 3.333 and e^x ≈ 28.0. Its heat-capacity fraction is C V/k ≈ x²e^x/(e^x − 1)² = (11.11)(28.0)/(27.0²) ≈ 0.427. The mode is partially activated: neither the frozen limit zero nor the high-T limit k is accurate. Its mean thermal excitation energy is hν/(e^x − 1) ≈ hν/27, in addition to zero-point energy.

Quick check

1. If a mode's frequency doubles at fixed T, what happens to x and its tendency to be thermally excited? Answer: θ vib and x both double. The first gap becomes larger relative to kT, so excited-level population and heat-capacity contribution generally decrease.

Exam focus

Write x = hν/(kT) or hcν̃/(kT) with correct frequency or wavenumber units. Distinguish total U from thermal energy above zero-point. Use C V = kx²e^x/(e^x − 1)² and check that it approaches zero at low T and k at high T. Sum modes only when separability is justified.

Advanced insight

For a polyatomic molecule, the number of normal vibrational modes is 3N atoms − 6 if nonlinear and 3N atoms − 5 if linear. Their frequencies are not equal, so a one-temperature activation picture is oversimplified. A measured heat-capacity curve can encode a spectrum of mode frequencies rather than one Einstein temperature.

Summary

One harmonic mode has q = e^(−x/2)/(1 − e^(−x)), U = hν/2 + hν/(e^x − 1) and C V = kx²e^x/(e^x − 1)², with x = θ vib/T. The mode's thermal heat capacity is suppressed when T ≪ θ vib and approaches k when T ≫ θ vib. Zero-point energy remains but adds no constant-frequency heat capacity.

Practice questions

1. Find the high-temperature vibrational C V for three independent harmonic modes per molecule. Answer: Each approaches k, so their combined contribution approaches 3k per molecule or 3R per mole, excluding translation and rotation. 2. Why does changing the energy zero by subtracting hν/2 not alter C V? Answer: It subtracts a temperature-independent constant from U. Differentiating U with respect to T removes that constant. 3. A mode has θ vib/T = 10. Is the classical limit appropriate? Answer: No. Excitation probability is strongly suppressed by e⁻¹⁰, and its heat-capacity contribution is far below k.