Quantum reactive scattering of H + hydrocarbon reactions

A practical quantum-dynamical method is described for predicting accurate rate constants for general chemical reactions. The ab initio potential energy surfaces for these reactions can be built from a minimal number of grid points (average of 50 points) and expressed in terms of analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimised using the MP2 method with a cc-pVTZ basis set. Single point energies are calculated on the optimised geometries at the CCSD(T) level of theory with the same basis set. The dynamics of these reactions occur on effective reduced dimensionality hyper-surfaces accounting for the zero-point energy of the optimised degrees of freedom. Bonds being broken and formed are treated with explicit hyperspherical time independent quantum dynamics. Application of the method to the H + CH(4)--> H(2)+ CH(3), H + C(2)H(6)--> H(2)+ C(2)H(5), H + C(3)H(8)--> H(2)+n-C(3)H(7)/H(2)+i-C(3)H(7) and H + CH(3)OH --> H(2)+ CH(3)O/H(2)+ CH(2)OH reactions illustrate the potential of the approach in predicting rate constants, kinetic isotope effects and branching ratios. All studied reactions exhibit large quantum tunneling in the rate constants at lower temperatures. These quantum calculations compare well with the experimental results.