Derivation of instanton rate theory from first principles

A Perspective on the Investigation of Spectroscopy and Kinetics of Complex Molecular Systems with Semiclassical Approaches

In this Perspective we show that semiclassical methods provide a rigorous hierarchical way to study the vibrational spectroscopy and kinetics of complex molecular systems. The time averaged approach to spectroscopy and the semiclassical transition state theory for kinetics, which have been first adopted and then further developed in our group, provide accurate quantum results on rigorous physical grounds and can be applied even when dealing with a large number of degrees of freedom. In spectroscopy, the multiple coherent, divide-and-conquer, and adiabatically switched semiclassical approaches have practically permitted overcoming issues related to the convergence of results. In this Perspective we demonstrate the possibility of studying the semiclassical vibrational spectroscopy of a molecule adsorbed on an anatase (101) surface, a system made of 51 atoms. In kinetics, the semiclassical transition state theory is able to account for anharmonicity and the coupling between the reactive and bound modes. Our group has developed this technique for practical applications involving the study of phenomena like kinetic isotope effect, heavy atom tunneling, and elusive conformer lifetimes. Here, we show that our multidimensional anharmonic quantum approach is able to tackle on-the-fly the thermal kinetic rate constant of a 135 degree-of-freedom system. Overall, semiclassical methods open up the possibility to describe at the quantum mechanical level systems characterized by hundreds of degrees of freedom leading to the accurate spectroscopic and kinetic description of biomolecules and complex molecular systems.

Extending non-adiabatic rate theory to strong electronic couplings in the Marcus inverted regime

Electron transfer reactions play an essential role in many chemical and biological processes. Fermi's golden rule (GR), which assumes that the coupling between electronic states is small, has formed the foundation of electron transfer rate theory; however, in short range electron/energy transfer reactions, this coupling can become very large, and, therefore, Fermi's GR fails to make even qualitatively accurate rate predictions. In this paper, I present a simple modified GR theory to describe electron transfer in the Marcus inverted regime at arbitrarily large electronic coupling strengths. This theory is based on an optimal global rotation of the diabatic states, which makes it compatible with existing methods for calculating GR rates that can account for nuclear quantum effects with anharmonic potentials. Furthermore, the optimal GR (OGR) theory can also be combined with analytic theories for non-adiabatic rates, such as Marcus theory and Marcus-Levich-Jortner theory, offering clear physical insights into strong electronic coupling effects in non-adiabatic processes. OGR theory is also tested on a large set of spin-boson models and an anharmonic model against exact quantum dynamics calculations, where it performs well, correctly predicting rate turnover at large coupling strengths. Finally, an example application to a boron-dipyrromethane-anthracene photosensitizer reveals that strong coupling effects inhibit excited state charge recombination in this system, reducing the rate of this process by a factor of 4. Overall, OGR theory offers a new approach to calculating electron transfer rates at strong couplings, offering new physical insights into a range of non-adiabatic processes.

A time averaged semiclassical approach to IR spectroscopy

We propose a new semiclassical approach to the calculation of molecular IR spectra. The method employs the time averaging technique of Kaledin and Miller upon symmetrization of the quantum dipole-dipole autocorrelation function. Spectra at high and low temperatures are investigated. In the first case, we are able to point out the possible presence of hot bands in the molecular absorption line shape. In the second case, we are able to reproduce accurate IR spectra as demonstrated by a calculation of the IR spectrum of the water molecule, which is within 4% of the exact intensity. Our time averaged IR spectra can be directly compared to time averaged semiclassical power spectra as shown in an application to the CO2 molecule, which points out the differences between IR and power spectra and demonstrates that our new approach can identify active IR transitions correctly. Overall, the method features excellent accuracy in calculating absorption intensities and provides estimates for the frequencies of vibrations in agreement with the corresponding power spectra. In perspective, this work opens up the possibility to interface the new method with the semiclassical techniques developed for power spectra, such as the divide-and-conquer one, to get accurate IR spectra of complex and high-dimensional molecular systems.

A personal perspective of the present status and future challenges facing thermal reaction rate theory

