Classical flux integrals in transition state theory: Generalized reaction coordinates

An Efficient and Accurate Formalism for the Treatment of Large Amplitude Intramolecular Motion

We propose a general approach to describe large amplitude motions (LAM) with multiple degrees of freedom (DOF) in molecules or reaction intermediates, which is useful for the computation of thermochemical or kinetic data. The kinetic part of the LAM Lagrangian is derived using a Z-matrix internal coordinate representation within a new numerical procedure. This derivation is exact for a classical system, and the uncertainties on the prediction of observable quantities largely arise from uncertainties on the LAM potential energy surface (PES) itself. In order to rigorously account for these uncertainties, we present an approach based on Bayesian theory to infer a parametrized physical model of the PES using ab initio calculations. This framework allows for quantification of uncertainties associated with a PES model as well as the forward propagation of these uncertainties to the quantity of interest. A selection and generalization of some treatments accounting for the coupling of the LAM with other internal or external DOF are also presented. Finally, we discuss and validate the approach with two applications: the calculation of the partition function of 1,3-butadiene and the calculation of the high-pressure reaction rate of the CH₃ + H → CH₄ recombination.

Predictive Theory for Hydrogen Atom−Hydrocarbon Radical Association Kinetics

Procedures for accurately predicting the kinetics of hydrogen atom associations with hydrocarbon radicals are described and applied to a series of reactions. The approach is based on CASPT2/cc-pvdz evaluations of the orientation-dependent interaction energies within variable reaction coordinate transition state theory. One-dimensional corrections to the interaction energies are estimated from CAS+1+2/aug-cc-pvtz evaluations for the H + CH3 reaction, and a dynamical correction factor of 0.9 is applied. This corrected CASPT2 approach yields results that are within 10% of those obtained with the full CAS+1+2/aug-cc-pvtz potential for the H + CH3, H + C2H5, H + C2H3, and H + C2H reactions. New predictions are made for the H + iso-C3H7, H + tert-C4H9, H + C6H5, and H + C10H7 reactions. For the H + CH3 and H + C2H3 reactions, where the experimental values appear to be the most well-determined, theory and experiment essentially agree to within their error bars. For the other reactions, the agreement is reasonably satisfactory given the often large dispersion in the experimental results. For the reactions with saturated alkyl radicals, the theory predicts that each additional CH3 group increases the steric factor by approximately a factor of 2. In contrast, for the unsaturated radicals, the H + C6H5 and H + C10H7 high-pressure association rate coefficients are nearly identical to that for H + C2H3.