1
|
Gantzer P, Staub R, Harabuchi Y, Maeda S, Varnek A. Chemography-guided analysis of a reaction path network for ethylene hydrogenation with a model Wilkinson's catalyst. Mol Inform 2024:e202400063. [PMID: 39121023 DOI: 10.1002/minf.202400063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 07/11/2024] [Accepted: 07/19/2024] [Indexed: 08/11/2024]
Abstract
Visualization and analysis of large chemical reaction networks become rather challenging when conventional graph-based approaches are used. As an alternative, we propose to use the chemical cartography ("chemography") approach, describing the data distribution on a 2-dimensional map. Here, the Generative Topographic Mapping (GTM) algorithm - an advanced chemography approach - has been applied to visualize the reaction path network of a simplified Wilkinson's catalyst-catalyzed hydrogenation containing some 105 structures generated with the help of the Artificial Force Induced Reaction (AFIR) method using either Density Functional Theory or Neural Network Potential (NNP) for potential energy surface calculations. Using new atoms permutation invariant 3D descriptors for structure encoding, we've demonstrated that GTM possesses the abilities to cluster structures that share the same 2D representation, to visualize potential energy surface, to provide an insight on the reaction path exploration as a function of time and to compare reaction path networks obtained with different methods of energy assessment.
Collapse
Affiliation(s)
- Philippe Gantzer
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido, 001-0021, Japan
| | - Ruben Staub
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido, 001-0021, Japan
| | - Yu Harabuchi
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido, 001-0021, Japan
| | - Satoshi Maeda
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido, 001-0021, Japan
| | - Alexandre Varnek
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido, 001-0021, Japan
- Laboratory of Chemoinformatics, UMR 7140, CNRS, University of Strasbourg, Strasbourg, 67081, France
| |
Collapse
|
2
|
Pliego JR. Hybrid Cluster-Continuum Method for Single-Ion Solvation Free Energy in Acetonitrile Solvent. J Phys Chem A 2024; 128:6440-6449. [PMID: 39052560 PMCID: PMC11317976 DOI: 10.1021/acs.jpca.4c03593] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2024] [Revised: 07/08/2024] [Accepted: 07/16/2024] [Indexed: 07/27/2024]
Abstract
A new hybrid discrete-continuum approach named the cluster-continuum static approximation (CCSA) has been proposed for acetonitrile solvent. The continuum part uses the conductor-like polarizable continuum model for electrostatic and a surface area-dependent term for nonelectrostatic solvation. The CCSA includes only one explicit acetonitrile solvent molecule and a damping function, which makes the CCSA method reduce to pure continuum solvation in the case of weaker potential of mean force for solute-solvent interaction. The performance of the model was tested for 22 anions and 22 cations, including challenge species that cannot be adequately described by pure continuum solvation. A comparison was done with the widely used solvent model density (SMD) model. For anions, the CCSA reduces to pure continuum solvation and the method has the same performance as the SMD model, with a standard deviation of the mean signed error (SD-MSE) of 2.7 kcal mol-1 for both models. However, the CCSA method for cations considerably outperforms the SMD model, with an SD-MSE of 3.3 kcal mol-1 for the former and 8.4 kcal mol-1 for the latter. The method can be automated, and the present study suggests that continuum solvation models could be parameterized taking into account the explicit solvation as proposed in this work.
Collapse
Affiliation(s)
- Josefredo R. Pliego
- Departamento de Ciências Naturais, Universidade Federal de São João del-Rei, São João del-Rei, Minas Gerais 36301-160, Brazil
| |
Collapse
|
3
|
Piskor T, Pinski P, Mast T, Rybkin V. Multi-Level Protocol for Mechanistic Reaction Studies Using Semi-Local Fitted Potential Energy Surfaces. Int J Mol Sci 2024; 25:8530. [PMID: 39126098 PMCID: PMC11312657 DOI: 10.3390/ijms25158530] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2024] [Revised: 07/18/2024] [Accepted: 07/31/2024] [Indexed: 08/12/2024] Open
Abstract
In this work, we propose a multi-level protocol for routine theoretical studies of chemical reaction mechanisms. The initial reaction paths of our investigated systems are sampled using the Nudged Elastic Band (NEB) method driven by a cheap electronic structure method. Forces recalculated at the more accurate electronic structure theory for a set of points on the path are fitted with a machine learning technique (in our case symmetric gradient domain machine learning or sGDML) to produce a semi-local reactive potential energy surface (PES), embracing reactants, products and transition state (TS) regions. This approach has been successfully applied to a unimolecular (Bergman cyclization of enediyne) and a bimolecular (SN2 substitution) reaction. In particular, we demonstrate that with only 50 to 150 energy-force evaluations with the accurate reference methods (here complete-active-space self-consistent field, CASSCF, and coupled-cluster singles and doubles, CCSD) it is possible to construct a semi-local PES giving qualitative agreement for stationary-point geometries, intrinsic reaction coordinates and barriers. Furthermore, we find a qualitative agreement in vibrational frequencies and reaction rate coefficients. The key aspect of the method's performance is its multi-level nature, which not only saves computational effort but also allows extracting meaningful information along the reaction path, characterized by zero gradients in all but one direction. Agnostic to the nature of the TS and computationally economic, the protocol can be readily automated and routinely used for mechanistic reaction studies.
Collapse
Affiliation(s)
- Tomislav Piskor
- HQS Quantum Simulations GmbH, Rintheimer Straße 23, 76131 Karlsruhe, Germany
- Theoretical Physics, Saarland University, 66123 Saarbrücken, Germany
| | - Peter Pinski
- HQS Quantum Simulations GmbH, Rintheimer Straße 23, 76131 Karlsruhe, Germany
| | - Thilo Mast
- HQS Quantum Simulations GmbH, Rintheimer Straße 23, 76131 Karlsruhe, Germany
| | - Vladimir Rybkin
- HQS Quantum Simulations GmbH, Rintheimer Straße 23, 76131 Karlsruhe, Germany
| |
Collapse
|
4
|
Shaik RA, Jasper AW, Lynch PT, Sivaramakrishnan R, Tranter RS. Initiation and Carbene Induced Radical Chain Reactions in CH 2F 2 Pyrolysis. Chemphyschem 2024; 25:e202400362. [PMID: 38714496 DOI: 10.1002/cphc.202400362] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2024] [Revised: 05/01/2024] [Accepted: 05/07/2024] [Indexed: 05/10/2024]
Abstract
High temperature dissociations of organic molecules typically involve a competition between radical and molecular processes. In this work, we use a modeling, experiment, theory (MET) framework to characterize the high temperature thermal dissociation of CH2F2, a flammable hydrofluorocarbon (HFC) that finds widespread use as a refrigerant. Initiation in CH2F2 proceeds via a molecular elimination channel; CH2F2→CHF+HF. Here we show that the subsequent self-reactions of the singlet carbene, CHF, are fast multichannel processes and a facile source of radicals that initiate rapid chain propagation reactions. These have a marked influence on the decomposition kinetics of CH2F2. The inclusion of these reactions brings the simulations into better agreement with the present and literature experiments. Additionally, flame simulations indicate that inclusion of the CHF+CHF multichannel reaction leads to a noticeable enhancement in predictions of laminar flame speeds, a key parameter that is used to determine the flammability of a refrigerant.
Collapse
Affiliation(s)
- Rizwan A Shaik
- Department of Mechanical & Industrial Engineering, University of Illinois at Chicago, 842 W Taylor Street, Chicago, IL-60607, USA
| | - Ahren W Jasper
- Chemical Dynamics Group, Chemical Sciences & Engineering Division, Argonne National Laboratory, 9700 S Cass Avenue, Lemont, IL-60439, USA
| | - Patrick T Lynch
- Department of Mechanical & Industrial Engineering, University of Illinois at Chicago, 842 W Taylor Street, Chicago, IL-60607, USA
| | - Raghu Sivaramakrishnan
- Chemical Dynamics Group, Chemical Sciences & Engineering Division, Argonne National Laboratory, 9700 S Cass Avenue, Lemont, IL-60439, USA
| | - Robert S Tranter
- Chemical Dynamics Group, Chemical Sciences & Engineering Division, Argonne National Laboratory, 9700 S Cass Avenue, Lemont, IL-60439, USA
| |
Collapse
|
5
|
Guo CM, Zhang H, Cheng XL. Exhaustive State-to-State Cross Sections and Rate Coefficients for Inelastic N 2-N 2 Collisions using QCT Combined with Neural Network Models. J Phys Chem A 2024; 128:5435-5444. [PMID: 38953499 DOI: 10.1021/acs.jpca.4c00590] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Using the quasi-classical trajectory method, we systematically studied the state-to-state vibrational relaxation process of N2(v1) + N2(v2) collisions over a wide temperature range (5000-30,000 K). Different temperature dependencies of the single- and multiquantum VV and VT events in various (v1,v2) collisions are captured, with the dominant channel being related to the initial vibrational energy levels (vmax = 50). At a specified relative translational energy, there is a monotonic relationship of the VT cross sections with the vibrational energy level, particularly in high-energy collisions. Additionally, we constructed well-trained neural network models (R-values reaching 0.99) using limited quasi-classical trajectory (QCT) data sets, which can be used to predict the state-to-state cross sections and rate coefficients of the VV processes N2(v1) + N2(v2) → N2(v1 - Δv) + N2(v2 + Δv) and VT processes N2(v1) + N2(v2) → N2(v1 - Δv) + N2(v2) (Δv = ±1, ±2, ±3) for collisions with arbitrary initial vibrational states. This work not only significantly reduces computational resources but also serves as a reference for the study of the state-to-state dynamics of all four-atom collision systems in hypersonic flows.
Collapse
Affiliation(s)
- Chang-Min Guo
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Hong Zhang
- College of Physics, Sichuan University, Chengdu 610065, China
| | - Xin-Lu Cheng
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| |
Collapse
|
6
|
Jiang W, Chen Y, Li Y. Reactions dynamics for X + H2 insertion reactions (X = C(1D), N(2D), O(1D), S(1D)) with Cayley propagator ring-polymer molecular dynamics. J Chem Phys 2024; 160:234107. [PMID: 38899683 DOI: 10.1063/5.0209143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2024] [Accepted: 05/31/2024] [Indexed: 06/21/2024] Open
Abstract
In this work, rate coefficients of four prototypical insertion reactions, X + H2 → H + XH (X = C(1D), N(2D), O(1D), S(1D)), and associated isotope reactions are calculated based on ring polymer molecular dynamics (RPMD) with Cayley propagator (Cayley-RPMD). The associated kinetic isotope effects are systematically studied too. The Cayley propagator used in this work increases the stability of numerical integration in RPMD calculations and also supports a larger evolution time interval, allowing us to reach both high accuracy and efficiency. So, our results do not only provide chemical kinetic data for the title reactions in an extended temperature range but also consist of experimental results, standard RPMD, and other theoretical methods. The results in this work also reflect that Cayley-RPMD has strong consistency and high reliability in its investigations of chemical dynamics for insertion reactions.
Collapse
Affiliation(s)
- Wenbin Jiang
- Department of Physics, International Center of Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
| | - Yuhao Chen
- Department of Physics, International Center of Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures, Shanghai University, Shanghai 200444, China
- Shanghai Key Laboratory of High Temperature Superconductors, Institute for Quantum Science and Technology, Shanghai University, Shanghai 200444, China
| |
Collapse
|
7
|
Gao Q, Shen C, Zhang H, Long B, Truhlar DG. Quantitative kinetics reveal that reactions of HO 2 are a significant sink for aldehydes in the atmosphere and may initiate the formation of highly oxygenated molecules via autoxidation. Phys Chem Chem Phys 2024; 26:16160-16174. [PMID: 38787752 DOI: 10.1039/d4cp00693c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2024]
Abstract
Large aldehydes are widespread in the atmosphere and their oxidation leads to secondary organic aerosols. The current understanding of their chemical transformation processes is limited to hydroxyl radical (OH) oxidation during daytime and nitrate radical (NO3) oxidation during nighttime. Here, we report quantitative kinetics calculations of the reactions of hexanal (C5H11CHO), pentanal (C4H9CHO), and butanal (C3H7CHO) with hydroperoxyl radical (HO2) at atmospheric temperatures and pressures. We find that neither tunneling nor multistructural torsion anharmonicity should be neglected in computing these rate constants; strong anharmonicity at the transition states is also important. We find rate constants for the three reactions in the range 3.2-7.7 × 10-14 cm3 molecule-1 s-1 at 298 K and 1 atm, showing that the HO2 reactions can be competitive with OH and NO3 oxidation under some conditions relevant to the atmosphere. Our findings reveal that HO2-initiated oxidation of large aldehydes may be responsible for the formation of highly oxygenated molecules via autoxidation. We augment the theoretic studies with laboratory flow-tube experiments using an iodide-adduct time-of-flight chemical ionization mass spectrometer to confirm the theoretical predictions of peroxy radicals and the autoxidation pathway. We find that the adduct from HO2 + C5H11CHO undergoes a fast unimolecular 1,7-hydrogen shift with a rate constant of 0.45 s-1. We suggest that the HO2 reactions make significant contributions to the sink of aldehydes.
Collapse
Affiliation(s)
- Qiao Gao
- School of Physics and Mechatronic Engineering, Guizhou Minzu University, Guiyang 550025, China.
| | - Chuanyang Shen
- Department of Chemistry, University of California, Riverside, California, 92507, USA.
| | - Haofei Zhang
- Department of Chemistry, University of California, Riverside, California, 92507, USA.
| | - Bo Long
- School of Physics and Mechatronic Engineering, Guizhou Minzu University, Guiyang 550025, China.
- College of Materials Science and Engineering, Guizhou Minzu university, Guiyang 550025, China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA.
| |
Collapse
|
8
|
Pham TV. Kinetics and mechanism of the gas-phase reaction of the hydroxyl radical with meta-aminotoluene compound. J Mol Model 2024; 30:169. [PMID: 38750281 DOI: 10.1007/s00894-024-05965-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2024] [Accepted: 05/06/2024] [Indexed: 06/07/2024]
Abstract
CONTEXT An ab initio investigation into the potential energy landscape of the meta-aminotoluene + •OH reaction has been conducted in this study. The calculated results reveal that the reaction channel leading to the product (NHC6H4CH3 + H2O) prevails under the 300-1700 K temperature range, while the reaction path forming the product (NH2C6H4CH2 + H2O) dominates in the higher-temperature region (T ≥ 1800 K). Within the specified temperature range, the product branching ratio for the former declines from 48 to 30%, while the latter shows an increase, reaching 29%. The overall second-order rate constants of the titled reaction obtained at the pressure 760 Torr (N2) can be illustrated by the modified Arrhenius expression of ktotal = 1.46 × 10-13 T0.58 exp[(-0.759 kcal.mol-1)/RT] cm3 molecule-1 s-1 and ktotal = 1.86 × 10-22 T3.24 exp[(-5.086 kcal.mol-1)/RT] cm3 molecule-1 s-1, covering the temperature range of T = 300-600 K and T > 600 K, respectively. The total rate constant at the ambient conditions in this work, 1.43 × 10-11 cm3 molecule-1 s-1, has been found to be roughly one order of magnitude lower than the available experimental data, ~ 1.2 × 10-10 cm3 molecule-1 s-1, measured by Atkinson et al., Rinke et al., and Witte et al., or the theoretical value, 4.4 × 10-10 cm3 molecule-1 s-1, and calculated by Abdel-Rahman and co-workers for the aniline + •OH reaction. METHODS The structures of reactants, transition states, intermediate states, and products of the meta-aminotoluene + •OH reaction are calculated with the aug-cc-pVTZ basis set and the methods DFT/B3LYP and CCSD(T). The rate constants and branching ratios in the 300-2000 K temperature range are calculated with the statistical theoretical TST and RRKM master equation computations including tunneling corrections, with potential energy surface constructed by the CCSD(T)//B3LYP/aug-cc-pVTZ approach.
Collapse
Affiliation(s)
- Tien V Pham
- School of Chemistry and Life Sciences, Hanoi University of Science and Technology, Hanoi, Vietnam.
| |
Collapse
|
9
|
Tang Y, Yang Y, Cao Q, Wang K, Ma D, Liu B, Zhu Q, Ma J. Theoretical Study of the Temperature- and Pressure-Dependent Rate Constants for Nine Reactions between CO n ( n = 0-4), O m ( m = 1-3), C 2O, and C 3O 2 during the Radiolysis of Carbon Dioxide. J Phys Chem A 2024. [PMID: 38691833 DOI: 10.1021/acs.jpca.4c00681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2024]
Abstract
We investigate the reaction pathways of nine important CO2-related reactions using the revDSD-PBEP86-D3(BJ)/jun-cc-pV(T+d)Z level and simultaneously employ an accurate composite method (jun-Cheap) based on coupled-cluster (CC) theory. Subsequently, the Rice-Ramsperger-Kassel-Marcus/master equation (RRKM/ME) is solved to calculate the temperature- and pressure-dependent rate constants. This work investigates reactions involving transition states that have been overlooked in previous literature, including the dissociation of singlet-state C3O2, the triple channel formation of C2O + CO to form C3O2, and the formation of O3 + CO. The results show that CO3 is highly prone to dissociation at high temperatures. Finally, the kinetic data show that over a wide temperature range, our calculations are consistent with previous experimental measurements. The majority of the reaction rate constants studied exhibit significant pressure dependence, while the O3 + CO reaction is pressure-independent at low temperatures. These results are instrumental in the development of detailed kinetic models for the CO2 radiolysis reaction network.
Collapse
Affiliation(s)
- Yuhong Tang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Yu Yang
- Reactor Operation and Application Research Sub-Institute, Nuclear Power Institute of China, Chengdu 610041, China
| | - Qi Cao
- Reactor Operation and Application Research Sub-Institute, Nuclear Power Institute of China, Chengdu 610041, China
| | - Keli Wang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Dandan Ma
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Beibei Liu
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Quan Zhu
- School of Chemical Engineering, Sichuan University, Chengdu 610065, PR China
- Engineering Research Center of Combustion and Cooling for Aerospace Power, Ministry of Education, Sichuan University, Chengdu 610065, PR China
| | - Jianyi Ma
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
- Engineering Research Center of Combustion and Cooling for Aerospace Power, Ministry of Education, Sichuan University, Chengdu 610065, PR China
| |
Collapse
|
10
|
Harabuchi Y, Yokoyama T, Matsuoka W, Oki T, Iwata S, Maeda S. Differentiating the Yield of Chemical Reactions Using Parameters in First-Order Kinetic Equations to Identify Elementary Steps That Control the Reactivity from Complicated Reaction Path Networks. J Phys Chem A 2024; 128:2883-2890. [PMID: 38564273 DOI: 10.1021/acs.jpca.4c00204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
The yield of a chemical reaction is obtained by solving its rate equation. This study introduces an approach for differentiating yields by utilizing the parameters of the rate equation, which is expressed as a first-order linear differential equation. The yield derivative for a specific pair of reactants and products is derived by mathematically expressing the rate constant matrix contraction method, which is a simple kinetic analysis method. The parameters of the rate equation are the Gibbs energies of the intermediates and transition states in the reaction path network used to formulate the rate equation. Thus, our approach for differentiating the yield allows a numerical evaluation of the contribution of energy variation to the yield for each intermediate and transition state in the reaction path network. In other words, a comparison of these values automatically extracts the factors affecting the yield from a complicated reaction path network consisting of numerous reaction paths and intermediates. This study verifies the behavior of the proposed approach through numerical experiments on the reaction path networks of a model system and the Rh-catalyzed hydroformylation reaction. Moreover, the possibility of using this approach for designing ligands in organometallic catalysts is discussed.
Collapse
Affiliation(s)
- Yu Harabuchi
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21, Nishi 10, Kita-ku, Sapporo, Hokkaido 001-0021, Japan
- JST, ERATO Maeda Artificial Intelligence in Chemical Reaction Design and Discovery Project, Kita 10, Nishi 8, Kita-ku, Sapporo, Hokkaido 060-0810, Japan
| | - Tomohiko Yokoyama
- Department of Mathematical Informatics, Graduate School of Information Science and Technology, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Wataru Matsuoka
- JST, ERATO Maeda Artificial Intelligence in Chemical Reaction Design and Discovery Project, Kita 10, Nishi 8, Kita-ku, Sapporo, Hokkaido 060-0810, Japan
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10, Nishi 8, Kita-ku, Sapporo, Hokkaido 060-0810, Japan
| | - Taihei Oki
- JST, ERATO Maeda Artificial Intelligence in Chemical Reaction Design and Discovery Project, Kita 10, Nishi 8, Kita-ku, Sapporo, Hokkaido 060-0810, Japan
- Department of Mathematical Informatics, Graduate School of Information Science and Technology, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Satoru Iwata
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21, Nishi 10, Kita-ku, Sapporo, Hokkaido 001-0021, Japan
- JST, ERATO Maeda Artificial Intelligence in Chemical Reaction Design and Discovery Project, Kita 10, Nishi 8, Kita-ku, Sapporo, Hokkaido 060-0810, Japan
- Department of Mathematical Informatics, Graduate School of Information Science and Technology, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Satoshi Maeda
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21, Nishi 10, Kita-ku, Sapporo, Hokkaido 001-0021, Japan
- JST, ERATO Maeda Artificial Intelligence in Chemical Reaction Design and Discovery Project, Kita 10, Nishi 8, Kita-ku, Sapporo, Hokkaido 060-0810, Japan
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10, Nishi 8, Kita-ku, Sapporo, Hokkaido 060-0810, Japan
- Research and Services Division of Materials Data and Integrated System (MaDIS), National Institute for Materials Science (NIMS), 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| |
Collapse
|
11
|
Rusdipoetra RA, Suwito H, Puspaningsih NNT, Haq KU. Theoretical insight of reactive oxygen species scavenging mechanism in lignin waste depolymerization products. RSC Adv 2024; 14:6310-6323. [PMID: 38380240 PMCID: PMC10877321 DOI: 10.1039/d3ra08346b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Accepted: 01/31/2024] [Indexed: 02/22/2024] Open
Abstract
Apart from natural products and synthesis, phenolic compounds can be produced from the depolymerization of lignin, a major waste in biofuel and paper production. This process yields a plethora of aryl propanoid phenolic derivatives with broad biological activities, especially antioxidant properties. Due to its versatility, our study focuses on investigating the antioxidant mechanisms of several phenolic compounds obtained from renewable and abundant resources, namely, syringol (Hs), 4-allylsyringol (HAs), 4-propenylsyringol (HPns), and 4-propylsyringol (HPs). Employing the density functional theory (DFT) approach in conjunction with the QM-ORSA protocol, we aim to explore the reactivity of these compounds in neutralizing hydroperoxyl radicals in physiological and non-polar media. Kinetic and thermodynamic parameter calculations on the antioxidant activity of these compounds were also included in this study. Additionally, our research utilizes the activation strain model (ASM) for the first time to explain the reactivity of the HT and RAF mechanisms in the peroxyl radical scavenging process. It is predicted that HPs has the best rate constant in both media (1.13 × 108 M-1 s-1 and 1.75 × 108 M-1 s-1, respectively). Through ASM analysis, it is observed that the increase in the interaction energy due to the formation of intermolecular hydrogen bonds during the reaction is an important feature for accelerating the hydrogen transfer process. Furthermore, by examining the physicochemical and toxicity parameters, only Hs is not suitable for further investigation as a therapeutic agent because of potential toxicity and mutagenicity. However, overall, all compounds are considered potent HOO˙ scavengers in lipid-rich environments compared to previously studied antioxidants.
Collapse
Affiliation(s)
- Rahmanto Aryabraga Rusdipoetra
- Bioinformatic Research Group, Research Centre of Bio-Molecule Engineering (BIOME), Airlangga University Jl. Ir. H. Soekarno Mulyorejo Surabaya Indonesia
- Department of Chemistry, Faculty of Science and Technology, Airlangga University Jl. Ir. H. Soekarno Mulyorejo Surabaya Indonesia
| | - Hery Suwito
- Department of Chemistry, Faculty of Science and Technology, Airlangga University Jl. Ir. H. Soekarno Mulyorejo Surabaya Indonesia
| | - Ni Nyoman Tri Puspaningsih
- Department of Chemistry, Faculty of Science and Technology, Airlangga University Jl. Ir. H. Soekarno Mulyorejo Surabaya Indonesia
- Proteomic Research Group, Research Centre of Bio-Molecule Engineering (BIOME), Airlangga University Jl. Ir. H. Soekarno Mulyorejo Surabaya Indonesia
| | - Kautsar Ul Haq
- Bioinformatic Research Group, Research Centre of Bio-Molecule Engineering (BIOME), Airlangga University Jl. Ir. H. Soekarno Mulyorejo Surabaya Indonesia
- Department of Chemistry, Faculty of Science and Technology, Airlangga University Jl. Ir. H. Soekarno Mulyorejo Surabaya Indonesia
| |
Collapse
|
12
|
Okada H, Maeda S. On Accelerating Substrate Optimization Using Computational Gibbs Energy Barriers: A Numerical Consideration Utilizing a Computational Data Set. ACS OMEGA 2024; 9:7123-7131. [PMID: 38371820 PMCID: PMC10870292 DOI: 10.1021/acsomega.3c09066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 01/05/2024] [Accepted: 01/16/2024] [Indexed: 02/20/2024]
Abstract
Substrate optimization is a time- and resource-consuming step in organic synthesis. Recent advances in chemo- and materials-informatics provide systematic and efficient procedures utilizing tools such as Bayesian optimization (BO). This study explores the possibility of reducing the required experiments further by utilizing computational Gibbs energy barriers. To thoroughly validate the impact of using computational Gibbs energy barriers in BO-assisted substrate optimization, this study employs a computational Gibbs energy barrier data set in the literature and performs an extensive numerical investigation virtually regarding the Gibbs energy barriers as virtual experimental results and those with systematic and random noises as virtual computational results. The present numerical investigation shows that even the computational reactivity affected by noises of as much as 20 kJ/mol helps reduce the number of required experiments.
Collapse
Affiliation(s)
- Hiroaki Okada
- Graduate
School of Chemical Sciences and Engineering, Hokkaido University, Sapporo, Hokkaido 060-8628, Japan
| | - Satoshi Maeda
- Department
of Chemistry, Graduate School of Science, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
- Institute
for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido 001-0021, Japan
- ERATO
Maeda Artificial Intelligence for Chemical Reaction Design and Discovery
Project, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
- Research
and Services Division of Materials Data and Integrated System (MaDIS), National Institute for Materials Science (NIMS), Tsukuba, Ibaraki 305-0044, Japan
| |
Collapse
|
13
|
Martínez AG, Siehl HU, de la Moya S, Gómez PC. Easy and accurate computation of energy barriers for carbocation solvation: an expeditious tool to face carbocation chemistry. Phys Chem Chem Phys 2023; 25:31012-31019. [PMID: 37938916 DOI: 10.1039/d3cp03544a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2023]
Abstract
An expeditious procedure for the challenging computation of the free energy barriers (ΔG≠) for the solvation of carbocations is presented. This procedure is based on Marcus Theory (MT) and the popular B3LYP/6-31G(d)//PCM method, and it allows the easy, accurate and inexpensive prediction of these barriers for carbocations of very different stability. This method was validated by the fair mean absolute error (ca. 1.5 kcal mol-1) achieved in the prediction of 19 known experimental barriers covering a range of ca. 50 kcal mol-1. Interestingly, the new procedure also uses an original method for the calculation of the required inner reorganization energy (Λi) and free energy of reaction (ΔG). This procedure should pave the way to face computationally the pivotal issue of carbocation chemistry and could be easily extended to any bimolecular organic reaction.
Collapse
Affiliation(s)
- Antonio G Martínez
- Departamento de Química Orgánica, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040, Madrid, Spain.
| | - Hans-Ulrich Siehl
- Abteilung Organische Chemie I, Universität Ulm, Albert Einstein Alee 11, 89069 Ulm, Germany
| | - Santiago de la Moya
- Departamento de Química Orgánica, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040, Madrid, Spain.
| | - Pedro C Gómez
- Departamento de Química Orgánica, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040, Madrid, Spain.
- Departamento de Química Física, Universidad Complutense de Madrid, Facultad de Ciencias Químicas, 28040 Madrid, Spain.
| |
Collapse
|
14
|
Li Y, Zhang RM, Xu X. Theoretical Kinetics studies of isoprene peroxy radical chemistry: The fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical. ECOTOXICOLOGY AND ENVIRONMENTAL SAFETY 2023; 266:115553. [PMID: 37839188 DOI: 10.1016/j.ecoenv.2023.115553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Revised: 09/19/2023] [Accepted: 10/05/2023] [Indexed: 10/17/2023]
Abstract
The OH radical recycling mechanism in isoprene oxidation is one of the most exciting topics in atmospheric chemistry, and the corresponding studies expand our understanding of oxidation mechanisms of volatile organic compounds in the troposphere and provide reliable evidence to improve and develop conventional atmospheric models. In this work, we performed a detailed theoretical kinetics study on the Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry, which is proposed as the heart of OH recycling in isoprene oxidation. With the full consideration of its accumulation and consumption channels, we studied and discussed the fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical by solving the energy-resolved master equation over a broad range of conditions, including not only room temperatures but also high temperatures of a forest fire or low temperatures and pressures of the upper troposphere. We found non-negligible pressure dependence of its fate at combustion temperatures (up to two orders of magnitude) and demonstrated the significance of both the multi-structural torsional anharmonicity and tunneling for accurately calculating kinetics of the studied system. More interestingly, the tunneling effect on the phenomenological rate constants of the H-shift reaction channel is also found to be pressure-dependent due to the competition with the O2 loss reaction. In addition, our time evolution calculations revealed a two-stage behavior of critical species in this reaction system and estimated the shortest half-lives for the Z-δ-(4-OH, 1-OO)-ISOPOO radical at various temperatures, pressures and altitudes. This detailed kinetics study of Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry offers a typical example to deeply understand the core mechanism of OH recycling pathways in isoprene oxidation, and provides valuable insights for promoting the development of relevant atmospheric models.
Collapse
Affiliation(s)
- Yan Li
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Rui Ming Zhang
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Xuefei Xu
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China.
| |
Collapse
|
15
|
Ćehić M, Nikšić-Franjić I. Reactions of a hydrogen atom with haloacetates in aqueous solutions: Computational evidence for proton-coupled electron transfer and competing mechanisms. J Comput Chem 2023; 44:2212-2222. [PMID: 37452597 DOI: 10.1002/jcc.27191] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Revised: 06/30/2023] [Accepted: 07/04/2023] [Indexed: 07/18/2023]
Abstract
A computational study of the mechanisms and kinetics of the aqueous reactions of a hydrogen atom with haloacetates is presented. Several mechanisms in the close competition are observed, such as proton-coupled electron transfer (PCET), hydrogen atom transfer (HAT), and halogen abstraction (XA). Computations predict that dechlorination takes place via PCET mechanisms and not via XA, as stated earlier, while XA is the fastest mechanism forIAc - . The reaction rate constants are reasonably well predicted within the theoretically most reliable canonical variational transition state theory with small curvature tunneling corrections and compared with the experimental ones. To reproduce the experimental rate constants of the debromination process it is necessary to include the PCET and XA cumulative values. Small curvature tunneling corrections to the rate constants are the highest for HAT and PCET mechanisms, up to 70 times larger than the Wigner, while variational effects for XA mechanisms are very small.
Collapse
Affiliation(s)
- Mirsada Ćehić
- Chair of Environmental Protection, University of Applied Sciences on Security and Safety, VSS, Zagreb, Croatia
| | | |
Collapse
|
16
|
Zankov D, Madzhidov T, Baskin I, Varnek A. Conjugated quantitative structure-property relationship models: Prediction of kinetic characteristics linked by the Arrhenius equation. Mol Inform 2023; 42:e2200275. [PMID: 37488968 DOI: 10.1002/minf.202200275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Revised: 07/08/2023] [Accepted: 07/24/2023] [Indexed: 07/26/2023]
Abstract
Conjugated QSPR models for reactions integrate fundamental chemical laws expressed by mathematical equations with machine learning algorithms. Herein we present a methodology for building conjugated QSPR models integrated with the Arrhenius equation. Conjugated QSPR models were used to predict kinetic characteristics of cycloaddition reactions related by the Arrhenius equation: rate constantl o g k ${{\rm l}{\rm o}{\rm g}k}$ , pre-exponential factorl o g A ${{\rm l}{\rm o}{\rm g}A}$ , and activation energyE a ${{E}_{{\rm a}}}$ . They were benchmarked against single-task (individual and equation-based models) and multi-task models. In individual models, all characteristics were modeled separately, while in multi-task modelsl o g k ${{\rm l}{\rm o}{\rm g}k}$ ,l o g A ${{\rm l}{\rm o}{\rm g}A}$ andE a ${{E}_{{\rm a}}}$ were treated cooperatively. An equation-based model assessedl o g k ${{\rm l}{\rm o}{\rm g}k}$ using the Arrhenius equation andl o g A ${{\rm l}{\rm o}{\rm g}A}$ andE a ${{E}_{{\rm a}}}$ values predicted by individual models. It has been demonstrated that the conjugated QSPR models can accurately predict the reaction rate constants at extreme temperatures, at which reaction rate constants hardly can be measured experimentally. Also, in the case of small training sets conjugated models are more robust than related single-task approaches.
Collapse
Affiliation(s)
- Dmitry Zankov
- Laboratory of Chemoinformatics, University of Strasbourg, France
| | - Timur Madzhidov
- Chemistry Solutions, Elsevier Ltd, Oxford, OX5 1GB, United Kingdom
| | - Igor Baskin
- Department of Materials Science and Engineering, Technion - Israel Institute of Technology, Israel
| | - Alexandre Varnek
- Laboratory of Chemoinformatics, University of Strasbourg, France
| |
Collapse
|
17
|
Ballotta B, Martínez-Núñez E, Rampino S, Barone V. New prebiotic molecules in the interstellar medium from the reaction between vinyl alcohol and CN radicals: unsupervised reaction mechanism discovery, accurate electronic structure calculations and kinetic simulations. Phys Chem Chem Phys 2023; 25:22840-22850. [PMID: 37584420 DOI: 10.1039/d3cp02571c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/17/2023]
Abstract
Vinyl alcohol (VyA) and cyanide (CN) radicals are relatively abundant in the interstellar medium (ISM). VyA is the enolic tautomer of acetaldehyde and has two low-lying conformers, characterized by the syn or anti placement of hydroxyl hydrogen with respect to the double bond. In this paper, we present a gas-phase model of the barrierless reactions of both VyA's conformers with CN employing accurate quantum chemical computations in the framework of a master equation approach based on the transition state theory. Our results indicate that both VyA conformers feature a similar reactivity with CN, starting with a barrierless addition to the double bond and followed by different isomerization, dissociation, and/or hydrogen elimination steps. The rate constants computed for temperatures up to 600 K show that several reaction channels are open even under the harsh conditions of the ISM, with the favoured one providing the first feasible formation route of a prebiotic molecule not yet detected in the ISM, namely cyanoacetaldehyde. This finding suggests looking for cyanoacetaldehyde in regions where both VyA and CN have already been detected, like, e.g., Sagittarius B2N or G+0.693-0.027.
Collapse
Affiliation(s)
- Bernardo Ballotta
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126, Pisa, Italy.
| | - Emilio Martínez-Núñez
- Departamento de Química Física, Facultade de Química, Campus Vida, Universidade de Santiago de Compostela, Avenida das Ciencias s/n, 15782, Santiago de Compostela, Spain
| | - Sergio Rampino
- Università degli Studi di Padova, Dipartimento di Scienze Chimiche, Via Marzolo 1, 35131 Padova, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126, Pisa, Italy.
| |
Collapse
|
18
|
Steffen J. Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package. J Chem Theory Comput 2023; 19:5334-5355. [PMID: 37555628 DOI: 10.1021/acs.jctc.3c00568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2023]
Abstract
A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.
Collapse
Affiliation(s)
- Julien Steffen
- Chair of Theoretical Chemistry, Friedrich-Alexander University Erlangen-Nürnberg, Egerlandstraße 3, 91058 Erlangen, Bavaria, Germany
| |
Collapse
|
19
|
Ahn W, Triana JF, Recabal F, Herrera F, Simpkins BS. Modification of ground-state chemical reactivity via light-matter coherence in infrared cavities. Science 2023; 380:1165-1168. [PMID: 37319215 DOI: 10.1126/science.ade7147] [Citation(s) in RCA: 34] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Accepted: 05/12/2023] [Indexed: 06/17/2023]
Abstract
Reaction-rate modifications for chemical processes due to strong coupling between reactant molecular vibrations and the cavity vacuum have been reported; however, no currently accepted mechanisms explain these observations. In this work, reaction-rate constants were extracted from evolving cavity transmission spectra, revealing resonant suppression of the intracavity reaction rate for alcoholysis of phenyl isocyanate with cyclohexanol. We observed up to an 80% suppression of the rate by tuning cavity modes to be resonant with the reactant isocyanate (NCO) stretch, the product carbonyl (CO) stretch, and cooperative reactant-solvent modes (CH). These results were interpreted using an open quantum system model that predicted resonant modifications of the vibrational distribution of reactants from canonical statistics as a result of light-matter quantum coherences, suggesting links to explore between chemistry and quantum science.
Collapse
Affiliation(s)
- Wonmi Ahn
- UNAM - National Nanotechnology Research Center and Institute of Materials Science and Nanotechnology, Bilkent University, Ankara, Turkey
| | - Johan F Triana
- Department of Physics, Universidad de Santiago de Chile, Santiago, Chile
| | - Felipe Recabal
- Department of Physics, Universidad de Santiago de Chile, Santiago, Chile
| | - Felipe Herrera
- Department of Physics, Universidad de Santiago de Chile, Santiago, Chile
- Millennium Institute for Research in Optics (MIRO), Concepción, Chile
| | - Blake S Simpkins
- Chemistry Division, US Naval Research Laboratory, Washington, DC, USA
| |
Collapse
|
20
|
Douroudgari H, Zarepour H, Vahedpour M, Jaberi M, Zarepour M. The atmospheric relevance of primary alcohols and imidogen reactions. Sci Rep 2023; 13:9150. [PMID: 37277419 DOI: 10.1038/s41598-023-35473-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2022] [Accepted: 05/18/2023] [Indexed: 06/07/2023] Open
Abstract
Organic alcohols as very volatile compounds play a crucial role in the air quality of the atmosphere. So, the removal processes of such compounds are an important atmospheric challenge. The main goal of this research is to discover the atmospheric relevance of degradation paths of linear alcohols by imidogen with the aid of simulation by quantum mechanical (QM) methods. To this end, we combine broad mechanistic and kinetic results to get more accurate information and to have a deeper insight into the behavior of the designed reactions. Thus, the main and necessary reaction pathways are explored by well-behaved QM methods for complete elucidation of the studying gaseous reactions. Moreover, the potential energy surfaces as a main factor are computed for easier judging of the most probable pathways in the simulated reactions. Our attempt to find the occurrence of the considered reactions in the atmospheric conditions is completed by precisely evaluating the rate constants of all elementary reactions. All of the computed bimolecular rate constants have a positive dependency on both temperature and pressure. The kinetic results show that H-abstraction from the α carbon is dominant relative to the other sites. Finally, by the results of this study, we conclude that at moderate temperatures and pressures primary alcohols can degrade with imidogen, so they can get atmospheric relevance.
Collapse
Affiliation(s)
- Hamed Douroudgari
- Department of Chemistry, University of Zanjan, PO Box 38791-45371, Zanjan, Iran.
| | - Hadi Zarepour
- Department of Chemistry, University of Zanjan, PO Box 38791-45371, Zanjan, Iran
| | - Morteza Vahedpour
- Department of Chemistry, University of Zanjan, PO Box 38791-45371, Zanjan, Iran.
| | - Mahdi Jaberi
- Department of Chemistry, University of Zanjan, PO Box 38791-45371, Zanjan, Iran
| | - Mahdi Zarepour
- Department of Chemistry, University of Zanjan, PO Box 38791-45371, Zanjan, Iran
| |
Collapse
|
21
|
Giese TJ, York DM. Estimation of frequency factors for the calculation of kinetic isotope effects from classical and path integral free energy simulations. J Chem Phys 2023; 158:174105. [PMID: 37125722 PMCID: PMC10154067 DOI: 10.1063/5.0147218] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Accepted: 04/17/2023] [Indexed: 05/02/2023] Open
Abstract
We use the modified Bigeleisen-Mayer equation to compute kinetic isotope effect values for non-enzymatic phosphoryl transfer reactions from classical and path integral molecular dynamics umbrella sampling. The modified form of the Bigeleisen-Mayer equation consists of a ratio of imaginary mode vibrational frequencies and a contribution arising from the isotopic substitution's effect on the activation free energy, which can be computed from path integral simulation. In the present study, we describe a practical method for estimating the frequency ratio correction directly from umbrella sampling in a manner that does not require normal mode analysis of many geometry optimized structures. Instead, the method relates the frequency ratio to the change in the mass weighted coordinate representation of the minimum free energy path at the transition state induced by isotopic substitution. The method is applied to the calculation of 16/18O and 32/34S primary kinetic isotope effect values for six non-enzymatic phosphoryl transfer reactions. We demonstrate that the results are consistent with the analysis of geometry optimized transition state ensembles using the traditional Bigeleisen-Mayer equation. The method thus presents a new practical tool to enable facile calculation of kinetic isotope effect values for complex chemical reactions in the condensed phase.
Collapse
Affiliation(s)
- Timothy J. Giese
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| |
Collapse
|
22
|
Alibakhshi A, Steffen J, Pinilla C, Hartke B. Comparison of Implicit and Explicit Solvent Approaches in Ab Initio Evaluation of Thermochemistry in Solution: Application in Studying Boron Isotope Fractionation in Water. J Phys Chem A 2023; 127:2503-2510. [PMID: 36917555 DOI: 10.1021/acs.jpca.3c00154] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
Evaluation of thermochemistry in solution plays a key role in numerous fields. For this task, the solvent effects are commonly included in theoretical computations based on either implicit or explicit solvent approaches. In the present study, we evaluate and compare the performance of some of the most widely applied methods based on these two approaches. For studying the solvent effect on thermochemistry with an explicit solvent, we demonstrate that partial normal mode analysis with frozen geometry of solvent molecules for multiple solute-solvent configurations can yield quite accurate and reliable results for a drastically reduced computational cost. As a case study, we consider the evaluation of the equilibrium constant for the boron isotope exchange between boric acid and borate (k3-4) in pure and saline water which is of high geochemical importance. Employing three different rigorous and high-precision theoretical approaches, we provide a reliable estimation of k3-4 which is a value within 1.028 to 1.030 for both pure and saline water which is in excellent agreement with experimental data.
Collapse
Affiliation(s)
- Amin Alibakhshi
- Theoretical Chemistry, Institute for Physical Chemistry, Christian-Albrechts-University, Olshausenstr. 40, 24118 Kiel, Germany
- Theoretical Chemistry, Ruhr-Universitaet Bochum, D-44780 Bochum, Germany
| | - Julien Steffen
- Theoretical Chemistry, Institute for Physical Chemistry, Christian-Albrechts-University, Olshausenstr. 40, 24118 Kiel, Germany
| | - Carlos Pinilla
- Departamento de Física y Geociencias, Universidad del Norte, Km 5 via Puerto Colombia, Barranquilla 080020, Colombia
- School of Chemistry, University of Bristol, Cantock's Close Road, BS8 1TS Bristol, U.K
| | - Bernd Hartke
- Theoretical Chemistry, Institute for Physical Chemistry, Christian-Albrechts-University, Olshausenstr. 40, 24118 Kiel, Germany
| |
Collapse
|
23
|
Su Y, Ren Q, Zhang WY, Chen F. Computational Studies on the Reactivity of Polycyclic Aromatic Hydrocarbons. Chemphyschem 2023; 24:e202200638. [PMID: 36409286 DOI: 10.1002/cphc.202200638] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2022] [Revised: 11/15/2022] [Accepted: 11/18/2022] [Indexed: 11/23/2022]
Abstract
Polycyclic aromatic hydrocarbons (PAHs) are widely present in the environment as toxic pollutants. In this study, quantum chemistry methods are used to study reactions of PAHs in both particle and gas phases. Seven theoretical methods are exploited to predict the reactive sites of 15 PAHs in the particle phase. Among these methods, the performance of the condensed Fukui function (CFF) is optimum. The gas-phase reactions of eight PAHs are also investigated. Except for fluorene, CFF predicts correctly the gas-phase mono-nitro products for seven systems. The products of fluorene predicted by CFF are 1-nitrofluorene and 3-nitrofluorene, which is however inconsistent with the experimental results. Transition state theory is then used to investigate the reaction mechanism of fluorene. Calculated rate constants for 3-nitrofluorene and 2-nitrofluorene formation are much bigger than that for 1-nitrofluorene formation, which is in agreement with the experimental results.
Collapse
Affiliation(s)
- Yingwei Su
- Department of Chemistry and Chemical Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing Key Laboratory for Science and Application of Functional Molecular and Crystalline Materials, Beijing, 100083, China
| | - Qing Ren
- Department of Chemistry and Chemical Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing Key Laboratory for Science and Application of Functional Molecular and Crystalline Materials, Beijing, 100083, China
| | - Wen-Yan Zhang
- Department of Chemistry and Chemical Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing Key Laboratory for Science and Application of Functional Molecular and Crystalline Materials, Beijing, 100083, China
| | - Feiwu Chen
- Department of Chemistry and Chemical Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing Key Laboratory for Science and Application of Functional Molecular and Crystalline Materials, Beijing, 100083, China
| |
Collapse
|
24
|
Fatahiyan L, Taghva Manesh A, Masnabadi N. Tautomerization reactions of thiobarbituric acid: A detailed kinetic study using combined canonical variational transition state theory and QTAIM approach. INT J CHEM KINET 2023. [DOI: 10.1002/kin.21632] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
Affiliation(s)
- Leila Fatahiyan
- Department of Chemistry, Central Tehran Branch, Islamic Azad University Tehran Iran
| | - Afshin Taghva Manesh
- Department of Chemistry, Central Tehran Branch, Islamic Azad University Tehran Iran
| | - Nasrin Masnabadi
- Department of Chemistry, Roudehen Branch, Islamic Azad University Roudehen Iran
| |
Collapse
|
25
|
Kee CW. Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities. Molecules 2023; 28:molecules28041715. [PMID: 36838703 PMCID: PMC9966076 DOI: 10.3390/molecules28041715] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2022] [Revised: 02/03/2023] [Accepted: 02/05/2023] [Indexed: 02/25/2023] Open
Abstract
Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.
Collapse
Affiliation(s)
- Choon Wee Kee
- Institute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR), 1 Pesek Road, Jurong Island, Singapore 627833, Republic of Singapore
| |
Collapse
|
26
|
Braunstein M, Bonnet L. An efficient algorithm for capturing quantum effects in classical reactive scattering: application to D + H+3 → H 2D + + H. Phys Chem Chem Phys 2023; 25:1602-1605. [PMID: 36541279 DOI: 10.1039/d2cp05108g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Motivated by a recent semiclassical analysis of chemical reaction thresholds [Bonnet et al., J. Chem. Phys., 2022, 157, 094114], we present an efficient algorithm for including zero-point energy (ZPE) effects in classical reactive scattering. The algorithm is an extension of the quasi-classical trajectory (QCT) Gaussian binning method. We apply it to the astrophysically important D + H+3 reaction, where there are significant quantum effects and where application of other methods is problematic [Braunstein et al., Phys. Chem. Chem. Phys., 2022, 24, 5489]. The rate constants computed with the new, general algorithm closely match recent Ring Polymer Molecular Dynamics (RPMD) [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] and experimentally derived [Bowen et al., J. Chem. Phys., 2021, 154, 084307] ones spanning ∼4 orders of magnitude from 70 to 1500 K.
Collapse
Affiliation(s)
- Matthew Braunstein
- Spectral Sciences Incorporated, 4 Fourth Avenue, Burlington, MA 01824, USA.
| | - Laurent Bonnet
- CNRS, Université de Bordeaux, ISM, UMR 5255, F-33400 Talence, France.
| |
Collapse
|
27
|
Chen J. Why Should the Reaction Order of a Bimolecular Reaction be 2.33 Instead of 2? J Phys Chem A 2022; 126:9719-9725. [PMID: 36520427 PMCID: PMC9805503 DOI: 10.1021/acs.jpca.2c07500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Predicting the reaction kinetics, that is, how fast a reaction can happen in a solution, is essential information for many processes, such as industrial chemical manufacturing, refining, synthesis and separation of petroleum products, environmental processes in air and water, biological reactions in cells, biosensing, and drug delivery. Collision theory was originally developed to explain the reaction kinetics of gas reactions with no dilution. For a reaction in a diluted inert gas solution or a diluted liquid solution, diffusion often dominates the collision process. Thus, it is necessary to include diffusion in such a calculation. Traditionally, the classical Smoluchowski rate is used as a starting point to predict the collision frequency of two molecules in a diluted solution. In this report, a different collision model is derived from the adsorption of molecules on a flat surface. A surprising result is obtained, showing that the reaction order for bimolecular reactions should be 2 and 1/3 instead of 2, following a fractal reaction kinetics.
Collapse
Affiliation(s)
- Jixin Chen
- Department of Chemistry and Biochemistry, Nanoscale and Quantum Phenomena Institute, Ohio University, Athens, OH 45701, USA
| |
Collapse
|
28
|
Couch DE, Mulvihill CR, Sivaramakrishnan R, Au K, Taatjes CA, Sheps L. Quantification of Key Peroxy and Hydroperoxide Intermediates in the Low-Temperature Oxidation of Dimethyl Ether. J Phys Chem A 2022; 126:9497-9509. [DOI: 10.1021/acs.jpca.2c06959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- David E. Couch
- Combustion Research Facility, Sandia National Laboratories, Livermore, California94551, United States
| | - Clayton R. Mulvihill
- Chemical Sciences and Engineering, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - Raghu Sivaramakrishnan
- Chemical Sciences and Engineering, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - Kendrew Au
- Combustion Research Facility, Sandia National Laboratories, Livermore, California94551, United States
| | - Craig A. Taatjes
- Combustion Research Facility, Sandia National Laboratories, Livermore, California94551, United States
| | - Leonid Sheps
- Combustion Research Facility, Sandia National Laboratories, Livermore, California94551, United States
| |
Collapse
|
29
|
Chi Y, Meng Q, He C, Zhang P. Metric-Based Assessment Method for MS-T Formalism with Small Subsets of Torsional Conformers. J Phys Chem A 2022; 126:8305-8314. [DOI: 10.1021/acs.jpca.2c04724] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yicheng Chi
- Department of Mechanical Engineering, The Hong Kong Polytechnic University, Hung Hom, Kowloon999077, Hong Kong
| | - Qinghui Meng
- Department of Mechanical Engineering, The Hong Kong Polytechnic University, Hung Hom, Kowloon999077, Hong Kong
- National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei230052, China
| | - Chengming He
- Wide Range Flight Engineering Science and Application Center, Institute of Mechanics, Chinese Academy of Sciences, Beijing100864, China
| | - Peng Zhang
- Department of Mechanical Engineering, The Hong Kong Polytechnic University, Hung Hom, Kowloon999077, Hong Kong
- Department of Mechanical Engineering, City University of Hong Kong, Kowloon Tong, Kowloon999077, Hong Kong
| |
Collapse
|
30
|
Long B, Xia Y, Truhlar DG. Quantitative Kinetics of HO 2 Reactions with Aldehydes in the Atmosphere: High-Order Dynamic Correlation, Anharmonicity, and Falloff Effects Are All Important. J Am Chem Soc 2022; 144:19910-19920. [PMID: 36264240 DOI: 10.1021/jacs.2c07994] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Kinetics provides the fundamental parameters for elucidating sources and sinks of key atmospheric species and for atmospheric modeling more generally. Obtaining quantitative kinetics in the laboratory for the full range of atmospheric temperatures and pressures is quite difficult. Here, we use computational chemistry to obtain quantitative rate constants for the reactions of HO2 with HCHO, CH3CHO, and CF3CHO. First, we calculate the high-pressure-limit rate constants by using a dual-level strategy that combines conventional transition state theory using a high level of electronic structure wave function theory with canonical variational transition state theory including small-curvature tunneling using density functional theory. The wave-function level is beyond-CCSD(T) for HCHO and CCSD(T)-F12a (Level-A) for XCHO (X = CH3, CF3), and the density functional (Level-B) is specifically validated for these reactions. Then, we calculate the pressure-dependent rate constants by using system-specific quantum RRK theory (SS-QRRK) and also by an energy-grained master equation. The two treatments of the pressure dependence agree well. We find that the Level-A//Level-B method gives good agreement with CCSDTQ(P)/CBS. We also find that anharmonicity is an important factor that increases the rate constants of all three reactions. We find that the HO2 + HCHO reaction has a significant dependence on pressure, but the HO2 + CF3CHO reaction is almost independent of pressure. Our findings show that the HO2 + HCHO reaction makes important contribution to the sink for HCHO, and the HO2 + CF3CHO reaction is the dominant sink for CF3CHO in the atmosphere.
Collapse
Affiliation(s)
- Bo Long
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Yu Xia
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
31
|
Bonnet L, Crespos C, Monnerville M. Chemical reaction thresholds according to classical-limit quantum dynamics. J Chem Phys 2022; 157:094114. [PMID: 36075739 DOI: 10.1063/5.0101311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Classical-limit quantum dynamics is used to explain the origin of the quantum thresholds of chemical reactions from their classical dynamics when these are vibrationally nonadiabatic across the interaction region. This study is performed within the framework of an elementary model of chemical reaction that mimics the passage from the free rotation of the reagents to the bending vibration at the transition state to the free rotation of the products.
Collapse
Affiliation(s)
- L Bonnet
- Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33400 Talence, France
| | - C Crespos
- Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33400 Talence, France
| | - M Monnerville
- Univ. Lille, CNRS, PhLAM, UMR 8523, 59655 Villeneuve d'Ascq, France
| |
Collapse
|
32
|
Qin J, Liu Y, Li J. Quantitative Dynamics of Paradigmatic SN2 reaction OH− + CH3F on Accurate Full-Dimensional Potential Energy Surface. J Chem Phys 2022; 157:124301. [DOI: 10.1063/5.0112228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The bimolecular reaction between OH− and CH3F is not just a prototypical SN2 process but also has three other product channels. Here, we develop an accurate full-dimensional potential energy surface (PES) based on 191 193 points calculated at the level CCSD(T)-F12a/aug-cc-pVTZ. A detailed dynamics and mechanism analysis were carried out on this PES by using the quasi-classical trajectory approach. It is verified that the trajectories do not follow the minimum energy path (MEP) but directly dissociate to F− and CH3OH. In addition, a new transition state for proton exchange and a new product complex CH2F−‧‧‧H2O for proton abstraction were discovered. The trajectories avoid the transition state or this complex, instead dissociate to H2O and CH2F− directly through the ridge regions of the MEP before the transition state. These non-MEP dynamics become more pronounced at high collision energies. Detailed dynamics simulations provide new insights into the atomic-level mechanisms of the title reaction thanks to the new chemically accurate PES with the aid of the machine learning.
Collapse
Affiliation(s)
- Jie Qin
- Chemistry and Chemical Engineering, Chongqing University Department of Chemical Engineering, China
| | | | - Jun Li
- School of Chemistry and Chemical Engineering, Chongqing University, China
| |
Collapse
|
33
|
Ventura E, Bezerra MG, Leitão EFV, do Monte SA. Acid Enriched with Methanol: Formation of a Prebiotic Cluster in the Interstellar Medium. Chemphyschem 2022; 23:e202200403. [PMID: 35962978 DOI: 10.1002/cphc.202200403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2022] [Revised: 08/12/2022] [Indexed: 11/12/2022]
Abstract
Organic molecules are a potential source of prebiotic chemistry in the interstellar medium (ISM). Methanol (MetOH) is a very important source of more complex molecules. H 3 O + (aq) and Cl - (aq) are fundamental to living organisms and can be generated in the ISM from the dissociation of HCl with just four water molecules, yielding the (H 3 O) + (H 2 O) 3 Cl - ion-pair. Here, a detailed mechanism, based on density functional theory (DFT) and ab-initio (2 nd order Mfller-Plesset perturbation theory, MP2) calculations, is suggested for the substitution reactions of these water molecules by MetOH. The time required for formation of an appreciable amount of the product ((H 3 O) + (MetOH) 3 Cl - ) can be only few years. Such reaction can take place in Sagittarius B2, where HCl, H 2 O and MetOH have already been identified and it can be an important source for the formation of more complex prebiotic structures.
Collapse
Affiliation(s)
- Elizete Ventura
- Universidade Federal da Paraiba, Departamento de Química, Centro de Ciências Exatas e da Natureza, 58-059-900, Joao Pessoa, BRAZIL
| | - Mariana Guedes Bezerra
- Universidade Federal da Paraiba, Departamento de Química, Centro de Ciências Exatas e da Natureza - Campus I, 58-059-900, Joao Pessoa, BRAZIL
| | - Ezequiel Fragoso Vieira Leitão
- Universidade Federal de Campina Grande, Unidade Acadêmica de Ciências Exatas e da Natureza, 58900-000, Cajazeiras, BRAZIL
| | - Silmar Andrade do Monte
- Universidade Federal da Paraiba, Quimica, Departamento de Quimica, CCEN,, Cidade Universitaria - Campus I, 58059-900, Joao Pessoa, BRAZIL
| |
Collapse
|
34
|
Haiwang Djefoulna VH, Abia D, Jules Fifen J, Nsangou M, Jaidane NE. Rotational thermodynamic parameters for asymmetric-top molecules: classical vs. quantum approaches and new analytical partition function. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2087565] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
| | - Daouda Abia
- Department of Chemistry, Faculty of Science, The University of Ngaoundere, Ngaoundere, Cameroon
| | - Jean Jules Fifen
- Quantum Theory and Aplications Unit, Department of Physics, Faculty of Science, The University of Ngaoundere, Ngaoundere, Cameroon
| | - Mama Nsangou
- Quantum Theory and Aplications Unit, Department of Physics, Faculty of Science, The University of Ngaoundere, Ngaoundere, Cameroon
- University of Maroua, Maroua, Cameroon
| | - Nejm-Eddine Jaidane
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications, Faculté des Sciences de Tunis, Université de Tunis El Manar, Tunis, Tunisie
| |
Collapse
|
35
|
Satta M, Catone D, Castrovilli MC, Bolognesi P, Avaldi L, Zema N, Cartoni A. Ion Chemistry of Carbon Dioxide in Nonthermal Reaction with Molecular Hydrogen. J Phys Chem A 2022; 126:3463-3471. [PMID: 35638704 PMCID: PMC9189832 DOI: 10.1021/acs.jpca.2c01695] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The exothermic hydrogen transfer from H2 to CO2·+ leading to H and HCO2+ is investigated in a combined experimental and theoretical work. The experimental mass/charge ratios of the ionic product (HCO2+) and the ionic reactant (CO2·+) are recorded as a function of the photoionization energy of the synchrotron radiation. Theoretical density functional calculations and variational transition state theory are employed and adapted to analyze the energetic and the kinetics of the reaction, which turns out to be barrierless and with nonthermal rate coefficients controlled by nonstatistical processes. This study aims to understand the mechanisms and energetics that drive the reactivity of the elementary reaction of CO2·+ with H2 in different processes.
Collapse
Affiliation(s)
- Mauro Satta
- Department of Chemistry, Institute of the Study of Nanostructured Materials-CNR (ISMN-CNR), Sapienza University of Rome, P. le Aldo Moro 5, Rome 00185, Italy
| | - Daniele Catone
- Institute of Structure of Matter-CNR (ISM-CNR), Area della Ricerca di Tor Vergata, Via del Fosso del Cavaliere, Rome 00133, Italy
| | - Mattea Carmen Castrovilli
- Institute of Structure of Matter-CNR (ISM-CNR), Area della Ricerca di Roma 1, Via Salaria km 29.300, Monterotondo 00015, Italy
| | - Paola Bolognesi
- Institute of Structure of Matter-CNR (ISM-CNR), Area della Ricerca di Roma 1, Via Salaria km 29.300, Monterotondo 00015, Italy
| | - Lorenzo Avaldi
- Institute of Structure of Matter-CNR (ISM-CNR), Area della Ricerca di Roma 1, Via Salaria km 29.300, Monterotondo 00015, Italy
| | - Nicola Zema
- Institute of Structure of Matter-CNR (ISM-CNR), Area della Ricerca di Tor Vergata, Via del Fosso del Cavaliere, Rome 00133, Italy
| | - Antonella Cartoni
- Department of Chemistry, Sapienza University of Rome, P. le Aldo Moro 5, Rome 00185, Italy
| |
Collapse
|
36
|
Rangel C, Espinosa-García J, Corchado JC. Full-dimensional potential energy surface for the H + CH 3OH reaction. Theoretical kinetics and dynamics study. Phys Chem Chem Phys 2022; 24:12501-12512. [PMID: 35578997 DOI: 10.1039/d2cp00864e] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dynamics and kinetics of the abstraction reactions of hydrogen atoms with methanol have been studied using quasi-classical trajectory calculations and variational transition state theory with tunnelling corrections, based on a new analytical potential energy surface (PES). The new PES is a valence-bond/molecular mechanics (VB/MM) expression that provides us with the potential energy for any set of Cartesian coordinates. Two reaction channels are considered: hydrogen abstraction from the methyl group (R1) and hydrogen abstraction from the alcohol group (R2), R1 being much more likely to occur in the wide temperature range under study (250-1000 K), as expected from the lower barrier height. Our dynamic calculations at a collision energy of 20 kcal mol-1 show that the H2 co-product is produced mainly in its vibrational ground-state and little rotation excitation is found. As for our kinetic results, they agree with those from previous theoretical studies as well as with those from kinetic experimental results (rate constants and kinetic isotopic effects), lending confidence to the analytical PES presented here. Thus, we expect this PES to be a simple yet powerful tool to understand such an important reaction in combustion chemistry at very high temperatures and interstellar chemistry at very low temperatures.
Collapse
Affiliation(s)
- Cipriano Rangel
- Área de Química Física, Facultad de Ciencias, and Instituto de Computación Científica Avanzada (ICCAEx). Universidad de Extremadura, Avenida de Elvas S/N, 06006 Badajoz, Spain.
| | - Joaquín Espinosa-García
- Área de Química Física, Facultad de Ciencias, and Instituto de Computación Científica Avanzada (ICCAEx). Universidad de Extremadura, Avenida de Elvas S/N, 06006 Badajoz, Spain.
| | - José C Corchado
- Área de Química Física, Facultad de Ciencias, and Instituto de Computación Científica Avanzada (ICCAEx). Universidad de Extremadura, Avenida de Elvas S/N, 06006 Badajoz, Spain.
| |
Collapse
|
37
|
Baiano C, Lupi J, Barone V, Tasinato N. Gliding on Ice in Search of Accurate and Cost-Effective Computational Methods for Astrochemistry on Grains: The Puzzling Case of the HCN Isomerization. J Chem Theory Comput 2022; 18:3111-3121. [PMID: 35446575 PMCID: PMC9097295 DOI: 10.1021/acs.jctc.1c01252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Indexed: 11/28/2022]
Abstract
The isomerization of hydrogen cyanide to hydrogen isocyanide on icy grain surfaces is investigated by an accurate composite method (jun-Cheap) rooted in the coupled cluster ansatz and by density functional approaches. After benchmarking density functional predictions of both geometries and reaction energies against jun-Cheap results for the relatively small model system HCN···(H2O)2, the best performing DFT methods are selected. A large cluster containing 20 water molecules is then employed within a QM/QM' approach to include a realistic environment mimicking the surface of icy grains. Our results indicate that four water molecules are directly involved in a proton relay mechanism, which strongly reduces the activation energy with respect to the direct hydrogen transfer occurring in the isolated molecule. Further extension of the size of the cluster up to 192 water molecules in the framework of a three-layer QM/QM'/MM model has a negligible effect on the energy barrier ruling the isomerization. Computation of reaction rates by the transition state theory indicates that on icy surfaces, the isomerization of HNC to HCN could occur quite easily even at low temperatures thanks to the reduced activation energy that can be effectively overcome by tunneling.
Collapse
Affiliation(s)
- Carmen Baiano
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, I-56126 Pisa, Italy
| | - Jacopo Lupi
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
38
|
Ma XR, Zhang J, Xiong YC, Zhou W. Revising the performance of the Landau–Zener surface hopping on some typical one-dimensional nonadiabatic models. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2051761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Xiang-Rui Ma
- Department of Material Physics, School of Mathematics, Physics and Optoelectronic Engineering, Hubei University of Automotive Technology, Shiyan, People's Republic of China
- Collaborative Innovation Center for Optoelectronic Technology, Hubei University of Automotive Technology, Shiyan, People's Republic of China
| | - Jun Zhang
- Department of Material Physics, School of Mathematics, Physics and Optoelectronic Engineering, Hubei University of Automotive Technology, Shiyan, People's Republic of China
- Collaborative Innovation Center for Optoelectronic Technology, Hubei University of Automotive Technology, Shiyan, People's Republic of China
| | - Yong-Chen Xiong
- Department of Material Physics, School of Mathematics, Physics and Optoelectronic Engineering, Hubei University of Automotive Technology, Shiyan, People's Republic of China
- Collaborative Innovation Center for Optoelectronic Technology, Hubei University of Automotive Technology, Shiyan, People's Republic of China
| | - Wanghuai Zhou
- Department of Material Physics, School of Mathematics, Physics and Optoelectronic Engineering, Hubei University of Automotive Technology, Shiyan, People's Republic of China
- Collaborative Innovation Center for Optoelectronic Technology, Hubei University of Automotive Technology, Shiyan, People's Republic of China
| |
Collapse
|
39
|
Xia Y, Long B, Lin S, Teng C, Bao JL, Truhlar DG. Large Pressure Effects Caused by Internal Rotation in the s-cis-syn-Acrolein Stabilized Criegee Intermediate at Tropospheric Temperature and Pressure. J Am Chem Soc 2022; 144:4828-4838. [PMID: 35262353 DOI: 10.1021/jacs.1c12324] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Criegee intermediates are important atmospheric oxidants, and quantitative kinetics for stabilized Criegee intermediates are key parameters for atmospheric modeling but are still limited. Here we report barriers and rate constants for unimolecular reactions of s-cis-syn-acrolein oxide (scsAO), in which the vinyl group makes it a prototype for Criegee intermediates produced in the ozonolysis of isoprene. We find that the MN15-L and M06-2X density functionals have CCSD(T)/CBS accuracy for the unimolecular cyclization and stereoisomerization of scsAO. We calculated high-pressure-limit rate constants by the dual-level strategy that combines (a) high-level wave function-based conventional transition-state theory (which includes coupled-cluster calculations with quasiperturbative inclusion of quadruple excitations because of the strongly multiconfigurational character of the electronic wave function) and (b) canonical variational transition-state theory with small-curvature tunneling based on a validated density functional. We calculated pressure-dependent rate constants both by system-specific quantum Rice-Ramsperger-Kassel theory and by solving the master equation. We report rate constants for unimolecular reactions of scsAO over the full range of atmospheric temperature and pressure. We found that the unimolecular reaction rates of this larger-than-previously studied Criegee intermediate depend significantly on pressure. Particularly, we found that falloff effects decrease the effective unimolecular cyclization rate constant of scsAO by about a factor of 3, but the unimolecular reaction is still the dominant atmospheric sink for scsAO at low altitudes. The large falloff caused by the inclusion of the stereoisomerization channel in the master equation calculations has broad implications for mechanistic analysis of reactions with competitive internal rotations that can produce stable rotamers.
Collapse
Affiliation(s)
- Yu Xia
- College of Mechanical and Electrical Engineering, Guizhou Minzu University, Guiyang 550025, China.,College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Bo Long
- College of Mechanical and Electrical Engineering, Guizhou Minzu University, Guiyang 550025, China.,College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Shiru Lin
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Chong Teng
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
40
|
A Theoretical Study on the Mechanism and Kinetics of the Dimeric TrimethylAluminum + O2 Reaction in the Gas Phase: A Potential Chain-initiation in the Hypergolic Combustion of TMA in Air. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113695] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
41
|
Liang N, Li Q, Ge M, Gu D, Liu Y. Mild Synthesis of a Dimethoxy‐Terminated Siloxane through a Ring‐Opening Reaction. ChemistrySelect 2022. [DOI: 10.1002/slct.202104617] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Nianjie Liang
- School of Chemistry Beihang University Beijing 100191 P. R. China
| | - Qiaosheng Li
- School of Chemistry Beihang University Beijing 100191 P. R. China
| | - Mengyuan Ge
- School of Chemistry Beihang University Beijing 100191 P. R. China
| | - Defa Gu
- School of Chemistry Beihang University Beijing 100191 P. R. China
| | - Yuzhou Liu
- School of Chemistry Beihang University Beijing 100191 P. R. China
- Beijing Advanced Innovation Center for Biomedical Engineering Beihang University Beijing 100191 P. R. China
- Beijing Shenyun Zhihe Technology Co. Ltd. Beijing 100094 P. R. China
| |
Collapse
|
42
|
Chen X, Liu M, Gao J. CARNOT: a Fragment-Based Direct Molecular Dynamics and Virtual-Reality Simulation Package for Reactive Systems. J Chem Theory Comput 2022; 18:1297-1313. [PMID: 35129348 DOI: 10.1021/acs.jctc.1c01032] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Traditionally, the study of reaction mechanisms of complex reaction systems such as combustion has been performed on an individual basis by optimizations of transition structure and minimum energy path or by reaction dynamics trajectory calculations for one elementary reaction at a time. It is effective, but time-consuming, whereas important and unexpected processes could have been missed. In this article, we present a direct molecular dynamics (DMD) approach and a virtual-reality simulation program, CARNOT, in which plausible chemical reactions are simulated simultaneously at finite temperature and pressure conditions. A key concept of the present ab initio molecular dynamics method is to partition a large, chemically reactive system into molecular fragments that can be adjusted on the fly of a DMD simulation. The theory represents an extension of the explicit polarization method to reactive events, called ReX-Pol. We propose a highest-and-lowest adapted-spin approximation to define the local spins of individual fragments, rather than treating the entire system by a delocalized wave function. Consequently, the present ab initio DMD can be applied to reactive systems consisting of an arbitrarily varying number of closed and open-shell fragments such as free radicals, zwitterions, and separate ions found in combustion and other reactions. A graph-data structure algorithm was incorporated in CARNOT for the analysis of reaction networks, suitable for reaction mechanism reduction. Employing the PW91 density functional theory and the 6-31+G(d) basis set, the capabilities of the CARNOT program were illustrated by a combustion reaction, consisting of 28 650 atoms, and by reaction network analysis that revealed a range of mechanistic and dynamical events. The method may be useful for applications to other types of complex reactions.
Collapse
Affiliation(s)
- Xin Chen
- Peking University Shenzhen Graduate School, Shenzhen, Guangdong 581055, China.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 581055, China
| | - Meiyi Liu
- Peking University Shenzhen Graduate School, Shenzhen, Guangdong 581055, China.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 581055, China
| | - Jiali Gao
- Peking University Shenzhen Graduate School, Shenzhen, Guangdong 581055, China.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 581055, China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
43
|
Ito T, Maeda S, Harabuchi Y. Kinetic Analysis of a Reaction Path Network Including Ambimodal Transition States: A Case Study of an Intramolecular Diels-Alder Reaction. J Chem Theory Comput 2022; 18:1663-1671. [PMID: 35099971 DOI: 10.1021/acs.jctc.1c01297] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This study proposes a methodology for the kinetic analysis of a reaction path network including ambimodal transition states (TSs), through which an ensemble of trajectories bifurcates to multiple minima in a phenomenon called dynamical bifurcation. The proposed methodology consists of three techniques: an automated reaction path search to construct a reaction path network including ambimodal TSs, an ab initio molecular dynamics simulation to evaluate the branching ratio, and the definition of rate constants incorporating this ratio. Applying the procedure to a Diels-Alder reaction, it was found that the inclusion of dynamical bifurcations is necessary to explain the experimental reaction yield of a byproduct. In addition, it was verified that the products take 1013 s to reach thermal equilibrium and that the experimental selectivity is determined by the dynamical bifurcations.
Collapse
Affiliation(s)
- Takuma Ito
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-8628, Japan
| | - Satoshi Maeda
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan.,Japan Science and Technology Agency (JST), Exploratory Research for Advanced Technology (ERATO) Maeda Artificial Intelligence for Chemical Reaction Design and Discovery Project, Sapporo 060-0810, Japan
| | - Yu Harabuchi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan.,Japan Science and Technology Agency (JST), Exploratory Research for Advanced Technology (ERATO) Maeda Artificial Intelligence for Chemical Reaction Design and Discovery Project, Sapporo 060-0810, Japan
| |
Collapse
|
44
|
Zhang RM, Chen W, Truhlar DG, Xu X. Master Equation Study of Hydrogen Abstraction from HCHO by OH Via a Chemically Activated Intermediate. Faraday Discuss 2022; 238:431-460. [DOI: 10.1039/d2fd00024e] [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
The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction...
Collapse
|
45
|
Zhou ZX, Zhang HX, Zheng QC. Predicting a Kind of Unusual Multiple-States Dimerization-Modes Transformation in Protein PD-L1 System by Computational Investigation and a Generalized Rate Theory. Front Chem 2021; 9:783444. [PMID: 34858950 PMCID: PMC8631179 DOI: 10.3389/fchem.2021.783444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2021] [Accepted: 10/14/2021] [Indexed: 11/20/2022] Open
Abstract
The new cancer immunotherapy has been carried out with an almost messianic zeal, but its molecular basis remains unclear due to the complexity of programmed death ligand 1 (PD-L1) dimerization. In this study, a new and integral multiple dimerization-modes transformation process of PD-L1s (with a new PD-L1 dimerization mode and a new transformation path discovered) and the corresponding mechanism are predicted using theoretical and computational methods. The results of the state analysis show that 5 stable binding states exist in system. A generalized inter-state transformation rate (GITR) theory is also proposed in such multiple-states self-assembly system to explore the kinetic characteristics of inter-state transformation. A “drug insertion” path was identified as the dominant path of the PD-L1 dimerization-modes transformation. Above results can provide supports for both the relative drug design and other multiple-states self-assembly system from the theoretical chemistry perspective.
Collapse
Affiliation(s)
- Zhong-Xing Zhou
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun, China
| | - Hong-Xing Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun, China
| | - Qing-Chuan Zheng
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun, China.,Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, College of Life Science, Jilin University, Changchun, China
| |
Collapse
|
46
|
Satta M, Casavola AR, Cartoni A, Castrovilli MC, Catone D, Chiarinelli J, Borocci S, Avaldi L, Bolognesi P. Ionization of 2- and 4(5)-Nitroimidazoles Radiosensitizers: A "Kinetic Competition" Between NO 2 and NO Losses. Chemphyschem 2021; 22:2387-2391. [PMID: 34597457 PMCID: PMC9293481 DOI: 10.1002/cphc.202100629] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2021] [Revised: 09/23/2021] [Indexed: 12/27/2022]
Abstract
Nitroimidazoles are a class of chemicals with a remarkable broad spectrum of applications from the production of explosives to the use as radiosensitizers in radiotherapy. The understanding of thedynamics of their fragmentation induced by ionizing sources is of fundamental interest. The goal of this work is to theoretically investigate the kinetic competition between the two most important decomposition channels of 2, 4 and 5‐Nitroimidazole cations: the NO and NO2 losses. The calculated rate constants of the two processes are in very good agreement with the experimental Photoelectron‐Photoion Coincidence (PEPICO) branching ratio. This study solves the intriguing and theoretically unexplained experimental observation that 2‐Nitroimidazole, at variance with the other two regio‐isomers is a source for only NO at low energies (<12.76 eV). This is a key point for biomedical application of the nitroimidazoles, because NO is the vasodilator that favors the reoxigenation of hypoxic tumor tissues.
Collapse
Affiliation(s)
- Mauro Satta
- Istituto per lo Studio dei Materiali Nanostrutturati (ISMN-CNR), Dipartimento di Chimica, Università degli studi di Roma La Sapienza, P.le Aldo Moro 5, 00185, Roma, Italy
| | - Anna Rita Casavola
- Istituto di Struttura della Materia (ISM-CNR), Area della Ricerca di Roma 1, Via Salaria Km 29,300, 00016, Monterotondo Scalo (RM), Italy
| | - Antonella Cartoni
- Dipartimento di Chimica, Università degli studi di Roma La Sapienza, Pl.e Aldo Moro 5, 00185, Roma, Italy
| | - Mattea Carmen Castrovilli
- Istituto di Struttura della Materia (ISM-CNR), Area della Ricerca di Roma 1, Via Salaria Km 29,300, 00016, Monterotondo Scalo (RM), Italy
| | - Daniele Catone
- Istituto di Struttura della Materia (ISM-CNR), Area della Ricerca di Roma 2, Via del Fosso del Cavaliere 10, 00133, Roma, Italy
| | - Jacopo Chiarinelli
- Istituto di Struttura della Materia (ISM-CNR), Area della Ricerca di Roma 1, Via Salaria Km 29,300, 00016, Monterotondo Scalo (RM), Italy
| | - Stefano Borocci
- Dipartimento per l'Innovazione nei Sistemi Biologici, Agroalimentari e Forestali (DIBAF), Università della Tuscia, Viterbo, Italy
| | - Lorenzo Avaldi
- Istituto di Struttura della Materia (ISM-CNR), Area della Ricerca di Roma 1, Via Salaria Km 29,300, 00016, Monterotondo Scalo (RM), Italy
| | - Paola Bolognesi
- Istituto di Struttura della Materia (ISM-CNR), Area della Ricerca di Roma 1, Via Salaria Km 29,300, 00016, Monterotondo Scalo (RM), Italy
| |
Collapse
|
47
|
Pelevkin AV, Loukhovitski BI, Sharipov AS. Reaction of the N Atom with Electronically Excited O 2 Revisited: A Theoretical Study. J Phys Chem A 2021; 125:8294-8312. [PMID: 34494840 DOI: 10.1021/acs.jpca.1c05733] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The kinetics of the reaction of N with electronically excited O2 (singlet a1Δg and b1Σg+ states), potentially relevant for NOx formation in nonthermal air plasma, is theoretically studied using the multireference second-order perturbation theory. The corresponding thermodynamically and kinetically favored reaction pathways together with possible intersystem crossings are identified. It has been revealed that the energy barrier for the N + O2(a1Δg) → NO + O reaction is approximately twice the barrier height for the counterpart process with O2(X3Σg-). The molecular oxygen in the b1Σg+ state, in turn, proved to be even less reactive to atomic nitrogen than O2(a1Δg). Appropriate thermal rate constants for specified reaction channels are calculated by the variational transition-state theory incorporating corrections for the tunneling effect, nonadiabatic transitions, and anharmonicity of vibrations for transition states and reactants. The corresponding three-parameter Arrhenius expressions for the broad temperature range (T = 300-4000 K) are reported. At last, post-transition-state molecular dynamics simulations indicate that the N + O2(a1Δg) reaction produces vibrationally much colder NO molecules than the N + O2(X3Σg-) process.
Collapse
Affiliation(s)
- Alexey V Pelevkin
- Central Institute of Aviation Motors, Aviamotornaya 2, Moscow 111116, Russia
| | - Boris I Loukhovitski
- Central Institute of Aviation Motors, Aviamotornaya 2, Moscow 111116, Russia.,Joint Institute for High Temperatures of the Russian Academy of Sciences, Izhorskaya 13 Bldg. 2, Moscow 125412, Russia
| | - Alexander S Sharipov
- Central Institute of Aviation Motors, Aviamotornaya 2, Moscow 111116, Russia.,Joint Institute for High Temperatures of the Russian Academy of Sciences, Izhorskaya 13 Bldg. 2, Moscow 125412, Russia
| |
Collapse
|
48
|
Wu X, Huang C, Chai J, Zhang F. Formation of Substituted Alkyls as Precursors of Peroxy Radicals with a Rapid H-Shift in the Atmosphere. J Phys Chem Lett 2021; 12:8790-8797. [PMID: 34491756 DOI: 10.1021/acs.jpclett.1c02503] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Long straight-chain alkyl peroxy (ROO) radicals substituted with C═C and oxo functional groups are expected to undergo a rapid hydrogen shift (H-shift), which is a critical step in the atmospheric autoxidation mechanism. The existence of a weak tertiary C-H bond plays a key role in the rapid H-shift. Here, the reaction kinetics between OH and two typical long straight-chain functionalized volatile organic compounds, 3-methyl-1-hexene (3-MH) and 2-methylpentanal (2-MP), was theoretically investigated to reveal the fate of the weak C-H bond. The results indicate that the most favored reaction pathways are direct consumption (H-abstraction of 2-MP) and indirect destruction (addition of OH to 3-MH) of the "weak" tertiary C-H bond. The yields of abstraction pathways producing precursors of ROO radicals that undergo rapid H-shifts are computed to be less than 10% for both 3-MH + OH and 2-MP + OH reactions.
Collapse
Affiliation(s)
- Xiaoqing Wu
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, P. R. China
- National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei, Anhui 230029, P. R. China
| | - Can Huang
- Chair of Technical Thermodynamics, RWTH Aachen University, 52062 Aachen, Germany
| | - Jiajue Chai
- Institute at Brown for Environment and Society, and Department of Earth, Environmental and Planetary Sciences, Brown University, 182 Hope Street, Providence, Rhode Island 02912, United States
| | - Feng Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, P. R. China
| |
Collapse
|
49
|
Sikes T, Banyon C, Schwind RA, Lynch PT, Comandini A, Sivaramakrishnan R, Tranter RS. Initiation reactions in the high temperature decomposition of styrene. Phys Chem Chem Phys 2021; 23:18432-18448. [PMID: 34612384 PMCID: PMC8409502 DOI: 10.1039/d1cp02437j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The thermal decomposition of styrene was investigated in a combined experimental, theory and modeling study with particular emphasis placed on the initial dissociation reactions. Two sets of shock tube/time-of-flight mass spectrometry (TOF-MS) experiments were performed to identify reaction products and their order of appearance. One set of experiments was conducted with a miniature high repetition rate shock tube at the Advanced Light Source at Lawrence Berkeley National Laboratory using synchrotron vacuum ultraviolet photoionization. The other set of experiments was performed in a diaphragmless shock tube (DFST) using electron impact ionization. The datasets span 1660–2260 K and 0.5–12 atm. The results show a marked transition from aromatic products at low temperatures to polyacetylenes, up to C8H2, at high temperatures. The TOF-MS experiments were complemented by DFST/LS (laser schlieren densitometry) experiments covering 1800–2250 K and 60–240 Torr. These were particularly sensitive to the initial dissociation reactions. These reactions were investigated theoretically and revealed the dissociation of styrene to be a complex multichannel process with strong pressure and temperature dependencies that were evaluated with multi-well master equation simulations. Simulations of the LS data with a mechanism developed in this work are in excellent agreement with the experimental data. From these simulations, rate coefficients for the dissociation of styrene were obtained that are in good agreement with the theoretical predictions. The simulation results also provide fair predictions of the temperature and pressure dependencies of the products observed in the TOF-MS studies. Prior experimental studies of styrene pyrolysis concluded that the main products were benzene and acetylene. In contrast, this study finds that the majority of styrene dissociates to create five styryl radical isomers. Of these, α-styryl accounts for about 50% with the other isomers consuming approximately 20%. It was also found that C–C bond scission to phenyl and vinyl radicals consumes up to 25% of styrene. Finally the dissociation of styrene to benzene and vinylidene accounts for roughly 5% of styrene consumption. Comments are made on the apparent differences between the results of this work and prior literature. A combined theoretical and experimental study showing styrene primarily decomposes to styryl radicals + H.![]()
Collapse
Affiliation(s)
- Travis Sikes
- Chemical Sciences and Engineering Division, Argonne National Laboratory, 9700 South Cass Ave., Lemont, IL 60439, USA.
| | | | | | | | | | | | | |
Collapse
|
50
|
Edeleva M, Van Steenberge PH, Sabbe MK, D’hooge DR. Connecting Gas-Phase Computational Chemistry to Condensed Phase Kinetic Modeling: The State-of-the-Art. Polymers (Basel) 2021; 13:3027. [PMID: 34577928 PMCID: PMC8467432 DOI: 10.3390/polym13183027] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Revised: 09/02/2021] [Accepted: 09/03/2021] [Indexed: 02/06/2023] Open
Abstract
In recent decades, quantum chemical calculations (QCC) have increased in accuracy, not only providing the ranking of chemical reactivities and energy barriers (e.g., for optimal selectivities) but also delivering more reliable equilibrium and (intrinsic/chemical) rate coefficients. This increased reliability of kinetic parameters is relevant to support the predictive character of kinetic modeling studies that are addressing actual concentration changes during chemical processes, taking into account competitive reactions and mixing heterogeneities. In the present contribution, guidelines are formulated on how to bridge the fields of computational chemistry and chemical kinetics. It is explained how condensed phase systems can be described based on conventional gas phase computational chemistry calculations. Case studies are included on polymerization kinetics, considering free and controlled radical polymerization, ionic polymerization, and polymer degradation. It is also illustrated how QCC can be directly linked to material properties.
Collapse
Affiliation(s)
- Mariya Edeleva
- Laboratory for Chemical Technology (LCT), Ghent University, Technologiepark 125, 9052 Zwijnaarde, Belgium; (P.H.M.V.S.); (M.K.S.)
| | - Paul H.M. Van Steenberge
- Laboratory for Chemical Technology (LCT), Ghent University, Technologiepark 125, 9052 Zwijnaarde, Belgium; (P.H.M.V.S.); (M.K.S.)
| | - Maarten K. Sabbe
- Laboratory for Chemical Technology (LCT), Ghent University, Technologiepark 125, 9052 Zwijnaarde, Belgium; (P.H.M.V.S.); (M.K.S.)
- Industrial Catalysis and Adsorption Technology (INCAT), Ghent University, Valentin Vaerwyckweg 1, 9000 Ghent, Belgium
| | - Dagmar R. D’hooge
- Laboratory for Chemical Technology (LCT), Ghent University, Technologiepark 125, 9052 Zwijnaarde, Belgium; (P.H.M.V.S.); (M.K.S.)
- Centre for Textile Science and Engineering (CTSE), Ghent University, Technologiepark 70a, 9052 Zwijnaarde, Belgium
| |
Collapse
|