1
|
Menéndez M, Veselinova A, Zanchet A, Jambrina PG, Aoiz FJ. Rate coefficients for the O + H 2 and O + D 2 reactions: how well ring polymer molecular dynamics accounts for tunelling. Phys Chem Chem Phys 2024; 26:20947-20961. [PMID: 39046374 DOI: 10.1039/d4cp01711k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2024]
Abstract
We present here extensive calculations of the O(3P) + H2 and O(3P) + D2 reaction dynamics spanning the temperature range from 200 K to 2500 K. The calculations have been carried out using fully converged time-independent quantum mechanics (TI QM), quasiclassical trajectories (QCT) and ring polymer molecular dynamics (RPMD) on the two lowest lying adiabatic potential energy surfaces (PESs), 13A' and 13A'', calculated by Zanchet et al. [J. Chem. Phys., 2019, 151, 094307]. TI QM rate coefficients were determined using the cumulative reaction probability formalism on each PES including all of the total angular momenta and the Coriolis coupling and can be considered to be essentially exact within the Born-Oppenheimer approximation. The agreement between the rate coefficients calculated by using QM and RPMD is excellent for the reaction with D2 in almost the whole temperature range. For the reaction with H2, although the agreement is very good above 500 K, the deviations are significant at lower temperatures. In contrast, the QCT calculations largely underestimate the rate coefficients for the two isotopic variants due to their inability to account for tunelling. The differences found in the disagreements between RPMD and QM rate coefficients for the reactions for both the isotopologues are indicative of the ability of the RPMD method to accurately describe systems where tunelling plays a relevant role. Considering that both reactions are dominated by tunelling below 500 K, the present results show that RPMD is a very powerful tool for determining rate coefficients. The present QM rate coefficients calculated on adiabatic PESs slightly underestimate the best global fits of the experimental measurements, which we attribute to the intersystem crossing with the singlet 11A' PES.
Collapse
Affiliation(s)
- Marta Menéndez
- Departamento de Química Física, Unidad Asociada CSIC, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| | - Anzhela Veselinova
- Departamento de Química Física, Universidad de Salamanca, 37008 Salamanca, Spain
| | - Alexandre Zanchet
- Instituto de Física Fundamental, CSIC, C/Serrano 121-123, 28006 Madrid, Spain
| | - Pablo G Jambrina
- Departamento de Química Física, Universidad de Salamanca, 37008 Salamanca, Spain
| | - F Javier Aoiz
- Departamento de Química Física, Unidad Asociada CSIC, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| |
Collapse
|
2
|
Espinosa-Garcia J, Bhowmick S. Kinetic study of the CN + C 2H 6 hydrogen abstraction reaction based on an analytical potential energy surface. Phys Chem Chem Phys 2024; 26:8344-8355. [PMID: 38391269 DOI: 10.1039/d3cp05930h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2024]
Abstract
The temperature dependence of the thermal rate constants and kinetic isotope effects (KIE) of the CN + C2H6 gas-phase hydrogen abstraction reaction was theoretically determined within the 25-1000 K temperature range, i.e., from very low- to high-temperature regimes. Based on a recently developed full-dimensional analytical potential energy surface fitted to highly accurate explicitly correlated ab initio calculations, three different kinetic theories were used: canonical variational transition state theory (CVT), quasiclassical trajectory theory (QCT), and ring polymer molecular dynamics (RPMD) method for the computation of rate constants. We found that the thermal rate constants obtained with the three theories show a V-shaped temperature dependence, with a pronounced minimum near 200 K, qualitatively reproducing the experimental measurements. Among the three methods used in this work, the QCT and RPMD methods have the best agreement with the experiment at low and high temperatures, respectively, while the CVT model shows the largest discrepancies. The significant increase in the rate constant at very low temperatures in this very exothermic and practically barrierless reaction could be attributed to the large value of the impact parameter, possibly ruling out the role of the tunneling effect and the intermediate complexes in the entrance channel. The theoretical H/D KIE depicted a "normal" behaviour, i.e., values greater than unity, emulating the experimental measurements and improving previous theoretical results. Finally, the discrepancies between theory and experiments were analysed as a function of several factors, such as limitations of the kinetics theories and the potential energy surface, as well as the uncertainties in the experimental measurements.
Collapse
Affiliation(s)
- Joaquin Espinosa-Garcia
- Departamento de Química Física and Instituto de Computación Científica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - Somnath Bhowmick
- Climate and Atmosphere Research Centre, The Cyprus Institute, Nicosia 2121, Cyprus.
| |
Collapse
|
3
|
Li C, Li Y, Jiang B. First-principles surface reaction rates by ring polymer molecular dynamics and neural network potential: role of anharmonicity and lattice motion. Chem Sci 2023; 14:5087-5098. [PMID: 37206404 PMCID: PMC10189860 DOI: 10.1039/d2sc06559b] [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: 11/29/2022] [Accepted: 04/05/2023] [Indexed: 08/04/2023] Open
Abstract
Elementary gas-surface processes are essential steps in heterogeneous catalysis. A predictive understanding of catalytic mechanisms remains challenging due largely to difficulties in accurately characterizing the kinetics of such steps. Experimentally, thermal rates for elementary surface reactions can now be measured using a novel velocity imaging technique, providing a stringent testing ground for ab initio rate theories. Here, we propose to combine ring polymer molecular dynamics (RPMD) rate theory with state-of-the-art first-principles-determined neural network potential to calculate surface reaction rates. Taking NO desorption from Pd(111) as an example, we show that the harmonic approximation and the neglect of lattice motion in the commonly-used transition state theory overestimates and underestimates the entropy change during the desorption process, respectively, leading to opposite errors in rate coefficient predictions and artificial error cancellations. Including anharmonicity and lattice motion, our results reveal a generally neglected surface entropy change due to significant local structural change during desorption and obtain the right answer for the right reasons. Although quantum effects are found to be less important in this system, the proposed approach establishes a more reliable theoretical benchmark for accurately predicting the kinetics of elementary gas-surface processes.
Collapse
Affiliation(s)
- Chen Li
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China Hefei Anhui 230026 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
| | - Bin Jiang
- Key Laboratory of Precision and Intelligent Chemistry, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China Hefei Anhui 230026 China
| |
Collapse
|
4
|
Angiolari F, Huppert S, Spezia R. Quantum versus classical unimolecular fragmentation rate constants and activation energies at finite temperature from direct dynamics simulations. Phys Chem Chem Phys 2022; 24:29357-29370. [PMID: 36448557 DOI: 10.1039/d2cp03809a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
In the present work, we investigate how nuclear quantum effects modify the temperature dependent rate constants and, consequently, the activation energies in unimolecular reactions. In the reactions under study, nuclear quantum effects mainly stem from the presence of a large zero point energy. Thus, we investigate the behavior of methods compatible with direct dynamics simulations, the quantum thermal bath (QTB) and ring polymer molecular dynamics (RPMD). To this end, we first compare them with quantum reaction theory for a model Morse potential before extending this comparison to molecular models. Our results show that, in particular in the temperature range comparable with or lower than the zero point energy of the system, the RPMD method is able to correctly capture nuclear quantum effects on rate constants and activation energies. On the other hand, although the QTB provides a good description of equilibrium properties including zero-point energy effects, it largely overestimates the rate constants. The origin of the different behaviours is in the different distance distributions provided by the two methods and in particular how they differently describe the tails of such distributions. The comparison with transition state theory shows that RPMD can be used to study fragmentation of complex systems for which it may be difficult to determine the multiple reaction pathways and associated transition states.
Collapse
Affiliation(s)
- Federica Angiolari
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 4 Place Jussieu, 75005 Paris, France.
| | - Simon Huppert
- Sorbonne Université, Institut de Nanosciences de Paris, UMR 7588 CNRS, 4 Place Jussieu, 75005 Paris, France
| | - Riccardo Spezia
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 4 Place Jussieu, 75005 Paris, France.
| |
Collapse
|
5
|
Gui X, Fan W, Sun J, Li Y. New Stable and Fast Ring-Polymer Molecular Dynamics for Calculating Bimolecular Rate Coefficients with an Example of OH + CH 4. J Chem Theory Comput 2022; 18:5203-5212. [PMID: 35983956 DOI: 10.1021/acs.jctc.2c00522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The accurate and efficient calculation of the rate coefficients of chemical reactions is a key issue in the research of chemical dynamics. In this work, by applying the dimension-free ultrastable Cayley propagator, the thermal rate coefficients of a prototypic high dimensional chemical reaction OH + CH4 → H2O + CH3 in the temperature range of 200 to 1500 K are investigated with ring polymer molecular dynamics (RPMD) on a highly accurate full-dimensional potential energy surface. Kinetic isotope effects (KIEs) for three isotopologues of the title reaction are also studied. The results demonstrate excellent agreement with experimental data, even in the deep tunneling region. Especially, the Cayley propagator shows a high calculation efficiency with little loss of accuracy. The present results confirmed the applicability of the RPMD method, particularly the speed-up using a Cayley propagator, in theoretical calculations of bimolecular reaction rates.
Collapse
Affiliation(s)
- Xiongfei Gui
- Department of Physics, International Center of Quantum and Molecular Structures, and Shanghai Frontiers Science Center of Quantum and Superconducting Matter States, Shanghai 200444, China
| | - Wenbin Fan
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Jiace Sun
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures, and Shanghai Frontiers Science Center of Quantum and Superconducting Matter States, Shanghai 200444, China
| |
Collapse
|
6
|
Peng Y, Zhang H. Mechanism and Kinetics of Methane Combustion. Part II: Potential Energy Surface for Hydrogen-Abstraction Reaction of CH 4 + O( 3P). J Phys Chem A 2022; 126:1946-1959. [PMID: 35298157 DOI: 10.1021/acs.jpca.1c10860] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Methane combustion plays an important role in various fields such as combustion chemistry and atmospheric chemistry of the stratosphere. Highly accurate study of its initial reaction remains a key challenge. Here, through extensive studies with a state-of-the-art ab initio and neural network method, we present a potential energy surface of the O(3P) + CH4 → OH + CH3 reaction on the ground state 13A and the first excited state 23A. In this work, the energies of 10 167 points covering all important regions are obtained with state-averaged complete active space self-consistent field calculations and then fitted using the Levenberg-Marquardt algorithm with a root-mean-square error of 0.391 and 0.442 kcal/mol for the 13A and 23A states, respectively. This study explores the characteristics of the radical van der Waals (VdW) complex and reveals a detailed mechanism of the methane combustion initial reaction. Within the scope of this mechanism, this surface gives a fairly accurate description of the regions around the saddle point, conical intersection, and vdW wells in the entrance for efficient computational simulations. As a theoretical study on a prototypical polyatomic reaction, it is hopeful that this work will modify our understanding of the primary process in hydrocarbon combustion.
Collapse
Affiliation(s)
- Ya Peng
- Department of Engineering Physics, Tsinghua University, Beijing 100084, P.R. China
| | - Hui Zhang
- Department of Engineering Physics, Tsinghua University, Beijing 100084, P.R. China
| |
Collapse
|
7
|
Marjollet A, Inhester L, Welsch R. Initial state-selected scattering for the reactions H + CH4/CHD3 and F + CHD3 employing ring polymer molecular dynamics. J Chem Phys 2022; 156:044101. [DOI: 10.1063/5.0076216] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- A. Marjollet
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
- Department of Physics, Universität Hamburg, Notkestr. 9-11, 22607 Hamburg, Germany
| | - L. Inhester
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - R. Welsch
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
8
|
Li C, Liu Q, Zhang L, Li Y, Jiang B. Ring polymer molecular dynamics in gas-surface reactions: tests on initial sampling and potential energy landscape. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1941367] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Chen Li
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Qinghua Liu
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Liang Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai, People’s Republic of China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| |
Collapse
|
9
|
Cao J, Wu Y, Ma H, Shen Z, Bian W. Dynamics and kinetics of the Si( 1D) + H 2/D 2 reactions on a new global ab initio potential energy surface. Phys Chem Chem Phys 2021; 23:6141-6153. [PMID: 33684184 DOI: 10.1039/d0cp05540a] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recent studies on the exothermic complex-forming reactions have improved our understanding of complex-forming reactions greatly, however, so far a similar level of study on endothermic ones has been rather limited. In this work, the endothermic complex-forming reaction Si(1D) + H2 → SiH + H and its deuterated isotopic variant are investigated by quantum dynamics (QD) and ring polymer molecular dynamics (RPMD) calculations on a new global ab initio potential energy surface (PES) for the ground electronic state, which is constructed based on 8996 symmetry unique points computed at the icMRCI+Q/aug-cc-pVQZ level. The PES reproduces our ab initio data very well in the dynamically important regions, on which the ro-vibrational energy levels of SiH2 are calculated and general good agreement with experiment is obtained. The integral cross sections and product angular and state distributions are computed in a wide range of collision energies, and interesting dynamics behaviors are revealed. The rate coefficients are also investigated, which display an exponential rise from 2.09 × 10-20 to 6.00 × 10-12 cm3 s-1 for the Si(1D) + H2 reaction as the temperature increases from 300 to 1500 K, in contrast to the nearly temperature-independent behavior of exothermic complex-forming reactions. In addition, the applicability of the RPMD approach is demonstrated, and the kinetic isotope effect is investigated, the ratio of which decreases from 7.89 (300 K) to 1.70 (1500 K). The effects of tunneling and initial rotational excitation are also discussed.
Collapse
Affiliation(s)
- Jianwei Cao
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.
| | - Yanan Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China. and School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Haitao Ma
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.
| | - Zhitao Shen
- Henan Key Laboratory of Photovoltaic Materials, Henan University, Kaifeng 475004, China
| | - Wensheng Bian
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China. and School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
10
|
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
|
11
|
Song Q, Zhang Q, Meng Q. Revisiting the Gaussian process regression for fitting high-dimensional potential energy surface and its application to the OH + HO2→O2+ H2O reaction. J Chem Phys 2020; 152:134309. [PMID: 32268765 DOI: 10.1063/1.5143544] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qingfei Song
- School of Chemistry and Chemical Engineering, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Key Laboratory of Special Functional and Smart Polymer Materials of Ministry of Industry and Information Technology, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| | - Qiuyu Zhang
- School of Chemistry and Chemical Engineering, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Key Laboratory of Special Functional and Smart Polymer Materials of Ministry of Industry and Information Technology, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| | - Qingyong Meng
- School of Chemistry and Chemical Engineering, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Key Laboratory of Special Functional and Smart Polymer Materials of Ministry of Industry and Information Technology, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| |
Collapse
|
12
|
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
|
13
|
Lawrence JE, Manolopoulos DE. Path integral methods for reaction rates in complex systems. Faraday Discuss 2020; 221:9-29. [DOI: 10.1039/c9fd00084d] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We shall use this introduction to the Faraday Discussion on quantum effects in complex systems to review the recent progress that has been made in using imaginary time path integral methods to calculate chemical reaction rates.
Collapse
Affiliation(s)
- Joseph E. Lawrence
- Department of Chemistry
- University of Oxford
- Physical and Theoretical Chemistry Laboratory
- Oxford OX1 3QZ
- UK
| | - David E. Manolopoulos
- Department of Chemistry
- University of Oxford
- Physical and Theoretical Chemistry Laboratory
- Oxford OX1 3QZ
- UK
| |
Collapse
|
14
|
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
|
15
|
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
|
16
|
Houston PL, Nandi A, Bowman JM. A Machine Learning Approach for Prediction of Rate Constants. J Phys Chem Lett 2019; 10:5250-5258. [PMID: 31423788 DOI: 10.1021/acs.jpclett.9b01810] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We report a machine learning approach to train and predict bimolecular thermal rate constants over a large temperature range. The approach uses Gaussian process (GP) regression to evaluate the difference between accurate quantum results and Eckart-corrected conventional transition state theory, mostly for collinear reactions. Training is done on a database of rate constants for 13 reaction/potential surface combinations, and testing is performed on a set of 39 reaction/potential surface combinations. Averaged over all test reactions, the GP method is within 80% of the accurate answer, whereas transition state theory (TST) is only within 330% and Eckart-corrected TST (ECK) is within 110%. In the tunneling region, GP is generally (with a few exceptions) more accurate and sometimes much more accurate. In the high-temperature recrossing region, GP is significantly more accurate than either TST or ECK, neither of which addresses the possibility of recrossing. The GP predictions for the 3D reactions O(3P) + H2, OH + H2, O(3P) + CH4, and H + CH4, for which accurate quantum results are available, provide further encouragement to the machine learning approach.
Collapse
Affiliation(s)
- Paul L Houston
- Department of Chemistry & Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientic Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientic Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
17
|
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]
|
18
|
Zhu Y, Lu Y, Song H. Thermal rate coefficients and kinetic isotope effects of the reaction HO + H2O → H2O + OH. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2495-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
19
|
Li J. Ring-polymer molecular dynamics studies of thermal rate coefficients for reaction F + H2O → HF + OH. CHINESE J CHEM PHYS 2019. [DOI: 10.1063/1674-0068/cjcp1808186] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Jun Li
- School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 400044, China
| |
Collapse
|
20
|
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
|
21
|
Liu R, Song HW, Yang MH. Understanding rotational mode specificity in the O(3P)+CHD3→ OH+CD3 reaction by simple reactant alignment pictures. CHINESE J CHEM PHYS 2019. [DOI: 10.1063/1674-0068/cjcp1810238] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Rui Liu
- Key Laboratory of Magnetic Resonance in Biological Systems, National Center for Magnetic Resonance in Wuhan, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
| | - Hong-wei Song
- Key Laboratory of Magnetic Resonance in Biological Systems, National Center for Magnetic Resonance in Wuhan, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
| | - Ming-hui Yang
- Key Laboratory of Magnetic Resonance in Biological Systems, National Center for Magnetic Resonance in Wuhan, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
| |
Collapse
|
22
|
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
|
23
|
Espinosa-Garcia J, Corchado JC. The hydrogen abstraction reaction H + C2H6 → H2(v,j) + C2H5. Part II. Theoretical kinetics and dynamics study. Phys Chem Chem Phys 2019; 21:13356-13367. [DOI: 10.1039/c9cp00701f] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Two important issues motivated the present study: the role of the tunnelling contribution at low temperatures and the role of the alkyl fragment in the dynamics.
Collapse
Affiliation(s)
- Joaquin Espinosa-Garcia
- Área de Química Física and Instituto de Computacion Cientifica Avanzada
- Universidad de Extremadura
- 06071 Badajoz
- Spain
| | - Jose C. Corchado
- Área de Química Física and Instituto de Computacion Cientifica Avanzada
- Universidad de Extremadura
- 06071 Badajoz
- Spain
| |
Collapse
|
24
|
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
|
25
|
Meng Q. Ring-Polymer Molecular Dynamics with Coarse-Grained Treatment of the Rate Coefficients of Chlorine Atom Reactions with Methane, Ethane, and Propane. J Phys Chem A 2018; 122:8320-8325. [PMID: 30281311 DOI: 10.1021/acs.jpca.8b08052] [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]
Affiliation(s)
- Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, 710072 Xi’an, China
- State Key Laboratory
of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
| |
Collapse
|
26
|
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]
|
27
|
de Carvalho EFV, Roberto-Neto O. Effects of multidimensional tunneling in the kinetics of hydrogen abstraction reactions of O ( 3 P) with CH 3 OCHO. J Comput Chem 2018; 39:1424-1432. [PMID: 29607518 DOI: 10.1002/jcc.25211] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2018] [Revised: 02/20/2018] [Accepted: 02/25/2018] [Indexed: 11/11/2022]
Abstract
Quantum tunneling paths are important in reactions when there is a significant component of hydrogenic motion along the potential energy surface. In this study, variational transition state with multidimensional tunneling corrections are employed in the calculations of the thermal rate constants for hydrogen abstraction from the cis-CH3 OCHO by O (3 P) giving CH3 OCO + OH (R1) and CH2 OCHO + OH (R2). The structures and electronic energies are computed with the M06-2X method. Benchmark calculations with the CBSD-T approach give an enthalpy of reaction at 0 K for R1 (-2.8 kcal/mol) and R2 (-2.5 kcal/mol) which are in good agreement with the experiment, i.e. -2.61 and -1.81 kcal/mol. At the low and intermediate values of temperatures, small- and large-curvature tunneling dominate the kinetics of R1, which is the dominant path over the range of temperature from 250 to 1200 K. This study shows the importance of multidimensional tunneling corrections for both R1 and R2, for which the total rate constant at 298 K calculated with the CVT/μOMT method is 8.2 × 10-15 cm3 molecule-1 s-1 which agrees well with experiment value of 9.3 × 10-15 cm3 molecule-1 s-1 (Mori, Bull. Inst. Chem. Res. 1981, 59, 116). © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- E F V de Carvalho
- Departamento de Física, Universidade Federal do Maranhão, São Luís, Maranhão, 65085-580, Brazil
| | - O Roberto-Neto
- Divisão de Aerotermodinâmica e Hipersônica, Instituto de Estudos Avançados, São José dos Campos, São Paulo, 12228-001, Brazil
| |
Collapse
|
28
|
Novikov IS, Suleimanov YV, Shapeev AV. Automated calculation of thermal rate coefficients using ring polymer molecular dynamics and machine-learning interatomic potentials with active learning. Phys Chem Chem Phys 2018; 20:29503-29512. [DOI: 10.1039/c8cp06037a] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We propose a methodology for the fully automated calculation of thermal rate coefficients of gas phase chemical reactions, which is based on combining ring polymer molecular dynamics (RPMD) and machine-learning interatomic potentials actively learning on-the-fly.
Collapse
Affiliation(s)
- I. S. Novikov
- Skolkovo Institute of Science and Technology
- Skolkovo Innovation Center
- Moscow 143026
- Russia
| | - Y. V. Suleimanov
- Computation-based Science and Technology Research Center
- Cyprus Institute
- Nicosia 2121
- Cyprus
- Department of Chemical Engineering
| | - A. V. Shapeev
- Skolkovo Institute of Science and Technology
- Skolkovo Innovation Center
- Moscow 143026
- Russia
| |
Collapse
|
29
|
Peng Y, Jiang Z, Chen J. Mechanism and Kinetics of Methane Combustion, Part I: Thermal Rate Constants for Hydrogen-Abstraction Reaction of CH4 + O(3P). J Phys Chem A 2017; 121:2209-2220. [DOI: 10.1021/acs.jpca.6b12125] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ya Peng
- School of Civil and Resource
Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Zhong’an Jiang
- School of Civil and Resource
Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Jushi Chen
- School of Civil and Resource
Engineering, University of Science and Technology Beijing, Beijing 100083, China
| |
Collapse
|
30
|
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
|
31
|
Meng Q, Chen J. Ring-polymer molecular dynamics studies on the rate coefficient of the abstraction channel of hydrogen plus ethane, propane, and dimethyl ether. J Chem Phys 2017; 146:024108. [PMID: 28088140 DOI: 10.1063/1.4973831] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [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
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
| | - Jun Chen
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
- iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Siming South Road 422, 361005 Xiamen, China
| |
Collapse
|
32
|
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
|
33
|
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
|
34
|
Bai M, Lu D, Li Y, Li J. Ring-polymer molecular dynamical calculations for the F + HCl → HF + Cl reaction on the ground 1 2A' potential energy surface. Phys Chem Chem Phys 2016; 18:32031-32041. [PMID: 27407007 DOI: 10.1039/c6cp03306g] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
The reaction kinetics of the heavy-light-heavy abstraction reaction F + HCl → HF + Cl on the ground electronic state potential energy surface (PES) is investigated theoretically by a recently developed ring polymer molecular dynamics (RPMD) approach. First, a new PES is developed by the permutation invariant polynomial neural network (PIP-NN) approach based on 30 620 points sampled over a large configuration space from the latest and most accurate Deskevich-Hayes-Takahashi-Skodje-Nesbitt (DHTSN) PES (J. Chem. Phys., 2006, 124, 224303). Excellent fitting performance was achieved with only 521 parameters. The PIP-NN PES is 11 times faster than the DHTSN PES. Besides, the first analytical derivatives with respect to the coordinates of the atoms have been obtained for the PIP-NN PES. The RPMD rate coefficients on the PIP-NN PES are calculated and compared with available theoretical and experimental results. It is found that the experimental rate coefficients are significantly larger than the theoretical results on the DHTSN PES, due to its overestimated reaction barrier. We conclude that a reliable PES for this important heavy-light-heavy reaction is highly desirable.
Collapse
Affiliation(s)
- Mengna Bai
- School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, China.
| | - Dandan Lu
- School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, China.
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structure, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China.
| | - Jun Li
- School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, China.
| |
Collapse
|
35
|
Welsch R, Song K, Shi Q, Althorpe SC, Miller TF. Non-equilibrium dynamics from RPMD and CMD. J Chem Phys 2016; 145:204118. [DOI: 10.1063/1.4967958] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Affiliation(s)
- Ralph Welsch
- Division of Chemistry and Chemical Engineering, California Institute of Technology,1200 E. California Blvd., Pasadena, California 91125, USA
| | - Kai Song
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology,1200 E. California Blvd., Pasadena, California 91125, USA
| |
Collapse
|
36
|
Zhao H, Wang W, Zhao Y. Thermal Rate Constants for the O( 3P) + CH 4 → OH + CH 3 Reaction: The Effects of Quantum Tunneling and Potential Energy Barrier Shape. J Phys Chem A 2016; 120:7589-7597. [PMID: 27640428 DOI: 10.1021/acs.jpca.6b07029] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The rate constants and kinetic isotope effects for the O(3P) + CH4 reaction have been investigated with the quantum instanton method in full dimensionality. The calculated rate constants are in good agreement with the experimental values above 400 K, below which the measured values are scattered. Compared to other theoretical approaches, the quantum instanton method predicts the largest quantum tunneling effect, so it gives the largest rate constants at low temperatures. The calculated kinetic isotope effects are always much larger than 1 and increase with decreasing temperature, due to the zero-point energy and quantum tunneling. Our calculations on different potential energy surfaces demonstrate that the potential energy barrier shape dominates the magnitude of quantum tunneling and has a great effect on the kinetic isotope effect.
Collapse
Affiliation(s)
- Huali Zhao
- College of Science, Northwest A&F University , Yangling, 712100 Shaanxi Province, P. R. China
| | - Wenji Wang
- College of Science, Northwest A&F University , Yangling, 712100 Shaanxi Province, P. R. China
| | - Yi Zhao
- State Key Laboratory for Physical Chemistry of Solid Surfaces and Fujian Provincial Key Lab of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University , Xiamen 361005, P. R. China
| |
Collapse
|
37
|
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
|
38
|
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
|
39
|
Hele TJH, Althorpe SC. An alternative derivation of ring-polymer molecular dynamics transition-state theory. J Chem Phys 2016; 144:174107. [DOI: 10.1063/1.4947589] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [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
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
40
|
Meng Q, Chen J, Zhang DH. Ring polymer molecular dynamics fast computation of rate coefficients on accurate potential energy surfaces in local configuration space: Application to the abstraction of hydrogen from methane. J Chem Phys 2016; 144:154312. [DOI: 10.1063/1.4947097] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qingyong Meng
- 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, 116023 Dalian, China
| | - Jun Chen
- 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, 116023 Dalian, 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, 116023 Dalian, China
| |
Collapse
|
41
|
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
|
42
|
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
|
43
|
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
|
44
|
Hele TJH, Willatt MJ, Muolo A, Althorpe SC. Boltzmann-conserving classical dynamics in quantum time-correlation functions: “Matsubara dynamics”. J Chem Phys 2015; 142:134103. [DOI: 10.1063/1.4916311] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.9] [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
| | - Michael J. Willatt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Andrea Muolo
- 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
|
45
|
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
|
46
|
Zhang B, Liu K, Czakó G. Correlated Dynamics of the O(3P) + CHD3(v=0) Reaction: A Joint Crossed-Beam and Quasiclassical Trajectory Study. J Phys Chem A 2014; 119:7190-6. [DOI: 10.1021/jp510377k] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Bailin Zhang
- Institute
of Atomic and Molecular Sciences (IAMS), Academia Sinica, P.O. Box 23-166, Taipei, Taiwan 10617
| | - Kopin Liu
- Institute
of Atomic and Molecular Sciences (IAMS), Academia Sinica, P.O. Box 23-166, Taipei, Taiwan 10617
| | - Gábor Czakó
- Laboratory
of Molecular Structure and Dynamics, Institute of Chemistry, Eötvös University, P.O. Box 32, H-1518 Budapest 112, Hungary
| |
Collapse
|
47
|
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
|
48
|
Zhang Y, Stecher T, Cvitaš MT, Althorpe SC. Which Is Better at Predicting Quantum-Tunneling Rates: Quantum Transition-State Theory or Free-Energy Instanton Theory? J Phys Chem Lett 2014; 5:3976-3980. [PMID: 26276480 DOI: 10.1021/jz501889v] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Quantum transition-state theory (QTST) and free-energy instanton theory (FEIT) are two closely related methods for estimating the quantum rate coefficient from the free-energy at the reaction barrier. In calculations on one-dimensional models, FEIT typically gives closer agreement than QTST with the exact quantum results at all temperatures below the crossover to deep tunneling, suggesting that FEIT is a better approximation than QTST in this regime. Here we show that this simple trend does not hold for systems of greater dimensionality. We report tests on several collinear and three-dimensional reactions, in which QTST outperforms FEIT over a range of temperatures below crossover, which can extend down to half the crossover temperature (below which FEIT outperforms QTST). This suggests that QTST-based methods such as ring-polymer molecular dynamics (RPMD) may often give closer agreement with the exact quantum results than FEIT.
Collapse
Affiliation(s)
- Yanchuan Zhang
- †Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Thomas Stecher
- †Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Marko T Cvitaš
- ‡Ruđer Bošković Institute, Department of Physical Chemistry, Bijenička Cesta 54, 10000 Zagreb, Croatia
| | - Stuart C Althorpe
- †Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
49
|
González-Lavado E, Corchado JC, Espinosa-Garcia J. The hydrogen abstraction reaction O(3P) + CH4: a new analytical potential energy surface based on fit to ab initio calculations. J Chem Phys 2014; 140:064310. [PMID: 24527918 DOI: 10.1063/1.4864358] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Based exclusively on high-level ab initio calculations, a new full-dimensional analytical potential energy surface (PES-2014) for the gas-phase reaction of hydrogen abstraction from methane by an oxygen atom is developed. The ab initio information employed in the fit includes properties (equilibrium geometries, relative energies, and vibrational frequencies) of the reactants, products, saddle point, points on the reaction path, and points on the reaction swath, taking especial caution respecting the location and characterization of the intermediate complexes in the entrance and exit channels. By comparing with the reference results we show that the resulting PES-2014 reproduces reasonably well the whole set of ab initio data used in the fitting, obtained at the CCSD(T) = FULL/aug-cc-pVQZ//CCSD(T) = FC/cc-pVTZ single point level, which represents a severe test of the new surface. As a first application, on this analytical surface we perform an extensive dynamics study using quasi-classical trajectory calculations, comparing the results with recent experimental and theoretical data. The excitation function increases with energy (concave-up) reproducing experimental and theoretical information, although our values are somewhat larger. The OH rotovibrational distribution is cold in agreement with experiment. Finally, our results reproduce experimental backward scattering distribution, associated to a rebound mechanism. These results lend confidence to the accuracy of the new surface, which substantially improves the results obtained with our previous surface (PES-2000) for the same system.
Collapse
Affiliation(s)
| | - Jose C Corchado
- Departamento de Química Física, Universidad de Extremadura, 06071 Badajoz, Spain
| | | |
Collapse
|
50
|
Wolf S, Curotto E. Ring polymer dynamics for rigid tops with an improved integrator. J Chem Phys 2014; 141:024116. [DOI: 10.1063/1.4887460] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|