1
|
Zhang L, Zuo J, Suleimanov YV, Guo H. Ring Polymer Molecular Dynamics Approach to Quantum Dissociative Chemisorption Rates. J Phys Chem Lett 2023; 14:7118-7125. [PMID: 37531595 DOI: 10.1021/acs.jpclett.3c01848] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2023]
Abstract
A ring polymer molecular dynamics (RPMD) method is proposed for the calculation of the dissociative chemisorption rate coefficient on surfaces. The RPMD rate theory is capable of handling quantum effects such as the zero-point energy and tunneling in dissociative chemisorption, while it relies on classical trajectories for the simulation. Applications to H2 dissociative chemisorption are demonstrated. For the highly activated process on Ag(111), strong deviations from Arrhenius behavior are found at low temperatures and attributed to tunneling. On Pt(111), where the dissociation has a barrierless pathway, the RPMD rate coefficient is found to agree with the experimentally derived thermal sticking coefficient within a factor of 2 over a large temperature range. Significant quantum effects are also found.
Collapse
Affiliation(s)
- Liang Zhang
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Junxiang Zuo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yury V Suleimanov
- American Association for the Advancement of Science, 1200 New York Ave NW, Washington, D.C. 20005, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
2
|
Bhowmick S, Maisser A, Suleimanov YV, Schmidt-Ott A, Biskos G. Electronic Structure, Stability, and Electrical Mobility of Cationic Silver Oxide Atomic Clusters. J Phys Chem A 2022; 126:6376-6386. [PMID: 36099558 DOI: 10.1021/acs.jpca.2c02809] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Silver oxide cluster cations (AgnOm+) can readily be produced by a number of methods including atmospheric-pressure spark ablation of pure silver electrodes when trace amounts of oxygen are present in the carrier gas. Here we determine the equilibrium geometries of AgnOm+ clusters (n = 1-4; m = 1-5) using accurate coupled cluster with singles and doubles (CCSD) method, while the stabilization energies are calculated with additional perturbative triples correction (CCSD(T)). Although a number of stable states have been identified, our results show that the AgnOm+ clusters with m = 1 are more stable than those with m ≥ 2 due to the absence of the terminally attached O2 molecule, corroborating recent observations by mass spectrometry. Using the computed structures, we calculate the electrical mobilities of the AgnOm+ clusters and label the values on a respective experimentally determined spectrum in an attempt to better interpret the occurrence of the peaks and troughs in the measurements.
Collapse
Affiliation(s)
- Somnath Bhowmick
- Climate & Atmosphere Research Centre, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus
| | - Anne Maisser
- Climate & Atmosphere Research Centre, 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
| | - Andreas Schmidt-Ott
- Climate & Atmosphere Research Centre, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.,Faculty of Applied Sciences, Delft University of Technology, Delft, 2629 HZ, The Netherlands
| | - George Biskos
- Climate & Atmosphere Research Centre, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.,Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft, 2628 CN, The Netherlands
| |
Collapse
|
3
|
Hickson KM, Bhowmick S, Suleimanov YV, Brandão J, Coelho DV. Experimental and theoretical studies of the gas-phase reactions of O( 1D) with H 2O and D 2O at low temperature. Phys Chem Chem Phys 2021; 23:25797-25806. [PMID: 34761769 DOI: 10.1039/d1cp04614d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Here we report the results of an experimental and theoretical study of the gas-phase reactions between O(1D) and H2O and O(1D) and D2O at room temperature and below. On the experimental side, the kinetics of these reactions have been investigated over the 50-127 K range using a continuous flow Laval nozzle apparatus, coupled with pulsed laser photolysis and pulsed laser induced fluorescence for the production and detection of O(1D) atoms respectively. Experiments were also performed at 296 K in the absence of a Laval nozzle. On the theoretical side, the existing full-dimensional ground X 1A potential energy surface for the H2O2 system involved in this process has been reinvestigated and enhanced to provide a better description of the barrierless H-atom abstraction pathway. Based on this enhanced potential energy surface, quasiclassical trajectory calculations and ring polymer molecular dynamics simulations have been performed to obtain low temperature rate constants. The measured and calculated rate constants display similar behaviour above 100 K, showing little or no variation as a function of temperature. Below 100 K, the experimental rate constants increase dramatically, in contrast to the essentially temperature independent theoretical values. The possible origins of the divergence between experiment and theory at low temperatures are discussed.
Collapse
Affiliation(s)
- Kevin M Hickson
- Université de Bordeaux, Institut des Sciences Moléculaires, F-33400 Talence, France. .,CNRS, Institut des Sciences Moléculaires, F-33400 Talence, France
| | - Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.,Climate & Atmosphere Research Centre, 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
| | - João Brandão
- Departamento de Química e Farmácia - FCT, Universidade do Algarve, Campus de Gambelas, 8005-139 Faro, Portugal
| | - Daniela V Coelho
- Departamento de Química e Farmácia - FCT, Universidade do Algarve, Campus de Gambelas, 8005-139 Faro, Portugal
| |
Collapse
|
4
|
Bhowmick S, Hernández MI, Campos-Martínez J, Suleimanov YV. Isotopic separation of helium through graphyne membranes: a ring polymer molecular dynamics study. Phys Chem Chem Phys 2021; 23:18547-18557. [PMID: 34612392 DOI: 10.1039/d1cp02121d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Microscopic-level understanding of the separation mechanism for two-dimensional (2D) membranes is an active area of research due to potential implications of this class of membranes for various technological processes. Helium (He) purification from the natural resources is of particular interest due to the shortfall in its production. In this work, we applied the ring polymer molecular dynamics (RPMD) method to graphdiyne (Gr2) and graphtriyne (Gr3) 2D membranes having variable pore sizes for the separation of He isotopes, and compare for the first time with rigorous quantum calculations. We found that the transmission rate through Gr3 is many orders of magnitude greater than Gr2. The selectivity of either isotope at low temperatures is a consequence of a delicate balance between the zero-point energy effect and tunneling of 4He and 3He. In particular, a remarkable tunneling effect is reported on the Gr2 membrane at 10 K, leading to a much larger permeation of the lighter species as compared to the heavier isotope. RPMD provides an efficient approach for studying the separation of He isotopes, taking into account quantum effects of light nuclei motions at low temperatures, which classical methods fail to capture.
Collapse
Affiliation(s)
- Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| | | | | | | |
Collapse
|
5
|
Espinosa-Garcia J, Fernandez-Ramos A, Suleimanov YV, Corchado JC. Correction to Theoretical Kinetics Study of the F( 2P) + NH 3 Hydrogen Abstraction Reaction. J Phys Chem A 2021; 125:5709-5710. [PMID: 34137626 DOI: 10.1021/acs.jpca.1c04811] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
6
|
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] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
7
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
8
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
9
|
González-Lezana T, Bossion D, Scribano Y, Bhowmick S, Suleimanov YV. Dynamics of H + HeH +( v = 0, j = 0) → H 2+ + He: Insight on the Possible Complex-Forming Behavior of the Reaction. J Phys Chem A 2019; 123:10480-10489. [PMID: 31725286 DOI: 10.1021/acs.jpca.9b06122] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The H + HeH+→ He + H2+ reaction has been studied by means of a combination of theoretical approaches: a statistical quantum method (SQM), ring polymer molecular dynamics (RPMD), and the quasiclassical trajectory (QCT) method. Cross sections and rate constants have been calculated in an attempt to investigate the dynamics of the process. The comparison with previous calculations and experimental results reveals that despite the fact that statistical predictions seem to reproduce some of the overall observed features, the analysis at a more detailed state-to-state level shows noticeable deviations from a complex-forming dynamics. We find some differences in cross sections and rate constants obtained in the QCT calculation with a Gaussian binning procedure with respect to previous works in which the standard histogram binning was employed.
Collapse
Affiliation(s)
| | - Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299 , Université de Montpellier , 34095 Montpellier Cedex , France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299 , Université de Montpellier , 34095 Montpellier Cedex , France
| | - 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
|
10
|
Nuñez-Reyes D, Hickson KM, Larrégaray P, Bonnet L, González-Lezana T, Bhowmick S, Suleimanov YV. Experimental and Theoretical Study of the O( 1D) + HD Reaction. J Phys Chem A 2019; 123:8089-8098. [PMID: 31464440 DOI: 10.1021/acs.jpca.9b06133] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This work addresses the kinetics and dynamics of the gas-phase reaction between O(1D) and HD molecules down to low temperature. Here, measurements were performed by using a supersonic flow (Laval nozzle) reactor coupled with pulsed laser photolysis for O(1D) production and pulsed-laser-induced fluorescence for O(1D) detection to obtain rate constants over the 50-300 K range. Additionally, temperature-dependent branching ratios (OD + H/OH + D) were obtained experimentally by comparison of the H/D atom atom yields with those of a reference reaction. In parallel, theoretical rate constants and branching ratios were calculated by using three different techniques; mean potential phase space theory (MPPST), the statistical quantum mechanical method (SQM), and ring polymer molecular dynamics (RPMD). Although the agreement between experimental and theoretical rate constants is reasonably good, with differences not exceeding 30% over the entire temperature range, the theoretical branching ratios derived by the MPPST and SQM methods are as much as 50% larger than the experimental ones. These results are presented in the context of earlier work, while the possible origins of the discrepancies between experiment and theory are discussed.
Collapse
Affiliation(s)
- Dianailys Nuñez-Reyes
- Université de Bordeaux, Institut des Sciences Moléculaires , F-33400 Talence , France.,CNRS, Institut des Sciences Moléculaires , UMR 5255, F-33400 Talence , France
| | - Kevin M Hickson
- Université de Bordeaux, Institut des Sciences Moléculaires , F-33400 Talence , France.,CNRS, Institut des Sciences Moléculaires , UMR 5255, F-33400 Talence , France
| | - Pascal Larrégaray
- Université de Bordeaux, Institut des Sciences Moléculaires , F-33400 Talence , France.,CNRS, Institut des Sciences Moléculaires , UMR 5255, F-33400 Talence , France
| | - Laurent Bonnet
- Université de Bordeaux, Institut des Sciences Moléculaires , F-33400 Talence , France.,CNRS, Institut des Sciences Moléculaires , UMR 5255, F-33400 Talence , France
| | - Tomás González-Lezana
- Instituto de Física Fundamental , CSIC , IFF-CSIC Serrano 123 , 28006 Madrid , 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.,Department of Chemical Engineering , Massachusetts Institute of Technology , 77 Massachusetts Avenue , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
11
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
12
|
Naumkin F, del Mazo-Sevillano P, Aguado A, Suleimanov YV, Roncero O. Zero- and high-pressure mechanisms in the complex forming reactions of OH with methanol and formaldehyde at low temperatures. ACS Earth Space Chem 2019; 3:1158-1169. [PMID: 31511842 PMCID: PMC6739233 DOI: 10.1021/acsearthspacechem.9b00051] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
A recent Ring Polymer Molecular Dynamics study of the reactions of OH with methanol and formaldehyde, at zero pressure and below 100 K, has shown the formation of long lived complexes, with long lifetimes, longer than 100 ns for the lower temperatures studied, 20-100 K (del Mazo-Sevillano et al., 2019). These long lifetimes support the existence of multi collision events with the He buffer-gas atoms under experimental conditions, as suggested by several transition state theory studies of these reactions. In this work we study these secondary collisions, as a dynamical approach to study pressure effects on these reactions. For this purpose, the potential energy surfaces of He with H2CO, OH, H2O and HCO are calculated at highly accurate ab initio level. The stability of some of the complexes is studied using Path Integral Molecular dynamics techniques, determining that OH-H2CO complexes can be formed up to 100 K or higher temperatures, while the weaker He-H2CO complexes dissociate at approximately 50 K. The predicted IR intensity spectra shows new features which could help the identification of the OH-H2CO complex. Finally, the He-H2CO + OH and OH-H2CO + He collisions are studied using quassi-classical trajectories, finding that the cross section to produce HCO + H2O products increases with decreasing collision energy, and that it is ten times higher in the He-H2CO + OH case.
Collapse
Affiliation(s)
| | - Pablo del Mazo-Sevillano
- Unidad Asociada UAM-IFF-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias Módulo 14, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | - Alfredo Aguado
- Unidad Asociada UAM-IFF-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias Módulo 14, Universidad Autónoma de Madrid, 28049, 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, 77 Massachusetts Ave., Cambridge, Massachusetts 02139, United States
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
13
|
del Mazo-Sevillano P, Aguado A, Jiménez E, Suleimanov YV, Roncero O. Quantum Roaming in the Complex-Forming Mechanism of the Reactions of OH with Formaldehyde and Methanol at Low Temperature and Zero Pressure: A Ring Polymer Molecular Dynamics Approach. J Phys Chem Lett 2019; 10:1900-1907. [PMID: 30939028 PMCID: PMC6534501 DOI: 10.1021/acs.jpclett.9b00555] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
The quantum dynamics of the title reactions are studied using the ring polymer molecular dynamics (RPMD) method from 20 to 1200 K using recently proposed full dimensional potential energy surfaces which include long-range dipole-dipole interactions. A V-shaped dependence of the reaction rate constants is found with a minimum at 200-300 K, in rather good agreement with the current experimental data. For temperatures above 300 K the reaction proceeds following a direct H-abstraction mechanism. However, below 100 K the reaction proceeds via organic-molecule···OH collision complexes, with very long lifetimes, longer than 10-7 s, associated with quantum roaming arising from the inclusion of quantum effects by the use of RPMD. The long lifetimes of these complexes are comparable to the time scale of the tunnelling to form reaction products. These complexes are formed at zero pressure because of quantum effects and not only at high pressure as suggested by transition state theory (TST) calculations for OH + methanol and other OH reactions. The zero-pressure rate constants reproduce quite well measured ones below 200 K, and this agreement opens the question of how important the pressure effects on the reaction rate constants are, as implied in TST-like formalisms. The zero-pressure mechanism is applicable only to very low gas density environments, such as the interstellar medium, which are not repeatable by experiments.
Collapse
Affiliation(s)
- Pablo del Mazo-Sevillano
- Unidad Asociada UAM-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | - Alfredo Aguado
- Unidad Asociada UAM-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | - Elena Jiménez
- Departamento de Química Física, Facultad de Ciencias y Tecnologías Químicas, Universidad de Castilla La Mancha, Avda. Camilo José Cela 1B, 13071 Ciudad Real, Spain
- Instituto de Investigación en Combustión y Contaminación Atmosférica, Universidad de Castilla La Mancha, Camino de Moledores s/n, 13071 Ciudad Real, 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, 77 Massachusetts Ave., Cambridge, Massachusetts 02139, United States
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
14
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
15
|
Kumar SS, Grussie F, Suleimanov YV, Guo H, Kreckel H. Low temperature rates for key steps of interstellar gas-phase water formation. Sci Adv 2018; 4:eaar3417. [PMID: 29942857 PMCID: PMC6014714 DOI: 10.1126/sciadv.aar3417] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/27/2017] [Accepted: 05/15/2018] [Indexed: 06/08/2023]
Abstract
The gas-phase formation of water molecules in the diffuse interstellar medium (ISM) proceeds mainly via a series of reactions involving the molecular ions OH+, H2O+, and H3O+ and molecular hydrogen. These reactions form the backbone for the chemistry leading to the formation of several complex molecular species in space. A comprehensive understanding of the mechanisms involved in these reactions in the ISM necessitates an accurate knowledge of the rate coefficients at the relevant temperatures (10 to 100 K). We present measurements of the rate coefficients for two key reactions below 100 K, which, in both cases, are significantly higher than the values used in astronomical models thus far. The experimental rate coefficients show excellent agreement with dedicated theoretical calculations using a novel ring-polymer molecular dynamics approach that offers a first-principles treatment of low-temperature barrierless gas-phase reactions, which are prevalent in interstellar chemical networks.
Collapse
Affiliation(s)
- Sunil S. Kumar
- Max-Planck-Institut für Kernphysik, 69117 Heidelberg, Germany
| | - Florian Grussie
- Max-Planck-Institut für Kernphysik, 69117 Heidelberg, Germany
| | - 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, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, NM 87131, USA
| | - Holger Kreckel
- Max-Planck-Institut für Kernphysik, 69117 Heidelberg, Germany
| |
Collapse
|
16
|
Suleimanov YV, Aguado A, Gómez-Carrasco S, Roncero O. A Ring Polymer Molecular Dynamics Approach to Study the Transition between Statistical and Direct Mechanisms in the H 2 + H 3+ → H 3+ + H 2 Reaction. J Phys Chem Lett 2018; 9:2133-2137. [PMID: 29633841 PMCID: PMC6031303 DOI: 10.1021/acs.jpclett.8b00783] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Because of its fundamental importance in astrochemistry, the H2 + H3+ → H3+ + H2 reaction has been studied experimentally in a wide temperature range. Theoretical studies of the title reaction significantly lag primarily because of the challenges associated with the proper treatment of the zero-point energy (ZPE). As a result, all previous theoretical estimates for the ratio between a direct proton-hop and indirect exchange (via the H5+ complex) channels deviate from the experiment, in particular, at lower temperatures where the quantum effects dominate. In this work, the ring polymer molecular dynamics (RPMD) method is applied to study this reaction, providing very good agreement with the experiment. RPMD is immune to the shortcomings associated with the ZPE leakage and is able to describe the transition from direct to indirect mechanisms below room temperature. We argue that RPMD represents a useful tool for further studies of numerous ZPE-sensitive chemical reactions that are of high interest in astrochemistry.
Collapse
Affiliation(s)
- 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 Ave., Cambridge, Massachusetts 02139, United States
| | - Alfredo Aguado
- Unidad Asociada UAM-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | | | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
17
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
18
|
Grambow CA, Jamal A, Li YP, Green WH, Zádor J, Suleimanov YV. Unimolecular Reaction Pathways of a γ-Ketohydroperoxide from Combined Application of Automated Reaction Discovery Methods. J Am Chem Soc 2018; 140:1035-1048. [DOI: 10.1021/jacs.7b11009] [Citation(s) in RCA: 59] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Colin A. Grambow
- Department
of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave., Cambridge, Massachusetts 02139, United States
| | - Adeel Jamal
- Department
of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave., Cambridge, Massachusetts 02139, United States
| | - Yi-Pei Li
- Department
of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave., Cambridge, Massachusetts 02139, United States
| | - William H. Green
- Department
of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave., Cambridge, Massachusetts 02139, United States
| | - Judit Zádor
- Combustion
Research Facility, Sandia National Laboratories, 7011 East Ave., Livermore, California 94551, United States
| | - Yury V. Suleimanov
- Computation-based
Science and Technology Research Center, Cyprus Institute, 20
Kavafi Str., Nicosia 2121, Cyprus
| |
Collapse
|
19
|
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
|
20
|
Mägi R, Suleimanov YV, Clarke GM, Kaakinen M, Fischer K, Prokopenko I, Morris AP. SCOPA and META-SCOPA: software for the analysis and aggregation of genome-wide association studies of multiple correlated phenotypes. BMC Bioinformatics 2017; 18:25. [PMID: 28077070 PMCID: PMC5225593 DOI: 10.1186/s12859-016-1437-3] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2016] [Accepted: 12/17/2016] [Indexed: 11/10/2022] Open
Abstract
Background Genome-wide association studies (GWAS) of single nucleotide polymorphisms (SNPs) have been successful in identifying loci contributing genetic effects to a wide range of complex human diseases and quantitative traits. The traditional approach to GWAS analysis is to consider each phenotype separately, despite the fact that many diseases and quantitative traits are correlated with each other, and often measured in the same sample of individuals. Multivariate analyses of correlated phenotypes have been demonstrated, by simulation, to increase power to detect association with SNPs, and thus may enable improved detection of novel loci contributing to diseases and quantitative traits. Results We have developed the SCOPA software to enable GWAS analysis of multiple correlated phenotypes. The software implements “reverse regression” methodology, which treats the genotype of an individual at a SNP as the outcome and the phenotypes as predictors in a general linear model. SCOPA can be applied to quantitative traits and categorical phenotypes, and can accommodate imputed genotypes under a dosage model. The accompanying META-SCOPA software enables meta-analysis of association summary statistics from SCOPA across GWAS. Application of SCOPA to two GWAS of high-and low-density lipoprotein cholesterol, triglycerides and body mass index, and subsequent meta-analysis with META-SCOPA, highlighted stronger association signals than univariate phenotype analysis at established lipid and obesity loci. The META-SCOPA meta-analysis also revealed a novel signal of association at genome-wide significance for triglycerides mapping to GPC5 (lead SNP rs71427535, p = 1.1x10−8), which has not been reported in previous large-scale GWAS of lipid traits. Conclusions The SCOPA and META-SCOPA software enable discovery and dissection of multiple phenotype association signals through implementation of a powerful reverse regression approach.
Collapse
Affiliation(s)
- Reedik Mägi
- Estonian Genome Center, University of Tartu, Tartu, Estonia
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center, Cyprus Institute, Nicosia, Cyprus.,Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Geraldine M Clarke
- Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, UK
| | | | - Krista Fischer
- Estonian Genome Center, University of Tartu, Tartu, Estonia
| | | | - Andrew P Morris
- Estonian Genome Center, University of Tartu, Tartu, Estonia. .,Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, UK. .,Department of Biostatistics, University of Liverpool, Liverpool, UK.
| |
Collapse
|
21
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
22
|
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
|
23
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
24
|
Suleimanov YV, Aoiz FJ, Guo H. ChemInform Abstract: Chemical Reaction Rate Coefficients from Ring Polymer Molecular Dynamics: Theory and Practical Applications. ACTA ACUST UNITED AC 2016. [DOI: 10.1002/chin.201650259] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
25
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
26
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
27
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
28
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
29
|
Kabeshov MA, Kysilka O, Rulíšek L, Suleimanov YV, Bella M, Malkov AV, Kočovský P. Cross-Aldol Reaction of Isatin with Acetone Catalyzed by Leucinol: A Mechanistic Investigation. Chemistry 2015; 21:12026-33. [PMID: 26147182 DOI: 10.1002/chem.201500536] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2015] [Indexed: 11/09/2022]
Abstract
Comprehensive mechanistic studies on the enantioselective aldol reaction between isatin (1 a) and acetone, catalyzed by L-leucinol (3 a), unraveled that isatin, apart from being a substrate, also plays an active catalytic role. Conversion of the intermediate oxazolidine 4 into the reactive syn-enamine 6, catalyzed by isatin, was identified as the rate-determining step by both the calculations (ΔG(≠) =26.1 kcal mol(-1) for the analogous L-alaninol, 3 b) and the kinetic isotope effect (kH /kD =2.7 observed for the reaction using [D6 ]acetone). The subsequent reaction of the syn-enamine 6 with isatin produces (S)-2 a (calculated ΔG(≠) =11.6 kcal mol(-1) ). The calculations suggest that the overall stereochemistry is controlled by two key events: 1) the isatin-catalyzed formation of the syn-enamine 6, which is thermodynamically favored over its anti-rotamer 7 by 2.3 kcal mol(-1) ; and 2) the high preference of the syn-enamine 6 to produce (S)-2 a on reaction with isatin (1 a) rather than its enantiomer (ΔΔG(≠) =2.6 kcal mol(-1) ).
Collapse
Affiliation(s)
- Mikhail A Kabeshov
- Department of Chemistry, Loughborough University, Loughborough, Leics LE11 3TU (UK). .,Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW (UK).
| | - Ondřej Kysilka
- FARMAK, Klášterní Hradisko, 77900 Olomouc (Czech Republic)
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Prague 6 (Czech Republic)
| | - 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 (USA)
| | - Marco Bella
- Department of Chemistry, "Sapienza" University of Rome, Piazzale Aldo Moro 5, 00185 Rome (Italy)
| | - Andrei V Malkov
- Department of Chemistry, Loughborough University, Loughborough, Leics LE11 3TU (UK).
| | - Pavel Kočovský
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Prague 6 (Czech Republic). .,Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE 10691 Stockholm (Sweden) & Department of Organic Chemistry, Charles University, Hlavova 8, 12843 Prague 2 (Czech Republic).
| |
Collapse
|
30
|
Suleimanov YV, Green WH. Automated Discovery of Elementary Chemical Reaction Steps Using Freezing String and Berny Optimization Methods. J Chem Theory Comput 2015; 11:4248-59. [DOI: 10.1021/acs.jctc.5b00407] [Citation(s) in RCA: 99] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [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
| | - William H. Green
- Department
of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
31
|
Kabeshov MA, Kysilka O, Rulíšek L, Suleimanov YV, Bella M, Malkov AV, Kočovský P. Inside Back Cover: Cross-Aldol Reaction of Isatin with Acetone Catalyzed by Leucinol: A Mechanistic Investigation (Chem. Eur. J. 34/2015). Chemistry 2015. [DOI: 10.1002/chem.201590155] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
32
|
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
|
33
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
34
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
35
|
Prozument K, Suleimanov YV, Buesser B, Oldham JM, Green WH, Suits AG, Field RW. A Signature of Roaming Dynamics in the Thermal Decomposition of Ethyl Nitrite: Chirped-Pulse Rotational Spectroscopy and Kinetic Modeling. J Phys Chem Lett 2014; 5:3641-3648. [PMID: 26278732 DOI: 10.1021/jz501758p] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Chirped-pulse (CP) Fourier transform rotational spectroscopy is uniquely suited for near-universal quantitative detection and structural characterization of mixtures that contain multiple molecular and radical species. In this work, we employ CP spectroscopy to measure product branching and extract information about the reaction mechanism, guided by kinetic modeling. Pyrolysis of ethyl nitrite, CH3CH2ONO, is studied in a Chen type flash pyrolysis reactor at temperatures of 1000-1800 K. The branching between HNO, CH2O, and CH3CHO products is measured and compared to the kinetic models generated by the Reaction Mechanism Generator software. We find that roaming CH3CH2ONO → CH3CHO + HNO plays an important role in the thermal decomposition of ethyl nitrite, with its rate, at 1000 K, comparable to that of the radical elimination channel CH3CH2ONO → CH3CH2O + NO. HNO is a signature of roaming in this system.
Collapse
Affiliation(s)
| | | | - Beat Buesser
- §IBM Research, Smarter Cities Technology Centre, Dublin 15, Ireland
| | - James M Oldham
- ∥Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | | | - Arthur G Suits
- ∥Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | | |
Collapse
|
36
|
Gonzalez-Lavado E, Corchado JC, Suleimanov YV, Green WH, Espinosa-Garcia J. Theoretical Kinetics Study of the O(3P) + CH4/CD4 Hydrogen Abstraction Reaction: The Role of Anharmonicity, Recrossing Effects, and Quantum Mechanical Tunneling. J Phys Chem A 2014; 118:3243-52. [DOI: 10.1021/jp5028965] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Jose C. Corchado
- Departamento
de Química Física, Universidad de Extremadura, 06071 Badajoz, Spain
| | - Yury V. Suleimanov
- Department of Chemical
Engineering, Massachusetts Institute of Technology, 77 Massachusetts
Avenue, Cambridge, Massachusetts 02139, United States
- Department of Mechanical and Aerospace Engineering, Combustion Energy
Frontier Research Center, Princeton University, Olden Street, Princeton, New Jersey 08544, United States
| | - William H. Green
- Department of Chemical
Engineering, Massachusetts Institute of Technology, 77 Massachusetts
Avenue, Cambridge, Massachusetts 02139, United States
| | | |
Collapse
|
37
|
Li Y, Suleimanov YV, Green WH, Guo H. Quantum rate coefficients and kinetic isotope effect for the reaction Cl + CH4 → HCl + CH3 from ring polymer molecular dynamics. J Phys Chem A 2014; 118:1989-96. [PMID: 24558961 DOI: 10.1021/jp501043z] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Thermal rate coefficients and kinetic isotope effect have been calculated for prototypical heavy-light-heavy polyatomic bimolecular reactions Cl + CH4/CD4 → HCl/DCl + CH3/CD3, using a recently proposed quantum dynamics approach: ring polymer molecular dynamics (RPMD). Agreement with experimental rate coefficients, which are quite scattered, is satisfactory. However, differences up to 50% have been found between the RPMD results and those obtained from the harmonic variational transition-state theory on one of the two full-dimensional potential energy surfaces used in the calculations. Possible reasons for such discrepancy are discussed. The present work is an important step in a series of benchmark studies aimed at assessing accuracy for RPMD for chemical reaction rates, which demonstrates that this novel method is a quite reliable alternative to previously developed techniques based on transition-state theory.
Collapse
Affiliation(s)
- Yongle Li
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
| | | | | | | |
Collapse
|
38
|
Li Y, Suleimanov YV, Guo H. Ring-Polymer Molecular Dynamics Rate Coefficient Calculations for Insertion Reactions: X + H2 → HX + H (X = N, O). J Phys Chem Lett 2014; 5:700-705. [PMID: 26270840 DOI: 10.1021/jz500062q] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The thermal rate constants of two prototypical insertion-type reactions, namely, N/O + H2 → NH/OH + H, are investigated with ring polymer molecular dynamics (RPMD) on full-dimensional potential energy surfaces using recently developed RPMDrate code. It is shown that the unique ability of the RPMD approach among the existing theoretical methods to capture the quantum effects, e.g., tunneling and zero-point energy, as well as recrossing dynamics quantum mechanically with ring-polymer trajectories leads to excellent agreement with rigorous quantum dynamics calculations. The present result is encouraging for future applications of the RPMD method and the RPMDrate code to complex-forming chemical reactions involving polyatomic reactants.
Collapse
Affiliation(s)
- Yongle Li
- †Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yury V Suleimanov
- ‡Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- §Department of Mechanical and Aerospace Engineering, Combustion Energy Frontier Research Center, Princeton University, Olden Street, Princeton, New Jersey 08544, United States
| | - Hua Guo
- †Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
39
|
Allen JW, Green WH, Li Y, Guo H, Suleimanov YV. Communication: Full dimensional quantum rate coefficients and kinetic isotope effects from ring polymer molecular dynamics for a seven-atom reaction OH + CH4 → CH3 + H2O. J Chem Phys 2013; 138:221103. [DOI: 10.1063/1.4811329] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
40
|
Li Y, Suleimanov YV, Li J, Green WH, Guo H. Rate coefficients and kinetic isotope effects of the X + CH4 → CH3 + HX (X = H, D, Mu) reactions from ring polymer molecular dynamics. J Chem Phys 2013; 138:094307. [DOI: 10.1063/1.4793394] [Citation(s) in RCA: 68] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Li Y, Suleimanov YV, Yang M, Green WH, Guo H. Ring Polymer Molecular Dynamics Calculations of Thermal Rate Constants for the O((3)P) + CH4 → OH + CH3 Reaction: Contributions of Quantum Effects. J Phys Chem Lett 2013; 4:48-52. [PMID: 26291210 DOI: 10.1021/jz3019513] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The thermal rate constant of the O((3)P) + CH4 → OH + CH3 reaction is investigated with ring polymer molecular dynamics on a full-dimensional potential energy surface. Good agreement with experimental and full-dimensional quantum multiconfiguration time-dependent Hartree results between 300 and 1500 K was obtained. It is shown that quantum effects, for example, tunneling and zero-point energy, can be effectively and efficiently included in this path-integral based approach implemented with classical trajectories. Convergence with respect to the number of beads is rapid, suggesting wide applicability for other reactions involving polyatomic molecules.
Collapse
Affiliation(s)
- Yongle Li
- †Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yury V Suleimanov
- ‡Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Minghui Yang
- §Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Centre for Magnetic Resonance, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan, 430071, China
| | - William H Green
- ‡Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Hua Guo
- †Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
42
|
Suleimanov YV, de Tudela RP, Jambrina PG, Castillo JF, Sáez-Rábanos V, Manolopoulos DE, Aoiz FJ. A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction. Phys Chem Chem Phys 2013; 15:3655-65. [DOI: 10.1039/c2cp44364c] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
43
|
Pérez de Tudela R, Aoiz FJ, Suleimanov YV, Manolopoulos DE. Chemical Reaction Rates from Ring Polymer Molecular Dynamics: Zero Point Energy Conservation in Mu + H2 → MuH + H. J Phys Chem Lett 2012; 3:493-497. [PMID: 26286053 DOI: 10.1021/jz201702q] [Citation(s) in RCA: 86] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
A fundamental issue in the field of reaction dynamics is the inclusion of the quantum mechanical (QM) effects such as zero point energy (ZPE) and tunneling in molecular dynamics simulations, and in particular in the calculation of chemical reaction rates. In this work we study the chemical reaction between a muonium atom and a hydrogen molecule. The recently developed ring polymer molecular dynamics (RPMD) technique is used, and the results are compared with those of other methods. For this reaction, the thermal rate coefficients calculated with RPMD are found to be in excellent agreement with the results of an accurate QM calculation. The very minor discrepancies are within the convergence error even at very low temperatures. This exceptionally good agreement can be attributed to the dominant role of ZPE in the reaction, which is accounted for extremely well by RPMD. Tunneling only plays a minor role in the reaction.
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
| | - F J Aoiz
- †Departamento de Química Física I, Facultad de CC. Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Yury V Suleimanov
- ‡Physical and Theoretical Chemistry Laboratory, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - David E Manolopoulos
- ‡Physical and Theoretical Chemistry Laboratory, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
44
|
Suleimanov YV, Collepardo-Guevara R, Manolopoulos DE. Bimolecular reaction rates from ring polymer molecular dynamics: application to H + CH4 → H2 + CH3. J Chem Phys 2011; 134:044131. [PMID: 21280711 DOI: 10.1063/1.3533275] [Citation(s) in RCA: 138] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In a recent paper, we have developed an efficient implementation of the ring polymer molecular dynamics (RPMD) method for calculating bimolecular chemical reaction rates in the gas phase, and illustrated it with applications to some benchmark atom-diatom reactions. In this paper, we show that the same methodology can readily be used to treat more complex polyatomic reactions in their full dimensionality, such as the hydrogen abstraction reaction from methane, H + CH(4) → H(2) + CH(3). The present calculations were carried out using a modified and recalibrated version of the Jordan-Gilbert potential energy surface. The thermal rate coefficients obtained between 200 and 2000 K are presented and compared with previous results for the same potential energy surface. Throughout the temperature range that is available for comparison, the RPMD approximation gives better agreement with accurate quantum mechanical (multiconfigurational time-dependent Hartree) calculations than do either the centroid density version of quantum transition state theory (QTST) or the quantum instanton (QI) model. The RPMD rate coefficients are within a factor of 2 of the exact quantum mechanical rate coefficients at temperatures in the deep tunneling regime. These results indicate that our previous assessment of the accuracy of the RPMD approximation for atom-diatom reactions remains valid for more complex polyatomic reactions. They also suggest that the sensitivity of the QTST and QI rate coefficients to the choice of the transition state dividing surface becomes more of an issue as the dimensionality of the reaction increases.
Collapse
Affiliation(s)
- Yury V Suleimanov
- Physical and Theoretical Chemistry Laboratory, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom.
| | | | | |
Collapse
|
45
|
Suleimanov YV, Krems RV. Efficient numerical method for locating Feshbach resonances of ultracold molecules in external fields. J Chem Phys 2011. [DOI: 10.1063/1.3512627] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
46
|
Collepardo-Guevara R, Suleimanov YV, Manolopoulos DE. Erratum: “Bimolecular reaction rates from ring polymer molecular dynamics” [J. Chem. Phys. 130, 174713 (2009)]. J Chem Phys 2010. [DOI: 10.1063/1.3464477] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
47
|
Collepardo-Guevara R, Suleimanov YV, Manolopoulos DE. Publisher’s Note: “Bimolecular reaction rates from ring polymer molecular dynamics” [J. Chem. Phys. 130, 174713 (2009)]. J Chem Phys 2009. [DOI: 10.1063/1.3147536] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
48
|
Collepardo-Guevara R, Suleimanov YV, Manolopoulos DE. Bimolecular reaction rates from ring polymer molecular dynamics. J Chem Phys 2009; 130:174713. [DOI: 10.1063/1.3127145] [Citation(s) in RCA: 111] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
49
|
Akopyan ME, Lukashov SS, Poretsky SA, Pravilov AM, Torgashkova AS, Buchachenko AA, Suleimanov YV. Dynamics and mechanism of the E-->D, D', beta, gamma, and delta nonadiabatic transitions induced in molecular iodine by collisions with CF4 and SF6 molecules. J Chem Phys 2008; 129:114309. [PMID: 19044962 DOI: 10.1063/1.2978170] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Nonadiabatic transitions among the first-tier ion-pair states of the iodine molecule in collisions with CF(4) and SF(6) partners are investigated by detecting the luminescence following the optical-optical double resonance excitation of the E0(g) (+)-state to the vibrational levels v(E)=8, 13, and 19. Total and partial rate constants, as well as vibrational product state distributions, are determined. It is found that electronic energy transfer in all channels is predominantly assisted by excitation of the dipole-allowed nu(3) and nu(4) modes of the partner. The measurements are accompanied by quantum scattering calculations that implement a close coupling treatment for the electronic and vibrational degrees of freedom and combine diatomics-in-molecule and long-range models for diabatic potential energy surfaces and coupling matrix elements. The analysis of experimental and theoretical data shows that the transitions without excitation of the partner are due to short-range couplings, whereas the vibrational excitation of the partner in the D0(u) (+) channel originates from the long-range coupling of two transition dipole moments: electronic of the iodine molecule and vibrational of the partner. Unexpectedly efficient excitations of the partner in the other ion-pair states, which are not coupled to the initial E0(g) (+)-state by the transition dipole, are interpreted within the postcollision mechanism. Qualitatively, this implies that during a single collision the long-range nonadiabatic transitions to D, nu(3) and D, nu(4) channels are followed by secondary short-range transitions without changing the state of the partner.
Collapse
Affiliation(s)
- Mikhail E Akopyan
- V. A. Fock Institute of Physics, Physical Faculty, St. Petersburg State University, Ul'yanovskaya 1, Staryj Peterhof, 198504 St. Petersburg, Russia
| | | | | | | | | | | | | |
Collapse
|
50
|
Hutchison JM, O'Hern RR, Stephenson TA, Suleimanov YV, Buchachenko AA. Theoretical and experimental studies of collision-induced electronic energy transfer from v=0-3 of the E(0g +) ion-pair state of Br2: collisions with He and Ar. J Chem Phys 2008; 128:184311. [PMID: 18532816 DOI: 10.1063/1.2912057] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Collisions of Br(2), prepared in the E(0(g)+) ion-pair (IP) electronic state, with He or Ar result in electronic energy transfer to the D, D', and beta IP states. These events have been examined in experimental and theoretical investigations. Experimentally, analysis of the wavelength resolved emission spectra reveals the distribution of population in the vibrational levels of the final electronic states and the relative efficiencies of He and Ar collisions in promoting a specific electronic energy transfer channel. Theoretically, semiempirical rare gas-Br(2) potential energy surfaces and diabatic couplings are used in quantum scattering calculations of the state-to-state rate constants for electronic energy transfer and distributions of population in the final electronic state vibrational levels. Agreement between theory and experiment is excellent. Comparison of the results with those obtained for similar processes in the IP excited I(2) molecule points to the general importance of Franck-Condon effects in determining vibrational populations, although this effect is more important for He collisions than for Ar collisions.
Collapse
Affiliation(s)
- J Matthew Hutchison
- Department of Chemistry and Biochemistry, Swarthmore College, Swarthmore, PA 19081, USA
| | | | | | | | | |
Collapse
|