Conventional Transition State Theory Research Articles

Systematic study on the hydrogen abstraction reactions from oxygenated compounds by H and HO2

AbstractTo extend the rule‐based approach for hydrogen abstraction reactions from oxygenated compounds, a systematic investigation was performed to examine the reactivity of gas‐phase hydrogen abstraction reactions from alkyl groups (methyl and ethyl groups) bound to oxygen atoms in five types of oxygenated compounds (alcohols, ethers, formate esters, acetate esters, and carbonate esters) by H atoms and HO2 radicals comprehensively considering rotational conformers. Quantum chemical calculations were conducted at the CBS‐QB3 level for stationary points. Rate constants were determined employing conventional transition state theory (TST). For hydrogen abstraction reactions by H, the rotational conformer distribution partition function was employed to approximate partition functions, owing to the similarity in vibrational energy‐level structures among conformers. In hydrogen abstraction reactions by HO2, the vibrational structures of transition‐state (TS) conformers varied significantly due to the hydrogen bonding, leading to an inappropriate evaluation of rate constants when using the lowest‐energy conformer as a representative. Therefore, the rate constants were calculated by the multi‐structural TST. It was revealed that the differences in functional groups containing O atoms mainly affect the bond dissociation energies of the C–H bonds and the activation energies of hydrogen abstraction reactions only when the C atoms are adjacent to the O atoms. Additionally, it was found that hydrogen bonds formed in the TSs show minor effect on rate parameters for the overall rate constants, apart from the reduction of the pre‐exponential factors for the H‐abstraction reactions from the methylene position of ethyl groups. The comparison with the rate constants from previous studies showed reasonable results, indicating that the rate constants in this study, which thoroughly consider rotational conformers, can be the current best estimates.

Computational Modeling of the Conformation-Dependent Atmospheric Reactivity of Criegee Intermediates.

The impacts of Criegee intermediates (CIs) on atmospheric chemistry depend significantly on the CI conformation. In this Perspective, I highlight examples of how electronic structure and statistical rate theory calculations, in conjunction with experiment, have revealed conformation-dependent details of both CI ground-state reactivity and electronic excitation. Calculations using single-reference electronic structure methods and conventional transition state theory have predicted that CIs with syn-alkyl or syn-vinyl substituents isomerize rapidly to vinyl hydroperoxides (VHPs) or dioxoles, both of which can decompose rapidly under atmospheric conditions. Ongoing computational research on hydroxyl radical (OH) roaming initiated by VHP dissociation requires the application of multireference electronic structure methods and variational transition state theory. CIs that lack both syn-alkyl and syn-vinyl substituents undergo either bimolecular reaction or π* ← π electronic excitation in the atmosphere. Accurate predictions of CI ultraviolet-visible spectra require multireference calculations with large active spaces and at least a second-order perturbative treatment of dynamic electron correlation. The extent to which electronic spectra can be diagnostic of the presence of specific CI conformers varies significantly with CI chemical identity.

Insight into the Kinetics of Isomerization and β-Scission Reactions Following H-Abstraction in 2-Pentanone Combustion.

2-Pentanone is a significant carbon-neutral fuel. To better understand 2-pentanone combustion, the CCSD(T)/CBS method was used to calculate the potential energy surfaces (PES) for H-abstraction, isomerization, and β-scission reactions of 2-pentanone. Rate constants for the above reactions were calculated by the MESS employing conventional transition state theory (CTST) at 300-2000 K. The findings reveal adherence to the Evans-Polanyi principle in the H-abstraction reactions of 2-pentanone by H atoms. Specifically, the 2-site shows more competitive kinetics due to having the lowest energy barrier of 7.8 kcal/mol. The 4-position has the highest energy barrier, making the reaction difficult to occur. In the subsequent reactions, the breaking of C-H bonds is most competitive at low temperatures. In particular, the H transfer from INT1 to INT3 has a potential well height of 24.8 kcal/mol. The β-scission reactions become the main pathway for the depletion of pentanone radicals with increasing temperature. This study significantly extends the kinetic parameters of 2-pentanone combustion over a wide range of temperatures and pressures, which is expected to contribute to the development of 2-pentanone combustion models.

