Divide-and-Conquer Method for Instanton Rate Theory

Nonadiabatic Tunneling in Chemical Reactions

Quantum tunneling can have a dramatic effect on chemical reaction rates. In nonadiabatic reactions such as electron transfers or spin crossovers, nuclear tunneling effects can be even stronger than for adiabatic proton transfers. Ring-polymer instanton theory enables molecular simulations of tunneling in full dimensionality and has been shown to be far more reliable than commonly used separable approximations. First-principles instanton calculations predict significant nonadiabatic tunneling of heavy atoms even at room temperature and give excellent agreement with experimental measurements for the intersystem crossing of two nitrenes in cryogenic matrix isolation, the spin-forbidden relaxation of photoexcited thiophosgene in the gas phase, and singlet oxygen deactivation in water at ambient conditions. Finally, an outlook of further theoretical developments is discussed.

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.

Instanton theory for Fermi's golden rule and beyond

Instanton theory provides a semiclassical approximation for computing quantum tunnelling effects in complex molecular systems. It is typically applied to proton-transfer reactions for which the Born-Oppenheimer approximation is valid. However, many processes in physics, chemistry and biology, such as electron transfers, are non-adiabatic and are correctly described instead using Fermi's golden rule. In this work, we discuss how instanton theory can be generalized to treat these reactions in the golden-rule limit. We then extend the theory to treat fourth-order processes such as bridge-mediated electron transfer and apply the method to simulate an electron moving through a model system of three coupled quantum dots. By comparison with benchmark quantum calculations, we demonstrate that the instanton results are much more reliable than alternative approximations based on superexchange-mediated effective coupling or a classical sequential mechanism. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

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.

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.

Instanton theory of ground-state tunneling splittings with general paths

We derive a multidimensional instanton theory for calculating ground-state tunneling splittings in Cartesian coordinates for general paths. It is an extension of the method by Mil'nikov and Nakamura [J. Chem. Phys. 115, 6881 (2001)] to include asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. The approach avoids multiple expensive matrix diagonalizations to converge the fluctuation prefactor in the ring-polymer instanton (RPI) method, and instead replaces them by an integration of a Riccati differential equation. When combined with the string method for locating instantons, we avoid the need to converge the calculation with respect to the imaginary time period of the semiclassical orbit, thereby reducing the number of convergence parameters of the optimized object to just one: the number of equally spaced system replicas used to represent the instanton path. The entirety of the numerical effort is thus concentrated in optimizing the shape of the path and evaluating hessians along the path, which is a dramatic improvement over RPI. In addition to the standard instanton approximations, we neglect the coupling of vibrational modes to external rotations. The method is tested on the model potential of malonaldehyde and on the water dimer and trimer, giving close agreement with RPI at a much-reduced cost.