Reaction rate theory has been at the center of physical chemistry for well over one hundred years. The evolution of the theory is not only of historical interest. Reliable and accurate computation of reaction rates remains a challenge to this very day, especially in view of the development of quantum chemistry methods, which predict the relevant force fields. It is still not possible to compute the numerically exact rate on the fly when the system has more than at most a few dozen anharmonic degrees of freedom, so one must consider various approximate methods, not only from the practical point of view of constructing numerical algorithms but also on conceptual and formal levels. In this Perspective, I present some of the recent analytical results concerning leading order terms in an ℏ2m series expansion of the exact rate and their implications on various approximate theories. A second aspect has to do with the crossover temperature between tunneling and thermal activation. Using a uniform semiclassical transmission probability rather than the "primitive" semiclassical theory leads to the conclusion that there is no divergence problem associated with a "crossover temperature." If one defines a semiclassical crossover temperature as the point at which the tunneling energy of the instanton equals the barrier height, then it is a factor of two higher than its previous estimate based on the "primitive" semiclassical approximation. In the low temperature tunneling regime, the uniform semiclassical theory as well as the "primitive" semiclassical theory were based on the classical Euclidean action of a periodic orbit on the inverted potential. The uniform semiclassical theory wrongly predicts that the "half-point," which is the energy at which the transmission probability equals 1/2, for any barrier potential, is always the barrier energy. We describe here how augmenting the Euclidean action with constant terms of order ℏ2 can significantly improve the accuracy of the semiclassical theory and correct this deficiency. This also leads to a deep connection with and improvement of vibrational perturbation theory. The uniform semiclassical theory also enables an extension of the quantum version of Kramers' turnover theory to temperatures below the "crossover temperature." The implications of these recent advances on various approximate methods used to date are discussed at length, leading to the conclusion that reaction rate theory will continue to challenge us both on conceptual and practical levels for years to come.

Comparison of Matsubara dynamics with exact quantum dynamics for an oscillator coupled to a dissipative bath

We report the first numerical calculations in which converged Matsubara dynamics is compared directly with exact quantum dynamics with no artificial damping of the time-correlation functions (TCFs). The system treated is a Morse oscillator coupled to a harmonic bath. We show that, when the system-bath coupling is sufficiently strong, the Matsubara calculations can be converged by explicitly including up to M = 200 Matsubara modes, with the remaining modes included as a harmonic "tail" correction. The resulting Matsubara TCFs are in near-perfect agreement with the exact quantum TCFs, for non-linear as well as linear operators, at a temperature at which the TCFs are dominated by quantum thermal fluctuations. These results provide compelling evidence that incoherent classical dynamics can arise in the condensed phase at temperatures at which the statistics are dominated by quantum (Boltzmann) effects, as a result of smoothing of imaginary-time Feynman paths. The techniques developed here may also lead to efficient methods for benchmarking system-bath dynamics in the overdamped regime.

Improved microcanonical instanton theory

Canonical (thermal) instanton theory is now routinely applicable to complex gas-phase reactions and allows for the accurate description of tunnelling in highly non-separable systems. Microcanonical instanton theory is by contrast far less well established. Here, we demonstrate that the best established microcanonical theory [S. Chapman, B. C. Garrett and W. H. Miller, J. Chem. Phys., 1975, 63, 2710-2716], fails to accurately describe the deep-tunnelling regime for systems where the frequencies of the orthogonal modes change rapidly along the instanton path. By taking a first principles approach to the derivation of microcanonical instanton theory, we obtain an improved method, which accurately recovers the thermal instanton rate when integrated over energy. The resulting theory also correctly recovers the separable limit and can be thought of as an instanton generalisation of Rice-Ramsperger-Kassel-Marcus (RRKM) theory. When combined with the density-of-states approach [W. Fang, P. Winter and J. O. Richardson, J. Chem. Theory Comput., 2021, 17, 40-55], this new method can be straightforwardly applied to real molecular systems.

Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. I. Theory

Reactions involving adsorbates on metallic surfaces and impurities in bulk metals are ubiquitous in a wide range of technological applications. The theoretical modeling of such reactions presents a formidable challenge for theory because nuclear quantum effects (NQEs) can play a prominent role and the coupling of the atomic motion with the electrons in the metal gives rise to important non-adiabatic effects (NAEs) that alter atomic dynamics. In this work, we derive a theoretical framework that captures both NQEs and NAEs and, due to its high efficiency, can be applied to first-principles calculations of reaction rates in high-dimensional realistic systems. More specifically, we develop a method that we coin ring polymer instanton with explicit friction (RPI-EF), starting from the ring polymer instanton formalism applied to a system-bath model. We derive general equations that incorporate the spatial and frequency dependence of the friction tensor and then combine this method with the ab initio electronic friction formalism for the calculation of thermal reaction rates. We show that the connection between RPI-EF and the form of the electronic friction tensor presented in this work does not require any further approximations, and it is expected to be valid as long as the approximations of both underlying theories remain valid.

