1
|
London N, Bu S, Johnson B, Ananth N. Mean-Field Ring Polymer Rates Using a Population Dividing Surface. J Phys Chem A 2024. [PMID: 38976564 DOI: 10.1021/acs.jpca.4c00005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Mean-field ring polymer molecular dynamics offers a computationally efficient method for the simulation of reaction rates in multilevel systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, it is necessary to ensure reactive trajectories form kinked ring polymer configurations at the dividing surface. Building on this idea, we introduce a population difference coordinate and a reactive flux expression modified to only include contributions from kinked configurations. We test the accuracy of the resulting mean-field rate theory on a series of linear vibronic coupling model systems. We demonstrate that this new kMF-RP rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.
Collapse
Affiliation(s)
- Nathan London
- Department of Chemistry, Cornell University, Ithaca, New York 14853, United States
| | - Siyu Bu
- Department of Chemistry, Cornell University, Ithaca, New York 14853, United States
| | - Britta Johnson
- Department of Chemistry, Cornell University, Ithaca, New York 14853, United States
| | - Nandini Ananth
- Department of Chemistry, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
2
|
Pollak E. A personal perspective of the present status and future challenges facing thermal reaction rate theory. J Chem Phys 2024; 160:150902. [PMID: 38639316 DOI: 10.1063/5.0199557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Accepted: 03/06/2024] [Indexed: 04/20/2024] Open
Abstract
Reaction rate theory has been at the center of physical chemistry for well over one hundred years. The evolution of the theory is not only of historical interest. Reliable and accurate computation of reaction rates remains a challenge to this very day, especially in view of the development of quantum chemistry methods, which predict the relevant force fields. It is still not possible to compute the numerically exact rate on the fly when the system has more than at most a few dozen anharmonic degrees of freedom, so one must consider various approximate methods, not only from the practical point of view of constructing numerical algorithms but also on conceptual and formal levels. In this Perspective, I present some of the recent analytical results concerning leading order terms in an ℏ2m series expansion of the exact rate and their implications on various approximate theories. A second aspect has to do with the crossover temperature between tunneling and thermal activation. Using a uniform semiclassical transmission probability rather than the "primitive" semiclassical theory leads to the conclusion that there is no divergence problem associated with a "crossover temperature." If one defines a semiclassical crossover temperature as the point at which the tunneling energy of the instanton equals the barrier height, then it is a factor of two higher than its previous estimate based on the "primitive" semiclassical approximation. In the low temperature tunneling regime, the uniform semiclassical theory as well as the "primitive" semiclassical theory were based on the classical Euclidean action of a periodic orbit on the inverted potential. The uniform semiclassical theory wrongly predicts that the "half-point," which is the energy at which the transmission probability equals 1/2, for any barrier potential, is always the barrier energy. We describe here how augmenting the Euclidean action with constant terms of order ℏ2 can significantly improve the accuracy of the semiclassical theory and correct this deficiency. This also leads to a deep connection with and improvement of vibrational perturbation theory. The uniform semiclassical theory also enables an extension of the quantum version of Kramers' turnover theory to temperatures below the "crossover temperature." The implications of these recent advances on various approximate methods used to date are discussed at length, leading to the conclusion that reaction rate theory will continue to challenge us both on conceptual and practical levels for years to come.
Collapse
Affiliation(s)
- Eli Pollak
- Chemical and Biological Physics Department, Weizmann Institute of Science, 76100 Rehovoth, Israel
| |
Collapse
|
3
|
Sadhasivam VG, Meuser L, Reichman DR, Althorpe SC. Instantons and the quantum bound to chaos. Proc Natl Acad Sci U S A 2023; 120:e2312378120. [PMID: 38032936 PMCID: PMC10710067 DOI: 10.1073/pnas.2312378120] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 10/26/2023] [Indexed: 12/02/2023] Open
Abstract
The rate at which information scrambles in a quantum system can be quantified using out-of-time-ordered correlators. A remarkable prediction is that the associated Lyapunov exponent [Formula: see text] that quantifies the scrambling rate in chaotic systems obeys a universal bound [Formula: see text]. Previous numerical and analytical studies have indicated that this bound has a quantum-statistical origin. Here, we use path-integral techniques to show that a minimal theory to reproduce this bound involves adding contributions from quantum thermal fluctuations (describing quantum tunneling and zero-point energy) to classical dynamics. By propagating a model quantum-Boltzmann-conserving classical dynamics for a system with a barrier, we show that the bound is controlled by the stability of thermal fluctuations around the barrier instanton (a delocalized structure which dominates the tunneling statistics). This stability requirement appears to be general, implying that there is a close relation between the formation of instantons, or related delocalized structures, and the imposition of the quantum-chaos bound.
Collapse
Affiliation(s)
- Vijay Ganesh Sadhasivam
- Yusuf Hamied Department of Chemistry, University of Cambridge, CambridgeCB2 1EW, United Kingdom
| | - Lars Meuser
- Yusuf Hamied Department of Chemistry, University of Cambridge, CambridgeCB2 1EW, United Kingdom
- Department of Chemistry and Applied Biosciences, ETH Zurich (Swiss Federal Institute of Technology), Zürich8093, Switzerland
| | | | - Stuart C. Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, CambridgeCB2 1EW, United Kingdom
| |
Collapse
|
4
|
Cook L, Runeson JE, Richardson JO, Hele TJH. Which Algorithm Best Propagates the Meyer-Miller-Stock-Thoss Mapping Hamiltonian for Non-Adiabatic Dynamics? J Chem Theory Comput 2023; 19:6109-6125. [PMID: 37704193 PMCID: PMC10536990 DOI: 10.1021/acs.jctc.3c00709] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Indexed: 09/15/2023]
Abstract
A common strategy to simulate mixed quantum-classical dynamics is by propagating classical trajectories with mapping variables, often using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian or the related spin-mapping approach. When mapping the quantum subsystem, the coupled dynamics reduce to a set of equations of motion to integrate. Several numerical algorithms have been proposed, but a thorough performance comparison appears to be lacking. Here, we compare three time-propagation algorithms for the MMST Hamiltonian: the Momentum Integral (MInt) (J. Chem. Phys., 2018, 148, 102326), the Split-Liouvillian (SL) (Chem. Phys., 2017, 482, 124-134), and the algorithm in J. Chem. Phys., 2012, 136, 084101 that we refer to as the Degenerate Eigenvalue (DE) algorithm due to the approximation required during derivation. We analyze the accuracy of individual trajectories, correlation functions, energy conservation, symplecticity, Liouville's theorem, and the computational cost. We find that the MInt algorithm is the only rigorously symplectic algorithm. However, comparable accuracy at a lower computational cost can be obtained with the SL algorithm. The approximation implicitly made within the DE algorithm conserves energy poorly, even for small timesteps, and thus leads to slightly different results. These results should guide future mapping-variable simulations.
Collapse
Affiliation(s)
- Lauren
E. Cook
- Department
of Chemistry, University College London, Christopher Ingold Building, London WC1H 0AJ, U.K.
| | - Johan E. Runeson
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, Zürich 8093, Switzerland
| | - Jeremy O. Richardson
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, Zürich 8093, Switzerland
| | - Timothy J. H. Hele
- Department
of Chemistry, University College London, Christopher Ingold Building, London WC1H 0AJ, U.K.
| |
Collapse
|
5
|
Hasegawa T. Nuclear Quantum Dynamics of Three-Dimensional Condensed-Phase Systems by Constant Uncertainty Molecular Dynamics. J Phys Chem Lett 2023; 14:8043-8049. [PMID: 37651720 DOI: 10.1021/acs.jpclett.3c01254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Constant uncertainty molecular dynamics (CUMD), which was developed to include nuclear quantum effects (NQEs) in molecular dynamics (MD) simulations, is extended to three-dimensional condensed-phase systems. Its applicability was verified via CUMD simulations of bulk water and ice Ih with the q-TIP4P/F potential model. The simulated radial distribution functions, proton momentum distributions, and infrared spectra were compared with those of previous studies in which the simulations were carried out on the basis of the imaginary-time path-integral formalism. We found that CUMD can function as an alternative trajectory-based simulation method for including static and dynamic NQEs in MD simulations of three-dimensional condensed-phase systems.
Collapse
Affiliation(s)
- Taisuke Hasegawa
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan
| |
Collapse
|
6
|
Fiechter M, Runeson JE, Lawrence JE, Richardson JO. How Quantum is the Resonance Behavior in Vibrational Polariton Chemistry? J Phys Chem Lett 2023; 14:8261-8267. [PMID: 37676159 PMCID: PMC10510439 DOI: 10.1021/acs.jpclett.3c01154] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2023] [Accepted: 07/28/2023] [Indexed: 09/08/2023]
Abstract
Recent experiments in polariton chemistry have demonstrated that reaction rates can be modified by vibrational strong coupling to an optical cavity mode. Importantly, this modification occurs only when the frequency of the cavity mode is tuned to closely match a molecular vibrational frequency. This sharp resonance behavior has proved to be difficult to capture theoretically. Only recently did Lindoy et al. [ Nat. Commun. 2023, 14, 2733] report the first instance of a sharp resonant effect in the cavity-modified rate simulated in a model system using exact quantum dynamics. We investigate the same model system with a different method, ring-polymer molecular dynamics (RPMD), which captures quantum statistics but treats dynamics classically. We find that RPMD does not reproduce this sharp resonant feature at the well frequency, and we discuss the implications of this finding for future studies of vibrational polariton chemistry.
Collapse
Affiliation(s)
- Marit
R. Fiechter
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, 8093 Zürich, Switzerland
| | - Johan E. Runeson
- Department
of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Joseph E. Lawrence
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, 8093 Zürich, Switzerland
| | - Jeremy O. Richardson
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, 8093 Zürich, Switzerland
| |
Collapse
|
7
|
Bocus M, Goeminne R, Lamaire A, Cools-Ceuppens M, Verstraelen T, Van Speybroeck V. Nuclear quantum effects on zeolite proton hopping kinetics explored with machine learning potentials and path integral molecular dynamics. Nat Commun 2023; 14:1008. [PMID: 36823162 PMCID: PMC9950054 DOI: 10.1038/s41467-023-36666-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Accepted: 02/10/2023] [Indexed: 02/25/2023] Open
Abstract
Proton hopping is a key reactive process within zeolite catalysis. However, the accurate determination of its kinetics poses major challenges both for theoreticians and experimentalists. Nuclear quantum effects (NQEs) are known to influence the structure and dynamics of protons, but their rigorous inclusion through the path integral molecular dynamics (PIMD) formalism was so far beyond reach for zeolite catalyzed processes due to the excessive computational cost of evaluating all forces and energies at the Density Functional Theory (DFT) level. Herein, we overcome this limitation by training first a reactive machine learning potential (MLP) that can reproduce with high fidelity the DFT potential energy surface of proton hopping around the first Al coordination sphere in the H-CHA zeolite. The MLP offers an immense computational speedup, enabling us to derive accurate reaction kinetics beyond standard transition state theory for the proton hopping reaction. Overall, more than 0.6 μs of simulation time was needed, which is far beyond reach of any standard DFT approach. NQEs are found to significantly impact the proton hopping kinetics up to ~473 K. Moreover, PIMD simulations with deuterium can be performed without any additional training to compute kinetic isotope effects over a broad range of temperatures.
Collapse
Affiliation(s)
- Massimo Bocus
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Ruben Goeminne
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Aran Lamaire
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Maarten Cools-Ceuppens
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | - Toon Verstraelen
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052, Zwijnaarde, Belgium
| | | |
Collapse
|
8
|
Anderson MC, Schile AJ, Limmer DT. Nonadiabatic transition paths from quantum jump trajectories. J Chem Phys 2022; 157:164105. [DOI: 10.1063/5.0102891] [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
We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. We detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.
Collapse
Affiliation(s)
- Michelle C. Anderson
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Addison J. Schile
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - David T. Limmer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Kavli Energy NanoSciences Institute, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
9
|
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
|
10
|
Effect of Substituents on the Homopolymerization Activity of Methyl Alkyl Diallyl Ammonium Chloride. Molecules 2022; 27:molecules27154677. [PMID: 35897851 PMCID: PMC9330561 DOI: 10.3390/molecules27154677] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Revised: 07/18/2022] [Accepted: 07/18/2022] [Indexed: 12/04/2022] Open
Abstract
Among nitrogen-containing cationic electrolytes, diallyl quaternary ammonium salt is a typical monomer with the highest positive charge density, which has attracted the most attention, especially in the research on homopolymers and copolymers of dimethyl diallyl ammonium chloride (DMDAAC), which occupy a very unique and important position. In order to improve the lipophilicity of substituted diallyl ammonium chloride monomers under the premise of high cationic charge density, the simplest, most direct, and most efficient structure design strategy was selected in this paper. Only one of the substituents on DMDAAC quaternary ammonium nitrogen was modified by alkyl; the substituents were propyl and amyl groups, and their corresponding monomers were methyl propyl diallyl ammonium chloride (MPDAAC) and methyl amyl diallyl ammonium chloride (MADAAC), respectively. The effect of substituent structure on the homopolymerization activity of methyl alkyl diallyl ammonium chloride was illustrated by quantum chemical calculation and homopolymerization rate determination experiments via ammonium persulfate (APS) as the initiator system. The results of quantum chemistry simulation showed that, with the finite increase in substituted alkyl chain length, the numerical values of the bond length and the charge distribution of methyl alkyl diallyl ammonium chloride monomer changed little, with the activation energy of the reactions in the following order: DMDAAC < MPDAAC < MADAAC. The polymerization activities measured by the dilatometer method were in the order DMDAAC > MPDAAC > MADAAC. The activation energies Ea of homopolymerization were 96.70 kJ/mol, 97.25 kJ/mol, and 100.23 kJ/mol, and the rate equation of homopolymerization of each monomer was obtained. After analyzing and comparing these results, it could be easily found that the electronic effect of substituent was not obvious, whereas the effect of the steric hindrance was dominant. The above studies have laid a good foundation for an understanding of the polymerization activity of methyl alkyl diallyl ammonium chloride monomers and the possibility of preparation and application of these polymers with high molecular weight.
Collapse
|
11
|
García-Alfonso E, Barranco M, Bonhommeau DA, Halberstadt N, Pi M, Calvo F. Clustering, collision, and relaxation dynamics in pure and doped helium nanoclusters: Density- vs particle-based approaches. J Chem Phys 2022; 157:014106. [DOI: 10.1063/5.0091942] [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 clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum–classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other.
Collapse
Affiliation(s)
- Ernesto García-Alfonso
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
| | - Manuel Barranco
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
- Department FQA, Facultat de Física, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain
- Institute of Nanoscience and Nanotechnology (IN2UB), Universitat de Barcelona, Barcelona, Spain
| | - David A. Bonhommeau
- Université de Reims Champagne Ardenne, CNRS, GSMA UMR 7331, 51100 Reims, France
| | - Nadine Halberstadt
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
| | - Martí Pi
- Department FQA, Facultat de Física, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain
- Institute of Nanoscience and Nanotechnology (IN2UB), Universitat de Barcelona, Barcelona, Spain
| | - Florent Calvo
- Université Grenoble Alpes, CNRS, LIPHY, F38000 Grenoble, France
| |
Collapse
|
12
|
Meng Q, Chen J, Ma J, Zhang X, Chen J. Adiabatic models for the quantum dynamics of surface scattering with lattice effects. Phys Chem Chem Phys 2022; 24:16415-16436. [PMID: 35766107 DOI: 10.1039/d2cp01560a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.
Collapse
Affiliation(s)
- Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Junbo Chen
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China. .,Xi'an Modern Chemistry Research Institute, China North Industries Group Corp., Ltd., East Zhangba Road 168, 710065 Xi'an, China
| | - Jianxing Ma
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Jun Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Yangqiao Road West 155, 350002 Fuzhou, China.,Fujian Science & Technology Innovation Laboratory for Optoelectronic Information of China, Optoelectronic Industry Base at High-tech Zone, 350108 Fuzhou, China
| |
Collapse
|
13
|
Arasaki Y, Takatsuka K. Nature of chemical bond and potential barrier in an invariant energy-orbital picture. J Chem Phys 2022; 156:234102. [PMID: 35732517 DOI: 10.1063/5.0088340] [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
Physical nature of the chemical bond and potential barrier is studied in terms of energy natural orbitals (ENOs), which are extracted from highly correlated electronic wavefunctions. ENO provides an objective one-electron picture about energy distribution in a molecule, just as the natural orbitals (NOs) represent one electron view about electronic charge distribution. ENO is invariant in the same sense as NO is, that is, ENOs converge to the exact ones as a series of approximate wavefunctions approach the exact one, no matter how the methods of approximation are adopted. Energy distribution analysis based on ENO can give novel insights about the nature of chemical bonding and formation of potential barriers, besides information based on the charge distribution alone. With ENOs extracted from full configuration interaction wavefunctions in a finite yet large enough basis set, we analyze the nature of chemical bonding of three low-lying electronic states of a hydrogen molecule, all being in different classes of the so-called covalent bond. The mechanism of energy lowering in bond formation, which gives a binding energy, is important, yet not the only concern for this small molecule. Another key notion in chemical bonding is whether a potential basin is well generated stiff enough to support a vibrational state(s) on it. Based on the virial theorem in the adiabatic approximation and in terms of the energy and force analyses with ENOs, we study the roles of the electronic kinetic energy and its nuclear derivative(s) on how they determine the curvature (or the force constant) of the potential basins. The same idea is applied to the potential barrier and, thereby, the transition states. The rate constant within the transition-state theory is formally shown to be described in terms of the electronic kinetic energy and the nuclear derivatives only. Thus, the chemical bonding and rate process are interconnected behind the scenes. Besides this aspect, we pay attention to (1) the effects of electron correlation that manifests itself not only in the orbital energy but also in the population of ENOs and (2) the role of nonadiabaticity (diabatic state mixing), resulting in double basins and a barrier on a single potential curve in bond formation. These factors differentiate a covalent bond into subclasses.
Collapse
Affiliation(s)
- Yasuki Arasaki
- Fukui Institute for Fundamental Chemistry, Kyoto University, 606-8103 Kyoto, Japan
| | - Kazuo Takatsuka
- Fukui Institute for Fundamental Chemistry, Kyoto University, 606-8103 Kyoto, Japan
| |
Collapse
|
14
|
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
|
15
|
Ananth N. Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics. Annu Rev Phys Chem 2022; 73:299-322. [PMID: 35081325 DOI: 10.1146/annurev-physchem-082620-021809] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This review focuses on a recent class of path-integral-based methods that simulate nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. The remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, USA;
| |
Collapse
|
16
|
Li X, Huo P. Investigating Tunneling-Controlled Chemical Reactions through Ab Initio Ring Polymer Molecular Dynamics. J Phys Chem Lett 2021; 12:6714-6721. [PMID: 34261316 DOI: 10.1021/acs.jpclett.1c01630] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We use the ab initio ring polymer molecular dynamics (RPMD) approach to investigate tunneling-controlled reactions in methylhydroxycarbene. Nuclear tunneling effects enable molecules to overcome the barriers which cannot be overcome classically. Under low-temperature conditions, intrinsic quantum tunneling effects can facilitate the chemical reaction in a pathway that is favored neither thermodynamically nor kinetically. This behavior is referred to as the tunneling-controlled chemical reaction and is regarded as the third paradigm of chemical reaction controls. In this work, we use the ab initio RPMD approach to incorporate the tunneling effects in our quantum dynamics simulations and investigate the reaction kinetics of two competitive reaction pathways at various temperatures. The reaction rate constants obtained here agree extremely well with the experimentally measured rates. We demonstrate the feasibility of using ab initio RPMD rate calculations in a realistic molecular system and provide an interesting and important example for future investigations of reaction mechanisms dominated by quantum tunneling effects.
Collapse
Affiliation(s)
- Xinyang Li
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| |
Collapse
|
17
|
Yang H, Fan WB, Fang JH, Song J, Li Y. Rate coefficients of roaming reaction H+MgH using ring polymer molecular dynamics. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2007121] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Hui Yang
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Wen-bin Fan
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Jun-hua Fang
- Department of Chemistry, Shanghai University, Shanghai 200444, China
| | - Jianing Song
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Yongle Li
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
- Devision of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena 91125, USA
| |
Collapse
|
18
|
Chen Q, Hu X, Guo H, Xie D. Theoretical H + O 3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties. Phys Chem Chem Phys 2021; 23:3300-3310. [PMID: 33506830 DOI: 10.1039/d0cp05771a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. The theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O3 → OD + O2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.
Collapse
Affiliation(s)
- Qixin Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Xixi Hu
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| |
Collapse
|
19
|
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
|
20
|
Fang W, Thapa MJ, Richardson JO. Nonadiabatic quantum transition-state theory in the golden-rule limit. II. Overcoming the pitfalls of the saddle-point and semiclassical approximations. J Chem Phys 2019; 151:214101. [DOI: 10.1063/1.5131092] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wei Fang
- Laboratory of Physical Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Manish J. Thapa
- Laboratory of Physical Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | | |
Collapse
|
21
|
Vaillant CL, Thapa MJ, Vaníček J, Richardson JO. Semiclassical analysis of the quantum instanton approximation. J Chem Phys 2019; 151:144111. [DOI: 10.1063/1.5123800] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Christophe L. Vaillant
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Manish J. Thapa
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | | |
Collapse
|
22
|
Menéndez M, Jambrina PG, Zanchet A, Verdasco E, Suleimanov YV, Aoiz FJ. New Stress Test for Ring Polymer Molecular Dynamics: Rate Coefficients of the O( 3P) + HCl Reaction and Comparison with Quantum Mechanical and Quasiclassical Trajectory Results. J Phys Chem A 2019; 123:7920-7931. [PMID: 31461272 DOI: 10.1021/acs.jpca.9b06695] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the past decade, ring polymer molecular dynamics (RPMD) has emerged as a very efficient method to determine thermal rate coefficients for a great variety of chemical reactions. This work presents the application of this methodology to study the O(3P) + HCl reaction, which constitutes a stringent test for any dynamical calculation due to rich resonant structure and other dynamical features. The rate coefficients, calculated on the 3A' and 3A″ potential energy surfaces (PESs) by Ramachandran and Peterson [ J. Chem. Phys. 2003 , 119 , 9590 ], using RPMD and quasiclassical trajectories (QCT) are compared with the existing experimental and the quantum mechanical (QM) results by Xie et al. [ J. Chem. Phys. 2005 122 , 014301 ]. The agreement is very good at T > 600 K, although RPMD underestimates rate coefficients by a factor between 4 and 2 in the 200-500 K interval. The origin of these discrepancies lies in the large contribution from tunneling on the 3A″ PES, which is enhanced by resonances due to quasibound states in the van der Waals wells. Although tunneling is fairly well accounted for by RPMD even below the crossover temperature, the effect of resonances, a long-time effect, is not included in the methodology. At the highest temperatures studied in this work, 2000-3300 K, the RPMD rate coefficients are somewhat larger than the QM ones, but this is shown to be due to limitations in the QM calculations and the RPMD are believed to be more reliable.
Collapse
Affiliation(s)
- M Menéndez
- Departamento de Química Física I, Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain
| | - P G Jambrina
- Departamento de Química Física, Facultad de Ciencias Químicas , Universidad de Salamanca , 37008 Salamanca , Spain
| | - A Zanchet
- Departamento de Química Física I, Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain
| | - E Verdasco
- Departamento de Química Física I, Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain
| | - Y V Suleimanov
- Computation-based Science and Technology Research Center , Cyprus Institute , 20 Kavafi Strasse , Nicosia 2121 , Cyprus.,Department of Chemical Engineering , Massachusetts Institute of Technology , 77 Massachusetts Avenue , Cambridge , Massachusetts 02139 , United States
| | - F J Aoiz
- Departamento de Química Física I, Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain
| |
Collapse
|
23
|
Trenins G, Willatt MJ, Althorpe SC. Path-integral dynamics of water using curvilinear centroids. J Chem Phys 2019. [DOI: 10.1063/1.5100587] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- George Trenins
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Michael J. Willatt
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
24
|
Lu X, Wang X, Fu B, Zhang D. Theoretical Investigations of Rate Coefficients of H + H2O2 → OH + H2O on a Full-Dimensional Potential Energy Surface. J Phys Chem A 2019; 123:3969-3976. [DOI: 10.1021/acs.jpca.9b02526] [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)
- Xiaoxiao Lu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - Xingan Wang
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| | - Donghui Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
Tao X, Shushkov P, Miller TF. Simple Flux-Side Formulation of State-Resolved Thermal Reaction Rates for Ring-Polymer Surface Hopping. J Phys Chem A 2019; 123:3013-3020. [DOI: 10.1021/acs.jpca.9b00877] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Philip Shushkov
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
28
|
Meng Q, Chen J. Ring-polymer molecular dynamics study on rate coefficient of the barrierless OH + CO system at low temperature. J Chem Phys 2019; 150:044307. [PMID: 30709288 DOI: 10.1063/1.5065657] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, 710072 Xi’an, China
| | - Jun Chen
- iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Siming South Road 422, 361005 Xiamen, China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
| |
Collapse
|
29
|
Lu X, Meng Q, Wang X, Fu B, Zhang DH. Rate coefficients of the H + H2O2→ H2+ HO2reaction on an accurate fundamental invariant-neural network potential energy surface. J Chem Phys 2018; 149:174303. [PMID: 30409010 DOI: 10.1063/1.5063613] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xiaoxiao Lu
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| | - Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, Xi’an 710072, China
| | - Xingan Wang
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| |
Collapse
|
30
|
Bhowmick S, Bossion D, Scribano Y, Suleimanov YV. The low temperature D + + H 2→ HD + H + reaction rate coefficient: a ring polymer molecular dynamics and quasi-classical trajectory study. Phys Chem Chem Phys 2018; 20:26752-26763. [PMID: 30324962 DOI: 10.1039/c8cp05398g] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The reaction between D+ and H2 plays an important role in astrochemistry at low temperatures and also serves as a prototype for a simple ion-molecule reaction. Its ground X[combining tilde]1A' state has a very small thermodynamic barrier (up to 1.8 × 10-2 eV) and the reaction proceeds through the formation of an intermediate complex lying within the potential well with a depth of at least 0.2 eV, thus representing a challenge for dynamical studies. In the present work, we analyze the title reaction within the temperature range of 20-100 K by means of ring polymer molecular dynamics (RPMD) and quasi-classical trajectory (QCT) methods over the full-dimensional global potential energy surface developed by Aguado et al. [A. Aguado, O. Roncero, C. Tablero, C. Sanz and M. Paniagua, J. Chem. Phys., 2000, 112, 1240]. The computed thermal RPMD and QCT rate coefficients are found to be almost independent of temperature and fall within the range of 1.34-2.01 × 10-9 cm3 s-1. They are also in very good agreement with previous time-independent quantum mechanical and statistical quantum method calculations. Furthermore, we observe that the choice of asymptotic separation distance between the reactants can markedly alter the rate coefficient in the low temperature regime (20-50 K). Therefore it is of utmost importance to correctly assign the value of this parameter for dynamical studies, particularly at very low temperatures of astrochemical importance. We finally conclude that the experimental rate measurements for the title reaction are highly desirable in future.
Collapse
Affiliation(s)
- Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| | - Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| |
Collapse
|
31
|
Meng Q. Ring-polymer molecular dynamics study on rate coefficients of hydrogen abstraction of methane: A reduced-dimensional model. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.06.042] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
32
|
Trenins G, Althorpe SC. Mean-field Matsubara dynamics: Analysis of path-integral curvature effects in rovibrational spectra. J Chem Phys 2018; 149:014102. [DOI: 10.1063/1.5038616] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- George Trenins
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
33
|
Jung KA, Videla PE, Batista VS. Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy. J Chem Phys 2018; 148:244105. [PMID: 29960352 DOI: 10.1063/1.5036768] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.
Collapse
Affiliation(s)
- Kenneth A Jung
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Pablo E Videla
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
34
|
Thompson KM, Gao Y, Marshall P, Wang H, Zhou L, Li Y, Guo H. Experimental and theoretical studies of the reactions of ground-state sulfur atoms with hydrogen and deuterium. J Chem Phys 2018; 147:134302. [PMID: 28987091 DOI: 10.1063/1.4991418] [Citation(s) in RCA: 6] [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 gas-phase kinetics of S(3P) atoms with H2 and D2 have been studied via the laser flash photolysis-resonance fluorescence technique. S atoms were generated by pulsed photolysis of CS2 at 193 nm and monitored by time-resolved fluorescence at 181 nm. The rate coefficients for H2 (k1) and D2 (k2), respectively, are summarized as k1(600-1110 K) = 3.0 × 10-9 exp-1.317×105-2.703×107K/T8.314 T/K cm3 molecule-1 s-1 and k2(770-1110 K) = 2.2 × 10-14 (T/298 K)3.55 exp(-5420 K/T) cm3 molecule-1 s-1. Error limits are discussed in the text. The rate coefficients for formation of SH(SD) + H(D) on a newly developed triplet potential energy surface were characterized via ring polymer molecular dynamics and canonical variational transition-state theory. There is excellent agreement above about 1000 K between theory and experiment. At lower temperatures, the experimental rate coefficient is substantially larger than the results computed for the adiabatic reaction, suggesting a significant role for intersystem crossing to the singlet potential energy surface at lower temperatures.
Collapse
Affiliation(s)
| | - Yide Gao
- Department of Chemistry, University of North Texas, Denton, Texas 76203, USA
| | - Paul Marshall
- Department of Chemistry, University of North Texas, Denton, Texas 76203, USA
| | - Han Wang
- Department of Physics, International Center of Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Linsen Zhou
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
35
|
Jang S, Voth GA. Non-uniqueness of quantum transition state theory and general dividing surfaces in the path integral space. J Chem Phys 2018; 146:174106. [PMID: 28477603 DOI: 10.1063/1.4982053] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Despite the fact that quantum mechanical principles do not allow the establishment of an exact quantum analogue of the classical transition state theory (TST), the development of a quantum TST (QTST) with a proper dynamical justification, while recovering the TST in the classical limit, has been a long standing theoretical challenge in chemical physics. One of the most recent efforts of this kind was put forth by Hele and Althorpe (HA) [J. Chem. Phys. 138, 084108 (2013)], which can be specified for any cyclically invariant dividing surface defined in the space of the imaginary time path integral. The present work revisits the issue of the non-uniqueness of QTST and provides a detailed theoretical analysis of HA-QTST for a general class of such path integral dividing surfaces. While we confirm that HA-QTST reproduces the result based on the ring polymer molecular dynamics (RPMD) rate theory for dividing surfaces containing only a quadratic form of low frequency Fourier modes, we find that it produces different results for those containing higher frequency imaginary time paths which accommodate greater quantum fluctuations. This result confirms the assessment made in our previous work [Jang and Voth, J. Chem. Phys. 144, 084110 (2016)] that HA-QTST does not provide a derivation of RPMD-TST in general and points to a new ambiguity of HA-QTST with respect to its justification for general cyclically invariant dividing surfaces defined in the space of imaginary time path integrals. Our analysis also offers new insights into similar path integral based QTST approaches.
Collapse
Affiliation(s)
- Seogjoo Jang
- Department of Chemistry and Biochemistry, Queens College, City University of New York, 65-30 Kissena Boulevard, Queens, New York 11367, USA and Ph.D. Programs in Chemistry and Physics, and Initiative for the Theoretical Sciences, Graduate Center, City University of New York, 365 Fifth Avenue, New York, New York 10016, USA
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, USA
| |
Collapse
|
36
|
|
37
|
|
38
|
Robertson C, Habershon S. Harmonic-phase path-integral approximation of thermal quantum correlation functions. J Chem Phys 2018; 148:102316. [PMID: 29544325 DOI: 10.1063/1.5002189] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an approximation to the thermal symmetric form of the quantum time-correlation function in the standard position path-integral representation. By transforming to a sum-and-difference position representation and then Taylor-expanding the potential energy surface of the system to second order, the resulting expression provides a harmonic weighting function that approximately recovers the contribution of the phase to the time-correlation function. This method is readily implemented in a Monte Carlo sampling scheme and provides exact results for harmonic potentials (for both linear and non-linear operators) and near-quantitative results for anharmonic systems for low temperatures and times that are likely to be relevant to condensed phase experiments. This article focuses on one-dimensional examples to provide insights into convergence and sampling properties, and we also discuss how this approximation method may be extended to many-dimensional systems.
Collapse
Affiliation(s)
- Christopher Robertson
- Department of Chemistry and Centre for Scientific Computing, University Of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University Of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
39
|
Ásgeirsson V, Arnaldsson A, Jónsson H. Efficient evaluation of atom tunneling combined with electronic structure calculations. J Chem Phys 2018; 148:102334. [DOI: 10.1063/1.5007180] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Vilhjálmur Ásgeirsson
- Science Institute and Faculty of Physical Sciences, University of Iceland VR-III, 107 Reykjavík, Iceland
| | - Andri Arnaldsson
- Science Institute and Faculty of Physical Sciences, University of Iceland VR-III, 107 Reykjavík, Iceland
- Vatnaskil, Sídumúli 28, 108 Reykjavík, Iceland
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland VR-III, 107 Reykjavík, Iceland
- Center for Nonlinear Studies, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
40
|
Mattiat J, Richardson JO. Effects of tunnelling and asymmetry for system-bath models of electron transfer. J Chem Phys 2018; 148:102311. [PMID: 29544261 DOI: 10.1063/1.5001116] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We apply the newly derived nonadiabatic golden-rule instanton theory to asymmetric models describing electron-transfer in solution. The models go beyond the usual spin-boson description and have anharmonic free-energy surfaces with different values for the reactant and product reorganization energies. The instanton method gives an excellent description of the behaviour of the rate constant with respect to asymmetry for the whole range studied. We derive a general formula for an asymmetric version of the Marcus theory based on the classical limit of the instanton and find that this gives significant corrections to the standard Marcus theory. A scheme is given to compute this rate based only on equilibrium simulations. We also compare the rate constants obtained by the instanton method with its classical limit to study the effect of tunnelling and other quantum nuclear effects. These quantum effects can increase the rate constant by orders of magnitude.
Collapse
Affiliation(s)
- Johann Mattiat
- Laboratory of Physical Chemistry, ETH Zurich, 8093 Zurich, Switzerland
| | | |
Collapse
|
41
|
Barrera O, Bombac D, Chen Y, Daff TD, Galindo-Nava E, Gong P, Haley D, Horton R, Katzarov I, Kermode JR, Liverani C, Stopher M, Sweeney F. Understanding and mitigating hydrogen embrittlement of steels: a review of experimental, modelling and design progress from atomistic to continuum. JOURNAL OF MATERIALS SCIENCE 2018; 53:6251-6290. [PMID: 31258179 PMCID: PMC6560796 DOI: 10.1007/s10853-017-1978-5] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/11/2017] [Accepted: 12/28/2017] [Indexed: 05/21/2023]
Abstract
Hydrogen embrittlement is a complex phenomenon, involving several length- and timescales, that affects a large class of metals. It can significantly reduce the ductility and load-bearing capacity and cause cracking and catastrophic brittle failures at stresses below the yield stress of susceptible materials. Despite a large research effort in attempting to understand the mechanisms of failure and in developing potential mitigating solutions, hydrogen embrittlement mechanisms are still not completely understood. There are controversial opinions in the literature regarding the underlying mechanisms and related experimental evidence supporting each of these theories. The aim of this paper is to provide a detailed review up to the current state of the art on the effect of hydrogen on the degradation of metals, with a particular focus on steels. Here, we describe the effect of hydrogen in steels from the atomistic to the continuum scale by reporting theoretical evidence supported by quantum calculation and modern experimental characterisation methods, macroscopic effects that influence the mechanical properties of steels and established damaging mechanisms for the embrittlement of steels. Furthermore, we give an insight into current approaches and new mitigation strategies used to design new steels resistant to hydrogen embrittlement.
Collapse
Affiliation(s)
- O. Barrera
- Oxford Brookes University, Wheatley Campus, Wheatley, Oxford, OX33 1HX UK
- Department of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ UK
| | - D. Bombac
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge, CB3 0FS UK
| | - Y. Chen
- Department of Materials, University of Oxford, Parks Road, Oxford, OX1 3PH UK
| | - T. D. Daff
- Engineering Laboratory, University of Cambridge, Trumpington Street, Cambridge, CB2 1PZ UK
| | - E. Galindo-Nava
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge, CB3 0FS UK
| | - P. Gong
- Department of Materials Science and Engineering, University of Sheffield, Mappin Street, Sheffield, S1 3JD UK
| | - D. Haley
- Department of Materials, University of Oxford, Parks Road, Oxford, OX1 3PH UK
| | - R. Horton
- Department of Physics, Imperial College London, Prince Consort Road, London, SW7 2BB UK
| | - I. Katzarov
- Department of Physics, King’s College London, Strand, London, WC2R 2LS UK
| | - J. R. Kermode
- Warwick Centre for Predictive Modelling, School of Engineering, University of Warwick, Coventry, CV4 7AL UK
| | - C. Liverani
- Department of Physics, Imperial College London, Prince Consort Road, London, SW7 2BB UK
| | - M. Stopher
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge, CB3 0FS UK
| | - F. Sweeney
- Department of Materials Science and Engineering, University of Sheffield, Mappin Street, Sheffield, S1 3JD UK
| |
Collapse
|
42
|
Pierre S, Duke JR, Hele TJH, Ananth N. A mapping variable ring polymer molecular dynamics study of condensed phase proton-coupled electron transfer. J Chem Phys 2017; 147:234103. [DOI: 10.1063/1.4986517] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Sadrach Pierre
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Jessica R. Duke
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Timothy J. H. Hele
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
43
|
Buxton SJ, Habershon S. Accelerated path-integral simulations using ring-polymer interpolation. J Chem Phys 2017; 147:224107. [DOI: 10.1063/1.5006465] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Samuel J. Buxton
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
44
|
Fang W, Richardson JO, Chen J, Li XZ, Michaelides A. Simultaneous Deep Tunneling and Classical Hopping for Hydrogen Diffusion on Metals. PHYSICAL REVIEW LETTERS 2017; 119:126001. [PMID: 29341641 DOI: 10.1103/physrevlett.119.126001] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/09/2017] [Indexed: 06/07/2023]
Abstract
Hydrogen diffusion on metals exhibits rich quantum behavior, which is not yet fully understood. Using simulations, we show that many hydrogen diffusion barriers can be categorized into those with parabolic tops and those with broad tops. With parabolic-top barriers, hydrogen diffusion evolves gradually from classical hopping, to shallow tunneling, to deep tunneling as the temperature (T) decreases, and noticeable quantum effects persist at moderate T. In contrast, with broad-top barriers quantum effects become important only at low T and the classical-to-quantum transition is sharp, at which classical hopping and deep tunneling both occur. This coexistence indicates that more than one mechanism contributes to the quantum reaction rate. The conventional definition of the classical-to-quantum crossover T is invalid for the broad tops, and we give a new definition. Extending this, we propose a model to predict the transition T for broad-top diffusion, providing a general guide for theory and experiment.
Collapse
Affiliation(s)
- Wei Fang
- Thomas Young Centre, London Centre for Nanotechnology and Department of Chemistry, University College London, London WC1E 6BT, United Kingdom
| | | | - Ji Chen
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Xin-Zheng Li
- School of Physics and the Collaborative Innovation Center of Quantum Matters, Peking University, Beijing 100871, People's Republic of China
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
45
|
Zuo J, Xie C, Guo H, Xie D. Accurate Determination of Tunneling-Affected Rate Coefficients: Theory Assessing Experiment. J Phys Chem Lett 2017; 8:3392-3397. [PMID: 28685568 DOI: 10.1021/acs.jpclett.7b01296] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The thermal rate coefficients of a prototypical bimolecular reaction are determined on an accurate ab initio potential energy surface (PES) using ring polymer molecular dynamics (RPMD). It is shown that quantum effects such as tunneling and zero-point energy (ZPE) are of critical importance for the HCl + OH reaction at low temperatures, while the heavier deuterium substitution renders tunneling less facile in the DCl + OH reaction. The calculated RPMD rate coefficients are in excellent agreement with experimental data for the HCl + OH reaction in the entire temperature range of 200-1000 K, confirming the accuracy of the PES. On the other hand, the RPMD rate coefficients for the DCl + OH reaction agree with some, but not all, experimental values. The self-consistency of the theoretical results thus allows a quality assessment of the experimental data.
Collapse
Affiliation(s)
- Junxiang Zuo
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University , Nanjing 210093, China
| | - Changjian Xie
- 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
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University , Nanjing 210093, China
- Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China , Hefei, Anhui 230026, China
| |
Collapse
|
46
|
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
|
47
|
Affiliation(s)
- Timothy J. H. Hele
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, USA
| |
Collapse
|
48
|
An analysis of nonadiabatic ring-polymer molecular dynamics and its application to vibronic spectra. Chem Phys 2017. [DOI: 10.1016/j.chemphys.2016.09.036] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
49
|
Castillo JF, Suleimanov YV. A ring polymer molecular dynamics study of the OH + H2(D2) reaction. Phys Chem Chem Phys 2017; 19:29170-29176. [DOI: 10.1039/c7cp05266a] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Using ring polymer molecular dynamics we have calculated the rate coefficients for the OH + H2 reaction.
Collapse
Affiliation(s)
- J. F. Castillo
- Departamento de Química Física I
- Facultad de CC. Químicas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| | - Y. V. Suleimanov
- Computation-based Science and Technology Research Center
- Cyprus Institute
- Nicosia 2121
- Cyprus
- Department of Chemical Engineering
| |
Collapse
|
50
|
Abstract
This paper summarizes the contributions to the Faraday Discussion on reaction rate theory. The topics range from contemporary usage of transition state theory, including rare event sampling, to instantons and non-adiabatic dynamics.
Collapse
Affiliation(s)
- David Chandler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - David E Manolopoulos
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford OX1 3QZ, UK
| |
Collapse
|