Theoretical kinetics study of key reactions between ammonia and fuel molecules, part II: H-atom abstraction from alcohols and ethers by ṄH2 radicals

Ammonia, as a carbon-free fuel, is expected to be a carbon-neutral alternative fuel to control pollution emissions. However, with lower reactivity, pure ammonia needs to be blended with highly reactive biofuels including alcohols and ethers, for a better combustion performance. Alcohols and ethers can be produced from various sources, including renewable materials and fossil fuels. Therefore, the combination of ammonia with alcohols/ethers is a promising choice for developing carbon-neutral energy solutions. As an important intermediate, ṄH2 radicals can easily be formed during the pyrolysis or oxidation processes for ammonia combustion. The H-atom abstraction reaction between ṄH2 radicals and fuel molecules is vital in determining the blending fuel reactivity. High-level ab initio calculations have been carried out in this study to perform a systematic theoretical chemical kinetic investigation of H-atom abstraction by ṄH2 radicals from C1–C5 alcohols, 2,3-dimethyl-2-butanol, C1–C4 ethers and methyl isobutyl ether with a total of 61 H-atom abstraction reaction channels. The QCISD(T)/CBS//M06–2X/6–311++G(d, p) level of theory was employed to investigate the potential energy surfaces of the title reactions involved. Electronic geometry optimization, frequency analysis, and zero-point energy correction were performed for reactants, transition states and products related to 61 reaction channels at the M06–2X/6–311++G(d, p) level of theory. Single-point energy calculations were performed at the QCISD(T)/cc-pVXZ (X = D, T) and MP2/cc-pVYZ (Y = D, T, Q) level of theory and then extrapolated to the complete basis set. Conventional transition state theory was used to calculate the rate constants for different channels in the temperature range from 500 to 2000 K. The calculated individual and average rate constants for H-atom abstraction from different sites of alcohols and ethers were also obtained to explore the kinetic influence from the functional group. This is the first systematic study providing the rate constants for the title reactions which can be used to develop the combustion kinetic models for ammonia/alcohols and ammonia/ethers blends.

Thermal decomposition study of 1,3-dioxolane as a representative of battery solvent

As industrial compounds and organic solvent, 1,3-dioxolane (C3H6O2) is used in lithium-ion batteries arousing extensive interest. The pyrolysis of C3H6O2 at 150 torr pressure was conducted in the temperature range of 743–1013 K. In order to quantify the major products, experimental measurements of the photoionization cross-section (PICS) and photoionization efficiency (PIE) for C3H6O2 were carried out by using synchrotron radiation photoionization mass spectrometry (SR-PIMS). Twenty-three intermediates were detected in pyrolysis of C3H6O2 including the newly observed ketene, acetaldehyde, ethenol, ethanol, dimethyl ether, diacetylene and 2-propen-1-ol which played important roles in developing mechanism. The H-abstraction of C3H6O2 by H, CH3 and OH radicals were studied via high level ab initio calculations. EStokTP code was used to find and calculate optimal geometries, vibration frequencies and intrinsic reaction coordinate (IRC) at the theory level of M06–2X/6–311+G(d,p). Conventional transition state theory (TST) was used and Rice-Ramsperger-Kassel-Marcus/Master Equation (RRKM/ME) was solved to obtain rate coefficients. For H-abstraction by H and CH3, the rates are increasing versus temperature, and rates of H-abstraction by OH show different tendency. These rate parameters were used to update the model of C3H6O2, which could well predict the ignition delay times (IDTs) in shock tube and oxidation data in jet stirred reactor (JSR) from literature. For the pyrolysis, H-abstraction reactions by H, CH3 and OH and unimolecular reactions contribute to solvent decay. Detailed pathways of C3H6O2 decomposition were also given and analyzed through this experiment and previous literature. The findings, including the experimental data and the calculations will assist to better understand pyrolysis characteristics of C3H6O2 as well.