Nonadiabatic instanton rate theory beyond the golden-rule limit

Fermi's golden rule describes the leading-order behaviour of the reaction rate as a function of the diabatic coupling. Its asymptotic (ℏ →0) limit is the semiclassical golden-rule instanton rate theory, which rigorously approximates nuclear quantum effects, lends itself to efficient numerical computation and gives physical insight into reaction mechanisms. However the golden rule by itself becomes insufficient as the strength of the diabatic coupling increases, so higher-order terms must be additionally considered. In this work we give a first-principles derivation of the next-order term beyond the golden rule, represented as a sum of three components. Two of them lead to new instanton pathways that extend the golden-rule case and, among other factors, account for the effects of recrossing on the full rate. The remaining component derives from the equilibrium partition function and accounts for changes in potential energy around the reactant and product wells due to diabatic coupling. The new semiclassical theory demands little computational effort beyond a golden-rule instanton calculation. It makes it possible to rigorously assess the accuracy of the golden-rule approximation and sets the stage for future work on general semiclassical nonadiabatic rate theories.

Determination of concerted or stepwise mechanism of hydrogen tunneling from isotope effects: Departure between experiment and theory

Isotope substitution is an important experimental technique that offers deep insight into reaction mechanisms, as the measured kinetic isotope effects (KIEs) can be directly compared with theory. For multiple proton transfer processes, there are two types of mechanisms: stepwise transfer and concerted transfer. The Bell-Limbach model provides a simple theory to determine whether the proton transfer mechanism is stepwise or concerted from KIEs. Recent STM experiments have studied the proton switching process in water tetramers on NaCl(001). Theoretical studies predict that this process occurs via a concerted mechanism, however, the experimental KIEs resemble the Bell-Limbach model for stepwise tunneling, raising question on the underlying mechanism or the validity of the model. We study this system using ab initio instanton theory, and in addition to thermal rates, we also considered microcanonical rates, as well as tunneling splittings. Instanton theory predicts a concerted mechanism, and the KIEs for tunneling rates (both thermal and microcanonical) upon deuteration are consistent with the Bell-Limbach model for concerted tunneling, but could not explain the experiments. For tunneling splittings, partial and full deuteration changes the size of it in a similar fashion to how it changes the rates. We further examined the Bell-Limbach model in another system, porphycene, which has both stepwise and concerted tunneling pathways. The KIEs predicted by instanton theory are again consistent with the Bell-Limbach model. This study highlights differences between KIEs in stepwise and concerted tunneling, and the discrepancy between theory and recent STM experiments. New theory/experiments are desired to settle this problem.

Semiclassical descriptions of rotational transitions in natural and shifted angles: Analysis of unexpected results

In the semiclassical theory of rotational transitions, S-matrix elements are expressed as integrals over initial and final angles of probability amplitudes calculated along the classical paths joining these angles, before final passage to an initial value representation [W. H. Miller, J. Phys. Chem. A 105, 2942 (2001)]. These angles can be either natural angles fixing the orientation of the rotor or angles shifted with respect to the previous ones so as to vary only within the interaction region causing the transitions. The two approaches, however, were recently shown to lead to different predictions. While the theory in natural angles lacks precision and exhibits unphysical behavior, the theory in shifted angles is much more accurate and physically well behaved [L. Bonnet, J. Chem. Phys. 153, 174102 (2020)]. The present work is devoted to the analysis of this unexpected finding.

Entanglement Effect and Angular Momentum Conservation in a Nonseparable Tunneling Treatment

