The interplay of membrane fluidity, acyl chain order and area per lipid on the partitioning of two antidepressants paroxetine and sertraline.

Selective serotonin reuptake inhibitors (SSRIs) are among the popular drugs for treating depression and mental disorders. Membrane fluidity has previously been considered as the main factor in modulating the membrane partitioning of SSRIs, while other biophysical properties, such as the acyl chain order and area per lipid, were often neglected. Varying the lipid membrane composition and temperature can significantly modify the physical phase and, in turn, affect its fluidity, acyl chain order and area per lipid. Here, we investigate the role of membrane fluidity, acyl chain order and area per lipid in the partitioning of two SSRIs, paroxetine (PAX) and sertraline (SER). The model membranes were either POPC : SM (1 : 1 mol ratio) or POPC : SM : Chol (1 : 1 : 1 mol ratio) and studied in the temperature range of 25-45 °C. The order parameters and area per lipid in the two lipid mixtures were calculated using molecular dynamics simulations. The membrane partitioning of PAX and SER was determined via second derivative spectrophotometry. In a lower temperature range (25-32 °C), membrane fluidity favors the SSRI partitioning into Lo/Ld POPC:SM:Chol. In a higher temperature range (37-45 °C), the interplay between membrane fluidity, acyl chain order and area per lipid favors drug partitioning into Ld POPC:SM. The findings offer indication for the inconsistent distribution of SSRIs in tissues as well as the possible interaction of SSRIs with lipid domains and membrane-bound proteins.

Mechanistic and Kinetic Insights into OH-Initiated Atmospheric Oxidation of Hymexazol: A Computational Study.

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.

OH-initiated oxidation of vinyl butyrate: ab initio insights.

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.

Ab initio kinetics of the CH3NH + NO2 reaction: formation of nitramines and N-alkyl nitroxides.

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).

On The Radical Scavenging and DNA Repairing Activities by Natural Oxygenated Diterpenoids: Theoretical Insights.

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.

Reaction kinetics of 1,4-cyclohexadienes with OH radicals: an experimental and theoretical study.

This work presents the OH-initiated oxidation kinetics of 1,4-cyclochexadiene (1,4-CHD). The temperature dependence of the reaction was investigated by utilizing a laser flash photolysis flow reactor and laser-induced fluorescence (LPFR/LIF) technique over the temperature range of 295-438 K and a pressure of ∼50 torr. The kinetics of the reaction was followed by measuring the LIF signal of OH radicals near 308 nm. The reaction of OH radicals with 1,4-CHD exhibited a clear negative temperature dependence. To discern the role of various channels, ab initio and RRKM-based ME calculations (RRKM-ME) were performed over temperatures of 200-2000 K and pressures of 0.76-7600 torr. The computed energy profile revealed that the reaction proceeds via the formation of a pre-reaction van der Waals complex at the entrance channel. The complex was found to be more stable than that usually seen in other alkenes + OH reactions. Both the addition channel and the abstraction reaction of allylic hydrogen were found to have negative energy barriers. Interestingly, the abstraction reaction exhibited a negative temperature dependence at low temperatures and contributed significantly (∼37%) to the total rate coefficients even under atmospheric conditions. At T ≥ 900 K, the reaction was found to proceed exclusively (>95%) via the abstraction channel. Due to the competing channels, the reaction of OH radicals with 1,4-CHD displays complicated kinetic behaviours, reflecting the salient features of the energy profile. The role of competing channels was fully characterized by our kinetic model. The calculated rate coefficients showed excellent agreement with the available experimental data.

New Mechanistic Insights into Atmospheric Oxidation of Aniline Initiated by OH Radicals.

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.

Multiple dehydrogenation of fluorene cation and neutral fluorene using the statistical molecular fragmentation model.

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.

Detailed kinetics of hydrogen abstraction from trans-decalin by OH radicals: the role of hindered internal rotation treatment.

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.

Comment on "Atmospheric oxidation reactions of imidazole initiated by hydroxyl radicals: kinetics and mechanism of reactions and atmospheric implications" by Safaei et al., Phys. Chem. Chem. Phys., 2019, 21, 8445.

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.

Ab initio kinetics of the C2H2 + NH2 reaction: a revisited study.

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.

Ab initio dynamics of hydrogen abstraction from N2H4 by OH radicals: an RRKM-based master equation study.

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.

Prediction Model for Antimalarial Activities of Hemozoin Inhibitors by Using Physicochemical Properties.

The rapid spread of strains of malaria parasites that are resistant to several drugs has threatened global malaria control. Hence, the aim of this study was to predict the antimalarial activity of chemical compounds that possess anti-hemozoin-formation activity as a new means of antimalarial drug discovery. After the initial in vitro anti-hemozoin-formation high-throughput screening (HTS) of 9,600 compounds, a total of 224 hit compounds were identified as hemozoin inhibitors. These 224 compounds were tested for in vitro erythrocytic antimalarial activity at 10 µM by using chloroquine-mefloquine-sensitive Plasmodium falciparum strain 3D7A. Two independent experiments were conducted. The physicochemical properties of the active compounds were extracted from the ChemSpider and SciFinder databases. We analyzed the extracted data by using Bayesian model averaging (BMA). Our findings revealed that lower numbers of S atoms; lower distribution coefficient (log D) values at pH 3, 4, and 5; and higher predicted distribution coefficient (ACD log D) values at pH 7.4 had significant associations with antimalarial activity among compounds that possess anti-hemozoin-formation activity. The BMA model revealed an accuracy of 91.23%. We report new prediction models containing physicochemical properties that shed light on effective chemical groups for synthetic antimalarial compounds and help with in silico screening for novel antimalarial drugs.

Detailed kinetics of tetrafluoroethene ozonolysis.

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.

Ab initio kinetics of the HOSO2 + 3O2 → SO3 + HO2 reaction.

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.

Global minimum profile error (GMPE) - a least-squares-based approach for extracting macroscopic rate coefficients for complex gas-phase chemical reactions.

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.

On-the-fly kinetics of hydrogen abstraction from polycyclic aromatic hydrocarbons by methyl/ethyl radicals.

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.

Low Temperature Oxidation Kinetics of Biodiesel Molecules: Rate Rules for Concerted HO2 Elimination from Alkyl Ester Peroxy Radicals.

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.

A combined high-temperature experimental and theoretical kinetic study of the reaction of dimethyl carbonate with OH radicals.

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.

Kinetics of Thermal Unimolecular Decomposition of Acetic Anhydride: An Integrated Deterministic and Stochastic Model.

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.