Radical scavenging activity of bromophenol analogs: analysis of kinetics and mechanisms.

This theoretical study explores the antioxidant activity of five bromophenol analogs, with a particular focus on their interaction with different solvent environments of varying polarities. Key findings include the correlation between increased solvent polarity and enhanced antioxidant activity of these analogs, comparable in some instances to ascorbic acid. Notably, compound 5, developed by our research team, demonstrates superior antioxidant activity in both lipid and aqueous solutions, surpassing that of ascorbic acid and other tested analogs. This research contributes to the understanding of bromophenol analogs, presenting the first known kinetic and chemical stability data such as rate constants, pKa values, and branching ratios for reactions with the methylperoxyl radical (CH3OO•). The computational analyses were conducted using the Gaussian 09 software suite at the M05-2X/6-31 + G(d) computational level. These analyses employed conventional transition state theory to account for various potential mechanisms and effects of solvent polarity on the antioxidant activities of bromophenol analogs. The study meticulously calculated enthalpy under standard conditions (298.15K and 1atm) with necessary thermodynamic corrections. Additionally, the Quantum Mechanics-based Test for Overall Radical Scavenging Activity (QMORSA) protocol guided the evaluation of radical scavenging activity, ensuring a comprehensive assessment of the antioxidant potential of the compounds.

Theoretical studies on the kinetics and mechanism of hydroxyl radical reaction with quinclorac and quinmerac herbicides in aqueous media

AbstractThe kinetic and mechanistic studies for the reaction of hydroxyl radical with two quinoline based herbicides, namely, quinclorac and quinmerac has been performed using various computational methods in aqueous media. Geometry optimizations were performed using Density Functional Theory (DFT) methods including water as the solvent. Local reactivity parameters of these herbicides towards the •OH radical are predicted using condensed Fukui function. Single point energies of various species were calculated using double hybrid method, namely, B2PLYP–D for better accuracy. The pKa values for these acid based herbicides allow them to exist in deprotonated form in aqueous condition. Hence, the calculations are also performed for the deprotonated or the anionic form apart from the neutral species. Individual rate coefficients for •OH radical addition reaction with each carbon atoms were evaluated using conventional transition state theory using one–dimensional tunneling corrections. The solvent effect on reaction is implemented through Collins–Kimball formulations. Both the approaches, namely, the Fukui index and individual rate constant determination confirms that the most reactive site for the •OH radical addition in these two herbicide is the carbon atom attached to the COOH group. The total rate constant for the •OH radical reaction with both neutral and anionic forms of these two herbicides are relatively high and equal to its diffusion‐limit value. Evaluation of the ecotoxicities of the parent herbicides and their OH adducts is estimated using the structure–activity relationship concept.

Theoretical Investigation on the H Atom Abstraction Reaction from C1-C4 Alkanes and Alkenes by ṄH2 Radicals.

The H atom abstraction reactions from alkanes and alkenes by ṄH2 are decisive in predicting the combustion characteristics of NH3/CxHy binary fuels. Theoretical investigation is carried out on the energy barriers of H atom abstraction reactions from C1-C4 alkanes/alkenes by ṄH2 radicals at the QCISD(T)/CBS//M06-2X/6-311++G(d,p) level of theory. Single-point energies of each species are computed using QCISD(T)/cc-pVDZ, TZ level of theories with basis set corrections from MP2/cc-pVDZ, TZ, and QZ methods. One-dimensional hindered rotor potentials are obtained by the M06-2X/cc-pVTZ method with 10° increment. Rate constants of each channel across temperatures of 298.15-2000 K are calculated by solving the RRKM/Master Equation with conventional transition state theory. For alkanes, rate constants order follows ktertiary> ksecondary> kprimary, while for alkenes the order follows kallylic> kprimary> kvinylic. Among the vinylic carbon sites within the same alkene species, the hydrogen atom sharing the same carbon with the allylic carbon on the C-C double bond is the preferred site for the H atom abstraction reaction. The branching ratio results indicate that the abstraction from tertiary or secondary carbon sites on alkanes and allylic carbon sites on alkenes are dominating during the investigated temperature range but become less important as the temperature increases. The data provided in this work are in good agreement with the literature data, but for the ṄH2+alkenes system, the literature data are scarce and further investigation is needed.

