1
|
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
|
2
|
Murakami T, Ibuki S, Hashimoto Y, Kikuma Y, Takayanagi T. Dynamics study of the post-transition-state-bifurcation process of the (HCOOH)H + → CO + H 3O +/HCO + + H 2O dissociation: application of machine-learning techniques. Phys Chem Chem Phys 2023; 25:14016-14027. [PMID: 37161528 DOI: 10.1039/d3cp00252g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
The process of protonated formic acid dissociating from the transition state was studied using ring-polymer molecular dynamics (RPMD), classical MD, and quasi-classical trajectory (QCT) simulations. Temperature had a strong influence on the branching fractions for the HCO+ + H2O and CO + H3O+ dissociation channels. The RPMD and classical MD simulations showed similar behavior, but the QCT dynamics were significantly different owing to the excess energies in the quasi-classical trajectories. Machine-learning analysis identified several key features in the phase information of the vibrational motions at the transition state. We found that the initial configuration and momentum of a hydrogen atom connected to a carbon atom and the shrinking coordinate of the CO bond at the transition state play a role in the dynamics of HCO+ + H2O production.
Collapse
Affiliation(s)
- Tatsuhiro Murakami
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
- Department of Materials & Life Sciences, Faculty of Science & Technology, Sophia University, 7-1 Kioicho, Chiyoda-ku, Tokyo, 102-8554, Japan
| | - Shunichi Ibuki
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| | - Yu Hashimoto
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| | - Yuya Kikuma
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| |
Collapse
|
3
|
Murakami T, Iida R, Hashimoto Y, Takahashi Y, Takahashi S, Takayanagi T. Ring-Polymer Molecular Dynamics and Kinetics for the H – + C 2H 2 → H 2 + C 2H – Reaction Using the Full-Dimensional Potential Energy Surface. J Phys Chem A 2022; 126:9244-9258. [DOI: 10.1021/acs.jpca.2c05851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Tatsuhiro Murakami
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
- Department of Materials & Life Sciences, Faculty of Science & Technology, Sophia University, 7-1 Kioicho, Chiyoda-ku, Tokyo102-8554, Japan
| | - Ryusei Iida
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Yu Hashimoto
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Yukinobu Takahashi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Soma Takahashi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| |
Collapse
|
4
|
Spiekermann KA, Pattanaik L, Green WH. Fast Predictions of Reaction Barrier Heights: Toward Coupled-Cluster Accuracy. J Phys Chem A 2022; 126:3976-3986. [PMID: 35727075 DOI: 10.1021/acs.jpca.2c02614] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Quantitative estimates of reaction barriers are essential for developing kinetic mechanisms and predicting reaction outcomes. However, the lack of experimental data and the steep scaling of accurate quantum calculations often hinder the ability to obtain reliable kinetic values. Here, we train a directed message passing neural network on nearly 24,000 diverse gas-phase reactions calculated at CCSD(T)-F12a/cc-pVDZ-F12//ωB97X-D3/def2-TZVP. Our model uses 75% fewer parameters than previous studies, an improved reaction representation, and proper data splits to accurately estimate performance on unseen reactions. Using information from only the reactant and product, our model quickly predicts barrier heights with a testing MAE of 2.6 kcal mol-1 relative to the coupled-cluster data, making it more accurate than a good density functional theory calculation. Furthermore, our results show that future modeling efforts to estimate reaction properties would significantly benefit from fine-tuning calibration using a transfer learning technique. We anticipate this model will accelerate and improve kinetic predictions for small molecule chemistry.
Collapse
Affiliation(s)
- Kevin A Spiekermann
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Lagnajit Pattanaik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - William H Green
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
5
|
Rangel C, Espinosa-Garcia J. Theoretical study of the O( 3P) + SiH 4 reaction: global potential energy surface, kinetics and dynamics study. Phys Chem Chem Phys 2022; 24:9735-9742. [PMID: 35412544 DOI: 10.1039/d2cp00524g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In order to understand the gas-phase hydrogen abstraction reaction between O(3P) and silane we began by developing the first full-dimensional analytical potential energy surface, named PES-2022. It is basically a valence bond function augmented with molecular mechanic terms describing in an intuitive way stretching and bending nuclei motions, and it is fitted to high level ab initio calculations. The surface presents continuous and smooth behaviour, with analytical first energy derivatives, on which the hydrogen atoms in silane are permutationally symmetric. Based on PES-2022, a kinetics study was performed using the variational transition-state theory with multidimensional tunnelling corrections in the temperature range of 300-1000 K. We observed that experimental and theoretical results show widely spread results, both in absolute value and temperature dependence, possibly because they include the reactivity from both O(3P) and O(1D) electronic states, which present different mechanisms and multiple channels. When the comparison is performed on the same footing, O(3P) + SiH4 → HO + SiH3, the present results agree with Ding and Marshall's experiments and with Zhang et al.'s theoretical rate constants. The kinetic isotope effects (KIEs) reproduced the only experimental value, improving previous theoretical results. Finally, a dynamics study was performed on PES-2022 using quasi-classical trajectory calculations under two different initial conditions, at fixed room temperature and at a fixed collision energy of 8.0 kcal mol-1. In the first case, the available energy deposited as HO(v) vibration was 47%, with population inversion, P(v = 0)/P(v = 1) = 11/89%, reproducing the experimental evidence. In the second case, the experimental product translational distribution was reasonably simulated, while the angular product distribution presented opposite behaviour, backward versus forward. On analysing this discrepancy, we found that while in the present work the O(3P) + SiH4 reaction was reported, in the experiment both O(3P) and O(1D) electronic states are reported. So, the comparison was not performed on the same footing. In sum, agreement of the present results with experiments permits us to be reasonably optimistic about the quality and accuracy of the new PES, and at the same time to highlight the fact that theory/experiment comparisons must be performed on the same footing.
Collapse
Affiliation(s)
- C Rangel
- Área de Química Física and Instituto de Computación Científica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - J Espinosa-Garcia
- Área de Química Física and Instituto de Computación Científica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| |
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
|
Saito K, Hashimoto Y, Takayanagi T. Ring-Polymer Molecular Dynamics Calculations of Thermal Rate Coefficients and Branching Ratios for the Interstellar H 3+ + CO → H 2 + HCO +/HOC + Reaction and Its Deuterated Analogue. J Phys Chem A 2021; 125:10750-10756. [PMID: 34918514 DOI: 10.1021/acs.jpca.1c09160] [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 reaction between H3+ and CO is important in understanding the H3+ destruction mechanism in the interstellar medium. In this work, thermal rate coefficients for the H3+ + CO and D3+ + CO reactions are calculated using ring-polymer molecular dynamics (RPMD) on a high-level machine-learning potential energy surface. The RPMD results agree well with the classical molecular dynamics results, where nuclear quantum effects are completely ignored, whereas the agreement between the RPMD results and the previous quasi-classical trajectory is good only at low temperatures. The calculated [HCO+]/[HOC+] product branching ratios decrease as the temperature increases, and the product branching is exclusively determined by the initial collisional orientation, which governs the formation of an ion-dipole complex, H3+···CO or H3+···OC, that dissociates into H2 + HCO+ or H2 + HOC+, respectively, via a direct mechanism. However, the contribution of the indirect mechanism via the rearrangement between H3+···CO and H3+···OC increases as the temperature increases, although its absolute fraction is small.
Collapse
Affiliation(s)
- Kohei Saito
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| | - Yu Hashimoto
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| |
Collapse
|
8
|
Talotta F, Morisset S, Rougeau N, Lauvergnat D, Agostini F. Electronic Structure and Excited States of the Collision Reaction O( 3P) + C 2H 4: A Multiconfigurational Perspective. J Phys Chem A 2021; 125:6075-6088. [PMID: 34259520 DOI: 10.1021/acs.jpca.1c02923] [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/30/2022]
Abstract
We present a study of the O(3P) + C2H4 scattering reaction, a process that takes place in the interstellar medium and is of relevance in atmospheric chemistry as well. A comprehensive investigation of the electronic properties of the system has been carried out based on multiconfigurational ab initio CASSCF/CASPT2 calculations, using a robust and consistent active space that can deliver accurate potential energy surfaces in the key regions visited by the system. The paper discloses detailed description of the primary reaction pathways and the relevant singlet and triplet excited states at the CASSCF and CASPT2 level, including an accurate description of the critical configurations, such as minima and transition states. The chosen active space and the CASSCF/CASPT2 computational protocol are assessed against coupled-cluster calculations to further check the stability and reliability of the entire multiconfigurational procedure.
Collapse
Affiliation(s)
- Francesco Talotta
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France.,Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, UMR8214, 91405 Orsay, France
| | - Sabine Morisset
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, UMR8214, 91405 Orsay, France
| | - Nathalie Rougeau
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, UMR8214, 91405 Orsay, France
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
9
|
Novikov IS, Suleimanov YV, Shapeev AV. Assessing parameters for ring polymer molecular dynamics simulations at low temperatures: DH + H chemical reaction. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
10
|
Kinetics theoretical study of the O(3P) + C2H6 reaction on an ab initio-based global potential energy surface. Theor Chem Acc 2020. [DOI: 10.1007/s00214-020-02695-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
11
|
Quiton SJ, Mitra U, Mallikarjun Sharada S. A matrix completion algorithm to recover modes orthogonal to the minimum energy path in chemical reactions. J Chem Phys 2020; 153:054122. [PMID: 32770897 DOI: 10.1063/5.0018326] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Structured statistical methods are promising for recovering or completing information from noisy and incomplete data with high fidelity. In particular, matrix completion exploits underlying structural properties such as rank or sparsity. Our objective is to employ matrix completion to reduce computational effort associated with the calculation of multiple quantum chemical Hessians, which are necessary for identification of temperature-dependent free energy maxima under canonical variational transition state theory (VTST). We demonstrate proof-of-principle of an algebraic variety-based matrix completion method for recovering missing elements in a matrix of transverse Hessian eigenvalues constituting the minimum energy path (MEP) of a reaction. The algorithm, named harmonic variety-based matrix completion (HVMC), utilizes the fact that the points lying on the MEP of a reaction step constitute an algebraic variety in the reaction path Hamiltonian representation. We demonstrate that, with as low as 30% random sampling of matrix elements for the largest system in our test set (46 atoms), the complete matrix of eigenvalues can be recovered. We further establish algorithm performance for VTST rate calculations by quantifying errors in zero-point energies and vibrational free energies. Motivated by this success, we outline next steps toward developing a practical HVMC algorithm, which utilizes a gradient-based sampling protocol for low-cost VTST rate computations.
Collapse
Affiliation(s)
- Stephen Jon Quiton
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, USA
| | - Urbashi Mitra
- Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California 90089, USA
| | - Shaama Mallikarjun Sharada
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
12
|
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
|
13
|
Espinosa-Garcia J, García-Chamorro M, Corchado JC. Rethinking the description of water product in polyatomic OH/OD + XH (X ≡ D, Br, NH2 and GeH3) reactions: theory/experimental comparison. Theor Chem Acc 2020. [DOI: 10.1007/s00214-020-2577-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
14
|
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
|
15
|
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
|
16
|
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
|
17
|
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
|
18
|
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
|
19
|
Nuñez-Reyes D, Hickson KM, Larrégaray P, Bonnet L, González-Lezana T, Suleimanov YV. A combined theoretical and experimental investigation of the kinetics and dynamics of the O( 1D) + D 2 reaction at low temperature. Phys Chem Chem Phys 2018; 20:4404-4414. [PMID: 29372194 DOI: 10.1039/c7cp07843a] [Citation(s) in RCA: 21] [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) + H2 reaction is a prototype for simple atom-diatom insertion type mechanisms considered to involve deep potential wells. While exact quantum mechanical methods can be applied to describe the dynamics, such calculations are challenging given the numerous bound quantum states involved. Consequently, efforts have been made to develop alternative theoretical strategies to portray accurately the reactive process. Here we report an experimental and theoretical investigation of the O(1D) + D2 reaction over the 50-296 K range. The calculations employ three conceptually different approaches - mean potential phase space theory, the statistical quantum mechanical method and ring polymer molecular dynamics. The calculated rate constants are in excellent agreement over the entire temperature range, exhibiting only weak temperature dependence. The agreement between experiment and theory is also very good, with discrepancies smaller than 26%. Taken together, the present and previous theoretical results validate the hypothesis that long-lived complex formation dominates the reaction dynamics at low temperature.
Collapse
Affiliation(s)
- Dianailys Nuñez-Reyes
- Université de Bordeaux, Institut des Sciences Moléculaires, F-33400 Talence, France.
| | | | | | | | | | | |
Collapse
|
20
|
Steffen J, Hartke B. Cheap but accurate calculation of chemical reaction rate constants from ab initio data, via system-specific, black-box force fields. J Chem Phys 2017; 147:161701. [DOI: 10.1063/1.4979712] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Vibrational distribution and dynamics study of the HCN(v
1, v
2, v
3) product in the CN + CH4 hydrogen abstraction reaction. Theor Chem Acc 2017. [DOI: 10.1007/s00214-017-2155-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
22
|
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
|
23
|
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
|
24
|
Espinosa-Garcia J, Rangel C, Suleimanov YV. Kinetics study of the CN + CH4 hydrogen abstraction reaction based on a new ab initio analytical full-dimensional potential energy surface. Phys Chem Chem Phys 2017; 19:19341-19351. [DOI: 10.1039/c7cp03499g] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have developed an analytical full-dimensional potential energy surface, named PES-2017, for the gas-phase hydrogen abstraction reaction between the cyano radical and methane.
Collapse
Affiliation(s)
- Joaquin Espinosa-Garcia
- Departamento de Química Física and Instituto de Computación Científica Avanzada
- Universidad de Extremadura
- Badajoz
- Spain
| | - Cipriano Rangel
- Departamento de Química Física and Instituto de Computación Científica Avanzada
- Universidad de Extremadura
- Badajoz
- Spain
| | - Yury V. Suleimanov
- Computation-based Science and Technology Research Center
- Cyprus Institute
- Nicosia 2121
- Cyprus
| |
Collapse
|
25
|
Bao JL, Truhlar DG. Variational transition state theory: theoretical framework and recent developments. Chem Soc Rev 2017; 46:7548-7596. [DOI: 10.1039/c7cs00602k] [Citation(s) in RCA: 207] [Impact Index Per Article: 29.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications.
Collapse
Affiliation(s)
- Junwei Lucas Bao
- Department of Chemistry
- Chemical Theory Center, and Minnesota Supercomputing Institute
- University of Minnesota
- Minneapolis
- USA
| | - Donald G. Truhlar
- Department of Chemistry
- Chemical Theory Center, and Minnesota Supercomputing Institute
- University of Minnesota
- Minneapolis
- USA
| |
Collapse
|
26
|
Fu B, Shan X, Zhang DH, Clary DC. Recent advances in quantum scattering calculations on polyatomic bimolecular reactions. Chem Soc Rev 2017; 46:7625-7649. [DOI: 10.1039/c7cs00526a] [Citation(s) in RCA: 63] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This review surveys quantum scattering calculations on chemical reactions of polyatomic molecules in the gas phase published in the last ten years.
Collapse
Affiliation(s)
- Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| | - Xiao Shan
- Physical and Theoretical Chemistry Laboratory
- Department of Chemistry
- University of Oxford
- Oxford
- UK
| | - 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
- Dalian 116023
- China
| | - David C. Clary
- Physical and Theoretical Chemistry Laboratory
- Department of Chemistry
- University of Oxford
- Oxford
- UK
| |
Collapse
|
27
|
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
|
28
|
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
|
29
|
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
|
30
|
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
|
31
|
Greene SM, Shan X, Clary DC. Rate constants of chemical reactions from semiclassical transition state theory in full and one dimension. J Chem Phys 2016; 144:244116. [DOI: 10.1063/1.4954840] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Samuel M. Greene
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Xiao Shan
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - David C. Clary
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
32
|
Minoshima Y, Seki Y, Takayanagi T, Shiga M. Effects of temperature and isotopic substitution on electron attachment dynamics of guanine–cytosine base pair: Ring-polymer and classical molecular dynamics simulations. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.02.019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
33
|
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
|
34
|
Greene SM, Shan X, Clary DC. An investigation of one- versus two-dimensional semiclassical transition state theory for H atom abstraction and exchange reactions. J Chem Phys 2016; 144:084113. [DOI: 10.1063/1.4942161] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Samuel M. Greene
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Xiao Shan
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - David C. Clary
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
35
|
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
|
36
|
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
|
37
|
Martínez R, Enríquez PA, Puyuelo MP, González M. Exploring the stereodynamics and microscopic mechanism of the O(3P) + CH4, CD4→ OH + CH3, OD + CD3 combustion reactions. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.08.019] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
38
|
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
|
39
|
|
40
|
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
|
41
|
Espinosa-Garcia J. Comment on "Correlated dynamics of the O(3P) + CHD3(v=0) reaction: a joint crossed-beam and quasiclassical trajectory study". J Phys Chem A 2015; 119:1439-41. [PMID: 25654278 DOI: 10.1021/acs.jpca.5b00538] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- J Espinosa-Garcia
- Departamento de Quimica Fisica, Universidad de Extremadura , 06071 Badajoz, Spain
| |
Collapse
|
42
|
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
|
43
|
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
|
44
|
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
|
45
|
Monge-Palacios M, González-Lavado E, Espinosa-Garcia J. Quasiclassical trajectory study of the effect of antisymmetric stretch mode excitation on the O(3P) + CH4(ν3 = 1) → OH + CH3 reaction on an analytical potential energy surface. Comparison with experiment. J Chem Phys 2014; 141:094307. [DOI: 10.1063/1.4893988] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
46
|
Mechanism and kinetics of the atmospheric degradation of perfluoropolymethylisopropyl ether by OH radical: a theoretical study. Struct Chem 2014. [DOI: 10.1007/s11224-014-0448-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|