14.2 Quantum chemical calculations for reaction rates

Quantum Chemical Calculations for Reaction Rates

Quantum chemical calculations predict reaction rates from first principles by solving the Schrödinger equation and mapping out potential energy surfaces. This lets you determine energy barriers and locate transition states without ever running an experiment. Combined with transition state theory, these calculations yield rate constants that complement and sometimes guide experimental kinetic studies.

Theoretical Basis of Quantum Calculations

Quantum chemistry applies quantum mechanics directly to chemical systems. The Schrödinger equation describes how electrons and nuclei behave, giving you wave functions and energy levels for a molecular system. In practice, solving this equation exactly is impossible for anything beyond a hydrogen atom, so approximations are essential.

The Born-Oppenheimer approximation is the most fundamental of these. It separates electronic and nuclear motion by assuming nuclei are essentially stationary compared to electrons. This lets you solve the electronic structure problem for a fixed set of nuclear positions, then repeat at different geometries.

• A potential energy surface (PES) maps the energy of a system as a function of nuclear coordinates
  • Minima on the PES correspond to stable species: reactants, products, and intermediates at their equilibrium geometries
  • Saddle points represent transition states, which sit at the highest energy point along the minimum-energy reaction path
• Reaction rates depend directly on the energy barrier separating reactants from the transition state
  • Higher barriers mean slower rates, since molecules need more energy to reach the transition state and proceed to products
  • Quantum tunneling allows molecules to pass through the barrier at energies below its classical height. This is especially important for light-atom transfers like H-atom transfer and proton-coupled electron transfer.

Application of Quantum Methods

Electronic structure methods solve the Schrödinger equation (approximately) for a given nuclear configuration. The main families differ in how they handle electron-electron interactions:

• Hartree-Fock (HF) is the simplest ab initio approach. It treats each electron as moving in the average field of all others, which means it neglects electron correlation. This often leads to overestimated barrier heights.
• Density functional theory (DFT) includes electron correlation effects through exchange-correlation functionals. It strikes a practical balance between accuracy and computational cost. Common functionals include B3LYP (a general-purpose hybrid functional) and M06-2X (better for barrier heights and noncovalent interactions).
• Post-HF methods systematically recover correlation energy at greater computational expense:
  • Møller-Plesset perturbation theory (MP2) adds a second-order correction to HF
  • Coupled cluster methods (CCSD, CCSD(T)) are often considered the "gold standard" for thermochemical accuracy, but scale steeply with system size

Once you've chosen a method, you need to locate the key structures on the PES. Here's the typical workflow:

1. Geometry optimization of reactants and products by minimizing the energy on the PES
2. Transition state optimization to locate the saddle point. This requires a reasonable initial guess for the TS geometry and uses algorithms that search for a first-order saddle point.
3. Frequency calculation at each stationary point. Reactants and products should have all real frequencies. A valid transition state is characterized by exactly one imaginary frequency, corresponding to motion along the reaction coordinate.
4. Intrinsic reaction coordinate (IRC) calculation starting from the transition state. This traces the minimum-energy path in both directions, confirming that the TS actually connects the intended reactants and products.

Rate Constants from Transition States

Transition state theory (TST) connects the quantum chemical results to a rate constant. The central equation is the Eyring equation:

$k = \frac{k_B T}{h} e^{-\Delta G^{\ddagger}/RT}$

where:

• $k_B$ = Boltzmann constant ($1.381 \times 10^{-23}$ J/K)
• $h$ = Planck's constant ($6.626 \times 10^{-34}$ J·s)
• $\Delta G^{\ddagger}$ = Gibbs free energy of activation (the free energy difference between the transition state and reactants)
• $T$ = temperature, $R$ = gas constant

Quantum chemical calculations supply the pieces you need to evaluate this equation:

• Electronic energies at the reactant and transition state geometries give the potential energy component of the barrier
• Vibrational frequencies and rotational constants from frequency calculations feed into statistical mechanics expressions for translational, rotational, and vibrational partition functions
• Thermochemical corrections convert the raw electronic energy difference into $\Delta G^{\ddagger}$ at a given temperature

For reactions involving light atoms (especially hydrogen), tunneling corrections improve the TST rate constant. Two common approaches:

• Wigner correction: $\kappa = 1 + \frac{1}{24}\left(\frac{h\nu^{\ddagger}}{k_BT}\right)^2$, where $\nu^{\ddagger}$ is the magnitude of the imaginary frequency at the transition state. This is a simple multiplicative factor but only works well for small tunneling contributions.
• Eckart barrier model handles asymmetric barriers (where reactant and product energies differ) and captures larger tunneling effects more reliably.

Accuracy of Quantum Predictions

Computational predictions are only useful if you understand their reliability. Validation against experiment is essential.

• Kinetic isotope effects (KIEs) are a particularly sensitive test. Replacing H with D changes vibrational frequencies and tunneling contributions, so comparing computed and measured KIEs reveals whether your model captures the right physics. Primary KIEs (bond breaking/forming) are typically larger than secondary KIEs (adjacent bonds).
• Temperature dependence of rate constants produces Arrhenius plots. Comparing computed and experimental activation energies and pre-exponential factors tests both the barrier height and the entropy of activation.

Accuracy depends heavily on the level of theory and basis set (the mathematical functions used to represent orbitals):

• Larger basis sets and higher-level methods generally improve results, but at significant computational cost
• Benchmark studies compare methods against high-accuracy reference data for specific reaction classes: organic reactions, enzymatic catalysis, atmospheric chemistry, and so on. These benchmarks guide method selection for new problems.

Several sources of error deserve attention:

• Anharmonicity in vibrational modes can matter for floppy molecules or hydrogen-bonded systems, where the harmonic approximation used in standard frequency calculations breaks down
• Non-Born-Oppenheimer effects become important near conical intersections, where two electronic states become degenerate
• Conformational sampling is critical for flexible molecules: if you miss a low-energy conformer of the reactant or TS, your barrier will be wrong
• Solvent effects can shift barriers substantially. Implicit solvation models (like PCM or SMD) are computationally cheap but approximate; explicit solvent molecules add accuracy at much greater cost