Theoretical Prediction of the Reaction Probabilities of H, O, and OH Radicals on the Polypropylene Surface.

To determine the H-abstraction reaction probabilities of H/O/OH radicals with a polypropylene (PP) surface, a first-principles calculation was performed based on the DLPNO-CCSD(T)/CBS//M06-2X-D3/def-TZVP theory level. The PP chain model used in this study was 2,4,6-trimethylheptane. The rate constants of the H/O/OH radicals with the isolated PP chain model were calculated based on the conventional transition-state theory. By comparing the experimental values and considering the error factors and their compensation, it was concluded that the orders of magnitude of the predicted rate constants were accurate. The resulting rate constants were converted to reaction probabilities between the H/O/OH radicals and the PP surface. The method used in this study is applicable for obtaining theoretical values of surface reaction probabilities based on first-principles calculations. The calculation at the DLPNO-CCSD(T)/CBS theory level has high accuracy but consumes a large amount of computational resources. The study also demonstrated that the double-hybrid functionals, wB97x-2-D3(BJ) and rev-DSD-PBEP86-D3(BJ), with a 3-ζ or 4-ζ basis set, could reproduce the electronic energy values obtained from DLPNO-CCSD(T)/CBS while using only approximately 1/100 of the computational resources required by the latter under our computer configuration.

Rapid Atmospheric Reactions between Criegee Intermediates and Hypochlorous Acid.

Hypochlorous acid (HOCl) is a paramount compound in the atmosphere due to its significant contribution to both tropospheric oxidation capacity and ozone depletion. The main removal routes for HOCl are photolysis and the reaction with OH during the daytime, while these processes are unimportant during the nighttime. Here, we report the rapid reactions of Criegee intermediates (CH2OO and anti/syn-CH3CHOO) with HOCl by using high-level quantum chemical methods as the benchmark; their accuracy is close to coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis limit by extrapolation [CCSDT(Q)/CBS]. Their rate constants have been calculated by using a dual-level strategy; this combines conventional transition state theory calculated at the benchmark level with variational transition state theory with small-curvature tunneling by a validated density functional method. We find that the introduction of the methyl group into Criegee intermediates not only affects their reactivities but also exerts a remarkable influence on anharmonicity. The calculated results uncover that anharmonicity increases the rate constants of CH2OO + HOCl by a factor of 18-5, while it is of minor importance in the anti/syn-CH3CHOO + HOCl reaction at 190-350 K. The present findings reveal that the loose transition state for anti-CH3CHOO and HOCl is a rate-determining step at 190-350 K. We also find that the reaction of Criegee intermediates with HOCl contributes significantly to the sink of HOCl during the nighttime in the atmosphere.

Theoretical Study of the NO Reduction Mechanism on Biochar Surfaces Modified by Li and Na Single Adsorption and OH Co-Adsorption

Theoretical and experimental investigations have shown that biochar, following KOH activation, enhances the efficiency of NO removal. Similarly, NaOH activation also improves NO removal efficiency, although the underlying mechanism remains unclear. In this study, zigzag configurations were employed as biochar models. Density functional theory (DFT) was utilized to examine how Li and Na single adsorption and OH co-adsorption affect the reaction pathways of NO reduction on the biochar surface. The rate constants for all reaction-determining steps (RDSs) within a temperature range of 200 to 1000 K were calculated using conventional transition state theory (TST). The results indicate a decrease in the activation energy for NO reduction reactions on biochar when activated by Li and Na adsorption, thus highlighting their beneficial role in NO reduction. Compared to the case with Na activation, Li-activated biochar exhibited superior performance in terms of the NO elimination rate. Furthermore, upon the adsorption of the OH functional group onto the Li-decorated and Na-decorated biochar models (LiOH-decorated and NaOH-decorated chars), the RDS energy barriers were higher than those of Li and Na single adsorption but easily overcome, suggesting effective NO reduction. In conclusion, Li-decorated biochar showed the highest reactivity due to its low RDS barrier and exothermic reaction on the surface.

