The Vibrational Partition Function

Introduction

Molecules vibrate, but at room temperature most of them are not vibrating any more vigorously than quantum mechanics demands. The vibrational partition function tells us exactly how many vibrational states are thermally accessible, and the answer is usually "about one". That single number explains why the heat capacities of N₂ and O₂ are close to 5/2 R rather than 7/2 R, and why heavy, floppy molecules such as I₂ behave differently. Because the harmonic oscillator has equally spaced levels, q V can be summed exactly in closed form.

Core explanation

The energy levels. In the harmonic approximation a diatomic molecule has vibrational levels

E v = (v + ½)hcν̃, v = 0, 1, 2, …

where ν̃ is the vibrational wavenumber. Each level is non-degenerate. The lowest level lies ½hcν̃ above the bottom of the potential well: this is the zero-point energy .

Choosing the zero of energy. Partition functions are usually written with energies measured from the lowest state, so that q → 1 as T → 0. Measuring from v = 0, the levels are ε v = vhcν̃, and

q V = Σ e^(−vβhcν̃) = 1 + x + x² + … with x = e^(−βhcν̃)

This is a geometric series with ratio less than one, so it sums exactly:

q V = 1/(1 − e^(−βhcν̃)) = 1/(1 − e^(−θ V/T))

If instead energies are measured from the bottom of the well, every term gains a factor e^(−βhcν̃/2), giving q V = e^(−θ V/2T)/(1 − e^(−θ V/T)). The two forms describe the same physics; what matters is using one energy zero consistently, especially when combining partition functions for reactions.

The vibrational temperature. θ V = hcν̃/k sets the scale. Since hc/k = 1.4388 cm K, θ V (in K) is 1.4388 × ν̃ (in cm⁻¹).

Molecule ν̃ / cm⁻¹ θ V / K --- --- --- H₂ 4401 6332 N₂ 2358 3393 Cl₂ 560 806 I₂ 214.5 309

Limits. When T ≪ θ V, e^(−θ V/T) is tiny and q V ≈ 1: only the ground level is populated. When T ≫ θ V, expanding the exponential gives q V ≈ T/θ V = kT/(hcν̃), the classical result, and many levels are occupied.

Populations. The fraction of molecules in level v is p v = e^(−vθ V/T)/q V. The ground-state fraction is simply p₀ = 1/q V = 1 − e^(−θ V/T).

Polyatomic molecules. A molecule with N atoms has 3N − 6 normal modes (3N − 5 if linear). In the harmonic approximation the modes are independent, so q V is the product of one factor per mode, with degenerate modes appearing as repeated factors. CO₂ has four modes: symmetric stretch, antisymmetric stretch, and a doubly degenerate bend at 667 cm⁻¹ that contributes a squared factor.

Formulae

q V = 1/(1 − e^(−θ V/T)) (zero at v = 0). θ V = hcν̃/k. High-T limit: q V ≈ T/θ V. Ground-state population: p₀ = 1 − e^(−θ V/T). Polyatomic: q V = Π over modes of 1/(1 − e^(−θ i/T)).

Step-by-step reasoning

To evaluate q V for a mode:

1. Convert ν̃ to θ V = 1.4388 cm K × ν̃. 2. Form the ratio θ V/T. 3. Compute e^(−θ V/T). 4. Evaluate q V = 1/(1 − e^(−θ V/T)). 5. Check the result: q V ≥ 1 always, close to 1 for stiff bonds and close to T/θ V for soft ones.

Visual explanation

Picture an evenly spaced ladder of vibrational levels inside a parabola. Draw a horizontal line at height kT above the lowest rung. For N₂ at 298 K the line barely rises a tenth of the way to the first rung; for I₂ it reaches almost to the first rung, so several rungs carry noticeable populations that decay geometrically upward.

Real-world analogy

A vending machine that only accepts coins of one large denomination sells nothing to someone holding small change. A vibration is similar: energy arrives in quanta of hcν̃, and if typical collisions offer only kT, much smaller than the quantum, the vibration simply cannot be bought.

Real-world example

Vibrational excitation matters in combustion and atmospheric science. In the hot gases of a flame, N₂ and CO₂ vibrations become significantly populated, raising heat capacities and affecting flame temperatures. The low-frequency CO₂ bend is also thermally active at atmospheric temperatures, which is relevant to how CO₂ absorbs and emits infrared radiation.

Why?

Why is q V so close to 1 for most small molecules at room temperature? Stretching vibrations of light, strongly bonded molecules have wavenumbers of 1000–4000 cm⁻¹, while kT/hc at 298 K is only about 207 cm⁻¹. The Boltzmann factor for the first excited level is therefore tiny.

Common misconception

"Because q V ≈ 1, molecules are not vibrating at room temperature." They are: the ground state has zero-point motion. q V ≈ 1 means only that vibrational excitation above v = 0 is rare, so vibrations contribute little to thermal energy.

Worked example

Question: Calculate q V for Cl₂ (ν̃ = 560 cm⁻¹) at 500 K, and the fraction of molecules in v = 0. (kT/hc = 347.5 cm⁻¹ at 500 K.)

Reasoning: θ V/T = 560/347.5 = 1.611. e^(−1.611) = 0.200. q V = 1/(1 − 0.200) = 1.25. p₀ = 1/q V = 0.80.

Answer: q V ≈ 1.25; about 80 % of Cl₂ molecules are in the vibrational ground state at 500 K.

Quick check

1. For N₂ at 298 K, θ V/T ≈ 11.4. Is q V closer to 1 or to 10? Answer: Very close to 1, since e^(−11.4) ≈ 1 × 10⁻⁵, so almost every molecule is in v = 0.

Exam focus

Examiners expect the geometric-series derivation, the high-temperature limit q V ≈ kT/(hcν̃), and careful statement of the energy zero. Show θ V/T explicitly in numerical work; it makes errors in unit conversion obvious.

Advanced insight

Real vibrations are anharmonic: levels converge toward the dissociation limit, so the true sum differs from the harmonic result at high temperature. For accurate thermochemistry, spectroscopic constants such as ωe and ωexe are used, or the vibrational sum is carried out directly over measured levels. Low-frequency torsions are often better treated as hindered rotors than as oscillators.

Summary

For a harmonic oscillator with energies measured from v = 0, q V = 1/(1 − e^(−θ V/T)), where θ V = hcν̃/k. Stiff, light molecules have large θ V and q V ≈ 1 at room temperature; soft modes approach the classical limit T/θ V. Measuring from the well bottom introduces the zero-point factor e^(−θ V/2T). Polyatomic molecules multiply one factor per normal mode.

Practice questions

1. Calculate θ V for I₂ (ν̃ = 214.5 cm⁻¹). Answer: θ V = 1.4388 × 214.5 ≈ 309 K. 2. Evaluate q V for I₂ at 298 K and the fraction in v = 0. Answer: θ V/T = 1.035, e^(−1.035) = 0.355, so q V = 1/0.645 ≈ 1.55 and p₀ ≈ 0.645. 3. Show that q V → T/θ V when T ≫ θ V. Answer: For small y = θ V/T, e^(−y) ≈ 1 − y, so 1 − e^(−y) ≈ y and q V ≈ 1/y = T/θ V. 4. How many vibrational factors appear in q V for H₂O and for CO₂? Answer: H₂O has 3N − 6 = 3 modes; CO₂, being linear, has 3N − 5 = 4 modes (the bend counted twice).