The important, and often dominant, role of tunneling in low temperature kinetics has resulted in numerous theoretical explorations into the methodology for predicting it. Nevertheless, there are still key aspects of the derivations that are lacking, particularly for nonseparable systems in the low temperature regime, and further explorations of the physical factors affecting the tunneling rate are warranted. In this work we obtain a closed-form rate expression for the tunneling rate constant that is a direct analog of the rigid-rotor-harmonic-oscillator expression. This expression introduces a novel "entanglement factor" that modulates the reaction rate. Furthermore, we are able to extend this expression, which is valid for nonseparable systems at low temperatures, to properly account for the conservation of angular momentum. In contrast, previous calculations have considered only vibrational transverse modes and so effectively employ a decoupled rotational partition function for the orientational modes. We also suggest a simple theoretical model to describe the tunneling effects in the vicinity of the crossover temperature (the temperature where tunneling becomes the dominating mechanism). This model allows one to naturally classify, interpret, and predict experimental data. Among other things, it quantitatively explains in simple terms the so-called "quantum bobsled" effect, also known as the negative centrifugal effect, which is related to curvature of the reaction path. Taken together, the expressions obtained here allow one to predict the thermal and E-resolved rate constants over broad ranges of temperatures and energies.

Microcanonical Tunneling Rates from Density-of-States Instanton Theory

Semiclassical instanton theory is a form of quantum transition-state theory which can be applied to the computation of thermal reaction rates in complex molecular systems including quantum tunneling effects. There have been a number of attempts to extend the theory to treat microcanonical rates. However, the previous formulations are either computationally unfeasible for large systems due to an explicit sum over states or they involve extra approximations, which make them less reliable. We propose a robust and practical microcanonical formulation called density-of-states instanton theory, which avoids the sum over states altogether. In line with the semiclassical approximations inherent to the instanton approach, we employ the stationary-phase approximation to the inverse Laplace transform to obtain the densities of states. This can be evaluated using only post-processing of the data available from a small set of instanton calculations, such that our approach remains computationally efficient. We show that the new formulation predicts results that agree well with quantum scattering theory for an atom-diatom reaction and with experiments for a photoexcited unimolecular hydrogen transfer in a Criegee intermediate. When the thermal rate is evaluated from a Boltzmann average over our new microcanonical formalism, it can overcome some problems of conventional instanton theory. In particular, it predicts a smooth transition at the crossover temperature and is able to describe bimolecular reactions with pre-reactive complexes such as CH₃OH + OH.

Calculations of quantum tunnelling rates for muonium reactions with methane, ethane and propane

Thermal rate constants for Mu + CH₄, Mu + C₂H₆ and Mu + C₃H₈ and their equivalent reactions with H were evaluated with ab initio instanton rate theory. The potential-energy surfaces are fitted using Gaussian process regression to high-level electronic-structure calculations evaluated around the tunnelling pathway. This method was able to successfully reproduce various experimental measurements for the rate constant of these reactions. However, it was not able to reproduce the faster-than-expected rate of Mu + C₃H₈ at 300 K reported by Fleming et al. [Phys. Chem. Chem. Phys., 2015, 17, 19901 and Phys. Chem. Chem. Phys., 2020, 22, 6326]. Analysis of our results indicates that the kinetic isotope effect at this temperature is not significantly influenced by quantum tunnelling. We consider many possible factors for the discrepancy between theory and experiment but conclude that in each case, the instanton approximation is unlikely to be the cause of the error. This is in part based on the good agreement we find between the instanton predictions and new multiconfigurational time-dependent Hartree (MCTDH) calculations for Mu + CH₄ using the same potential-energy surface. Further experiments will therefore be needed to resolve this issue.

Advances and challenges for experiment and theory for multi-electron multi-proton transfer at electrified solid–liquid interfaces

Understanding microscopic mechanism of multi-electron multi-proton transfer reactions at complexed systems is important for advancing electrochemistry-oriented science in the 21st century.

Revisiting nuclear tunnelling in the aqueous ferrous–ferric electron transfer

We find that golden-rule quantum transition-state theory predicts nearly an order of magnitude less tunnelling than some of the previous estimates. This may indicate that the spin-boson model of electron transfer is not valid in the quantum regime.

Parallel Implementation of Semiclassical Transition State Theory