Theoretical kinetic study of key reactions between ammonia and fuel molecules, part I: Hydrogen atom abstraction from alkanes by ṄH2 radical

Cross reactions between ṄH2 radical and fuel molecules are essential in determining the reactivity for ammonia/fuel blends. A systematic theoretical chemical kinetic investigation of hydrogen atom abstraction by ṄH2 radical from C1–C6 alkanes, n-heptane, and isooctane has been performed in this study. The geometry optimization, frequency, and zero-point energy calculations for all species related to the title reactions were calculated at the M06-2X/6-311++G(d,p) level of theory. One-dimensional hindered rotor treatment has also been performed at the same level of theory to obtain the potential energy surface for the low frequency torsional modes. Single-point energy was obtained at the QCISD(T)/cc-pVDZ, QCISD(T)/cc-pVTZ, MP2/cc-pVDZ, MP2/cc-pVTZ and MP2/cc-pVQZ levels of theory and then extrapolated to the complete basis set. Benchmark calculations for single point energy have also been carried out at the CCSD(T)-F12/cc-pVTZ, CCSD(T)-F12/cc-pVQZ, CCSD(T)-F12/aug-cc-pVTZ and CCSD(T)-F12/aug-cc-pVQZ levels of theory based on the geometries obtained at the CCSD(T)/cc-pVDZ level of theory for the reaction systems of methane + ṄH2 and ethane + ṄH2. The rate constants of 41 reactions in the temperature range 298.15–2000 K were calculated by using the Master Equation System Solver (MESS) with conventional transition state theory. Good agreement has been gained between our calculated results and the available theoretical and experimental data in the literature. Rate rules for different reaction sites have also been derived, and the average rate constants for primary, secondary, and tertiary sites were also reported on a per hydrogen atom basis. The impact of calculation data for the ammonia/alkane blends model has also been investigated. The calculated rate constants are valuable for the combustion kinetic model development for ammonia blending fuels.

Theoretical and kinetic study of the H-atom abstraction reactions by Ḣ atom from alkyl cyclohexanes.

Reaction kinetics of hydrogen atom abstraction from six alkyl cyclohexanes, methyl cyclohexane (MCH), ethyl cyclohexane (ECH), n-propyl cyclohexane (nPCH), iso-propyl cyclohexane (iPCH), sec-butyl cyclohexane (sBCH) and iso-butyl cyclohexane (iBCH), by the Ḣ atom are systematically studied in this work. The M06-2X method combined with the 6-311++G(d,p) basis set is used to perform geometry optimization, frequency analysis and zero-point energy calculations for all species. The intrinsic reaction coordinate (IRC) calculations are performed to confirm the transition states connecting the reactants and products correctly. One-dimensional hindered rotors are used to treat the low frequency torsional models with potentials scanned at the M06-2X/6-31G level of theory. Electronic single-point energy calculations for all reactants, transition states, and products are performed at the QCISD(T)/CBS level of theory. High-pressure limiting rate constants of 39 reaction channels are obtained using conventional transition state theory with asymmetric Eckart tunneling corrections in the temperature range 298.15-2000 K. Reaction rate rules for H-atom abstraction by the Ḣ atom from alkyl cyclohexanes on primary, secondary and tertiary carbon sites on both the side chain and ring are provided. The obtained rate constants are given by the Arrhenius expression in the temperature range 500-2000 K, which can be used for the combustion kinetics model development for alkyl cyclohexanes.

