1
|
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
|
2
|
Fang JH, Fan WB, Yang H, Song JN, Li YL. Rate coefficients and kinetic isotope effects of Cl+XCl→XCl+Cl (X=H, D, Mu) reactions from ring polymer molecular dynamics. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2007117] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Jun-hua Fang
- Department of Chemistry, 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
| | - 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
| | - Jia-ning Song
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Yong-le Li
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| |
Collapse
|
3
|
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
|
4
|
Espinosa-Garcia J, Garcia-Chamorro M, Corchado JC, Bhowmick S, Suleimanov YV. VTST and RPMD kinetics study of the nine-body X + C 2H 6 (X ≡ H, Cl, F) reactions based on analytical potential energy surfaces. Phys Chem Chem Phys 2020; 22:13790-13801. [PMID: 32538410 DOI: 10.1039/d0cp02238a] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate constants of nine-atom hydrogen abstraction reactions, X + C2H6 → HX + C2H5 (X ≡ H, Cl, F) with qualitatively different reaction paths, have been investigated using two kinetics approaches - variational transition state theory with multidimensional tunnelling (VTST/MT) and ring polymer molecular dynamics (RPMD) - and full dimensional analytical potential energy surfaces. For the H + C2H6 reaction, which proceeds through a noticeable barrier height of 11.62 kcal mol-1, kinetics approaches showed excellent agreement between them (with differences less than 30%) and with the experiment (with differences less than 60%) in the wide temperature range of 200-2000 K. For X = Cl and F, however, the situation is very different. The barrier height is either low or very low, 2.44 and 0.23 kcal mol-1, respectively, and the presence of van der Waals complexes in the entrance channel leads to a very flat topography and, consequently, imposes theoretical challenges. For the Cl(2P) reaction, VTST/MT underestimates the experimental rate constants (with differences less than 86%), and RPMD demonstrates better agreement (with differences less than 47%), although the temperature dependence is opposite to the experiment at low temperatures. Finally, for the F(2P) reaction, available experimental information shows discrepancies, both in the absolute values of the rate constants and also in the temperature dependence. Unfortunately, kinetics theories did not resolve this discrepancy. Different possible causes of these theory/experiment discrepancies were analyzed.
Collapse
Affiliation(s)
- Joaquin Espinosa-Garcia
- Departamento de Quimica Fisica and Instituto de Computacion Cientifica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - Moises Garcia-Chamorro
- Departamento de Quimica Fisica and Instituto de Computacion Cientifica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - Jose C Corchado
- Departamento de Quimica Fisica and Instituto de Computacion Cientifica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| |
Collapse
|
5
|
Liu Y, Li J. An accurate potential energy surface and ring polymer molecular dynamics study of the Cl + CH4→ HCl + CH3reaction. Phys Chem Chem Phys 2020; 22:344-353. [DOI: 10.1039/c9cp05693a] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate coefficients for the Cl + CH4/CD4reactions were studied on a new full-dimensional accurate potential energy surface with the spin–orbit corrections considered in the entrance channel.
Collapse
Affiliation(s)
- Yang Liu
- School of Chemistry and Chemical Engineering
- Chongqing University
- Chongqing 401331
- China
| | - Jun Li
- School of Chemistry and Chemical Engineering
- Chongqing University
- Chongqing 401331
- China
| |
Collapse
|
6
|
Novikov IS, Shapeev AV, Suleimanov YV. Ring polymer molecular dynamics and active learning of moment tensor potential for gas-phase barrierless reactions: Application to S + H2. J Chem Phys 2019; 151:224105. [DOI: 10.1063/1.5127561] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ivan S. Novikov
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Nobel St. 3, Moscow 143026, Russia
| | - Alexander V. Shapeev
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Nobel St. 3, Moscow 143026, Russia
| | - Yury V. Suleimanov
- Computation-Based Science and Technology Research Center, Cyprus Institute, 20 Kavafi Street, Nicosia 2121, Cyprus
| |
Collapse
|
7
|
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
|
8
|
Wang H, Fang J, Yang H, Song J, Li Y. Ring-polymer molecular dynamical benchmarks for X + H2 insertion reactions. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.06.012] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
9
|
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
|
10
|
Hickson KM, Suleimanov YV. Low-Temperature Experimental and Theoretical Rate Constants for the O(1D) + H2 Reaction. J Phys Chem A 2017; 121:1916-1923. [DOI: 10.1021/acs.jpca.7b00722] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kevin M. Hickson
- Institut
des Sciences Moléculaires, Université de Bordeaux, F-33400 Talence, France
- Institut
des Sciences Moléculaires, CNRS, F-33400 Talence, France
| | - Yury 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
| |
Collapse
|
11
|
Hickson KM, Suleimanov YV. An experimental and theoretical investigation of the C(1D) + D2 reaction. Phys Chem Chem Phys 2017; 19:480-486. [DOI: 10.1039/c6cp07381f] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Rate constants derived from ring polymer molecular dynamics calculations confirm the validity of this method for studying low-temperature complex-forming reactions
Collapse
Affiliation(s)
- Kevin M. Hickson
- Université de Bordeaux
- Institut des Sciences Moléculaires
- F-33400 Talence
- France
- CNRS
| | - Yury V. Suleimanov
- Computation-based Science and Technology Research Center
- Cyprus Institute
- Nicosia 2121
- Cyprus
- Department of Chemical Engineering
| |
Collapse
|
12
|
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
|
13
|
Rampino S, Suleimanov YV. Thermal Rate Coefficients for the Astrochemical Process C + CH+ → C2+ + H by Ring Polymer Molecular Dynamics. J Phys Chem A 2016; 120:9887-9893. [DOI: 10.1021/acs.jpca.6b10592] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sergio Rampino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italia
| | - Yury V. Suleimanov
- Computation-based
Science and Technology Research Center, Cyprus Institute, 20
Kavafi Street, Nicosia 2121, Cyprus
- Department
of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
14
|
Suleimanov YV, Aoiz FJ, Guo H. Chemical Reaction Rate Coefficients from Ring Polymer Molecular Dynamics: Theory and Practical Applications. J Phys Chem A 2016; 120:8488-8502. [DOI: 10.1021/acs.jpca.6b07140] [Citation(s) in RCA: 92] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yury V. Suleimanov
- Computation-based Science
and Technology Research Center, Cyprus Institute, 20 Kavafi Street, Nicosia 2121, Cyprus
- Department
of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - F. Javier Aoiz
- Departamento de Química
Física I, Facultad de CC. Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Hua Guo
- Department of Chemistry and
Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
15
|
Zuo J, Li Y, Guo H, Xie D. Rate Coefficients of the HCl + OH → Cl + H2O Reaction from Ring Polymer Molecular Dynamics. J Phys Chem A 2016; 120:3433-40. [DOI: 10.1021/acs.jpca.6b03488] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
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
| | - Yongle Li
- Department
of Physics, International Center of Quantum and Molecular Structures, 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, 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
|
16
|
Meng Q, Hickson KM, Shao K, Loison JC, Zhang DH. Theoretical and experimental investigations of rate coefficients of O(1D) + CH4 at low temperature. Phys Chem Chem Phys 2016; 18:29286-29292. [DOI: 10.1039/c6cp05517f] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The O(1D) + CH4 rates from RPMD and VUV-LIF are close to each other in the temperature range of the atmosphere.
Collapse
Affiliation(s)
- Qingyong Meng
- Department of Applied Chemistry
- Northwestern Polytechnical University
- 710072 Xi'an
- China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
| | - Kevin M. Hickson
- CNRS UMR 5255
- Institut des Sciences Moléculaires
- Université de Bordeaux
- F-33400 Talence
- France
| | - Kejie Shao
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- 116023 Dalian
- China
| | - Jean-Christophe Loison
- CNRS UMR 5255
- Institut des Sciences Moléculaires
- Université de Bordeaux
- F-33400 Talence
- France
| | - 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
- 116023 Dalian
- China
| |
Collapse
|
17
|
Hickson KM, Loison JC, Guo H, Suleimanov YV. Ring-Polymer Molecular Dynamics for the Prediction of Low-Temperature Rates: An Investigation of the C((1)D) + H2 Reaction. J Phys Chem Lett 2015; 6:4194-4199. [PMID: 26538033 DOI: 10.1021/acs.jpclett.5b02060] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Quantum mechanical calculations are important tools for predicting the rates of elementary reactions, particularly for those involving hydrogen and at low temperatures where quantum effects become increasingly important. These approaches are computationally expensive, however, particularly when applied to complex polyatomic systems or processes characterized by deep potential wells. While several approximate techniques exist, many of these have issues with reliability. The ring-polymer molecular dynamics method was recently proposed as an accurate and efficient alternative. Here, we test this technique at low temperatures (300-50 K) by analyzing the behavior of the barrierless C((1)D) + H2 reaction over the two lowest singlet potential energy surfaces. To validate the theory, rate coefficients were measured using a supersonic flow reactor down to 50 K. The experimental and theoretical rates are in excellent agreement, supporting the future application of this method for determining the kinetics and dynamics of a wide range of low-temperature reactions.
Collapse
Affiliation(s)
- Kevin M Hickson
- Institut des Sciences Moléculaires, Université de Bordeaux , F-33400 Talence, France
- CNRS, Institut des Sciences Moléculaires , F-33400 Talence, France
| | - Jean-Christophe Loison
- Institut des Sciences Moléculaires, Université de Bordeaux , F-33400 Talence, France
- CNRS, Institut des Sciences Moléculaires , F-33400 Talence, France
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center, Cyprus Institute , 20 Kavafi Str., Nicosia 2121, Cyprus
- Department of Chemical Engineering, Massachusetts Institute of Technology , 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
18
|
Arseneau DJ, Fleming DG, Li Y, Li J, Suleimanov YV, Guo H. Rate Coefficient for the 4Heμ + CH4 Reaction at 500 K: Comparison between Theory and Experiment. J Phys Chem B 2015; 120:1641-8. [DOI: 10.1021/acs.jpcb.5b08368] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Donald J. Arseneau
- TRIUMF
and Department of Chemistry, University of British Columbia, Vancouver, BC V6T 2Z1, Canada
| | | | | | | | | | | |
Collapse
|
19
|
Hele TJH, Suleimanov YV. Should thermostatted ring polymer molecular dynamics be used to calculate thermal reaction rates? J Chem Phys 2015; 143:074107. [DOI: 10.1063/1.4928599] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Timothy J. H. Hele
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Yury V. Suleimanov
- Computation-based Science and Technology Research Center, Cyprus Institute, 20 Kavafi St., Nicosia 2121, Cyprus
- Department of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave., Cambridge, Massachusetts 02139, USA
| |
Collapse
|
20
|
Suleimanov YV, Espinosa-Garcia J. Recrossing and Tunneling in the Kinetics Study of the OH + CH4 → H2O + CH3 Reaction. J Phys Chem B 2015; 120:1418-28. [DOI: 10.1021/acs.jpcb.5b02103] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yury V. Suleimanov
- Department
of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts
Avenue, Cambridge, Massachusetts 02139, United States
- Computation-based Science and Technology
Research Center, Cyprus Institute, 20 Kavafi Street, Nicosia 2121, Cyprus
| | - J. Espinosa-Garcia
- Departamento de Química
Física, Universidad de Extremadura, 06071 Badajoz, Spain
| |
Collapse
|
21
|
Menéndez M, Castillo JF, Martínez-Haya B, Aoiz FJ. The Cl + O3 reaction: a detailed QCT simulation of molecular beam experiments. Phys Chem Chem Phys 2015; 17:25471-82. [DOI: 10.1039/c5cp04323a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
QCT calculations have been carried out to determine angle–velocity differential cross-sections to simulate the results of molecular beam experiments.
Collapse
Affiliation(s)
- M. Menéndez
- Departamento de Qumíca Física I
- Facultad de CC. Qumícas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| | - J. F. Castillo
- Departamento de Qumíca Física I
- Facultad de CC. Qumícas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| | - B. Martínez-Haya
- Department of Physical
- Chemical and Natural Systems
- Universidad Pablo de Olavide
- ES-41013 Seville
- Spain
| | - F. J. Aoiz
- Departamento de Qumíca Física I
- Facultad de CC. Qumícas
- Universidad Complutense de Madrid
- 28040 Madrid
- Spain
| |
Collapse
|
22
|
Suleimanov YV, Kong WJ, Guo H, Green WH. Ring-polymer molecular dynamics: Rate coefficient calculations for energetically symmetric (near thermoneutral) insertion reactions (X + H2) → HX + H(X = C(1D), S(1D)). J Chem Phys 2014; 141:244103. [DOI: 10.1063/1.4904080] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yury V. Suleimanov
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
- Computation-based Science and Technology Research Center, Cyprus Institute, 20 Kavafi Street, Nicosia 2121, Cyprus
| | - Wendi J. Kong
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Hua Guo
- Department of Chemical and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - William H. Green
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
23
|
Pérez de Tudela R, Suleimanov YV, Richardson JO, Sáez Rábanos V, Green WH, Aoiz FJ. Stress Test for Quantum Dynamics Approximations: Deep Tunneling in the Muonium Exchange Reaction D + HMu → DMu + H. J Phys Chem Lett 2014; 5:4219-4224. [PMID: 26278957 DOI: 10.1021/jz502216g] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Quantum effects play a crucial role in chemical reactions involving light atoms at low temperatures, especially when a light particle is exchanged between two heavier partners. Different theoretical methodologies have been developed in the last decades attempting to describe zero-point energy and tunneling effects without abandoning a classical or semiclassical framework. In this work, we have chosen the D + HMu → DMu + H reaction as a stress test system for three well-established methods: two representative versions of transition state theory (TST), canonical variational theory and semiclassical instanton, and ring polymer molecular dynamics (RPMD). These calculations will be compared with accurate quantum mechanical results. Despite its apparent simplicity, the exchange of the extremely light muonium atom (0.114 u) becomes a most challenging reaction for conventional methods. The main result of this work is that RPMD provides an overall better performance than TST-based methods for such a demanding reaction. RPMD might well turn out to be a useful tool beyond TST applicability.
Collapse
Affiliation(s)
- Ricardo Pérez de Tudela
- †Departamento de Quı́mica Fı́sica I, Facultad de CC. Quı́micas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Yury V Suleimanov
- ‡Computation-based Science and Technology Research Center, Cyprus Institute, 20 Kavafi Str., Nicosia 2121, Cyprus
- ∥Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Jeremy O Richardson
- ¶Institut für Theoretische Physik und Interdisziplinäres Zentrum für Molekulare Materialien, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), 91058 Erlangen, Germany
| | - Vicente Sáez Rábanos
- §Departamento de Quı́mica y Bioquı́mica, ETS Ingenieros de Montes, Universidad Politécnica, 28040 Madrid, Spain
| | - William H Green
- ∥Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - F J Aoiz
- †Departamento de Quı́mica Fı́sica I, Facultad de CC. Quı́micas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
24
|
Gonzalez-Lavado E, Corchado JC, Suleimanov YV, Green WH, Espinosa-Garcia J. Theoretical Kinetics Study of the O(3P) + CH4/CD4 Hydrogen Abstraction Reaction: The Role of Anharmonicity, Recrossing Effects, and Quantum Mechanical Tunneling. J Phys Chem A 2014; 118:3243-52. [DOI: 10.1021/jp5028965] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Jose C. Corchado
- Departamento
de Química Física, Universidad de Extremadura, 06071 Badajoz, Spain
| | - Yury V. Suleimanov
- Department of Chemical
Engineering, Massachusetts Institute of Technology, 77 Massachusetts
Avenue, Cambridge, Massachusetts 02139, United States
- Department of Mechanical and Aerospace Engineering, Combustion Energy
Frontier Research Center, Princeton University, Olden Street, Princeton, New Jersey 08544, United States
| | - William H. Green
- Department of Chemical
Engineering, Massachusetts Institute of Technology, 77 Massachusetts
Avenue, Cambridge, Massachusetts 02139, United States
| | | |
Collapse
|
25
|
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: 50] [Impact Index Per Article: 5.0] [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
|
26
|
Li Y, Suleimanov YV, Guo H. Ring-Polymer Molecular Dynamics Rate Coefficient Calculations for Insertion Reactions: X + H2 → HX + H (X = N, O). J Phys Chem Lett 2014; 5:700-705. [PMID: 26270840 DOI: 10.1021/jz500062q] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The thermal rate constants of two prototypical insertion-type reactions, namely, N/O + H2 → NH/OH + H, are investigated with ring polymer molecular dynamics (RPMD) on full-dimensional potential energy surfaces using recently developed RPMDrate code. It is shown that the unique ability of the RPMD approach among the existing theoretical methods to capture the quantum effects, e.g., tunneling and zero-point energy, as well as recrossing dynamics quantum mechanically with ring-polymer trajectories leads to excellent agreement with rigorous quantum dynamics calculations. The present result is encouraging for future applications of the RPMD method and the RPMDrate code to complex-forming chemical reactions involving polyatomic reactants.
Collapse
Affiliation(s)
- Yongle Li
- †Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yury V Suleimanov
- ‡Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- §Department of Mechanical and Aerospace Engineering, Combustion Energy Frontier Research Center, Princeton University, Olden Street, Princeton, New Jersey 08544, United States
| | - Hua Guo
- †Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|