1
|
Welz O, Pfeifle M, Plehiers PM, Sure R, Deglmann P. Reaction of OH with Aliphatic and Aromatic Isocyanates. J Phys Chem A 2022; 126:9333-9352. [DOI: 10.1021/acs.jpca.2c06011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Oliver Welz
- BASF SE, Scientific Modelling − Quantum Chemistry, Group Research, Carl-Bosch-Straße 38, Ludwigshafen am Rhein67056, Germany
| | - Mark Pfeifle
- BASF SE, Scientific Modelling − Quantum Chemistry, Group Research, Carl-Bosch-Straße 38, Ludwigshafen am Rhein67056, Germany
| | - Patrick M. Plehiers
- International Isocyanate Institute Inc. (III), 333 Route 46 West, Suite. 206, Mountain Lakes, New Jersey07046, United States
| | - Rebecca Sure
- BASF SE, Scientific Modelling − Quantum Chemistry, Group Research, Carl-Bosch-Straße 38, Ludwigshafen am Rhein67056, Germany
| | - Peter Deglmann
- BASF SE, Scientific Modelling − Quantum Chemistry, Group Research, Carl-Bosch-Straße 38, Ludwigshafen am Rhein67056, Germany
| |
Collapse
|
2
|
Klippenstein SJ. Spiers Memorial Lecture: theory of unimolecular reactions. Faraday Discuss 2022; 238:11-67. [DOI: 10.1039/d2fd00125j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our...
Collapse
|
3
|
Rath G, Kopp WA, Leonhard K. Coupled Anharmonic Thermochemistry from Stratified Monte Carlo Integration. J Chem Inf Model 2021; 61:5853-5870. [PMID: 34874733 DOI: 10.1021/acs.jcim.1c00668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This study presents configuration integral Monte Carlo integration (CIMCI), a new semiclassical method for handling fully coupled anharmonicity in gas-phase thermodynamics that promises to be black boxable, to be applicable to all kinds of anharmonicity, and to scale better at higher dimensionality than other methods for handling gas-phase molecular anharmonicity. The method does so using automatically and recursively stratified, simultaneous Monte Carlo (MC) integration of multiple functions, following a modified version of the standard MISER scheme that converges at a rate of about the square of naïve MC integration. For the small systems analyzed by this study where proper reference data is available (H2O and H2O2), the method's anharmonic entropy corrections match reference data better than those of other black box anharmonic methods, e.g., vibrational perturbation theory (VPT2) and the McClurg hindered rotor model used with automatic detection of rotors; for H2O2 and NH2OH, the method is also in general agreement with one-dimensional hindered rotor treatments at low temperatures. This holds even when sampling with CIMCI is done with primitive force fields, e.g., UFF, while the competing methods are used with proper, comprehensive potentials, e.g., the M06-2X metahybrid density-functional theory (DFT) functional.
Collapse
Affiliation(s)
- Gabriel Rath
- Institute of Technical Thermodynamics, RWTH Aachen University, 52062 Aachen, Germany.,Software for Chemistry & Materials, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Wassja A Kopp
- Institute of Technical Thermodynamics, RWTH Aachen University, 52062 Aachen, Germany
| | - Kai Leonhard
- Institute of Technical Thermodynamics, RWTH Aachen University, 52062 Aachen, Germany
| |
Collapse
|
4
|
Rousso AC, Jasper AW, Ju Y, Hansen N. Extreme Low-Temperature Combustion Chemistry: Ozone-Initiated Oxidation of Methyl Hexanoate. J Phys Chem A 2020; 124:9897-9914. [PMID: 33174431 DOI: 10.1021/acs.jpca.0c07584] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The accelerating chemical effect of ozone addition on the oxidation chemistry of methyl hexanoate [CH3(CH2)4C(═O)OCH3] was investigated over a temperature range from 460 to 940 K. Using an externally heated jet-stirred reactor at p = 700 Torr (residence time τ = 1.3 s, stoichiometry φ = 0.5, 80% argon dilution), we explored the relevant chemical pathways by employing molecular-beam mass spectrometry with electron and single-photon ionization to trace the temperature dependencies of key intermediates, including many hydroperoxides. In the absence of ozone, reactivity is observed in the so-called low-temperature chemistry (LTC) regime between 550 and 700 K, which is governed by hydroperoxides formed from sequential O2 addition and isomerization reactions. At temperatures above 700 K, we observed the negative temperature coefficient (NTC) regime, in which the reactivity decreases with increasing temperatures, until near 800 K, where the reactivity increases again. Upon addition of ozone (1000 ppm), the overall reactivity of the system is dramatically changed due to the time scale of ozone decomposition in comparison to fuel oxidation time scales of the mixtures at different temperatures. While the LTC regime seems to be only slightly affected by the addition of ozone with respect to the identity and quantity of the observed intermediates, we observed an increased reactivity in the intermediate NTC temperature range. Furthermore, we observed experimental evidence for an additional oxidation regime in the range near 500 K, herein referred to as the extreme low-temperature chemistry (ELTC) regime. Experimental evidence and theoretical rate constant calculations indicate that this ELTC regime is likely to be initiated by H abstraction from methyl hexanoate via O atoms, which originate from thermal O3 decomposition. The theoretical calculations show that the rate constants for methyl ester initiation via abstraction by O atoms increase dramatically with the size of the methyl ester, suggesting that ELTC is likely not important for the smaller methyl esters. Experimental evidence is provided indicating that, similar to the LTC regime, the chemistry in the ELTC regime is dominated by hydroperoxide chemistry. However, mass spectra recorded at various reactor temperatures and at different photon energies provide experimental evidence of some differences in chemical species between the ELTC and the LTC temperature ranges.
Collapse
Affiliation(s)
- Aric C Rousso
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ahren W Jasper
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Yiguang Ju
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Nils Hansen
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, United States
| |
Collapse
|
5
|
Sun Y, Zhou CW, Somers KP, Curran HJ. An ab Initio/Transition State Theory Study of the Reactions of Ċ 5H 9 Species of Relevance to 1,3-Pentadiene, Part II: Pressure Dependent Rate Constants and Implications for Combustion Modeling. J Phys Chem A 2020; 124:4605-4631. [PMID: 32396376 DOI: 10.1021/acs.jpca.0c02244] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The temperature- and pressure-dependence of rate constants for several radicals and unsaturated hydrocarbons reactions (1,3-C5H8/1,4-C5H8/cyC5H8 + Ḣ, C2H4 + Ċ3H5-a, C3H6 + Ċ2H3) are analyzed in this paper. The abstraction reactions of these systems are also calculated and compared with available literature data. Ċ5H9 radicals can be produced via Ḣ atom addition reactions to the pentadiene isomers and cyclopentene, and also by H-atom abstraction reactions from 1- and 2-pentene and cyclopentane. Comprehensive Ċ5H9 potential energy surface (PES) analyses and high-pressure limiting rate constants for related reactions have been explored in part I of this work ( J. Phys. Chem. A 2019, 123 (22), 9019-9052). In this work, a chemical kinetic model is constructed based on the computed thermochemistry and high-pressure limiting rate constants from part I, to further understand the chemistry of different C5H8 molecules. The most important channels for these addition reactions are discussed in the present work based on reaction pathway analyses. The dominant reaction pathways for these five systems are combined together to generate a simplified Ċ5H9 PES including nine reactants, 25 transition states (TSs), and nine products. Spin-restricted single point energies are calculated for radicals and TSs on the simplified PES at the ROCCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). Temperature- and pressure-dependent rate constants are calculated using RRKM theory with a Master Equation analysis, with restricted energies used for minima on the simplified Ċ5H9 PES and unrestricted energies for other species, over a temperature range of 300-2000 K and in the pressure range 0.01-100 atm. The rate constants calculated are in good agreement with those in the literature. The chemical kinetic model is updated with pressure-dependent rate constants and is used to simulate the species concentration profiles for Ḣ atom addition to cyclopentane and cyclopentene. Through detailed analyses and comparisons, this model can reproduce the experimental measurements of species qualitatively and quantitatively with reasonably good agreement.
Collapse
Affiliation(s)
- Yanjin Sun
- Combustion Chemistry Centre, National University of Ireland, Galway, Ireland
| | - Chong-Wen Zhou
- Combustion Chemistry Centre, National University of Ireland, Galway, Ireland.,School of Energy and Power Engineering, Beihang University, Beijing 100191, P. R. China
| | - Kieran P Somers
- Combustion Chemistry Centre, National University of Ireland, Galway, Ireland
| | - Henry J Curran
- Combustion Chemistry Centre, National University of Ireland, Galway, Ireland
| |
Collapse
|
6
|
Schmalz F, Kopp WA, Kröger LC, Leonhard K. Correcting Rate Constants from Anharmonic Molecular Dynamics for Quantum Effects. ACS OMEGA 2020; 5:2242-2253. [PMID: 32064385 PMCID: PMC7016917 DOI: 10.1021/acsomega.9b03383] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Accepted: 12/30/2019] [Indexed: 06/10/2023]
Abstract
Anharmonicity can greatly affect rate constants. One or even several orders of magnitude of deviation are found for obtaining rate constants using the standard rigid-rotor harmonic-oscillator model. In turn, reactive molecular dynamics (MD) simulations are a powerful way to explore chemical reaction networks and calculate rate constants from the fully anharmonic potential energy surface. However, the classical nature of the dynamics and the required numerical efficiency of the force field limit the accuracy of the resulting kinetics. We combine the best of both worlds by presenting an approximation that pairs anharmonic information intrinsic to classical MD with high-accuracy energies and frequencies from quantum-mechanical electronic structure calculations. The proposed scheme is applied to hydrogen abstractions in the methane system, which allows for the benchmarking of rate constants corrected by our approach against experimental rate constants. This comparison reveals a standard deviation of factor 2.6. Two archetypes of possible failure are identified in the course of a detailed investigation of the CH3 • + H• → CH2 2• + H2 reaction. From this follows the application range of the method, within which the method shows a standard deviation of factor 2.1. The computational efficiency and beneficial scaling of the method allow for application to larger systems, as shown for hydrogen abstraction from 2-butanone by HO2 •.
Collapse
|
7
|
Wu J, Gao LG, Ren W, Truhlar DG. Anharmonic kinetics of the cyclopentane reaction with hydroxyl radical. Chem Sci 2020; 11:2511-2523. [PMID: 34084417 PMCID: PMC8157450 DOI: 10.1039/c9sc05632g] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2019] [Accepted: 01/14/2020] [Indexed: 12/02/2022] Open
Abstract
Cyclopentane is one of the major constituents of transportation fuels, especially jet fuel and diesel, and also is a volatile organic compound with a significant presence in the atmosphere. Hydrogen abstraction from cyclopentane by hydroxyl radical plays a significant role in combustion and atmospheric chemistry. In this work we study the kinetics of this reaction at 200-2000 K using direct dynamics calculations in which the potential energy surface is obtained by quantum mechanical electronic structure calculations. The forward and reverse barrier heights and reaction energies obtained by the CCSD(T)-F12a/jun-cc-pVTZ coupled cluster calculations are used as a benchmark to select an accurate electronic structure method among 36 combinations of exchange-correlation functional and basis set. The selected M06-2X/MG3S method shows the best performance with a mean unsigned deviation from the benchmark of only 0.22 kcal mol-1 for reaction energies and barrier heights. A quadratic-quartic function is adopted to describe the ring bending potential of cyclopentane, and the quartic anharmonicity in the bending mode is treated by a one-dimensional Schrödinger equation using both Wentzel-Kramers-Brillouin (WKB) and Fourier Grid Hamiltonian (FGH) methods. The torsional anharmonicity in the transition state is treated in turn by the free rotor approximation, the one-dimensional hindered rotor approximation, and the multi-structural torsional anharmonicity method. Rate constants of the title reaction are computed by canonical variational transition state theory including tunneling by the multi-dimensional small-curvature tunneling approximation (CVT/SCT). The final rate constants include the quasiharmonic, quadratic-quartic, and torsional anharmonicity. Our calculations are in excellent agreement with all the experimental data available at both combustion and atmospheric temperatures with a deviation of less than 30%.
Collapse
Affiliation(s)
- Junjun Wu
- Department of Mechanical and Automation Engineering, Shenzhen Research Institute, The Chinese University of Hong Kong New Territories Hong Kong SAR China
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota Minneapolis USA
| | - Lu Gem Gao
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota Minneapolis USA
- Center for Combustion Energy, Department of Energy and Power Engineering, Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University Beijing China
| | - Wei Ren
- Department of Mechanical and Automation Engineering, Shenzhen Research Institute, The Chinese University of Hong Kong New Territories Hong Kong SAR China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota Minneapolis USA
| |
Collapse
|
8
|
Jasper AW. Microcanonical Rate Constants for Unimolecular Reactions in the Low-Pressure Limit. J Phys Chem A 2020; 124:1205-1226. [DOI: 10.1021/acs.jpca.9b10693] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ahren W. Jasper
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
9
|
Jasper AW, Harding LB, Knight C, Georgievskii Y. Anharmonic Rovibrational Partition Functions at High Temperatures: Tests of Reduced-Dimensional Models for Systems with up to Three Fluxional Modes. J Phys Chem A 2019; 123:6210-6228. [DOI: 10.1021/acs.jpca.9b03592] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
10
|
Bross DH, Yu HG, Harding LB, Ruscic B. Active Thermochemical Tables: The Partition Function of Hydroxymethyl (CH2OH) Revisited. J Phys Chem A 2019; 123:4212-4231. [DOI: 10.1021/acs.jpca.9b02295] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David H. Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Hua-Gen Yu
- Division of Chemistry, Department of Energy and Photon Sciences, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Lawrence B. Harding
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Consortium for Advanced Science and Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
11
|
|