Theoretical study on the H-atom abstraction reactions of pentanol + HȮ2, part I: five branched pentanol isomers.

The chemical kinetic studies of hydrogen atom (H-atom) abstraction reactions by hydroperoxyl (HȮ2) radicals from five branched pentanol isomers, including 3-methyl-1-butanol, 2-methyl-1-butanol, 1,1-dimethyl-1-propanol, 1,2-dimethyl-1-propanol, and 2,2-dimethyl-1-propanol were investigated systematically through high-level ab initio calculations. Geometry optimization, frequency analysis, and zero-point energy (ZPE) corrections were performed for six reactants, twenty-three transition states (TSs), and twenty-four products at the M06-2X/6-311++G(d,p) level of theory. The intrinsic reaction coordinate calculation was performed at the same level of theory to confirm the transition state connection. The one-dimensional hindered rotor treatment for low-frequency torsional modes was also treated at the M06-2X/6-311++G(d,p) level of theory. The QCISD(T)/CBS level of theory was used to calculate the single-point energies for the species whose T1 diagnostic value was lower than 0.035. At the same time, the CASPT2/CBS level of theory was used to calculate the single-point energies for the channel in which the T1 diagnostic value of transition states was greater than 0.035. Rate constants for the H-atom abstraction reactions from the five branched pentanol isomers by HȮ2 radicals were calculated by using conventional transition state theory with asymmetric Eckart tunneling corrections in the temperature range of 500-2000 K. Rate constants and branching ratios for the title reactions and the rate rules for ten different H-atom abstraction types were investigated. Temperature-dependent thermochemistry properties for all reactants and products were calculated by the composite methods of G3/G4/CBS-QB3/CBS-APNO, which were in good agreement with the data available in the literature. Rate constants for the H-atom abstraction reactions by HȮ2 radical from branched pentanol isomers were investigated in this work as part I, and those for linear pentanol isomers will be analyzed in part II. All the calculated kinetics and thermochemistry data can be utilized in the model development for branched pentanol isomers oxidation.

Oxidation pathways and kinetics of the 1,1,2,3-tetrafluoropropene (CF2CF-CH2F) reaction with Cl-atoms and subsequent aerial degradation of its product radicals in the presence of NO.

To give a comprehensive account of the environmental acceptability of 1,1,2,3-tetrafluoropropene (CF2CF-CH2F) in the troposphere, we have examined the oxidation reaction pathways and kinetics of CF2CF-CH2F initiated by Cl-atoms using the second-order Møller-Plesset perturbation (MP2) theory along with the 6-31+G(d,p) basis set. We also performed single-point energy calculations to further refine the energies at the CCSD(T) level along with the basis sets 6-31+G(d,p) and 6-311++G(d,p). The estimation of the relative energies and thermodynamic parameters of the CF2CF-CH2F + Cl reaction clearly shows that Cl-atom addition reaction pathways are more dominant compared to H-abstraction reaction pathways. The value of the rate coefficient for each reaction channel is calculated using the conventional transition state theory (TST) over the temperature range of 200-1000 K at 1 atm. The estimated overall rate coefficients for the title reaction are found to be 1.10 × 10-12, 1.21 × 10-10, and 1.13 × 10-8 cm3 per molecule per s via the respective calculation methods viz. MP2/6-31+G(d,p), CCSD(T)//MP2/6-31+G(d,p), and CCSD(T)/6-311++G(d,p)//MP2/6-31+G(d,p), at 298.15 K. Moreover, the calculated rate coefficients and percentage branching ratio values suggest that the Cl-atom addition reaction at the β-carbon atom is more preferable to that of the α-carbon addition to CF2CF-CH2F. Based on the rate coefficient values calculated by the three different methods, the atmospheric lifetime for the title reaction at 298.15 K is estimated. The radiative efficiency (RE) and Global Warming Potential (GWP) results of the title molecule show that its GWP would be negligible. Further, we have explored the degradation of its product radicals in the presence of O2 and NO. From the degradation results, we have found that CF2(Cl)COF, FCOCH2F, FCFO and FCOCl are formed as stable end products along with various radicals such as ˙CF2Cl and ˙CH2F. Therefore, these findings of kinetic and mechanistic data can be applied to the development and implementation of a novel CFC replacement.