This paper presents the parsctst code, an efficient parallel implementation of the semiclassical transition state theory (SCTST) for reaction rate constant calculations. Parsctst is developed starting from a previously presented approach for the computation of the vibrational density of states of fully coupled anharmonic molecules ( Nguyen et al. Chem. Phys. Lett. 2010 , 499 , 915 ). The parallel implementation makes it practical to tackle reactions involving more than 100 fully coupled anharmonic vibrational degrees of freedom and also includes multidimensional tunneling effects. After describing the pseudocode and demonstrating its computational efficiency, we apply the new code for estimating the rate constant of the proton transfer isomerization reaction of the 2,4,6-tri- tert-butylphenyl to 3,5-di- tert-butylneophyl. Comparison with both theoretical and experimental results is presented. Parsctst code is user-friendly and provides a significant computational time saving compared to serial calculations. We believe that parsctst can boost the application of SCTST as an alternative to the basic transition state theory for accurate kinetics modeling not only in combustion or atmospheric chemistry, but also in organic synthesis, where bigger reactive systems are encountered.

Instanton rate constant calculations using interpolated potential energy surfaces in nonredundant, rotationally and translationally invariant coordinates

A trivial flaw in the utilization of artificial neural networks in interpolating chemical potential energy surfaces (PES) whose descriptors are Cartesian coordinates is their dependence on simple translations and rotations of the molecule under consideration. A different set of descriptors can be chosen to circumvent this problem, internuclear distances, inverse internuclear distances or z-matrix coordinates are three such descriptors. The objective is to use an interpolated PES in instanton rate constant calculations, hence information on the energy, gradient, and Hessian is required at coordinates in the vicinity of the tunneling path. Instanton theory relies on smoothly fitted Hessians, therefore we use energy, gradients, and Hessians in the training procedure. A major challenge is presented in the proper back-transformation of the output gradients and Hessians from internal coordinates to Cartesian coordinates. We perform comparisons between our method, a previous approach and on-the-fly rate constant calcuations on the hydrogen abstraction from methanol and on the hydrogen addition to isocyanic acid. © 2018Wiley Periodicals, Inc.

Path Integral Energy Landscapes for Water Clusters

The energy landscapes for a discretized path integral representation of the water dimer, trimer and pentamer are characterized in terms of the localized (classical) and delocalized minima and transition states. The transition states are finite-temperature approximations to the exact instanton path, and they are typically used to calculate the tunneling splittings or reaction rates. The features of the path integral landscape are explored, thus elucidating procedures that could usefully be automated when searching for instantons in larger systems. Our work not only clarifies the role of minima and transition states in path integral calculations but also enables us to analyze the quantum-to-classical transition.

An Extended Computational Study of Criegee Intermediate-Alcohol Reactions

