Statistical Mechanics for Chemistry
Anchor (Master): McQuarrie, Statistical Mechanics; Hill, An Introduction to Statistical Thermodynamics
Intuition [Beginner]
Thermodynamics tells you what happens at the macroscopic level -- pressure, temperature, volume, energy. But it is silent on why. Statistical mechanics bridges the gap by counting the microscopic configurations available to a system and showing how that count determines every thermodynamic property you already know.
The core idea is simple: a mole of gas contains roughly 10^23 molecules, each with its own energy. The measurable internal energy U is just the average energy across all those molecules, weighted by how likely each energy level is to be occupied. The probability of occupying a given level is given by the Boltzmann distribution, which states that lower-energy states are more populated than higher-energy states by a factor that depends on temperature.
The single most important quantity in statistical mechanics is the partition function q. Think of it as a bookkeeping device -- it sums up information about every accessible energy level. Once you know q, you can derive U, S, H, G, and the heat capacities directly, without ever measuring a calorimeter.
Visual [Beginner]
Energy levels
|
E4 ------|------ (rarely occupied at low T)
E3 ------|------
E2 ------|------ (moderately populated)
E1 ------|------ (heavily populated)
E0 ------|------ (ground state, most populated)
|
Population of each level:
Low T: High T:
E4 * E4 ****
E3 ** E3 *****
E2 **** E2 ******
E1 ******* E1 ******
E0 *************** E0 ******
As T increases, populations spread toward higher levels.
Partition function:
q = g0*exp(-E0/kT) + g1*exp(-E1/kT) + g2*exp(-E2/kT) + ...
\____________________________________________________/
sum over ALL accessible microstates
The partition function encodes how many states are thermally accessible. At low T, only the ground state contributes significantly (q is small). At high T, many states contribute (q is large). The larger q is, the more dispersed the energy, and the greater the entropy.
Worked example [Beginner]
Problem: A molecule has three energy levels at 0, 400 cm^-1, and 800 cm^-1, each with degeneracy 1. Calculate the molecular partition function at T = 300 K.
Solution:
First, convert energy units. The Boltzmann factor uses E/kT, and at 300 K:
kT = (1.381 x 10^-23 J/K)(300 K) = 4.14 x 10^-21 J
Convert wavenumbers to joules. E (J) = E (cm^-1) x hc, where hc = 1.986 x 10^-23 J cm.
E0 = 0 J E1 = (400)(1.986 x 10^-23) = 7.94 x 10^-21 J E2 = (800)(1.986 x 10^-23) = 1.59 x 10^-20 J
Now compute Boltzmann factors:
exp(-E0/kT) = exp(0) = 1.000 exp(-E1/kT) = exp(-7.94/4.14) = exp(-1.918) = 0.147 exp(-E2/kT) = exp(-15.87/4.14) = exp(-3.834) = 0.0217
The molecular partition function:
q = 1.000 + 0.147 + 0.0217 = 1.169
Interpretation: at 300 K, only the ground state is significantly populated. The "effective number of accessible states" is barely more than 1.
Check your understanding [Beginner]
Formal definition [Intermediate+]
Canonical partition function. For a system of N distinguishable, non-interacting molecules at temperature T in contact with a thermal reservoir, the canonical partition function is:
Q = q^N (distinguishable particles)
For indistinguishable particles (an ideal gas):
Q = q^N / N!
where q is the molecular partition function:
q = sum over i of g_i * exp(-epsilon_i / k_B T)
Here epsilon_i is the energy of level i, g_i is its degeneracy, and k_B is the Boltzmann constant.
Thermodynamic quantities from Q. All macroscopic thermodynamic properties follow from Q and its temperature derivative:
- Internal energy: U = k_B T^2 (d ln Q / dT)_V
- Entropy: S = k_B ln Q + U/T = k_B ln Q + k_B T (d ln Q / dT)_V
- Helmholtz free energy: A = -k_B T ln Q
- Pressure: P = k_B T (d ln Q / dV)_T
- Enthalpy: H = U + PV
- Gibbs free energy: G = A + PV
- Heat capacity: C_V = (dU/dT)_V
Molecular partition functions by degree of freedom. For an ideal gas molecule, the total molecular partition function factors:
q = q_trans * q_rot * q_vib * q_elec
Translational: q_trans = (2 pi m k_B T / h^2)^(3/2) * V
Rotational (linear molecule): q_rot = 8 pi^2 I k_B T / (sigma h^2)
where I is the moment of inertia and sigma is the symmetry number.
Rotational (nonlinear molecule): q_rot = (pi^(1/2) / sigma) * (8 pi^2 k_B T / h^2)^(3/2) * (I_A I_B I_C)^(1/2)
Vibrational (each normal mode): q_vib = 1 / (1 - exp(-h nu / k_B T))
For a molecule with s vibrational modes, q_vib is the product over all modes.
Electronic: q_elec = sum over electronic states of g_elec,i * exp(-epsilon_elec,i / k_B T)
For most molecules at ordinary temperatures, only the ground electronic state contributes and q_elec = g_elec,0.
Equipartition theorem. In the classical limit, each quadratic degree of freedom contributes (1/2)k_B T to the average energy per molecule, or (1/2)R per mole. A monatomic ideal gas has 3 translational degrees of freedom, giving U = (3/2)RT. A linear molecule adds 2 rotational degrees for (5/2)RT, and a nonlinear molecule adds 3 for 3RT.
Key results [Intermediate+]
Boltzmann distribution. The probability that a molecule occupies state i is:
p_i = (g_i exp(-epsilon_i / k_B T)) / q
This is the single most important equation in statistical thermodynamics. Every population ratio follows from it.
Sackur-Tetrode equation. The translational entropy of a monatomic ideal gas:
S = N k_B [ln(q_trans / N) + 5/2]
This is one of the few cases where absolute entropy can be calculated exactly from first principles, and it agrees with the third-law entropies measured calorimetrically.
Heat capacity from partition functions. The constant-volume heat capacity:
C_V = (partial U / partial T)_V = k_B T^2 (d^2 ln Q / dT^2)_V + 2 k_B T (d ln Q / dT)_V
At low T, vibrational contributions to C_V freeze out (they follow the Einstein or Debye model). Rotational and translational contributions reach their classical equipartition values at moderate temperatures.
Equilibrium constants from partition functions. For a reaction aA + bB -> cC + dD:
K = (q_C^c * q_D^d / q_A^a * q_B^b) * (1/V)^(Delta n) * exp(-Delta epsilon_0 / k_B T)
where Delta epsilon_0 is the difference in ground-state energies. This provides a purely molecular route to equilibrium constants.
Advanced treatment [Master]
Ensemble theory. The canonical partition function Q(N, V, T) is the Laplace transform of the microcanonical density of states Omega(E). The connection:
Q(N, V, T) = integral from 0 to infinity of Omega(E) exp(-E / k_B T) dE
This integral is dominated by the saddle point at the most probable energy, which is why the microcanonical and canonical ensembles give identical thermodynamics in the thermodynamic limit.
Ideal gas corrections beyond classical behavior. At high density or low temperature, quantum statistics become important. For fermions (e.g., electrons in a metal), the Fermi-Dirac distribution replaces the Boltzmann distribution. For bosons (e.g., He-4), the Bose-Einstein distribution applies and permits Bose-Einstein condensation. The Maxwell-Boltzmann distribution used in introductory treatments is the classical limit valid when the thermal de Broglie wavelength lambda = h / (2 pi m k_B T)^(1/2) is much smaller than the mean interparticle spacing.
Real gases and the configuration integral. For interacting molecules, the partition function includes a configuration integral:
Q = (q_int^N / N!) * (1 / V^N) * integral ... integral exp(-U(r_1,...,r_N) / k_B T) dr_1 ... dr_N
where q_int is the internal (rot + vib + elec) partition function and U(r_1,...,r_N) is the potential energy of the N-particle configuration. Expanding this integral in powers of density yields the virial equation of state. The second virial coefficient B_2(T) is directly related to the pair potential through:
B_2(T) = -2 pi integral from 0 to infinity [exp(-u(r)/k_B T) - 1] r^2 dr
Quantum statistical mechanics of vibration. The Einstein model treats each vibrational mode as an independent quantum harmonic oscillator. The Debye model replaces the discrete spectrum with a continuous distribution of modes up to a cutoff frequency nu_D. The Debye heat capacity at low T follows C_V proportional to T^3, in agreement with experiment, whereas the Einstein model predicts an exponential freeze-out that is too rapid.
Path integral formulations. For systems where quantum effects are important but perturbation theory fails, the Feynman path integral provides a formally exact approach. Each quantum particle is represented as a classical ring polymer of P beads, and the partition function becomes a classical configuration integral over all bead positions. Path integral molecular dynamics (PIMD) exploits this isomorphism to compute quantum statistical properties using classical sampling.
Connections [Master]
Statistical mechanics is the theoretical backbone that unifies thermodynamics, kinetics, and quantum chemistry:
Thermodynamics: Every thermodynamic potential (U, S, H, G, A) is derivable from the appropriate partition function. The laws of thermodynamics emerge as consequences of counting microstates.
Chemical kinetics: Transition state theory is built directly on the partition function formalism. The Eyring equation expresses the rate constant in terms of the partition function of the transition state relative to the reactants. The pre-exponential factor in the Arrhenius equation is interpretable as a ratio of partition functions.
Quantum chemistry: The energy levels epsilon_i that enter the partition function are computed by solving the Schrodinger equation. Ab initio electronic structure methods (Hartree-Fock, DFT, coupled cluster) provide the input for statistical mechanical calculations of thermochemistry.
Spectroscopy: Population distributions among energy levels, governed by the Boltzmann distribution, determine spectral line intensities. Rotational and vibrational partition functions predict temperature-dependent spectral envelopes.
Materials science: The statistical mechanics of lattice models (Ising model for magnetism, lattice gas for adsorption) connects microscopic interactions to macroscopic phase behavior and critical phenomena.
Biochemistry: Protein folding, ligand binding, and allosteric regulation are all equilibrium processes governed by free energy differences, which are computable from the partition functions of folded, unfolded, and bound states.
Bibliography [Master]
McQuarrie, D. A. Statistical Mechanics. University Science Books, 2000. The standard graduate text for chemistry students. Chapters 1--6 cover ensembles and ideal gases; Chapters 7--10 cover interacting systems.
Hill, T. L. An Introduction to Statistical Thermodynamics. Dover, 1986. Concise and rigorous. Particularly strong on the configuration integral and real gases.
Atkins, P. W. and de Paula, J. Physical Chemistry, 11th ed. Oxford University Press, 2018. Chapters 12--13 provide an intermediate-level treatment with many chemical applications.
Chandler, D. Introduction to Modern Statistical Mechanics. Oxford University Press, 1987. Elegant presentation of ensemble theory and the connection to thermodynamics. The chapter on correlation functions is essential reading.
Pathria, R. K. and Beale, P. D. Statistical Mechanics, 4th ed. Academic Press, 2021. Comprehensive reference for formal theory, including quantum statistics and phase transitions.