Experimental and kinetic study on 2-methoxyethanol part I: Molecular chemistry and pyrolysis investigation

Glycol ethers have been widely recognized as an alternative biofuels. 2-Methoxyethanol (ME2, HOCH2CH2OCH3) is the simplest and most representative glycol ether. Hence, the combustion chemistry of ME2 has aroused extensive interest. The ME2 pyrolysis has been investigated experimentally in a jet-stirred reactor (JSR) at atmospheric pressure within the temperature range of 650–1075 K for the first time. FTIR, GC–MS, and GC-FID were used to identify and quantify species profiles in ME2 pyrolysis. Sixteen intermediates were newly detected by the present work including acetaldehyde, hydroxyacetaldehyde, and methyl vinyl ether which are essential for mechanism development. The H abstraction reactions of ME2 by H, CH3, HO2, OH, O2, and the unimolecular reactions of four fuel-related radicals were first studied via high level ab initio calculations. The geometries optimization, vibrational frequencies, and intrinsic reaction coordinate (IRC) calculations were carried out at the M06–2X/6–311++G(d,p) level of theory. Conventional transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus/Master Equation (RRKM/ME) were used to calculate the temperature- and pressure-dependent rate constants from 300 to 2000 K and 0.01 to 100 atm for the first time. It is found that the H abstraction reactions at the two middle carbon atoms are the most favored pathways, and b-scissions for fuel-related radicals have the lowest energy barriers. The thermochemical properties of ME2 and its derived radicals have been obtained via the combination of composite methods. A mechanism involving 596 species and 3123 reactions was developed for the first time to predict the ME2 pyrolysis characteristics based on experimental and calculated results obtained by the present work. Good agreements were achieved between the measured and simulated results. Rate of production (ROP), as well as sensitivity analysis, were performed to get a deeper insight into the ME2 pyrolysis. It has been revealed that the major ME2 decomposition channels are the H abstractions by H, CH3, and OH radicals, and unimolecular reactions have scarce contributions to fuel decay. The results, in terms of the measured data and of the calculated results, and the newly developed kinetic model as well as findings obtained within this study will assist and contribute to a better understanding of combustion characteristics of ME2 and of glycol ethers as well.

Insight into the fate of nitrogen during char thermal conversion and the influence mechanism of potassium: A theoretical research

Nitrogen oxides (NOₓ) are primary pollutants produced during biomass combustion. During the devolatilization stage, char nitrogen (char(N)) is formed. In the subsequent char combustion stage, char(N) can decompose directly into NOx precursors or engage in heterogeneous reactions with O2 or NO to form NO and N2. Nonetheless, a comprehensive understanding of the reaction mechanisms and competitiveness of char(N) migration, especially the influence of the alkali metal potassium (K) present in biomass, remains incomplete. Building on the Zigzag char(N) models, the present study delves into the migration reactions of char(N), assessing their competitive dynamics through the integration of density functional theory, electronic structure analysis, and conventional transition state theory. Furthermore, it examines the impact of K on char(N) conversion. The competitiveness of the heterogeneous reactions follows the sequence: heterogeneous reduction of NO to N2 > heterogeneous oxidation of char(N) to NO > decomposition of char(N) to NOx precursors. Moreover, the formation of HCN is more favorable than NH₃ production. The successive conversion from char(N) to NO and then to N2 is the predominant migration route for char(N), with NO generation as the pivotal step. The less preferred char(N) migration route involves decomposition to NH3/HCN, followed by oxidation to NOx within the main combustion zone, which cannot be mitigated by char. K can accelerate NO generation and sustain the primacy of the heterogeneous NO reduction, consequently enhancing the oxidation-reduction process of char(N). As a result, K plays a constructive role in managing NOx emissions during the thermal conversion of char.