High-level ab initio calculations (DF-LCCSD(T)-F12a//B3LYP/aug-cc-pVTZ) are performed on a range of stabilized Criegee intermediate (sCI)-alcohol reactions, computing reaction coordinate energies, leading to the formation of α-alkoxyalkyl hydroperoxides (AAAHs). These potential energy surfaces are used to model bimolecular reaction kinetics over a range of temperatures. The calculations performed in this work reproduce the complicated temperature-dependent reaction rates of CH₂OO and (CH₃)₂COO with methanol, which have previously been experimentally determined. This methodology is then extended to compute reaction rates of 22 different Criegee intermediates with methanol, including several intermediates derived from isoprene ozonolysis. In some cases, sCI-alcohol reaction rates approach those of sCI-(H₂O)₂. This suggests that in regions with elevated alcohol concentrations, such as urban Brazil, these reactions may generate significant quantities of AAAHs and may begin to compete with sCI reactions with other trace tropospheric pollutants such as SO₂. This work also demonstrates the ability of alcohols to catalyze the 1,4-H transfer unimolecular decomposition of α-methyl substituted sCIs.

Atom tunnelling in the reaction NH3+ + H2 → NH4+ + H and its astrochemical relevance

The title reaction is involved in the formation of ammonia in the interstellar medium. We have calculated thermal rates including atom tunnelling using different rate theories. Canonical variational theory with microcanonically optimised multidimensional tunnelling was used for bimolecular rates, modelling the gas-phase reaction and also a surface-catalysed reaction of the Eley-Rideal type. Instanton theory provided unimolecular rates, which model the Langmuir-Hinshelwood type surface reaction. The potential energy was calculated on the CCSD(T)-F12 level of theory on the fly. We report thermal rates and H/D kinetic isotope effects. The latter have implications for observed H/D fractionation in molecular clouds. Tunnelling causes rate constants to be sufficient for the reaction to play a role in interstellar chemistry even at cryogenic temperature. We also discuss intricacies and limitations of the different tunnelling approximations to treat this reaction, including its pre-reactive minimum.

Approximating Matsubara dynamics using the planetary model: Tests on liquid water and ice

Matsubara dynamics is the quantum-Boltzmann-conserving classical dynamics which remains when real-time coherences are taken out of the exact quantum Liouvillian [T. J. H. Hele et al., J. Chem. Phys. 142, 134103 (2015)]; because of a phase-term, it cannot be used as a practical method without further approximation. Recently, Smith et al. [J. Chem. Phys. 142, 244112 (2015)] developed a "planetary" model dynamics which conserves the Feynman-Kleinert (FK) approximation to the quantum-Boltzmann distribution. Here, we show that for moderately anharmonic potentials, the planetary dynamics gives a good approximation to Matsubara trajectories on the FK potential surface by decoupling the centroid trajectory from the locally harmonic Matsubara fluctuations, which reduce to a single phase-less fluctuation particle (the "planet"). We also show that the FK effective frequency can be approximated by a direct integral over these fluctuations, obviating the need to solve iterative equations. This modification, together with use of thermostatted ring-polymer molecular dynamics, allows us to test the planetary model on water (gas-phase, liquid, and ice) using the q-TIP4P/F potential surface. The "planetary" fluctuations give a poor approximation to the rotational/librational bands in the infrared spectrum, but a good approximation to the bend and stretch bands, where the fluctuation lineshape is found to be motionally narrowed by the vibrations of the centroid.

Effects of tunnelling and asymmetry for system-bath models of electron transfer

We apply the newly derived nonadiabatic golden-rule instanton theory to asymmetric models describing electron-transfer in solution. The models go beyond the usual spin-boson description and have anharmonic free-energy surfaces with different values for the reactant and product reorganization energies. The instanton method gives an excellent description of the behaviour of the rate constant with respect to asymmetry for the whole range studied. We derive a general formula for an asymmetric version of the Marcus theory based on the classical limit of the instanton and find that this gives significant corrections to the standard Marcus theory. A scheme is given to compute this rate based only on equilibrium simulations. We also compare the rate constants obtained by the instanton method with its classical limit to study the effect of tunnelling and other quantum nuclear effects. These quantum effects can increase the rate constant by orders of magnitude.

Quantum tunneling during interstellar surface-catalyzed formation of water: the reaction H + H2O2 → H2O + OH

The final step of the water formation network on interstellar grain surfaces starting from the H + O₂ route is the reaction between H and H₂O₂. This reaction is known to have a high activation energy and therefore at low temperatures it can only proceed via tunneling. To date, however, no rate constants are available at temperatures below 200 K. In this work, we use instanton theory to compute rate constants for the title reaction with and without isotopic substitutions down to temperatures of 50 K. The calculations are based on density functional theory, with additional benchmarks for the activation energy using unrestricted single-reference and multireference coupled-cluster single-point energies. Gas-phase bimolecular rate constants are calculated and compared with available experimental data not only for H + H₂O₂ → H₂O + OH, but also for H + H₂O₂ → H₂ + HO₂. We find a branching ratio where the title reaction is favored by at least two orders of magnitude at 114 K. In the interstellar medium this reaction predominantly occurs on water surfaces, which increases the probability that the two reactants meet. To mimic this, one, two, or three spectator H₂O molecules are added to the system. Eley-Rideal bimolecular and Langmuir-Hinshelwood unimolecular rate constants are presented here. The kinetic isotope effects for the various cases are compared to experimental data as well as to expressions commonly used in astrochemical models. Both the rectangular barrier and the Eckart approximations lead to errors of about an order of magnitude. Finally, fits of the rate constants are provided as input for astrochemical models.

Ab initio instanton rate theory made efficient using Gaussian process regression

Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into the calculations of chemical reaction rates based on a potential-energy 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 + CH4 results and evaluate new low-temperature rates of H + C2H6 in full dimensionality at a UCCSD(T)-F12b/cc-pVTZ-F12 level.

Kinetic isotope effects and how to describe them

We review several methods for computing kinetic isotope effects in chemical reactions including semiclassical and quantum instanton theory. These methods describe both the quantization of vibrational modes as well as tunneling and are applied to the ⋅H + H2 and ⋅H + CH4 reactions. The absolute rate constants computed with the semiclassical instanton method both using on-the-fly electronic structure calculations and fitted potential-energy surfaces are also compared directly with exact quantum dynamics results. The error inherent in the instanton approximation is found to be relatively small and similar in magnitude to that introduced by using fitted surfaces. The kinetic isotope effect computed by the quantum instanton is even more accurate, and although it is computationally more expensive, the efficiency can be improved by path-integral acceleration techniques. We also test a simple approach for designing potential-energy surfaces for the example of proton transfer in malonaldehyde. The tunneling splittings are computed, and although they are found to deviate from experimental results, the ratio of the splitting to that of an isotopically substituted form is in much better agreement. We discuss the strengths and limitations of the potential-energy surface and based on our findings suggest ways in which it can be improved.

Simultaneous Deep Tunneling and Classical Hopping for Hydrogen Diffusion on Metals

Hydrogen diffusion on metals exhibits rich quantum behavior, which is not yet fully understood. Using simulations, we show that many hydrogen diffusion barriers can be categorized into those with parabolic tops and those with broad tops. With parabolic-top barriers, hydrogen diffusion evolves gradually from classical hopping, to shallow tunneling, to deep tunneling as the temperature (T) decreases, and noticeable quantum effects persist at moderate T. In contrast, with broad-top barriers quantum effects become important only at low T and the classical-to-quantum transition is sharp, at which classical hopping and deep tunneling both occur. This coexistence indicates that more than one mechanism contributes to the quantum reaction rate. The conventional definition of the classical-to-quantum crossover T is invalid for the broad tops, and we give a new definition. Extending this, we propose a model to predict the transition T for broad-top diffusion, providing a general guide for theory and experiment.

The dissociation and recombination rates of CH4 through the Ni(111) surface: The effect of lattice motion

Methane dissociation is a prototypical system for the study of surface reaction dynamics. The dissociation and recombination rates of CH₄ through the Ni(111) surface are calculated by using the quantum instanton method with an analytical potential energy surface. The Ni(111) lattice is treated rigidly, classically, and quantum mechanically so as to reveal the effect of lattice motion. The results demonstrate that it is the lateral displacements rather than the upward and downward movements of the surface nickel atoms that affect the rates a lot. Compared with the rigid lattice, the classical relaxation of the lattice can increase the rates by lowering the free energy barriers. For instance, at 300 K, the dissociation and recombination rates with the classical lattice exceed the ones with the rigid lattice by 6 and 10 orders of magnitude, respectively. Compared with the classical lattice, the quantum delocalization rather than the zero-point energy of the Ni atoms further enhances the rates by widening the reaction path. For instance, the dissociation rate with the quantum lattice is about 10 times larger than that with the classical lattice at 300 K. On the rigid lattice, due to the zero-point energy difference between CH₄ and CD₄, the kinetic isotope effects are larger than 1 for the dissociation process, while they are smaller than 1 for the recombination process. The increasing kinetic isotope effect with decreasing temperature demonstrates that the quantum tunneling effect is remarkable for the dissociation process.

Comparison of classical reaction paths and tunneling paths studied with the semiclassical instanton theory

Comparison of classical reaction paths and semiclassical instanton paths for a proton transfer reaction mechanism.

A quantum chemical study on ˙Cl-initiated atmospheric degradation of acrylonitrile

Degradation of acrylonitrile (CH₂CHCN) by reaction with atomic chlorine was studied using quantum chemical methods.

Quantum Tunneling Rates of Gas-Phase Reactions from On-the-Fly Instanton Calculations

The instanton method obtains approximate tunneling rates from the minimum-action path (known as the instanton) linking reactants to the products at a given temperature. An efficient way to find the instanton is to search for saddle-points on the ring-polymer potential surface, which is obtained by expressing the quantum Boltzmann operator as a discrete path-integral. Here we report a practical implementation of this ring-polymer form of instanton theory into the Molpro electronic-structure package, which allows the rates to be computed on-the-fly, without the need for a fitted analytic potential-energy surface. As a test case, we compute tunneling rates for the benchmark H + CH₄ reaction, showing how the efficiency of the instanton method allows the user systematically to converge the tunneling rate with respect to the level of electronic-structure theory.