Unit 9: Statistical Thermodynamics of Chemical Equilibrium & Ideal Gas Reactions

The chemical potential \(\mu\) is the fundamental driving force for chemical reactions and phase equilibria.

Microscopic Expression for Chemical Potential

In the canonical ensemble, the chemical potential of component \(j\) is:

For an ideal gas mixture of indistinguishable particles:

Differentiating with respect to \(N_j\):

Therefore, the chemical potential per molecule is:

On a molar basis (\(\mu_j^{\text{molar}} = N_A \mu_j\)):

Standard Chemical Potential \(\mu^\circ\)

Recall that \(q_j = q_{j, \text{trans}} q_{j, \text{int}}\) where \(q_{j, \text{trans}} = \frac{V}{\Lambda_j^3}\). Using the ideal gas equation \(V = \frac{N_j k_B T}{P_j}\):

where \(q_j^\circ\) is the standard molar partition function evaluated at standard pressure \(P^\circ = 1\text{ bar}\). Substituting into the chemical potential:

This precisely reproduces the classical thermodynamic pressure-dependent chemical potential equation, providing its exact microscopic definition:

Consider a general reversible gas-phase chemical reaction:

where \(\nu_j\) are the stoichiometric coefficients (positive for products, negative for reactants).

Condition of Chemical Equilibrium

At thermodynamic equilibrium at constant \(T\) and \(P\), the Gibbs free energy of the reaction is at a minimum:

Substitute the microscopic expression for chemical potential \(\mu_j = -k_B T \ln(q_j / N_j)\):

Combining terms using logarithm laws:

Rearranging into terms of number densities \(\rho_j = N_j / V\):

The left-hand side is the equilibrium concentration quotient in molecules per unit volume:

Standard Pressure Equilibrium Constant \(K_p(T)\)

Using partial pressures \(P_j = \rho_j k_B T\):

where \(\Delta \nu = \sum_j \nu_j\) is the change in moles of gas. In terms of the standard partition functions \(q_j^\circ = q_j(V = V^\circ)\) where \(V^\circ = \frac{N_A k_B T}{P^\circ}\):

where \(\Delta \varepsilon_0\) is the difference in zero-point ground state energies between products and reactants. This remarkable formula enables the absolute theoretical calculation of chemical equilibrium constants directly from spectroscopic molecular constants with zero adjustable parameters.

In calculating partition functions for isolated molecules, each species usually references its energy to its own ground state (\(\varepsilon_{0, j} = 0\)). However, in a chemical reaction, atoms rearrange, making it imperative that all participating species share a common universal energy zero.

Common Energy Reference: Dissociated Atoms

Let the universal energy zero be the completely separated neutral atoms at infinite distance. For molecule \(j\), its ground state lies at energy \(-\mathcal{D}_{0, j}\) below the atomic dissociation limit, where \(\mathcal{D}_{0, j}\) is the spectroscopic dissociation energy (including zero-point vibrational energy). When all molecular energies are referenced to their individual ground states, the ground state of species \(j\) has energy \(\varepsilon_{0, j}\). The reaction energy difference at absolute zero is:

The Zero-Point Shift Factor

If \(q_j'\) represents the partition function referenced to its own ground state (\(q_j' = \sum e^{-\beta(\varepsilon - \varepsilon_0)}\)), then the partition function referenced to the universal zero is:

Substituting into the equilibrium constant formula:

