1
|
Richardson JO. Nonadiabatic Tunneling in Chemical Reactions. J Phys Chem Lett 2024; 15:7387-7397. [PMID: 38995660 DOI: 10.1021/acs.jpclett.4c01098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
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.
Collapse
Affiliation(s)
- Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland
| |
Collapse
|
2
|
Lawrence JE, Richardson JO. Improved microcanonical instanton theory. Faraday Discuss 2022; 238:204-235. [PMID: 35929848 DOI: 10.1039/d2fd00063f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
|
3
|
Trenins G, Richardson JO. Nonadiabatic instanton rate theory beyond the golden-rule limit. J Chem Phys 2022; 156:174115. [DOI: 10.1063/5.0088518] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- George Trenins
- ETH Zurich Department of Chemistry and Applied Biosciences, Switzerland
| | | |
Collapse
|
4
|
Ásgeirsson V, Birgisson BO, Bjornsson R, Becker U, Neese F, Riplinger C, Jónsson H. Nudged Elastic Band Method for Molecular Reactions Using Energy-Weighted Springs Combined with Eigenvector Following. J Chem Theory Comput 2021; 17:4929-4945. [PMID: 34275279 DOI: 10.1021/acs.jctc.1c00462] [Citation(s) in RCA: 90] [Impact Index Per Article: 30.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The climbing image nudged elastic band method (CI-NEB) is used to identify reaction coordinates and to find saddle points representing transition states of reactions. It can make efficient use of parallel computing as the calculations of the discretization points, the so-called images, can be carried out simultaneously. In typical implementations, the images are distributed evenly along the path by connecting adjacent images with equally stiff springs. However, for systems with a high degree of flexibility, this can lead to poor resolution near the saddle point. By making the spring constants increase with energy, the resolution near the saddle point is improved. To assess the performance of this energy-weighted CI-NEB method, calculations are carried out for a benchmark set of 121 molecular reactions. The performance of the method is analyzed with respect to the input parameters. Energy-weighted springs are found to greatly improve performance and result in successful location of the saddle points in less than a thousand energy and force evaluations on average (about a hundred per image) using the same set of parameter values for all of the reactions. Even better performance is obtained by stopping the calculation before full convergence and complete the saddle point search using an eigenvector following method starting from the location of the climbing image. This combination of methods, referred to as NEB-TS, turns out to be robust and highly efficient as it reduces the average number of energy and force evaluations down to a third, to 305. An efficient and flexible implementation of these methods has been made available in the ORCA software.
Collapse
Affiliation(s)
- Vilhjálmur Ásgeirsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Benedikt Orri Birgisson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Ragnar Bjornsson
- Max-Planck-Institute for Chemical Energy Conversion, Mülheim an der Ruhr 45470, Germany
| | - Ute Becker
- Max-Planck-Institute for Kohlenforschung, Mülheim an der Ruhr 45470, Germany
| | - Frank Neese
- Max-Planck-Institute for Kohlenforschung, Mülheim an der Ruhr 45470, Germany
| | | | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| |
Collapse
|
5
|
Fang W, Winter P, Richardson JO. Microcanonical Tunneling Rates from Density-of-States Instanton Theory. J Chem Theory Comput 2020; 17:40-55. [PMID: 33351621 DOI: 10.1021/acs.jctc.0c01118] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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 CH3OH + OH.
Collapse
Affiliation(s)
- Wei Fang
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Pierre Winter
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
6
|
Laude G, Calderini D, Welsch R, Richardson JO. Calculations of quantum tunnelling rates for muonium reactions with methane, ethane and propane. Phys Chem Chem Phys 2020; 22:16843-16854. [PMID: 32666960 DOI: 10.1039/d0cp01346c] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate constants for Mu + CH4, Mu + C2H6 and Mu + C3H8 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 + C3H8 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 + CH4 using the same potential-energy surface. Further experiments will therefore be needed to resolve this issue.
Collapse
Affiliation(s)
- Gabriel Laude
- Laboratory of Physical Chemistry, ETH Zürich, Switzerland.
| | | | | | | |
Collapse
|
7
|
Abstract
The diffusion of water molecules and clusters across the surfaces of materials is important to a wide range of processes. Interestingly, experiments have shown that on certain substrates, water dimers can diffuse more rapidly than water monomers. Whilst explanations for anomalously fast diffusion have been presented for specific systems, the general underlying physical principles are not yet established. We investigate this through a systematic ab initio study of water monomer and dimer diffusion on a range of surfaces. Calculations reveal different mechanisms for fast water dimer diffusion, which is found to be more widespread than previously anticipated. The key factors affecting diffusion are the balance of water-water versus water-surface bonding and the ease with which hydrogen-bond exchange can occur (either through a classical over-the-barrier process or through quantum-mechanical tunnelling). We anticipate that the insights gained will be useful for understanding future experiments on the diffusion and clustering of hydrogen-bonded adsorbates. The experimental observation that water dimers diffuse more rapidly than monomers across materials’ surfaces is yet to be clarified. Here the authors show by ab initio calculations classical and quantum mechanical mechanisms for faster water dimer diffusion on a broad range of metal and non-metal surfaces.
Collapse
|
8
|
Heller ER, Richardson JO. Instanton formulation of Fermi’s golden rule in the Marcus inverted regime. J Chem Phys 2020; 152:034106. [DOI: 10.1063/1.5137823] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Eric R. Heller
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
9
|
Sakaushi K, Kumeda T, Hammes-Schiffer S, Melander MM, Sugino O. Advances and challenges for experiment and theory for multi-electron multi-proton transfer at electrified solid–liquid interfaces. Phys Chem Chem Phys 2020; 22:19401-19442. [DOI: 10.1039/d0cp02741c] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Understanding microscopic mechanism of multi-electron multi-proton transfer reactions at complexed systems is important for advancing electrochemistry-oriented science in the 21st century.
Collapse
Affiliation(s)
- Ken Sakaushi
- Center for Green Research on Energy and Environmental Materials
- National Institute for Materials Science
- Ibaraki 305-0044
- Japan
| | - Tomoaki Kumeda
- Center for Green Research on Energy and Environmental Materials
- National Institute for Materials Science
- Ibaraki 305-0044
- Japan
| | | | - Marko M. Melander
- Nanoscience Center
- Department of Chemistry
- University of Jyväskylä
- Jyväskylä
- Finland
| | - Osamu Sugino
- The Institute of Solid State Physics
- the University of Tokyo
- Chiba 277-8581
- Japan
| |
Collapse
|
10
|
Winter P, Richardson JO. Divide-and-Conquer Method for Instanton Rate Theory. J Chem Theory Comput 2019; 15:2816-2825. [DOI: 10.1021/acs.jctc.8b01267] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Pierre Winter
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
11
|
Thapa MJ, Fang W, Richardson JO. Nonadiabatic quantum transition-state theory in the golden-rule limit. I. Theory and application to model systems. J Chem Phys 2019; 150:104107. [DOI: 10.1063/1.5081108] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Manish J. Thapa
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Wei Fang
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
12
|
|
13
|
|
14
|
Laude G, Calderini D, Tew DP, Richardson JO. Ab initio instanton rate theory made efficient using Gaussian process regression. Faraday Discuss 2018; 212:237-258. [PMID: 30230495 DOI: 10.1039/c8fd00085a] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
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.
Collapse
Affiliation(s)
- Gabriel Laude
- Laboratory of Physical Chemistry, ETH Zurich, Switzerland. and On exchange from School of Chemistry, University of Edinburgh, UK
| | | | - David P Tew
- Max-Planck-Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | | |
Collapse
|