Statistical Basis of Reaction Equilibrium

Introduction

An equilibrium constant is often calculated from tabulated Gibbs energies. Statistical thermodynamics asks where those Gibbs energies come from. Molecular ground energies favour some products, while rotational, vibrational, electronic and translational states contribute entropy and temperature dependence. A partition-function calculation combines those effects to predict a standard chemical-potential difference and thus K.

Core explanation

For a reaction Σ iν i A i = 0, the thermodynamic relation is Δ rG° = Σ iν iμ i° and K = exp(−Δ rG°/RT). Each standard μ i° can be derived from a molecular partition function together with a specified standard state. For an ideal gas, molecular mass affects translation, moments of inertia affect rotation, vibrational frequencies affect vibrational excitation and zero-point energy, and electronic degeneracy affects state counting. A reaction's K therefore depends on both energy differences and the number of accessible states.

Write a molecular partition function as an energy-reference factor times a thermal state sum: q i = exp(−ε i,0/kT)q i,thermal. When reaction stoichiometry combines species, ground-state energies contribute exp(−Δε₀/kT) in a simple per-particle reaction, while thermal state sums contribute products raised to stoichiometric powers. Translational factors carry volume or pressure dependence, so a standard pressure or concentration must be built into the final dimensionless K. Neglecting these standard-state factors can produce an apparently dimensioned “equilibrium constant” or an incorrect dependence on an arbitrary container volume.

A particularly clean example is ideal isomerisation A ⇌ B, where each molecule has the same elemental composition and mass. At equilibrium μ A = μ B. If their standard states are defined identically and their translational factors cancel, an illustrative ideal-gas ratio has K = q B,int/q A,int × exp[−(ε B,0 − ε A,0)/kT], with internal thermal partition functions referenced to their respective ground energies. A lower B ground energy favours B through the exponential; more B internal states favour B through q B,int. If rotational symmetry differs, its correction belongs in those internal factors.

The exact expression for a general reaction requires attention to stoichiometric molecule-number changes, standard pressure, zero-point and electronic energy conventions, and indistinguishable-particle counting. The chemical-potential route is safer than trying to guess a ratio from molecular q alone. Real mixtures require activities or fugacities; the molecular ideal-gas expression predicts a standard-state K, while actual composition obeys Q = K.

Step-by-step reasoning

Balance the reaction and choose standard states. For each species, list ground energy, mass, rotational constants, vibrational modes and electronic degeneracy. Build consistent thermal partition factors and the standard chemical potential, then sum ν iμ i° to find Δ rG°. Calculate K = exp(−Δ rG°/RT). For a simple isomerisation, identify which factors cancel before forming the remaining ratio.

Visual explanation

Draw two molecular energy ladders A and B with different ground heights and densities of excited levels. An arrow from the ground-energy difference points to an exponential energetic factor, and an arrow from the number of accessible upper levels points to a partition-function ratio. The two feed K, illustrating energetic preference versus state-counting preference.

Real-world analogy

Two destinations can differ both in travel cost and in the number of available routes. A cheaper destination is attractive, but many route options can also make one destination more frequently reached. Reaction equilibrium balances molecular energy and state multiplicity; unlike travel, the statistical weights are precise Boltzmann factors.

Real-world example

Conformational or structural isomers can have different relative populations at equilibrium. A lower-energy isomer may dominate at low temperature, while an isomer with more accessible conformations can gain relative population as temperature rises. Spectroscopy can supply energy gaps and rotational or vibrational data used to model such ratios.

Why?

At equilibrium, reaction extent no longer lowers Gibbs energy, so the stoichiometric sum of chemical potentials is zero. Partition functions determine those potentials from molecular state probabilities. Ground-state energy differences and entropy from accessible states therefore appear together in the macroscopic K.

Common misconception

The lowest-energy molecule need not dominate at every temperature; degeneracy and excited-state populations can matter. Another error is to use a raw q B/q A ratio for a reaction changing molecule number without considering translational standard-state factors. K must be dimensionless under a stated activity convention.

Worked example

Suppose ideal isomers A and B have the same mass and compatible standards. B's ground energy is higher by Δε₀ = kT ln 2 at the temperature of interest, but B has three times A's internal thermal partition factor: q B,int/q A,int = 3. Then K = 3e^(−ln 2) = 1.5. B is energetically disfavoured by a factor 1/2 but favoured by a factor 3 in accessible-state weighting, giving an equilibrium B/A activity ratio of 1.5 under this simple model.

Quick check

1. If A and B have equal internal state sums and B's ground energy exceeds A's by a positive Δε, which isomer is favoured at low T? Answer: A. K for A ⇌ B contains exp(−Δε/kT), which becomes very small as T falls when Δε > 0.

Exam focus

State whether partition functions include zero-point and electronic ground energies. Keep k with per-molecule energies and R with molar energies. Explain standard-state factors for reactions changing gas molecule count. Use Δ rG° = −RT ln K as the final thermodynamic check and distinguish K from the current quotient Q.

Advanced insight

High-quality statistical K calculations may include anharmonic vibrational levels, hindered internal rotation and electronically excited states. These details can matter greatly for flexible molecules even if their ground electronic energies are accurate. Comparing predicted K(T) with measured equilibrium data tests both energetic and entropic parts of the molecular model.

Summary

Partition functions connect molecular spectra and state counts to standard chemical potentials. Combining them stoichiometrically gives Δ rG° and dimensionless K = exp(−Δ rG°/RT). Ground-energy differences supply Boltzmann preference, while accessible states supply entropy. Standard states and energy references must remain consistent, especially beyond simple isomerisation.

Practice questions

1. In a simple A ⇌ B model with equal thermal internal factors, B is lower by 2kT in ground energy. What is the illustrative K? Answer: Δε₀ = ε B − ε A = −2kT, so K = exp(2) ≈ 7.39 when translational and standard-state factors cancel. 2. Why does changing the standard pressure convention alter μ i° but not a correctly calculated physical equilibrium composition? Answer: The activity definition changes consistently with the standard potentials and K. Physical chemical potentials and the condition Q = K remain invariant when the convention is transformed coherently. 3. Name two spectroscopic inputs that can affect a calculated K besides electronic ground-state energies. Answer: Rotational constants or moments of inertia and vibrational frequencies affect state sums; electronic degeneracies and low excited levels can matter too.