Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into the calculations of chemical reaction rates based on a potentialenergy surface computed on the fly from electronic-structure theory. This approach is necessary to extend conventional transition-state theory into the deep-tunnelling regime, but it is also more computationally expensive as it requires many more ab initio calculations. We propose an approach which uses Gaussian process regression to fit the potential-energy surface locally around the dominant tunnelling pathway. The method can be converged to give the same result as from an on-the-fly ab initio instanton calculation but it requires far fewer electronic-structure calculations. This makes it a practical approach for obtaining accurate rate constants based on high-level electronic-structure methods. We show fast convergence to reproduce benchmark H + CH 4 results and evaluate new low-temperature rates of H + C 2 H 6 in full dimensionality at a UCCSD(T)-F12b/cc-pVTZ-F12 level.

We describe an integrated computational strategy aimed at providing reliable thermochemical and kinetic information on the formation processes of astrochemical complex organic molecules. The approach involves state-of-the-art quantum-mechanical computations, second-order vibrational perturbation theory, and kinetic models based on capture and transition state theory together with the master equation approach. Notably, tunneling, quantum reflection, and leading anharmonic contributions are accounted for in our model. Formamide has been selected as a case study in view of its interest as a precursor in the abiotic amino acid synthesis. After validation of the level of theory chosen for describing the potential energy surface, we have investigated several pathways of the OH + CHNH and NH + HCO reaction channels. Our results show that both reaction channels are essentially barrierless (in the sense that all relevant transition states lie below or only marginally above the reactants) and once tunneling is taken into the proper account indicate that the reaction can occur under the low temperature conditions of interstellar environments.

We present a method for calculating partition functions taking into account anharmonic contributions for systems involving both small-amplitude vibrations and hindered rotations. The Wang-Landau scheme is used in the first case, while two alternative schemes are used for hindered rotation based on imaginary time propagation and fitting of the exact energy levels as a function of quantum number. These two schemes are shown to be complementary in their ranges of applicability (in terms of the torsional rotational constant and the relevant potential). Partition functions for four different molecules are calculated and compared to simpler ones obtained using a harmonic model.

No abstract