• The prefactor \(\prod (q_j'^\circ / N_A)^{\nu_j}\) reflects the ratio of available quantum phase space (entropy factor).
• The exponential Boltzmann factor \(e^{-\Delta_r \varepsilon_0 / k_B T}\) reflects the electronic energy change (enthalpy factor).

At low temperatures, the exponential factor dominates, favoring species with the lowest ground-state energy (strongest chemical bonds). At high temperatures, the partition function ratio dominates, favoring species with higher multiplicity and greater density of states.

Isotopic exchange reactions provide an exceptionally pure test of statistical thermodynamics because the electronic potential energy surface is identical under the Born-Oppenheimer approximation. Consider the classic hydrogen-deuterium exchange equilibrium:

Cancellation of Electronic and Potential Factors

Because isotopes share identical electronic structures:

The electronic partition functions cancel identically (\(q_{\text{elec}} = 1\)). The only energy difference arises from the zero-point vibrational energies (ZPE):

Factorization of the Equilibrium Constant

The equilibrium constant is:

1. Translational factor:

2. Rotational factor:

Notice the symmetry numbers: \(\sigma_{\text{H}_2} = 2, \sigma_{\text{D}_2} = 2, \sigma_{\text{HD}} = 1\). Therefore, \(\frac{\sigma_{\text{H}_2} \sigma_{\text{D}_2}}{\sigma_{\text{HD}}^2} = \frac{2 \times 2}{1^2} = 4\). Because \(I = \mu R_e^2\) and \(R_e\) is identical, \(\frac{I_{\text{HD}}^2}{I_{\text{H}_2} I_{\text{D}_2}} = \frac{\mu_{\text{HD}}^2}{\mu_{\text{H}_2} \mu_{\text{D}_2}} = \frac{(2/3)^2}{(1/2)(1)} = \frac{4/9}{1/2} = \frac{8}{9}\). Multiplying:

3. Product of translation and rotation:

At high temperature where \(\Delta\varepsilon_0 / k_B T \rightarrow 0\) and vibrational functions approach unity, the mass factors cancel via the Teller-Redlich product rule, leaving:

The factor of 4 is purely entropic (statistical multiplicity).

Consider the thermal dissociation of a homonuclear diatomic gas:

Examples include \(\text{I}_2 \rightleftharpoons 2\text{I}\), \(\text{Br}_2 \rightleftharpoons 2\text{Br}\), and \(\text{H}_2 \rightleftharpoons 2\text{H}\).

General Statistical Formula

The equilibrium constant \(K_p\) is:

where \(D_0\) is the ground-state dissociation energy of the molecule. Evaluating each partition function:

1. Atomic partition function \(q_{\text{X}}\):

2. Diatomic partition function \(q_{\text{X}_2}\):

Ratio of Translational Partition Functions

Because \(m_{\text{X}_2} = 2 m_{\text{X}}\):

Substituting \(V^\circ = \frac{N_A k_B T}{P^\circ}\) and assembling all components:

Notice the physical features:

• As \(T\) increases, the rotational factor \(T\) in the denominator is overwhelmed by the \(T^{5/2}\) translational prefactor and the exponential \(e^{-D_0 / k_B T}\), driving complete dissociation at high temperatures.
• The factor of 2 in the numerator arises from the symmetry number \(\sigma = 2\) of the homonuclear reactant.

Chemical kinetics can be derived from equilibrium statistical thermodynamics via Henry Eyring, Michael Polanyi, and Eugene Wigner's Transition State Theory (TST).

The Activated Complex Hypothesis

Consider an elementary bimolecular reaction:

TST posits that reactants \(\text{A}\) and \(\text{B}\) exist in quasi-equilibrium with a short-lived activated complex (transition state \([\text{AB}]^\ddagger\)) residing at the saddle point of the potential energy surface. The reaction rate equals the concentration of transition states multiplied by their rate of passage across the barrier:

Factorization of the Reaction Coordinate

The transition state possesses \(3N_{\text{at}} - 1\) standard bound degrees of freedom plus one special degree of freedom: translation along the reaction coordinate with frequency \(\nu^\ddagger\). In the harmonic limit, the partition function for this motion is:

Factoring this out of the transition state partition function:

where \(q^\ddagger\) contains the remaining \(3N_{\text{at}} - 1\) degrees of freedom. The quasi-equilibrium concentration is:

The Eyring Rate Constant Equation

Substituting into the rate equation \(v = \nu^\ddagger [\text{AB}^\ddagger] = k(T) [\text{A}][\text{B}]\), the crossing frequency \(\nu^\ddagger\) cancels identically:

where \(\kappa\) is the transmission coefficient (typically \(\approx 1\)). In thermodynamic formulation:

where \(\Delta S^{\ddagger\circ}\) and \(\Delta H^{\ddagger\circ}\) are the standard entropy and enthalpy of activation. This establishes that reaction rates depend not only on barrier height (\(\Delta H^\ddagger\)) but also on the activation entropy (\(\Delta S^\ddagger\))—the geometric tightness or looseness of the transition state relative to reactants.

The classical van 't Hoff equation governs how the equilibrium constant varies with temperature:

Statistical Mechanics Derivation

From the statistical formula:

Differentiating with respect to \(T\):

Recall the statistical relation for molar internal energy:

Thus:

Substituting into the derivative:

Now incorporate the \(P V\) work term for ideal gases where \(H_j = U_j + R T\):

Taking into account the volume derivative in \(q_j^\circ\) under constant pressure yields precisely:

Statistical Behavior of \(\Delta_r H^\circ(T)\)

Because the heat capacity varies with temperature due to the quantum activation of vibrational modes (\(\Delta_r C_P^\circ(T) = \sum \nu_j C_{P, j}(T)\)), the reaction enthalpy is not constant:

Statistical mechanics provides the exact non-linear van 't Hoff curve over thousands of Kelvins without requiring empirical polynomial fits.

Real gases deviate from the ideal gas law \(P V = N k_B T\) due to intermolecular attraction (van der Waals dispersion) and short-range Pauli repulsive cores.

The Virial Equation of State

Heike Kamerlingh Onnes represented real gas behavior using a power series in molar density \(\rho = n / V\):

where \(B_2(T)\) is the Second Virial Coefficient (governing pairwise intermolecular interactions) and \(B_3(T)\) is the Third Virial Coefficient (three-body interactions).

The Classical Configuration Integral

In statistical mechanics, the canonical partition function for \(N\) interacting particles of mass \(m\) with pairwise potential \(U(\mathbf{r}_1, \dots, \mathbf{r}_N) = \sum_{i < j} u(r_{i j})\) is:

where \(Z_N\) is the Configuration Integral:

Mayer Cluster Expansion and the Second Virial Coefficient

Joseph Mayer introduced the Mayer \(f\)-function:

• When particles are far apart (\(r_{i j} \rightarrow \infty\)), \(u(r_{i j}) \rightarrow 0 \implies f_{i j} \rightarrow 0\).
• At short distances, \(f_{i j}\) is non-zero, restricting cluster integrals to overlapping particles.

Expanding \(Z_N\) in products of \(f_{i j}\) yields the exact statistical expression for the second virial coefficient:

Physical Regimes of \(B_2(T)\)

1. Low Temperatures: Attractive potential dominates (\(u(r) < 0 \implies e^{-u/k_B T} > 1\)), making the integrand positive, so \(B_2(T) < 0\). Intermolecular attraction reduces gas pressure below ideal values.

2. High Temperatures: Repulsive core dominates (\(u(r) \gg 0 \implies e^{-u/k_B T} \approx 0\)), making the integrand negative, so \(B_2(T) > 0\). Volume exclusion increases pressure above ideal values.

3. The Boyle Temperature \(T_B\):

The unique temperature where attraction and repulsion balance exactly:

At \(T = T_B\), real gases obey the ideal gas law across an extended pressure range.

## Advanced Mathematical Supplement: Quantum Tunneling Rate Corrections in TST

Non-Classical Reaction Coordinates: Quantum Tunneling Corrections in TST

Classical Transition State Theory assumes that all reacting trajectories pass over the potential energy barrier. For light atoms (particularly protons \(\text{H}^+\), hydrogen atoms \(\text{H}^\bullet\), and hydride ions \(\text{H}^-\)), quantum mechanical tunneling through the barrier significantly accelerates reaction rates at ambient and sub-ambient temperatures.

The Wigner Semiclassical Tunneling Correction

In 1932, Eugene Wigner expanded the quantum transmission coefficient in powers of Planck's constant \(\hbar\). For an inverted parabolic barrier of imaginary barrier frequency \(\nu^\ddagger = \frac{1}{2\pi}\sqrt{\frac{k^\ddagger}{\mu^\ddagger}}\):

When \(\frac{h \nu^\ddagger}{k_B T} < 1\), the Wigner formula provides a rapid, accurate tunneling correction:

The Eckart Barrier Model for Deep Tunneling

For high, narrow barriers where \(h \nu^\ddagger > k_B T\), tunneling from levels below the barrier apex dominates. Carl Eckart modeled the reaction profile with an asymmetric continuous potential:

The quantum transmission probability \(P(E)\) has an exact analytical solution in terms of hypergeometric functions:

where \(\alpha = \frac{1}{2}\sqrt{E/C}\), \(\beta = \frac{1}{2}\sqrt{(E - A)/C}\), and \(\delta = \frac{1}{2}\sqrt{(B - C)/C}\) with \(C = \frac{h^2}{8 m L^2}\). The thermally averaged transmission coefficient is:

In enzyme-catalyzed hydrogen transfer reactions, \(\kappa_{\text{Eckart}}\) can exceed \(100\), explaining immense kinetic isotope effects (\(k_H / k_D \sim 10 - 80\)) that classical Arrhenius kinetics fails to capture.

## Research Monograph: Non-Adiabatic Ultrafast Dynamics & Trajectory Surface Hopping

Breakdown of the Born-Oppenheimer Approximation

The Born-Oppenheimer approximation fails catastrophically whenever multiple electronic potential energy surfaces approach each other closely (\(\Delta E \le \hbar \omega_{\text{vib}}\)), especially near conical intersections. In these regions, non-adiabatic derivative coupling vectors:

diverge, mediating ultrafast radiationless transitions (internal conversion and intersystem crossing) on femtosecond timescales (\(10 - 100\text{ fs}\)).

Tully's Fewest Switches Surface Hopping (FSSH)

John Tully formulated the most widely used semiclassical method for non-adiabatic chemical dynamics:

1. Nuclear Motion: Classical nuclei move on a single active electronic potential surface \(E_k(\mathbf{R})\) governed by Newton's equations \(\mathbf{F} = -\nabla E_k\).

2. Electronic Amplitudes: The electronic wavefunction \(\Psi(t) = \sum c_j(t) \psi_j\) is integrated simultaneously along the trajectory via the time-dependent Schrödinger equation:

3. Stochastic Hopping: At each time step, the probability of hopping from current state \(k\) to state \(j\) is:

4. Energy Conservation: When a hop occurs, the nuclear momentum is rescaled along the direction of the non-adiabatic coupling vector \(\mathbf{d}_{j k}\) to conserve total energy.

FSSH calculations explain how human retinal isomerizes in \(200\text{ fs}\) during vision and how DNA dissipates UV radiation safely as heat within \(100\text{ fs}\) to prevent photochemical carcinogenesis.