Utility of (MgO)12 nanocage as a chemical sensor for recognition of amphetamine drug: A computational inspection

DFT calculations on sensor-drug interactions are necessary for understanding binding mechanisms, predicting sensor performance, evaluating stability and reactivity, and rational design of sensor materials. We scrutinized the adsorption of amphetamine (AFE) on the pure magnesium oxide nano-cage (MgONC) by applying density functional theory. All geometries and single point energy computations were optimized at M06–2X/6–311 G (d, p). Furthermore, we performed an analysis of the natural bond orbital (NBO) and evaluated the values of partial natural charges. Additionally, we investigated donor-acceptor (D-A) interactions and examined the Wiberg bond index (WBI) in greater depth. The MgONC was capable of adsorbing AFE with greater strength with the energy of adsorption (Eads) of −48.19 kcal/mol (for stable configurations). Moreover, the NBO method demonstrated more effective D-A interactions between AFE and the MgONC. Based on the computations, for the most stable configuration, there was a substantial alteration in the HOMO-LUMO gap of the MgONC following the drug adsorption, thus increasing the electrical conductance (EC) of the MgONC. The sensing mechanism is related to the gap difference, which depends on the change in the EC. We adopted the conventional transition state theory for the prediction of recovery time. The computations indicated that the MgONC+ AFE configuration had a short recovery time for the desorption of AFE. Finally, based on our findings, we could conclude that the MgONC is an appropriate choice for the improvement of effective AFE sensors. DFT study of drug sensors will focus on enhancing sensitivity, selectivity, and stability while exploring novel materials and optimizing performance through theoretical simulations and analysis.

Theoretical study of H-atom abstraction by CH3OȮ radicals from aldehydes and alcohols: Ab initio and comprehensive kinetic modeling

Reaction kinetics of H-atom abstraction from C1C3 alcohols, including CH3OH, C2H5OH, NC3H7OH, IC3H7OH, and C1C2 aldehydes, including CH2O, CH3CHO, by CH3OȮ radicals is investigated in this work through high-level ab initio calculations. Electronic structure optimizations are performed for all the stationary points with M06–2X/6−311++g(d,p) method. Quadratic configuration interaction method QCISD(T)/cc-pVXZ (where X = D and T) and Møller–Plesset perturbation theory MP2/cc-pVXZ (where X = D, T, and Q) are used to calculate single point energies. Subsequently, rate coefficients for all H-atom abstraction channels are determined using conventional transition state theory with unsymmetric tunneling corrections. These calculated results are updated to a recently updated model and then compared against expansive experimental datasets to investigate their influence on the model prediction performance. The updated model shows obvious discrepancies with the original model, where the updated model is less reactive for CH3CHO, C2H5OH, NC3H7OH, and IC3H7OH, while showing negligible differences for CH3OH and CH2O due to the lack of CH3OȮ radical formation. Sensitivity and flux analyses are further conducted, through which the difference between the studied species in their oxidation pathways and the relative importance of H-atom abstraction by CH3OȮ radicals is highlighted. With the updated rate parameters, the branching ratios of H-atom abstractions from CH2O and CH3CHO are significantly altered for CH3CHO, C2H5OH, NC3H7OH and IC3H7OH, with an obvious shift away from the abstraction channel by CH3OȮ. Further analyses highlight the critical roles of CH2O and CH3CHO chemistry as core chemistries for large hydrocarbons (e.g., carbon number > 1), which have been inadequately described in existing chemistry models. Systematic efforts are urgently needed to improve the existing CH2O and CH3CHO chemistries.

Polyrate 2023: A computer program for the calculation of chemical reaction rates for polyatomics. New version announcement

Polyrate 2023: A computer program for the calculation of chemical reaction rates for polyatomics. New version announcement