1
|
Bridge O, Lazzaroni P, Martinazzo R, Rossi M, Althorpe SC, Litman Y. Quantum rates in dissipative systems with spatially varying friction. J Chem Phys 2024; 161:024110. [PMID: 38984959 DOI: 10.1063/5.0216823] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2024] [Accepted: 06/20/2024] [Indexed: 07/11/2024] Open
Abstract
We investigate whether making the friction spatially dependent on the reaction coordinate introduces quantum effects into the thermal reaction rates for dissipative reactions. Quantum rates are calculated using the numerically exact multi-configuration time-dependent Hartree method, as well as the approximate ring-polymer molecular dynamics (RPMD), ring-polymer instanton methods, and classical molecular dynamics. By conducting simulations across a wide range of temperatures and friction strengths, we can identify the various regimes that govern the reactive dynamics. At high temperatures, in addition to the spatial-diffusion and energy-diffusion regimes predicted by Kramer's rate theory, a (coherent) tunneling-dominated regime is identified at low friction. At low temperatures, incoherent tunneling dominates most of Kramer's curve, except at very low friction, when coherent tunneling becomes dominant. Unlike in classical mechanics, the bath's influence changes the equilibrium time-independent properties of the system, leading to a complex interplay between spatially dependent friction and nuclear quantum effects even at high temperatures. More specifically, a realistic friction profile can lead to an increase (or decrease) of the quantum (classical) rates with friction within the spatial-diffusion regime, showing that classical and quantum rates display qualitatively different behaviors. Except at very low frictions, we find that RPMD captures most of the quantum effects in the thermal reaction rates.
Collapse
Affiliation(s)
- Oliver Bridge
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Paolo Lazzaroni
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Rocco Martinazzo
- Department of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Mariana Rossi
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Stuart C Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Yair Litman
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
2
|
Fang W, Zhu YC, Cheng Y, Hao YP, Richardson JO. Robust Gaussian Process Regression Method for Efficient Tunneling Pathway Optimization: Application to Surface Processes. J Chem Theory Comput 2024; 20:3766-3778. [PMID: 38708859 PMCID: PMC11099967 DOI: 10.1021/acs.jctc.4c00158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2024] [Revised: 04/17/2024] [Accepted: 04/17/2024] [Indexed: 05/07/2024]
Abstract
Simulation of surface processes is a key part of computational chemistry that offers atomic-scale insights into mechanisms of heterogeneous catalysis, diffusion dynamics, and quantum tunneling phenomena. The most common theoretical approaches involve optimization of reaction pathways, including semiclassical tunneling pathways (called instantons). The computational effort can be demanding, especially for instanton optimizations with an ab initio electronic structure. Recently, machine learning has been applied to accelerate reaction-pathway optimization, showing great potential for a wide range of applications. However, previous methods still suffer from numerical and efficiency issues and were not designed for condensed-phase reactions. We propose an improved framework based on Gaussian process regression for general transformed coordinates, which has improved efficiency and numerical stability, and we propose a descriptor that combines internal and Cartesian coordinates suitable for modeling surface processes. We demonstrate with 11 instanton optimizations in three representative systems that the improved approach makes ab initio instanton optimization significantly cheaper, such that it becomes not much more expensive than a classical transition-state theory rate calculation.
Collapse
Affiliation(s)
- Wei Fang
- Department
of Chemistry, Shanghai Key Laboratory of Molecular Catalysis and Innovative
Materials, Fudan University, Shanghai 200438, P. R. China
- Laboratory
of Physical Chemistry, ETH Zürich, Zürich 8093, Switzerland
- State
Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical
Computational Chemistry, Dalian Institute
of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
| | - Yu-Cheng Zhu
- State
Key Laboratory for Artificial Microstructure and Mesoscopic Physics,
Frontier Science Center for Nano-optoelectronics and School of Physics, Peking University, Beijing 100871, China
| | - Yihan Cheng
- State
Key Laboratory for Artificial Microstructure and Mesoscopic Physics,
Frontier Science Center for Nano-optoelectronics and School of Physics, Peking University, Beijing 100871, China
| | - Yi-Ping Hao
- State
Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical
Computational Chemistry, Dalian Institute
of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
| | - Jeremy O. Richardson
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, Zürich 8093, Switzerland
| |
Collapse
|
3
|
Litman Y, Pós ES, Box CL, Martinazzo R, Maurer RJ, Rossi M. Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. I. Theory. J Chem Phys 2022; 156:194106. [PMID: 35597633 DOI: 10.1063/5.0088399] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Y Litman
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - E S Pós
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - C L Box
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - R Martinazzo
- Department of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - R J Maurer
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - M Rossi
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
4
|
Cheng YH, Zhu YC, Kang W, Li X, Fang W. Determination of concerted or stepwise mechanism of hydrogen tunneling from isotope effects: Departure between experiment and theory. J Chem Phys 2022; 156:124304. [DOI: 10.1063/5.0085010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | | | - Wei Kang
- Center for Applied Physics and Technology, Peking University, China
| | | | - Wei Fang
- Dalian Institute of Chemical Physics, China
| |
Collapse
|
5
|
Roque JPL, Nunes CM, Viegas LP, Pereira NAM, Pinho E Melo TMVD, Schreiner PR, Fausto R. Switching on H-Tunneling through Conformational Control. J Am Chem Soc 2021; 143:8266-8271. [PMID: 34048232 DOI: 10.1021/jacs.1c04329] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
H-tunneling is a ubiquitous phenomenon, relevant to fields from biochemistry to materials science, but harnessing it for mastering the manipulation of chemical structures still remains nearly illusory. Here, we demonstrate how to switch on H-tunneling by conformational control using external radiation. This is outlined with a triplet 2-hydroxyphenylnitrene generated in an N2 matrix at 10 K by UV-irradiation of an azide precursor. The anti-orientation of the nitrene's OH moiety was converted to syn by selective vibrational excitation at the 2ν(OH) frequency, thereby moving the H atom closer to the vicinal nitrene center. This triggers spontaneous H-tunneling to a singlet 6-imino-2,4-cyclohexadienone. Computations reveal that such fast H-tunneling occurs through crossing the triplet-to-singlet potential energy surfaces. Our experimental realization provides an exciting novel strategy to attain control over tunneling, opening new avenues for directing chemical transformations.
Collapse
Affiliation(s)
- José P L Roque
- University of Coimbra, CQC, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - Cláudio M Nunes
- University of Coimbra, CQC, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - Luís P Viegas
- University of Coimbra, CQC, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - Nelson A M Pereira
- University of Coimbra, CQC, Department of Chemistry, 3004-535 Coimbra, Portugal
| | | | - Peter R Schreiner
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany
| | - Rui Fausto
- University of Coimbra, CQC, Department of Chemistry, 3004-535 Coimbra, Portugal
| |
Collapse
|
6
|
Rossi M. Progress and challenges in ab initio simulations of quantum nuclei in weakly bonded systems. J Chem Phys 2021; 154:170902. [PMID: 34241065 DOI: 10.1063/5.0042572] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Atomistic simulations based on the first-principles of quantum mechanics are reaching unprecedented length scales. This progress is due to the growth in computational power allied with the development of new methodologies that allow the treatment of electrons and nuclei as quantum particles. In the realm of materials science, where the quest for desirable emergent properties relies increasingly on soft weakly bonded materials, such methods have become indispensable. In this Perspective, an overview of simulation methods that are applicable for large system sizes and that can capture the quantum nature of electrons and nuclei in the adiabatic approximation is given. In addition, the remaining challenges are discussed, especially regarding the inclusion of nuclear quantum effects (NQEs) beyond a harmonic or perturbative treatment, the impact of NQEs on electronic properties of weakly bonded systems, and how different first-principles potential energy surfaces can change the impact of NQEs on the atomic structure and dynamics of weakly bonded systems.
Collapse
Affiliation(s)
- Mariana Rossi
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|