RESUMO
This study provides a detailed understanding of how the reaction between CH3NH, one of the primary products of the CH3NH2 + OH/Cl reactions, and NOx occurs in the atmosphere since the reaction is expected to be a dominant sink for the tropospheric CH3NH radical. First, we focus on the reaction of the aminyl radical CH3NH with NO2, complementing the known reaction between CH3NH and NO, to provide the overall picture of the CH3NH + NOx system. The reaction was meticulously examined across the extended range of temperature (298-2000 K) and pressure (0.76-76 000 torr) using quantum chemistry calculations and kinetic modeling based on the framework of the Rice-Ramsperger-Kassel-Marcus (RRKM)-based master equation. Highly correlated electronic structure calculations unveil that the intricate reaction mechanism of the CH3NH + NO2 reaction, which can proceed through O-addition or N-addition to form NO2, encompasses numerous steps, channels, and various intermediates and products. The temperature-/pressure-dependent kinetic behaviors and product distribution of the CH3NH + NO2 reaction are revealed under atmospheric and combustion conditions. The main products under atmospheric conditions are found to be CH3NHO and NO, as well as CH3NHNO2, while under combustion conditions, the primary products are only CH3NHO and NO. Given its stability under ambient conditions, CH3NHNO2, a nitramine, is believed to have the potential to induce DNA damage, which can ultimately result in severe cancers. Secondly, by building upon prior research on the CH3NH + NO system, this study shows that the reaction of CH3NH with NOx holds greater importance in urban areas with elevated NOx emissions than other oxidants like O2. Furthermore, this reaction occurs swiftly and results in the creation of various compounds, such as the carcinogenic nitrosamine (CH3NHNO), carcinogenic nitramine (CH3NHNO2), CH3NNOH, (CH3NN + H2O) and (CH3NHO + NO).
RESUMO
Hymexazol is a volatile fungicide widely used in agriculture, causing its abundance in the atmosphere; thus, its atmospheric fate and conversion are of great importance when assessing its environmental impacts. Herein, we report a theoretical kinetic mechanism for the oxidation of hymexazol by OH radicals, as well as the subsequent reactions of its main products with O2 and then with NO by using the Rice-Ramsperger-Kassel-Marcus-based Master equation kinetic model on the potential energy surface explored at the ROCBS-QB3//M06-2X/aug-cc-pVTZ level. The predicted total rate constants ktotal(T, P) for the reaction between hymexazol and OH radicals show excellent agreement with scarcely available experimental values (e.g., 3.6 × 10-12 vs (4.4 ± 0.8) × 10-12 cm3/molecule/s at T = 300 K and P = 760 Torr); thus, the calculated kinetic parameters can be confidently used for modeling/simulation of N-heterocycle-related applications under atmospheric and even combustion conditions. The model shows that 3,4-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-5-yl (IM2), 3,5-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-4-yl (IM3), and (3-hydroxy-1,2-oxazol-5-yl)methyl (P8) are the main primary intermediates, which form the main secondary species of (3,4-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-5-yl)dioxidanyl (IM4), (3,5-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-4-yl)dioxidanyl (IM7), and ([(3-hydroxy-1,2-oxazol-5-yl)methyl]dioxidanyl (IM11), respectively, through the reactions with O2. The main secondary species then can react with NO to form the main tertiary species, namely, (3,4-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-5-yl)oxidanyl (P19), (3,5-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-4-yl)oxidanyl (P21), and [(3-hydroxy-1,2-oxazol-5-yl)methyl]oxidanyl (P23), respectively, together with NO2. Besides, hymexazol could be a persistent organic pollutant in the troposphere due to its calculated half-life τ1/2 of 13.7-68.1 h, depending on the altitude.
Assuntos
Atmosfera , Cinética , OxirreduçãoRESUMO
The widespread use of vinyl butyrate (CH2îCHOC(O)CH2CH2CH3 or VB) in the polymer industry and daily-life materials inevitably results in its emission into the atmosphere. Therefore, understanding the mechanism and kinetics of the VB conversion is critical for evaluating its fate and environmental impacts. Herein, we theoretically investigate the chemical transformation of VB initiated by OH radicals in the atmosphere using the stochastic Rice-Ramsperger-Kassel-Marcus (RRKM)-based master equation kinetic model on the potential energy surface explored at the M06-2X/aug-cc-pVTZ level of theory. Showing excellent agreement with limited experimental kinetic data, the VB + OH kinetic model reveals that H-abstraction from Cß (i.e., -CßH2CH3) prevails over the OH-addition to the double bond (CîC), even at low temperatures. The detailed analyses, including those of the time-resolved species profiles, reaction rate, and reaction flux, reveal the reaction mechanism shift with temperature (causing the U-shaped temperature dependence of k(T, P)) and the noticeable pressure dependence of k(T, P) at low temperatures. The secondary chemistry under atmospheric conditions (namely, the reaction of the main product with O2 and its subsequent reactions with NO) was then characterized within the same framework to reveal the detailed kinetic mechanism (e.g., [4-(ethenyloxy)-4-oxobutan-2-yl]oxidanyl (IM12) + NO2 is the dominant channel under atmospheric conditions), suggesting VB is not a persistent organic pollutant and a new environmental concern regarding the formed NO2. Also, the kinetic behaviors of vinyl butyrate and its oxidation products were extended from atmospheric to combustion conditions for further applications. Moreover, through TD-DFT calculations, it is shown that several related important species (i.e., 1-(ethenyloxy)-1-oxobutan-2-yl (P4), [4-(ethenyloxy)-4-oxobutan-2-yl]dioxidanyl (IM7), and IM12) can potentially undergo photolysis in the atmosphere.
RESUMO
In this study, a series of 14 Cu (II), Zn (II), Ni (II) and Ag (I) complexes containing bis-benzimidazole derivatives were successfully designed and synthesized from 2-(1H-benzimidazole-2-yl)-phenol derivatives and corresponding metal salt solutions. The compound structures were identified by FT-IR, 1H-NMR, powder X-ray diffraction and ESI-MS analyses, and the presence of the metal in the complexes was confirmed by ultraviolet-visible spectroscopy and ICP optical emission spectrometry. Electronic structure calculations were also carried out to describe the detailed structures in addition to the electronic absorption spectra of the ligands. The cytotoxic activity of the complexes was evaluated against three human cancer cell lines: lung (A549), breast (MDA-MB-231) and prostate (PC3) cancer cells. All complexes inhibited anti-proliferative cancer cells better than free ligands, especially Zn (II) and Ag (I) complexes, which are most sensitive to MDA-MB-231 cells. In addition, showing the growth inhibition of three cancer cell lines with IC50 < 10.4 µM, complexes C1 , C3 and C14 could be considered potential multi-targeted anti-cancer agents.
RESUMO
Diterpenoids are abundant and important compounds in Euphorbia species owing to their structural diversity; therefore, in this study, we investigate the modern-concept antioxidant activities, including free-radical scavenging and oxidative DNA damage repairing, of highly oxygenated diterpenoids originating from the aerial part of Euphorbia helioscopia. Four compounds with structural types of ent-abietane, containing a fused furan ring in their structures, including euphelionolide A (1), euphelionolide D (2), euphelionolide I (3), and euphelionolide L (4) are selected. First, the radical-scavenging activity of these compounds was evaluated with two typical radicals HOO⢠and HO⢠in water and pentyl ethanoate (PEA, to mimic lipid environment) via three main mechanisms, namely, hydrogen atom transfer (HAT), radical adduct formation (RAF), and single electron transfer. It is found that the studied compounds are able to scavenge free radicals at multiple reactive sites favorably via HAT and RAF mechanisms, in which the former dominates in the case with HOO⢠while both mechanisms are competitive in the reaction with HOâ¢. Second, chemical repairing of DNA damage is modeled with the H-atom and single electron being transferred from the studied molecules to damaged 2'-deoxyguanosine (2dG) (i.e., 2dG⢠radicals and 2dGâ¢+ radical cation). Among the four compounds, euphelionolide A is shown as the most effective radical scavenger and also the highest potential species for chemical repairing of radical-damaged DNA in both water and PEA.
Assuntos
Diterpenos , Antioxidantes , DNA , Diterpenos/química , Sequestradores de Radicais Livres/química , Sequestradores de Radicais Livres/farmacologia , Radicais Livres , Hidrogênio , ÁguaRESUMO
This work reports a detailed mechanism of the initial thermal pyrolysis of isopropyl propionate, (C2H5C(=O)OCH(CH3)2), an important biodiesel additive/surrogate, for a wide range of T = 500-2000 K and P = 7.6-76 000 Torr. The detailed kinetic behaviors of the title reaction on the potential energy surface constructed at the CBS-QB3 level were investigated using the RRKM-based master equation (RRKM-ME) rate model, including hindered internal rotation (HIR) and tunneling corrections. It is revealed that the C3H6 elimination occurring via a six-centered retro-ene transition state is dominant at low temperatures, while the homolytic fission of the C-C bonds becomes more competitive at higher temperatures. The tunneling treatment is found to slightly increase the rate constant at low temperatures (e.g., â¼1.59 times at 563 K), while the HIR treatment, being important at high temperatures, decreases the rate (e.g., by 5.9 times at 2000 K). Showing a good agreement with experiments in low-temperature kinetics, the kinetic model reveals that the pressure effect should be taken into account at high temperatures. Finally, the temperature- and pressure-dependent kinetic mechanism, consisting of the calculated thermodynamic and kinetic data, is provided for further modeling and simulation of any related systems.
RESUMO
This study theoretically reports the comprehensive kinetic mechanism of the aniline + OH reaction in the range of 200-2000 K and 0.76-7600 Torr. The temperature- and pressure-dependent behaviors, including time-resolved species profiles and rate coefficients, were studied within the stochastic RRKM-based master equation framework with the reaction energy profile, together with molecular properties of the species involved, characterized at the M06-2X/aug-cc-pVTZ level. Hindered internal rotation and Eckart tunneling treatments were included. The H-abstraction from the -NH2 moiety (to form C6H5NH (P1)) is found to prevail over the OH-addition on the C atom at the ortho site of aniline (to form 6-hydroxy-1-methylcyclohexa-2,4-dien-1-yl (I2)) with the atmospheric rate expressions (in cm3/molecule/s) as kabstraction(P1) = 3.41 × 101 × T-4.56 × exp (-255.2 K/T) for 200-2000 K and kaddition(I2) = 3.68 × 109 × T-7.39 × exp (-1163.9 K/T) for 200-800 K. The U-shaped temperature-dependent characteristics and weakly positive pressure dependence at low temperatures (e.g., T ≤ 800 K and P = 760 Torr) of ktotal(T) are also observed. The disagreement in ktotal(T) between the previous calculations and experimental studies is also resolved, and atmospheric aniline is found to be primarily removed by OH radicals (τOH â¼ 1.1 h) in the daytime. Also, via TD-DFT simulations, it is recommended to include P1 and I2 in any atmospheric photolysis-related model.
Assuntos
Compostos de Anilina , Radical Hidroxila , Cinética , Oxirredução , TemperaturaRESUMO
The statistical molecular fragmentation (SMF) model was used to analyze the 306 fragmentation channels (containing 611 different species) that result from the fluorene (C13H10+) cation losing up to three hydrogen atoms (neutral radicals and/or a proton). Breakdown curves from such analysis permit one to extract experimentally inaccessible information about the fragmentation of the fluorene cation, such as the locations of the lost hydrogen atoms (or proton), yields of the neutral fragments, electronic states of the residues, and quantification of very low probability channels that would be difficult to detect. Charge localization during the fragmentation pathways was studied to provide a qualitative understanding of the fragmentation process. Breakdown curves for both the fluorene cation and neutral fluorene were compared. The SMF results match the rise and fall of the one hydrogen loss yield experimentally measured by imaging photoelectron-photoion coincidence spectroscopy using a VUV synchrotron.
RESUMO
The comprehensive kinetic mechanism of the OH-initiated gas-phase oxidation of pyrrole is first theoretically reported in a broad range of conditions (T = 200-2000 K &P = 1-7600 Torr). On the potential energy surface constructed at the M06-2X/aug-cc-pVTZ level, the temperature- and pressure-dependent behaviors of the title reaction were characterized using the stochastic Rice-Ramsperger-Kassel-Marcus based Master Equation (RRKM-ME) rate model. The corrections of the hindered internal rotation and quantum tunneling treatments were included. The calculated results reveal the competition between the two distinct pathways: OH-addition and direct H-abstraction. The former channels are found favorable at low-temperature and high-pressure range (e.g., T < 900 K and P = 760 Torr) where non-Arrhenius and positive pressure-dependent behaviors of the rate constants are noticeably observed, while the latter predominate at temperatures higher than 900 K at atmospheric pressure and no pressure dependence on the rate constant is found. The predicted global rate constants are in excellent agreement with laboratory values; thus, the derived kinetic parameters are recommended for modeling/simulation of N-heterocycle-related applications in atmospheric and even in combustion conditions. Besides, pyrrole should not be considered as a persistent organic pollutant owing to its short atmospheric lifetime (â¼1 h) towards OH radicals. The secondary mechanisms of the subsequent reactions of two OH-pyrrole adducts (namely, I1 and I2) with two abundant species, O2/NO, which are relevant to the atmospheric degradation process, were also investigated. It is also revealed by TD-DFT calculations that two OH-pyrrole adducts (I1 &I2), nine intermediates, Ii (i = 3-11) and four products (P1, P2, P3 and P6) can undergo photodissociation under the sunlight.
Assuntos
Atmosfera , Radical Hidroxila , Cinética , Oxirredução , PirróisRESUMO
In this study, the detailed kinetic mechanism of the trans-decalin + OH reaction is firstly investigated for a wide range of conditions (i.e., T = 200-2000 K & P = 0.76-76 000 Torr) using the M06-2X/aug-cc-pVTZ level and stochastic Rice-Ramsperger-Kassel-Marcus based master equation (RRKM-ME) rate model, which includes corrections of the hindered internal rotor (HIR) and tunneling effects. Our predicted global rate constant excellently matches with the scarce experimental measurement (R. Atkinson, et al. Int. J. Chem. Kinet., 1983, 15, 37-50). The H-abstraction channel from Cα of trans-decalin is found to be dominant at low temperatures. A U-shaped temperature-dependent behavior and slightly positive pressure-dependence at low temperatures (e.g., T ≤ 400 K & P = 760 Torr) of the total rate constants are also observed. Detailed analysis reveals that the HIR treatment is essential to capture the kinetic behavior while the tunneling correction only plays a minor role.
RESUMO
The detailed reaction mechanism of the N2H4 + OH reaction is comprehensively reported for a wide range of conditions (i.e., T = 200-3000 K & P = 1-7600 Torr) using the CCSD(T)/CBS//M06-2X/6-311++G(3df,2p) level and the master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) rate model, which includes corrections of the hindered internal rotor (HIR) and tunneling treatments. Our calculated rate constants are found in excellent agreement with the latest experimental data (G. L. Vaghjiani, Int. J. Chem. Kinet., 2001, 33, 354-362), which helps to resolve the discrepancy between the previous experimental and theoretical studies. The reaction mechanism is revealed as: (i) the H-abstraction channel is more thermodynamically favorable than the OH-substitution mechanism; (ii) non-Arrhenius behaviors and slightly positive pressure-dependence at low temperature (T≤ 500 K) of the rate coefficients are observed and (iii) the HIR treatment plays a substantial role in obtaining the reliable rate constants. Moreover, the performance of several molecular electronic structure methods (i.e., M06-2X, B3LYP, BH&HLYP and MP2) on the rate coefficient calculations is also discussed thoroughly in this work.
RESUMO
In this short communication, we resolve the discrepancy in chemical kinetics between the recent theoretical work by Safaei et al. (Phys. Chem. Chem. Phys., 2019, 21, 8445) and the experimental measurement for the reaction of imidazole initiated by OH radicals. In particular, using a more comprehensive potential energy surface (PES) explored with the same electronic structure method used by Safaei and coworkers (i.e., M06-2X/aug-cc-pVTZ) and a more rigorous stochastic master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) rate model which includes corrections for the hindered internal rotation and tunneling treatment, we reported the calculated rate constants which are in excellent agreement with the experimental data. Furthermore, it is suggested that imidazole should not be considered as a persistent organic pollutant due to its short atmospheric lifetime of â¼3.7 hours towards the removal of OH radicals.
RESUMO
This work provides a rigorous detailed kinetic study on the C2H2 + NH2 reaction in a wide range of conditions (T = 250-2000 K & P = 1-76000 Torr). In particular, the composite method W1U was used to construct the potential energy surface on which the kinetic behaviors were characterized within the state-of-the-art master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) framework. Corrections of the hindered internal rotation (HIR) treatment and quantum tunneling effect were included. A clear reaction mechanism shift with respect to both temperature and pressure was revealed via detailed kinetic and species analyses. In particular, bimolecular products (i.e., CH2[double bond, length as m-dash]C[double bond, length as m-dash]NH + H, CH[triple bond, length as m-dash]CNH2 + H, CH3CN + H, CH[triple bond, length as m-dash]C· + NH3 in the decreasing mole fraction order) can be formed directly from the reactants at high temperature and/or low pressure while they can be produced indirectly via intermediates (e.g., ·CH[double bond, length as m-dash]CHNH2(cis), ·CH[double bond, length as m-dash]CHNH2(trans), CH2[double bond, length as m-dash]C·NH2, ) at low temperature and/or high pressure. The calculated rate constants are in good agreement with the literature data from ab initio calculations without any adjustment; thus, the proposed temperature- and pressure-dependent rate constants, together with the thermodynamic data of the species involved, can be confidently used for modeling NH2-related systems under atmospheric and combustion conditions.
RESUMO
The C2F4 + O3 reaction plays an important role in the oxidation process of perfluoroalkenes in the atmosphere. The detailed reaction mechanism was explored using the accurate electronic structure method, CCSD(T)/CBS//B3LYP/aug-cc-pVTZ. The 1,3-cycloaddition of O3 with C2F4 to form the primary ozonide was found to be the rate-determining step of the oxidation process with a small barrier (i.e., 7.3 kcal mol-1 at 0 K). The temperature- and pressure-dependent behaviors of the title reaction were characterized in the range of 200-1000 K & 0.1-760 Torr using the integrated deterministic and stochastic master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) rate model with the inclusion of the corrections for anharmonicity and tunneling treatments. It is found that the anharmonic effect plays a role in the kinetic behaviors (e.g., lower the rate by a factor of â¼ two at 298 K) while the tunneling correction is insignificant. The total rate constants were found to be pressure-independent under the considered conditions, shown as ktot(T) = 4.80 × 10-23 × T2.69 × exp(-2983.4 K/T) (cm3 per molecule per s), which confirms the latest experimental data by Acerboni et al. (G. Acerboni, N. R. Jensen, B. Rindone and J. Hjorth, Chem. Phys. Lett., 1999, 309, 364-368); thus this study helps to resolve a long-term controversy among the previous measurements. The sensitivity analyses on the derived rate coefficients and time-resolved species mole fraction with respect to the ab initio input parameters were also performed to further understand as well as quantify the kinetic behaviors for the title reaction.
RESUMO
In an attempt to construct detailed kinetic mechanisms for biodiesel fuels on the fly, high-pressure rate rules for the concerted HO2 elimination reaction class were derived using a comprehensive training reaction set of more than 60 reactions involving different alkyl methyl/ethyl ester peroxy radicals (RCOOR')-OOâ¢. Using the composite electronic structure method CBS-QB3 in combination with classical statistical mechanics and the transition state theory (TST) rate model, high-pressure rate constants for the reactions in the training set as well as thermodynamic properties for the species involved were calculated. The corrections from Eckart tunneling and hindered internal rotation (HIR) treatments were also included in the calculations. The results reveal that the ester group (-COO-) selectively promotes the reaction when compared with the traditional hydrocarbon fuels; thus it is recommended that the seven derived rate rules for the title reaction class (including the thermodynamic data of the species involved in the NASA format) should be used for construction of detailed kinetic mechanisms for real biodiesel molecules.
RESUMO
This work provides a rigorous procedure, within the framework of the Reaction Class Transition State Theory (RC-TST) and the Structure-Activity Relationship (SAR), for predicting reliable thermal rate constants on-the-fly for hydrogen abstraction reactions by methyl/ethyl radicals from Polycyclic Aromatic Hydrocarbons (PAHs) in a temperature range of 300-3000 K. All necessary RC-TST parameters were derived from ab initio calculations for a representative set of 36 reactions on which different error analyses and comparisons with available literature data were carried out. In addition to the good agreement between the RC-TST rate constants and the literature data, the detailed error analyses show that RC-TST/SAR, utilizing either the Linear Energy Relationship (LER) where only the reaction energy is needed or Barrier Height Grouping (BHG) where no additional data is needed, can predict the thermal rate constants for any reaction in the title reaction class with an average systematic error of less than 50% when compared to the explicit rate calculations. Therefore, the constructed RC-TST procedure can be confidently used to obtain reliable rate constants on the fly in an attempt to effectively construct detailed kinetic mechanisms for PAH-related fuels.
RESUMO
The detailed kinetic mechanism of the HOSO2 + 3O2 reaction, which plays a pivotal role in the atmospheric oxidation of SO2, was investigated using accurate electronic structure calculations and novel statistical thermodynamic/kinetic models. Explored using the accurate composite method W1U, the detailed potential energy surface (PES) revealed that the addition of O2 to a HOSO2 radical to form the adduct (HOSO4) proceeds via a transition state with a slightly positive barrier (i.e., 0.7 kcal mol-1 at 0 K). Such a finding compromises a long-term hypothesis about this channel of being a barrierless process. Moreover, the overall reaction was found to be slightly exothermic by 1.7 kcal mol-1 at 0 K, which is in good agreement with recent studies. On the newly-constructed PES, the temperature- and pressure-dependent behaviors of the title reaction were characterized in a wide range of conditions (T = 200-1000 K & P = 10-760 Torr) using the integrated deterministic and stochastic master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) rate model in which corrections for hindered internal rotation (HIR) and tunneling treatments were included. The calculated numbers were found to be in excellent agreement with literature data. The sensitivity analyses on the derived rate coefficients with respect to the ab initio input parameters (i.e., barrier height and energy transfer) were also performed to further understand the kinetic behaviors of the title reaction. The detailed kinetic mechanism, consisting of thermodynamic and kinetic data (in NASA polynomial and modified Arrhenius formats, respectively), was also provided at different T & P for further use in the modeling/simulation of any related systems.
RESUMO
Master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) has shown to be a powerful framework for modeling kinetic and dynamic behaviors of a complex gas-phase chemical system on a complicated multiple-species and multiple-channel potential energy surface (PES) for a wide range of temperatures and pressures. Derived from the ME time-resolved species profiles, the macroscopic or phenomenological rate coefficients are essential for many reaction engineering applications including those in combustion and atmospheric chemistry. Therefore, in this study, a least-squares-based approach named Global Minimum Profile Error (GMPE) was proposed and implemented in the MultiSpecies-MultiChannel (MSMC) code (Int. J. Chem. Kinet., 2015, 47, 564) to extract macroscopic rate coefficients for such a complicated system. The capability and limitations of the new approach were discussed in several well-defined test cases.
RESUMO
An integrated deterministic and stochastic model within the master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) framework was first used to characterize temperature- and pressure-dependent behaviors of thermal decomposition of acetic anhydride in a wide range of conditions (i.e., 300-1500 K and 0.001-100 atm). Particularly, using potential energy surface and molecular properties obtained from high-level electronic structure calculations at CCSD(T)/CBS, macroscopic thermodynamic properties and rate coefficients of the title reaction were derived with corrections for hindered internal rotation and tunneling treatments. Being in excellent agreement with the scattered experimental data, the results from deterministic and stochastic frameworks confirmed and complemented each other to reveal that the main decomposition pathway proceeds via a 6-membered-ring transition state with the 0 K barrier of 35.2 kcal·mol-1. This observation was further understood and confirmed by the sensitivity analysis on the time-resolved species profiles and the derived rate coefficients with respect to the ab initio barriers. Such an agreement suggests the integrated model can be confidently used for a wide range of conditions as a powerful postfacto and predictive tool in detailed chemical kinetic modeling and simulation for the title reaction and thus can be extended to complex chemical reactions.
RESUMO
The reaction kinetics of dimethyl carbonate (DMC) and OH radicals were investigated behind reflected shock waves over the temperature range of 872-1295 K and at pressures near 1.5 atm. Reaction progress was monitored by detecting OH radicals at 306.69 nm using a UV laser absorption technique. The rate coefficients for the reaction of DMC with OH radicals were extracted using a detailed kinetic model developed by Glaude et al. (Proc. Combust. Inst. 2005, 30(1), 1111-1118). The experimental rate coefficients can be expressed in Arrhenius form as: kexpt'l = 5.15 × 1013 exp(-2710.2/T) cm3 mol-1 s-1. To explore the detailed chemistry of the DMC + OH reaction system, theoretical kinetic analyses were performed using high-level ab initio and master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) calculations. Geometry optimization and frequency calculations were carried out at the second-order Møller-Plesset (MP2) perturbation level of theory using Dunning's augmented correlation consistent-polarized valence double-ζ basis set (aug-cc-pVDZ). The energy was extrapolated to the complete basis set using single point calculations performed at the CCSD(T)/cc-pVXZ (where X = D, T) level of theory. For comparison purposes, additional ab initio calculations were also carried out using composite methods such as CBS-QB3, CBS-APNO, G3 and G4. Our calculations revealed that the H-abstraction reaction of DMC by OH radicals proceeds via an addition elimination mechanism in an overall exothermic process, eventually forming dimethyl carbonate radicals and H2O. Theoretical rate coefficients were found to be in excellent agreement with those determined experimentally. Rate coefficients for the DMC + OH reaction were combined with literature rate coefficients of four straight chain methyl ester + OH reactions to extract site-specific rates of H-abstraction from methyl esters by OH radicals.