The molecular partition function Z = sum_i exp(-epsilon_i / k_BT) factorizes into independent contributions -- translational, rotational, vibrational, and electronic -- when these modes are approximately separable: Z_total = Z_trans * Z_rot * Z_vib * Z_elec. Each factor connects molecular parameters to bulk thermodynamic quantities through exact statistical mechanical relations: U = k_BT^2 * d(ln Z)/dT, C_v = dU/dT, S = k_B*ln Z + U/T. The translational partition function depends on mass and volume; rotational on moments of inertia and symmetry number; vibrational on normal mode frequencies. At high temperature each quadratic degree of freedom contributes (1/2)k_BT to energy (equipartition), but at low temperature quantum effects freeze out rotational and especially vibrational modes, explaining the temperature dependence of heat capacities that classical physics could not account for.
Calculate the partition function contributions and heat capacity for a diatomic molecule like HCl at several temperatures (100 K, 300 K, 1000 K, 5000 K). Show how C_v rises from (3/2)R (translation only) toward (7/2)R as rotational and vibrational modes become thermally accessible, reproducing the experimental Cv(T) curve.
You already know that the molecular partition function Z sums Boltzmann weights over all energy levels and that it factorizes into translational, rotational, vibrational, and electronic contributions when those modes are approximately independent. The power of this factorization is that each factor has a closed-form expression built from molecular constants you can look up or measure — mass, bond length, vibrational frequency, symmetry — and once you have Z, every equilibrium thermodynamic quantity follows from differentiation or algebraic manipulation.
The translational partition function depends on the particle mass m, temperature T, and container volume V. For any molecule in a macroscopic box, Z_trans is enormous (on the order of 10^30), reflecting the vast number of thermally accessible translational states. The rotational partition function depends on moments of inertia and a symmetry number σ that prevents overcounting indistinguishable orientations — σ = 1 for heteronuclear diatomics like HCl, σ = 2 for homonuclear ones like O₂. At room temperature most molecules have fully activated rotation, but light molecules like H₂ at cryogenic temperatures reveal discrete rotational level spacing. The vibrational partition function depends on normal mode frequencies and is the most temperature-sensitive factor because vibrational energy gaps are typically large compared to k_BT at ordinary temperatures.
The bridge from partition functions to thermodynamics is a set of exact relations. Internal energy U = k_BT² ∂(ln Z)/∂T, which extracts the average energy from the statistical distribution. Heat capacity C_v = ∂U/∂T tells you how that average energy changes with temperature. Entropy S = k_B ln Z + U/T combines the counting of accessible states with their energy content. Because Z factorizes, ln Z is additive, and each mode contributes independently to U, C_v, and S.
Consider a diatomic molecule like HCl as a concrete example. At very low temperature, only translation is active and C_v = (3/2)R — three translational degrees of freedom each contributing (1/2)R, which is the equipartition theorem prediction for quadratic energy terms. As temperature rises past about 50 K, rotation switches on and C_v climbs to (5/2)R. Vibrational modes, with characteristic temperatures often above 2000 K, only contribute significantly at high T, eventually pushing C_v toward (7/2)R. This stepwise activation is purely a quantum effect: classical equipartition would predict (7/2)R at all temperatures, which contradicts experiment. The partition function formalism naturally captures the freezing out of high-energy modes at low temperature because exp(−hν/k_BT) becomes negligibly small when hν ≫ k_BT.
This framework extends directly to polyatomic molecules by including all 3N−6 (or 3N−5 for linear molecules) vibrational normal modes and the appropriate rotational constants for symmetric, spherical, or asymmetric tops. Each vibrational mode has its own characteristic temperature, so different modes activate at different temperatures — you can predict which specific vibrations contribute to the heat capacity at any given temperature simply by comparing hν_i to k_BT. This is how statistical mechanics replaces the empirical curve-fitting of classical thermodynamics with first-principles prediction from molecular structure.