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
|
Steffen J. Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package. J Chem Theory Comput 2023; 19:5334-5355. [PMID: 37555628 DOI: 10.1021/acs.jctc.3c00568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2023]
Abstract
A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.
Collapse
Affiliation(s)
- Julien Steffen
- Chair of Theoretical Chemistry, Friedrich-Alexander University Erlangen-Nürnberg, Egerlandstraße 3, 91058 Erlangen, Bavaria, Germany
| |
Collapse
|
3
|
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
|
4
|
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. II. Benchmarks and applications. J Chem Phys 2022; 156:194107. [PMID: 35597654 DOI: 10.1063/5.0088400] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], we presented the ring-polymer instanton with explicit friction (RPI-EF) method and showed how it can be connected to the ab initio electronic friction formalism. This framework allows for the calculation of tunneling reaction rates that incorporate the quantum nature of the nuclei and certain types of non-adiabatic effects (NAEs) present in metals. In this paper, we analyze the performance of RPI-EF on model potentials and apply it to realistic systems. For a 1D double-well model, we benchmark the method against numerically exact results obtained from multi-layer multi-configuration time-dependent Hartree calculations. We demonstrate that RPI-EF is accurate for medium and high friction strengths and less accurate for extremely low friction values. We also show quantitatively how the inclusion of NAEs lowers the crossover temperature into the deep tunneling regime, reduces the tunneling rates, and, in certain regimes, steers the quantum dynamics by modifying the tunneling pathways. As a showcase of the efficiency of this method, we present a study of hydrogen and deuterium hopping between neighboring interstitial sites in selected bulk metals. The results show that multidimensional vibrational coupling and nuclear quantum effects have a larger impact than NAEs on the tunneling rates of diffusion in metals. Together with Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], these results advance the calculations of dissipative tunneling rates from first principles.
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
|
5
|
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
|
6
|
NAST: Nonadiabatic Statistical Theory Package for Predicting Kinetics of Spin-Dependent Processes. Top Curr Chem (Cham) 2022; 380:15. [PMID: 35201520 DOI: 10.1007/s41061-022-00366-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 01/15/2022] [Indexed: 10/19/2022]
Abstract
We present a nonadiabatic statistical theory (NAST) package for predicting kinetics of spin-dependent processes, such as intersystem crossings, spin-forbidden unimolecular reactions, and spin crossovers. The NAST package can calculate the probabilities and rates of transitions between the electronic states of different spin multiplicities. Both the microcanonical (energy-dependent) and canonical (temperature-dependent) rate constants can be obtained. Quantum effects, including tunneling, zero-point vibrational energy, and reaction path interference, can be accounted for. In the limit of an adiabatic unimolecular reaction proceeding on a single electronic state, NAST reduces to the traditional transition state theory. Because NAST requires molecular properties at only a few points on potential energy surfaces, it can be applied to large molecular systems, used with accurate high-level electronic structure methods, and employed to study slow nonadiabatic processes. The essential NAST input data include the nuclear Hessian at the reactant minimum, as well as the nuclear Hessians, energy gradients, and spin-orbit coupling at the minimum energy crossing point (MECP) between two states. The additional computational tools included in the NAST package can be used to extract the required input data from the output files of electronic structure packages, calculate the effective Hessian at the MECP, and fit the reaction coordinate for more advanced NAST calculations. We describe the theory, its implementation, and three examples of application to different molecular systems.
Collapse
|
7
|
Mandelli G, Aieta C, Ceotto M. Heavy Atom Tunneling in Organic Reactions at Coupled Cluster Potential Accuracy with a Parallel Implementation of Anharmonic Constant Calculations and Semiclassical Transition State Theory. J Chem Theory Comput 2022; 18:623-637. [PMID: 34995057 PMCID: PMC8830048 DOI: 10.1021/acs.jctc.1c01143] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
![]()
We describe and test
on some organic reactions a parallel implementation
strategy to compute anharmonic constants, which are employed in semiclassical
transition state theory reaction rate calculations. Our software can
interface with any quantum chemistry code capable of a single point
energy estimate, and it is suitable for both minimum and transition
state geometry calculations. After testing the accuracy and comparing
the efficiency of our implementation against other software, we use
it to estimate the semiclassical transition state theory (SCTST) rate
constant of three reactions of increasing dimensionality, known as
examples of heavy atom tunneling. We show how our method is improved
in efficiency with respect to other existing implementations. In conclusion,
our approach allows SCTST rates and heavy atom tunneling at a high
level of electronic structure theory (up to CCSD(T)) to be evaluated.
This work shows how crucial the possibility to perform high level
ab initio rate evaluations can be.
Collapse
Affiliation(s)
- Giacomo Mandelli
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
8
|
Nagahata Y, Hernandez R, Komatsuzaki T. Phase space geometry of isolated to condensed chemical reactions. J Chem Phys 2021; 155:210901. [PMID: 34879678 DOI: 10.1063/5.0059618] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The complexity of gas and condensed phase chemical reactions has generally been uncovered either approximately through transition state theories or exactly through (analytic or computational) integration of trajectories. These approaches can be improved by recognizing that the dynamics and associated geometric structures exist in phase space, ensuring that the propagator is symplectic as in velocity-Verlet integrators and by extending the space of dividing surfaces to optimize the rate variationally, respectively. The dividing surface can be analytically or variationally optimized in phase space, not just over configuration space, to obtain more accurate rates. Thus, a phase space perspective is of primary importance in creating a deeper understanding of the geometric structure of chemical reactions. A key contribution from dynamical systems theory is the generalization of the transition state (TS) in terms of the normally hyperbolic invariant manifold (NHIM) whose geometric phase-space structure persists under perturbation. The NHIM can be regarded as an anchor of a dividing surface in phase space and it gives rise to an exact non-recrossing TS theory rate in reactions that are dominated by a single bottleneck. Here, we review recent advances of phase space geometrical structures of particular relevance to chemical reactions in the condensed phase. We also provide conjectures on the promise of these techniques toward the design and control of chemical reactions.
Collapse
Affiliation(s)
- Yutaka Nagahata
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Rigoberto Hernandez
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Tamiki Komatsuzaki
- Research Center of Mathematics for Social Creativity, Research Institute for Electronic Science, Hokkaido University, Kita 20 Nishi 10, Kita-ku, Sapporo 001-0 020, Japan
| |
Collapse
|
9
|
Wang C, Ulusoy IS, Aebersold LE, Wilson AK. Multi-configuration electron-nuclear dynamics: An open-shell approach. J Chem Phys 2021; 155:154103. [PMID: 34686063 DOI: 10.1063/5.0063478] [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
The multi-configuration electron-nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.
Collapse
Affiliation(s)
- Cong Wang
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Inga S Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Lucas E Aebersold
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| |
Collapse
|
10
|
Preston RJ, Gelin MF, Kosov DS. First-passage time theory of activated rate chemical processes in electronic molecular junctions. J Chem Phys 2021; 154:114108. [PMID: 33752339 DOI: 10.1063/5.0045652] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Confined nanoscale spaces, electric fields, and tunneling currents make the molecular electronic junction an experimental device for the discovery of new out-of-equilibrium chemical reactions. Reaction-rate theory for current-activated chemical reactions is developed by combining the Keldysh nonequilibrium Green's function treatment of electrons, Fokker-Planck description of the reaction coordinate, and Kramers first-passage time calculations. The nonequilibrium Green's functions (NEGF) provide an adiabatic potential as well as a diffusion coefficient and temperature with local dependence on the reaction coordinate. Van Kampen's Fokker-Planck equation, which describes a Brownian particle moving in an external potential in an inhomogeneous medium with a position-dependent friction and diffusion coefficient, is used to obtain an analytic expression for the first-passage time. The theory is applied to several transport scenarios: a molecular junction with a single reaction coordinate dependent molecular orbital and a model diatomic molecular junction. We demonstrate the natural emergence of Landauer's blowtorch effect as a result of the interplay between the configuration dependent viscosity and diffusion coefficients. The resultant localized heating in conjunction with the bond-deformation due to current-induced forces is shown to be the determining factors when considering chemical reaction rates, each of which results from highly tunable parameters within the system.
Collapse
Affiliation(s)
- Riley J Preston
- College of Science and Engineering, James Cook University, Townsville, QLD 4811, Australia
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, 310018 Hangzhou, China
| | - Daniel S Kosov
- College of Science and Engineering, James Cook University, Townsville, QLD 4811, Australia
| |
Collapse
|
11
|
Bardakcioglu R, Reiff J, Feldmaier M, Main J, Hernandez R. Thermal decay rates of an activated complex in a driven model chemical reaction. Phys Rev E 2020; 102:062204. [PMID: 33466091 DOI: 10.1103/physreve.102.062204] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Accepted: 11/15/2020] [Indexed: 11/07/2022]
Abstract
Recent work has shown that in a nonthermal, multidimensional system, the trajectories in the activated complex possess different instantaneous and time-averaged reactant decay rates. Under dissipative dynamics, it is known that these trajectories, which are bound on the normally hyperbolic invariant manifold (NHIM), converge to a single trajectory over time. By subjecting these dissipative systems to thermal noise, we find fluctuations in the saddle-bound trajectories and their instantaneous decay rates. Averaging over these instantaneous rates results in the decay rate of the activated complex in a thermal system. We find that the temperature dependence of the activated complex decay in a thermal system can be linked to the distribution of the phase space resolved decay rates on the NHIM in the nondissipative case. By adjusting the external driving of the reaction, we show that it is possible to influence how the decay rate of the activated complex changes with rising temperature.
Collapse
Affiliation(s)
- Robin Bardakcioglu
- Institut für Theoretische Physik I, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Johannes Reiff
- Institut für Theoretische Physik I, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Matthias Feldmaier
- Institut für Theoretische Physik I, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Jörg Main
- Institut für Theoretische Physik I, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Rigoberto Hernandez
- Department of Chemistry, Department of Chemical and Biomolecular Engineering, and Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
12
|
Debnath J, Parrinello M. Gaussian Mixture-Based Enhanced Sampling for Statics and Dynamics. J Phys Chem Lett 2020; 11:5076-5080. [PMID: 32510225 DOI: 10.1021/acs.jpclett.0c01125] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
We introduce an enhanced sampling method that is based on constructing a model probability density from which a bias potential is derived. The model relies on the fact that in a physical system most of the configurations visited can be grouped into isolated metastable islands. With each island we associate a distribution that is fitted to a Gaussian mixture. The different distributions are linearly combined together with coefficients that are computed self-consistently. This leads to an integrated procedure for discovering new metastable states, exploring reaction pathways, computing free energy differences, and estimating reaction rates.
Collapse
Affiliation(s)
- Jayashrita Debnath
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Universitá della Svizzera Italiana (USI), Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zürich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali, Universitá della Svizzera Italiana (USI), Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
- Instituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| |
Collapse
|
13
|
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
|
14
|
Fang W, Zarotiadis RA, Richardson JO. Revisiting nuclear tunnelling in the aqueous ferrous–ferric electron transfer. Phys Chem Chem Phys 2020; 22:10687-10698. [DOI: 10.1039/c9cp06841d] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Wei Fang
- Laboratory of Physical Chemistry
- ETH Zürich
- 8093 Zürich
- Switzerland
| | | | | |
Collapse
|
15
|
Steffen J. A new class of reaction path based potential energy surfaces enabling accurate black box chemical rate constant calculations. J Chem Phys 2019; 150:154105. [DOI: 10.1063/1.5092589] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Julien Steffen
- Institut für Physikalische Chemie, Christian-Albrechts-Universität, Olshausenstraße 40, D–24098 Kiel, Germany
| |
Collapse
|
16
|
Aieta C, Gabas F, Ceotto M. Parallel Implementation of Semiclassical Transition State Theory. J Chem Theory Comput 2019; 15:2142-2153. [PMID: 30822385 DOI: 10.1021/acs.jctc.8b01286] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Chiara Aieta
- Dipartimento di Chimica , Università degli Studi di Milano , via C. Golgi 19 , 20133 Milano , Italy
| | - Fabio Gabas
- Dipartimento di Chimica , Università degli Studi di Milano , via C. Golgi 19 , 20133 Milano , Italy
| | - Michele Ceotto
- Dipartimento di Chimica , Università degli Studi di Milano , via C. Golgi 19 , 20133 Milano , Italy
| |
Collapse
|
17
|
Fábri C, Marquardt R, Császár AG, Quack M. Controlling tunneling in ammonia isotopomers. J Chem Phys 2019; 150:014102. [DOI: 10.1063/1.5063470] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Csaba Fábri
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
- MTA-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary
- Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Roberto Marquardt
- Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Laboratoire de Chimie Quantique, Institut de Chimie UMR 7177 CNRS/Université de Strasbourg, 4, Rue Blaise Pascal CS 90032, 67081 Strasbourg Cedex, France
| | - Attila G. Császár
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
- MTA-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary
| | - Martin Quack
- Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| |
Collapse
|
18
|
Quantum scattering theory for collisional energy transfer. ACTA ACUST UNITED AC 2019. [DOI: 10.1016/b978-0-444-64207-3.00002-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
19
|
Löhle A, Kästner J. Calculation of Reaction Rate Constants in the Canonical and Microcanonical Ensemble. J Chem Theory Comput 2018; 14:5489-5498. [DOI: 10.1021/acs.jctc.8b00565] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Andreas Löhle
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart,Germany
| | - Johannes Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart,Germany
| |
Collapse
|
20
|
Steffen J, Hartke B. Cheap but accurate calculation of chemical reaction rate constants from ab initio data, via system-specific, black-box force fields. J Chem Phys 2017; 147:161701. [DOI: 10.1063/1.4979712] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Buchholz M, Grossmann F, Ceotto M. Application of the mixed time-averaging semiclassical initial value representation method to complex molecular spectra. J Chem Phys 2017; 147:164110. [DOI: 10.1063/1.4998510] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Max Buchholz
- Institut für Theoretische Physik, Technische Universität Dresden, 01062 Dresden, Germany
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Frank Grossmann
- Institut für Theoretische Physik, Technische Universität Dresden, 01062 Dresden, Germany
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
22
|
Aieta C, Ceotto M. A quantum method for thermal rate constant calculations from stationary phase approximation of the thermal flux-flux correlation function integral. J Chem Phys 2017; 146:214115. [DOI: 10.1063/1.4984099] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
23
|
Montoya-Castillo A, Reichman DR. Approximate but accurate quantum dynamics from the Mori formalism. II. Equilibrium time correlation functions. J Chem Phys 2017; 146:084110. [DOI: 10.1063/1.4975388] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - David R. Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
24
|
McConnell SR, Löhle A, Kästner J. Rate constants from instanton theory via a microcanonical approach. J Chem Phys 2017; 146:074105. [DOI: 10.1063/1.4976129] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sean R. McConnell
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Andreas Löhle
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Johannes Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
25
|
Zhao B, Guo H. State‐to‐state quantum reactive scattering in four‐atom systems. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1301] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Bin Zhao
- Department of Chemistry and Chemical Biology University of New Mexico Albuquerque NM USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology University of New Mexico Albuquerque NM USA
| |
Collapse
|
26
|
Manivannan V, Padmanaban R. Quantum wavepacket dynamics of the N( 4S) + NO( X2Π) reaction and its isotopic variants: Integral cross sections and thermal rate constants. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.06.007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
27
|
Stanton JF. Semiclassical Transition-State Theory Based on Fourth-Order Vibrational Perturbation Theory: The Symmetrical Eckart Barrier. J Phys Chem Lett 2016; 7:2708-2713. [PMID: 27358083 DOI: 10.1021/acs.jpclett.6b01239] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Semiclassical transition-state theory based on fourth-order vibrational perturbation theory (VPT4-SCTST) is applied to compute the barrier transmission coefficient for the symmetric Eckart potential. For a barrier parametrized to mimic the H2 + H exchange reaction, the results obtained are in excellent agreement with exact quantum calculations over a range of energy that extends down to roughly 1% of the barrier height, V0, where tunneling is negligible. The VPT2-SCTST treatment, which is commonly used in chemical kinetics studies, also performs quite well but already shows an error of a few percent at ca. 0.8 V0 where tunneling is still important. This suggests that VPT4-SCTST could offer an improvement over VPT2-SCTST in applications studies. However, the computational effort for VPT4-SCTST treatments of molecules is excessive, and any improvement gained is unlikely to warrant the increased effort. Nevertheless, the treatment of the symmetric Eckart barrier problem here suggests a simple modification of the usual VPT2-SCTST protocol that warrants further investigation.
Collapse
Affiliation(s)
- John F Stanton
- Institute for Theoretical Chemistry, Department of Chemistry and Biochemistry, The University of Texas at Austin , Austin, Texas 78712, United States
| |
Collapse
|
28
|
The role of torsional motion on the properties of propiolic acid and its H/D isotopic analogs: A density functional study using SCTST and a full anharmonic VPT2 model. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.01.015] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
29
|
Zhao B, Sun Z, Guo H. State-to-State Mode Specificity: Energy Sequestration and Flow Gated by Transition State. J Am Chem Soc 2015; 137:15964-70. [DOI: 10.1021/jacs.5b11404] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bin Zhao
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Zhigang Sun
- State
Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical
and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
30
|
Coughtrie DJ, Tew DP. A Gaussian wave packet phase-space representation of quantum canonical statistics. J Chem Phys 2015; 143:044102. [DOI: 10.1063/1.4927235] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David J. Coughtrie
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - David P. Tew
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
31
|
Zhao B, Guo H. Modulations of Transition-State Control of State-to-State Dynamics in the F + H2O → HF + OH Reaction. J Phys Chem Lett 2015; 6:676-680. [PMID: 26262484 DOI: 10.1021/acs.jpclett.5b00071] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The full-dimensional quantum dynamics of the F + H2O → HF + OH reaction is investigated at the state-to-state level for the first time using a transition-state wave packet method on an accurate global potential energy surface. It is found that the H2O rotation enhances the reactivity and the product-state distribution is dominated by HF vibrational excitation while the OH moiety serves effectively as a spectator. These observations underscore the transition-state control of the reaction dynamics, as both the H2O rotational and HF vibrational modes are strongly coupled to the reaction coordinate at the transition state. It is also shown that the transition-state dominance of the reaction dynamics is modulated by other features on the potential energy surface, such as the prereaction well.
Collapse
Affiliation(s)
- Bin Zhao
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
32
|
Bredtmann T, Diestler DJ, Li SD, Manz J, Pérez-Torres JF, Tian WJ, Wu YB, Yang Y, Zhai HJ. Quantum theory of concerted electronic and nuclear fluxes associated with adiabatic intramolecular processes. Phys Chem Chem Phys 2015; 17:29421-64. [DOI: 10.1039/c5cp03982g] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Example of concerted electronic (right) and nuclear (left) fluxes: isomerization of B4.
Collapse
Affiliation(s)
- Timm Bredtmann
- State Key Laboratory of Quantum Optics and Quantum Optics Devices
- Institute of Laser Spectroscopy
- Shanxi University
- Taiyuan 030006
- China
| | - Dennis J. Diestler
- Freie Universität Berlin
- Institut für Chemie und Biochemie
- 14195 Berlin
- Germany
- University of Nebraska-Lincoln
| | - Si-Dian Li
- Nanocluster Laboratory
- Institute of Molecular Science
- Shanxi University
- Taiyuan 030006
- China
| | - Jörn Manz
- State Key Laboratory of Quantum Optics and Quantum Optics Devices
- Institute of Laser Spectroscopy
- Shanxi University
- Taiyuan 030006
- China
| | | | - Wen-Juan Tian
- Nanocluster Laboratory
- Institute of Molecular Science
- Shanxi University
- Taiyuan 030006
- China
| | - Yan-Bo Wu
- Nanocluster Laboratory
- Institute of Molecular Science
- Shanxi University
- Taiyuan 030006
- China
| | - Yonggang Yang
- State Key Laboratory of Quantum Optics and Quantum Optics Devices
- Institute of Laser Spectroscopy
- Shanxi University
- Taiyuan 030006
- China
| | - Hua-Jin Zhai
- Nanocluster Laboratory
- Institute of Molecular Science
- Shanxi University
- Taiyuan 030006
- China
| |
Collapse
|
33
|
|
34
|
Vazquez FX, Talapatra S, Sension RJ, Geva E. The Entropic Origin of Solvent Effects on the Single Bond cZt-tZt Isomerization Rate Constant of 1,3,5-cis-Hexatriene in Alkane and Alcohol Solvents: A Molecular Dynamics Study. J Phys Chem B 2014; 118:7869-77. [DOI: 10.1021/jp5002644] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Francisco X. Vazquez
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, United States
| | - Surma Talapatra
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, United States
| | - Roseanne J. Sension
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, United States
| |
Collapse
|
35
|
Li Y, Suleimanov YV, Green WH, Guo H. Quantum rate coefficients and kinetic isotope effect for the reaction Cl + CH4 → HCl + CH3 from ring polymer molecular dynamics. J Phys Chem A 2014; 118:1989-96. [PMID: 24558961 DOI: 10.1021/jp501043z] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Thermal rate coefficients and kinetic isotope effect have been calculated for prototypical heavy-light-heavy polyatomic bimolecular reactions Cl + CH4/CD4 → HCl/DCl + CH3/CD3, using a recently proposed quantum dynamics approach: ring polymer molecular dynamics (RPMD). Agreement with experimental rate coefficients, which are quite scattered, is satisfactory. However, differences up to 50% have been found between the RPMD results and those obtained from the harmonic variational transition-state theory on one of the two full-dimensional potential energy surfaces used in the calculations. Possible reasons for such discrepancy are discussed. The present work is an important step in a series of benchmark studies aimed at assessing accuracy for RPMD for chemical reaction rates, which demonstrates that this novel method is a quite reliable alternative to previously developed techniques based on transition-state theory.
Collapse
Affiliation(s)
- Yongle Li
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
| | | | | | | |
Collapse
|
36
|
Ulusoy IS, Stanton JF, Hernandez R. Effects of Roaming Trajectories on the Transition State Theory Rates of a Reduced-Dimensional Model of Ketene Isomerization. J Phys Chem A 2013; 117:7553-60. [DOI: 10.1021/jp402322h] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Inga S. Ulusoy
- Center for Computational and
Molecular Science and Technology, School of Chemistry
and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - John F. Stanton
- Department of Chemistry
and Biochemistry, The University of Texas at Austin, Austin, Texas 78712-0165, United States
| | - Rigoberto Hernandez
- Center for Computational and
Molecular Science and Technology, School of Chemistry
and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
37
|
Mátyus E. On the Calculation of Resonances in Pre-Born–Oppenheimer Molecular Structure Theory. J Phys Chem A 2013; 117:7195-206. [DOI: 10.1021/jp4010696] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Edit Mátyus
- Laboratory of Molecular
Structure and Dynamics, Institute
of Chemistry, Eötvös University, P.O. Box 32, H-1518, Budapest 112, Hungary
| |
Collapse
|
38
|
Zhuang Y, Siebert MR, Hase WL, Kay KG, Ceotto M. Evaluating the Accuracy of Hessian Approximations for Direct Dynamics Simulations. J Chem Theory Comput 2012; 9:54-64. [DOI: 10.1021/ct300573h] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yu Zhuang
- Department of Computer Science, Texas Tech University, Lubbock, Texas 79409-3104, United States
| | - Matthew R. Siebert
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
| | - William L. Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
| | - Kenneth G. Kay
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
39
|
Tao G, Miller WH. Time-dependent importance sampling in semiclassical initial value representation calculations for time correlation functions. II. A simplified implementation. J Chem Phys 2012; 137:124105. [DOI: 10.1063/1.4752206] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Barth I. Translational Effects on Electronic and Nuclear Ring Currents. J Phys Chem A 2012; 116:11283-303. [DOI: 10.1021/jp305318s] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ingo Barth
- Max-Born-Institut, Max-Born-Strasse 2A, 12489 Berlin, Germany
| |
Collapse
|
41
|
Reinisch G, Miki K, Vignoles GL, Wong BM, Simmons CS. An Efficient and Accurate Formalism for the Treatment of Large Amplitude Intramolecular Motion. J Chem Theory Comput 2012; 8:2713-2724. [PMID: 22904694 PMCID: PMC3419465 DOI: 10.1021/ct300278x] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2012] [Indexed: 11/29/2022]
Abstract
We propose a general approach to describe large amplitude
motions
(LAM) with multiple degrees of freedom (DOF) in molecules or reaction
intermediates, which is useful for the computation of thermochemical
or kinetic data. The kinetic part of the LAM Lagrangian is derived
using a Z-matrix internal coordinate representation
within a new numerical procedure. This derivation is exact for a classical
system, and the uncertainties on the prediction of observable quantities
largely arise from uncertainties on the LAM potential energy surface
(PES) itself. In order to rigorously account for these uncertainties,
we present an approach based on Bayesian theory to infer a parametrized
physical model of the PES using ab initio calculations. This framework
allows for quantification of uncertainties associated with a PES model
as well as the forward propagation of these uncertainties to the quantity
of interest. A selection and generalization of some treatments accounting
for the coupling of the LAM with other internal or external DOF are
also presented. Finally, we discuss and validate the approach with
two applications: the calculation of the partition function of 1,3-butadiene
and the calculation of the high-pressure reaction rate of the CH3 + H → CH4 recombination.
Collapse
|
42
|
Manz J, Yamamoto K. A selection rule for the directions of electronic fluxes during unimolecular pericyclic reactions in the electronic ground state. Mol Phys 2012. [DOI: 10.1080/00268976.2011.648666] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
43
|
Head-Gordon M. Publications. Mol Phys 2012. [DOI: 10.1080/00268976.2012.693252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
44
|
Roy T, Mahapatra S. Quantum dynamics of H + LiH reaction and its isotopic variants. J Chem Phys 2012; 136:174313. [DOI: 10.1063/1.4707144] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
GRAY STEPHENK. CHEMICAL REACTION DYNAMICS WITH REAL WAVE PACKETS. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2012. [DOI: 10.1142/s0219633602000178] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
An approach to carrying out accurate quantum dynamics simulations of chemical reactions, termed the real wave packet (RWP) method, is outlined. The method focuses on propagation of just the real part of a complex-valued wave packet, halving computational memory and effort requirements in comparison with comparable high accuracy quantum propagation methods. Applications to 3-atom and 4-atom chemical reactions are reviewed. Potential future directions are indicated.
Collapse
Affiliation(s)
- STEPHEN K. GRAY
- Chemistry Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
46
|
Bredtmann T, Hupf E, Paulus B. Electronic fluxes during large amplitude vibrations of single, double and triple bonds. Phys Chem Chem Phys 2012; 14:15494-501. [DOI: 10.1039/c2cp42494k] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
47
|
A comparison between reduced dimensionality and mixed quantum/classical evaluations of k(T) using the flux–flux correlation function formalism. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.07.069] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
48
|
Craig IR, Thoss M, Wang H. Accurate quantum-mechanical rate constants for a linear response Azzouz-Borgis proton transfer model employing the multilayer multiconfiguration time-dependent Hartree approach. J Chem Phys 2011; 135:064504. [DOI: 10.1063/1.3624342] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
49
|
Shi Q, Zhu L, Chen L. Quantum rate dynamics for proton transfer reaction in a model system: Effect of the rate promoting vibrational mode. J Chem Phys 2011; 135:044505. [DOI: 10.1063/1.3611050] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
50
|
Manthe U. Accurate calculations of reaction rates: predictive theory based on a rigorous quantum transition state concept. Mol Phys 2011. [DOI: 10.1080/00268976.2011.564594] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|