Carbon-Centered Radical Addition and β-Scission Reactions: Modeling of Activation Energies and Pre-exponential Factors

Machine learning applications for thermochemical and kinetic property prediction

Detailed kinetic models play a crucial role in comprehending and enhancing chemical processes. A cornerstone of these models is accurate thermodynamic and kinetic properties, ensuring fundamental insights into the processes they describe. The prediction of these thermochemical and kinetic properties presents an opportunity for machine learning, given the challenges associated with their experimental or quantum chemical determination. This study reviews recent advancements in predicting thermochemical and kinetic properties for gas-phase, liquid-phase, and catalytic processes within kinetic modeling. We assess the state-of-the-art of machine learning in property prediction, focusing on three core aspects: data, representation, and model. Moreover, emphasis is placed on machine learning techniques to efficiently utilize available data, thereby enhancing model performance. Finally, we pinpoint the lack of high-quality data as a key obstacle in applying machine learning to detailed kinetic models. Accordingly, the generation of large new datasets and further development of data-efficient machine learning techniques are identified as pivotal steps in advancing machine learning's role in kinetic modeling.

Bright Nanocomposites based on Quantum Dot-Initiated Photocatalysis

Integrating quantum dots (QDs) into polymer matrix to form nanocomposites without compromising the QD photoluminescence (PL) is crucial to emerging QD light-emitting and solar energy conversion fields. However, the most widely-used bulk polymerization technique, where monomers serve as the QD solvent, usually leads to QD PL quenching caused by radical initiators. Here we demonstrate high-brightness nanocomposites with near-unity PL quantum yield (QY), through a novel QDs-catalyzed (-initiated) bulk polymerization without using any radical initiators. Different from previous reports where QDs were designed as photo-sensitizers/catalysts (always with cocatalysts) and hence non-emissive in catalytic conditions, our QDs combine high brightness with highly effective catalysis, a combination that was previously considered to be hardly possible. In our case, apart from emitting light (at a large probability), the photoexcited QDs act as 'overall reaction' catalysts by simultaneously employing photoexcited electrons and holes to produce active radicals without the need of any sacrificial agents. These active radicals, though with a small amount, are sufficient to initiate effective chain reaction-dominated bulk polymerization, eliminating the requirement of extra radical initiators. This study provides new insights for understanding and development of QDs for energy applications.

Modeling the kinetics of hydrogen abstraction reactions in nitrogen-containing compounds via group additivity

New group additivity values are presented to enable the modeling of a broad range of intermolecular hydrogen abstraction reactions involving nitrogen-containing compounds. From a dataset of 316 reaction rate coefficients calculated at the CBS-QB3 level of theory in the high-pressure limit, 76 group additivity values and 14 resonance corrections have been estimated. The influence of substituents on both the attacked hydrogen and attacking radical, being a carbon or nitrogen atom, has been investigated systematically. The new group additivity models can be applied to approximate the Arrhenius parameters of hydrogen abstraction reactions of nitrogen-containing compounds by hydrogen atoms, carbon-centered and nitrogen-centered radicals in the 300-1800 K temperature range. Complementary to the group additivity model, correlations for the tunneling coefficients, which depend on both the temperature and the activation energy of the reaction in the exothermic direction, have been generated. The good performance of the new group additivity schemes has been demonstrated using a test set of reactions. At 1000 K, the rate coefficients for all test set reactions are approximated on average within a factor of 1.45, 1.47 and 1.34, for the hydrogen abstractions with a reactive center of the type H-H-N, N-H-N and C-H-N respectively.

Chemical upcycling of polyethylene, polypropylene, and mixtures to high-value surfactants

Conversion of plastic wastes to fatty acids is an attractive means to supplement the sourcing of these high-value, high-volume chemicals. We report a method for transforming polyethylene (PE) and polypropylene (PP) at ~80% conversion to fatty acids with number-average molar masses of up to ~700 and 670 daltons, respectively. The process is applicable to municipal PE and PP wastes and their mixtures. Temperature-gradient thermolysis is the key to controllably degrading PE and PP into waxes and inhibiting the production of small molecules. The waxes are upcycled to fatty acids by oxidation over manganese stearate and subsequent processing. PP ꞵ-scission produces more olefin wax and yields higher acid-number fatty acids than does PE ꞵ-scission. We further convert the fatty acids to high-value, large-market-volume surfactants. Industrial-scale technoeconomic analysis suggests economic viability without the need for subsidies.

The merit of pressure dependent kinetic modelling in steam cracking

Renewable cracking feedstocks from plastic waste and the need for novel reactor designs related to electrification of steam crackers drives the development of accurate and fundamental kinetic models for this process, despite its large scale implementation for more than half a century. Pressure dependent kinetics have mostly been omitted in fundamental steam cracking models, while they are crucial in combustion models. Therefore, we have assessed the importance of pressure dependent kinetics for steam cracking via in-depth modelling and experimental studies. In particular we have studied the influence of considering fall-off on the product yields for ethane and propane steam cracking. A high-pressure limit fundamental kinetic model is generated, based on quantum chemical data and group additive values, and supplemented with literature values for pressure dependent kinetic parameters for β-scission reactions and homolytic bond scissions of C₂ and C₃ species. Model simulations with high-pressure limit rate coefficients and pressure dependent kinetics are compared to new experimental measurements. Steam cracking experiments for pure ethane and propane feeds are performed on a tubular bench-scale reactor at 0.17 MPa and temperatures ranging from 1058 to 1178 K. All important product species are identified using a comprehensive GC × GC-FID/q-MS. For homolytic bond scissions, the inclusion of pressure dependent kinetics has a significant effect on the conversion profile for ethane steam cracking. On the other hand, pressure dependence of C₂ β-scissions significantly influences conversion and product species profiles for both ethane and propane steam cracking. The C₃ β-scissions pressure dependence has a negligible effect in ethane steam cracking, while for propane steam cracking the effect is non-negligible on the product species profiles.

Electrochemically-mediated C–H functionalization of allenes and 1,3-dicarbonyl compounds to construct tetrasubstituted furans

We reported an electrocatalytic C–H activation method to construct novel highly functionalized tetrasubstituted furan derivatives, which uses allenes and 1,3-dicarbonyl compounds as substrates.

High-Pressure-Limit and Pressure-Dependent Rate Rules for Unimolecular Reactions Related to Hydroperoxy Alkyl Radicals in Normal Alkyl Cyclohexane Combustion. 1. Concerted HO2 Elimination Reaction Class and β-Scission Reaction Class

The reactions of the concerted HO₂ elimination from alkyl peroxy radicals and the β-scission of the C-OOH bond from hydroperoxy alkyl radicals, which lead to the formation of olefins and HO₂ radicals, are two important reaction classes that compete with the second O₂ addition step of hydroperoxy alkyl radicals, which are responsible for the chain branching in the low-temperature oxidation of normal alkyl cycloalkanes. These two reaction classes are also believed to be responsible for the negative temperature coefficient behavior due to the formation of the relatively unreactive HO₂ radical, which has the potential to inhibit ignition of normal alkyl cycloalkanes. In this work, the kinetics of the above two reaction classes in normal alkyl cycloalkanes are studied, where reactions in the concerted elimination class are divided into subclasses depending upon the types of carbons from which the H atom is eliminated and the positions of the reaction center (on the alkyl side chain or on the cycle), and the reactions in the β-scission reaction class are divided into subclasses depending upon the types of the carbons on which the radical is located and the positions of the reaction center. Energy barriers by using quantum chemical methods at the CBS-QB3 level, high-pressure-limit rate constants by using canonical transition state theory, and pressure-dependent rate constants at pressures from 0.01 to 100 atm by using Rice-Ramsberger-Kassel-Marcus/Master Equation theory are calculated for a representative set of reactions from methyl cyclohexane to n-butyl cyclohexane in each subclass, from which high-pressure-limit rate rules and pressure-dependent rate rules for each subclass are derived from the average rate constants of reactions within each subclass. A comparison of the rate constants for the reactions in the two reaction classes calculated in this work is made with the rate constants of the same reactions from available mechanisms published in the literature, where most of the rate constants are approximately estimated from analogous reactions in alkanes or small alkyl cyclohexanes, and it is found that a large difference may exist between them, indicating that the present work, which provides more accurate kinetic parameters and reasonable rate rules for these reaction classes, can be helpful to construct higher-accuracy mechanism models for normal alkyl cyclohexane combustion.

Electronic structure-based rate rules for ipso addition-elimination reactions on mono-aromatic hydrocarbons with single and double OH/CH3/OCH3/CHO/C2H5 substituents: a systematic theoretical investigation

The recent interest in bio-oils combustion and the key role of mono-aromatic hydrocarbons (MAHs) in existing kinetic frameworks, both in terms of poly-aromatic hydrocarbons growth and surrogate fuels formulation, motivates the current systematic theoretical investigation of one of the relevant reaction classes in MAHs pyrolysis and oxidation: ipso substitution by hydrogen. State-of-the-art theoretical methods and protocols implemented in automatized computational routines allowed to investigate 14 different potential energy surfaces involving MAHs with hydroxy and methyl single (phenol and toluene) and double (o-,m-,p-C6H4(OH)2, o-,m-,p-CH3C6H4OH, and o-,m-,p-C6H4(CH3)2) substituents, providing rate constants for direct implementation in existing kinetic models. The accuracy of the adopted theoretical method was validated by comparison of the computed rate constants with the available literature data. Systematic trends in energy barriers, pre-exponential factors, and temperature dependence of the Arrhenius parameters were found, encouraging the formulation of rate rules for ipso substitutions on MAHs. The rules here proposed allow to extrapolate from a reference system the necessary activation energy and pre-exponential factor corrections for a large number of reactions from a limited set of electronic structure calculations. We were able to estimate rate constants for other 63 ipso addition-elimination reactions on di-substituted MAHs, reporting in total 75 rate constants for ipso substitution reactions o-,m-,p-R'C6H4R + → C6H5R + ', with R,R' = OH/CH3/OCH3/CHO/C2H5, in the 300-2000 K range. Additional calculations performed for validation showed that the proposed rate rules are in excellent agreement with the rate constants calculated using the full computational protocol in the 500-2000 K range, generally with errors below 20%, increasing up to 40% in a few cases. The main results of this work are the successful application of automatized electronic structure calculations for the derivation of accurate rate constants for ipso substitution reactions on MAHs, and an efficient and innovative approach for rate rules formulation for this reaction class.

Kinetic modeling of methyl pentanoate pyrolysis based on ab initio calculations

Recently, methyl pentanoate (MP) was proposed as a viable biodiesel surrogate to petroleum-based fuels. To better understand the pyrolysis chemistry of MP, the unimolecular decomposition kinetics of MP is theoretically investigated on the basis of ab initio calculations; ten primary channels, including four intramolecular H-shifts and six C-C and C-O bond fissions, are identified. The geometries are optimized at the M06-2X/cc-pVTZ level of theory, and accurate barrier heights are determined using the DLPNO-CCSD(T)/CBS(T-Q) method, which shows a good performance against the CCSD(T)/CBS(T-Q) method with an uncertainty of 0.5 kcal mol-1 for small methyl esters. The atomization enthalpy method is adopted to obtain the thermodynamics of involved species. The Rice-Ramsperger-Kassel-Marcus/master equation theory coupled with one-dimensional hindered rotor approximation is employed to calculate the phenomenological rate constants at 500-2000 K and 0.01-100 atm. The branching ratio analysis indicates that two reactions, MP ↔ CH3OC([double bond, length as m-dash]O)CH3 + CH2CHCH3 and MP ↔ CH3OC([double bond, length as m-dash]O)CH2 + CH2CH2CH3, are the dominant channels at low and high temperatures, respectively. The model from Diévart et al. [Proc. Combust. Inst., 2013, 34(1), 821-829] is updated with our calculations, and the modified model can yield a better prediction in reproducing the ignition delay times of MP at high temperatures. This work provides a comprehensive investigation of MP unimolecular decomposition, and can serve as a prototype for understanding the pyrolysis of larger alkyl esters.

An ab Initio/Transition State Theory Study of the Reactions of Ċ5H9 Species of Relevance to 1,3-Pentadiene, Part II: Pressure Dependent Rate Constants and Implications for Combustion Modeling

The temperature- and pressure-dependence of rate constants for several radicals and unsaturated hydrocarbons reactions (1,3-C₅H₈/1,4-C₅H₈/cyC₅H₈ + Ḣ, C₂H₄ + Ċ₃H₅-a, C₃H₆ + Ċ₂H₃) are analyzed in this paper. The abstraction reactions of these systems are also calculated and compared with available literature data. Ċ₅H₉ radicals can be produced via Ḣ atom addition reactions to the pentadiene isomers and cyclopentene, and also by H-atom abstraction reactions from 1- and 2-pentene and cyclopentane. Comprehensive Ċ₅H₉ potential energy surface (PES) analyses and high-pressure limiting rate constants for related reactions have been explored in part I of this work ( J. Phys. Chem. A 2019, 123 (22), 9019-9052). In this work, a chemical kinetic model is constructed based on the computed thermochemistry and high-pressure limiting rate constants from part I, to further understand the chemistry of different C₅H₈ molecules. The most important channels for these addition reactions are discussed in the present work based on reaction pathway analyses. The dominant reaction pathways for these five systems are combined together to generate a simplified Ċ₅H₉ PES including nine reactants, 25 transition states (TSs), and nine products. Spin-restricted single point energies are calculated for radicals and TSs on the simplified PES at the ROCCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). Temperature- and pressure-dependent rate constants are calculated using RRKM theory with a Master Equation analysis, with restricted energies used for minima on the simplified Ċ₅H₉ PES and unrestricted energies for other species, over a temperature range of 300-2000 K and in the pressure range 0.01-100 atm. The rate constants calculated are in good agreement with those in the literature. The chemical kinetic model is updated with pressure-dependent rate constants and is used to simulate the species concentration profiles for Ḣ atom addition to cyclopentane and cyclopentene. Through detailed analyses and comparisons, this model can reproduce the experimental measurements of species qualitatively and quantitatively with reasonably good agreement.

Pressure-dependent rate rules for cycloaddition, intramolecular H-shift, and concerted elimination reactions of alkenyl peroxy radicals at low temperature

The reactions of cycloaddition, intramolecular H-shift and concerted elimination of alkenyl peroxy radicals are three kinds of important reactions in the low temperature combustion of alkenes. In this study, the cycloaddition reactions are divided into classes considering endo-cycloaddition, exo-cycloaddition and the size of the transition states; the intramolecular H-shift reactions are divided into classes depending upon the ring size of the transition states and the type of C-H bonds from which the hydrogen atom is transferred; the concerted elimination reactions are divided into classes according to the type of H-C_βC_αOO bond that is broken. All geometry optimizations are performed at the B3LYP/6-31G(2df,p) level. With the electronic structure calculations being performed using the composite Gaussian-4 (G4) method, high pressure limit rate constants and pressure-dependent rate constants at pressures varying from 0.01 to 100 atm are calculated by using canonical transition state theory and the Rice-Ramsberger-Kassel-Marcus/master equation method, respectively. All rate constants are given in the form of the modified Arrhenius expression. The high pressure limit rate rules and the pressure-dependent rate rules are derived by averaging the rate constants of a representative set of reactions in each class. The results show that the rate rules for these three classes of reactions have a large uncertainty and the impact of the pressure on the rate constants increases as temperature increases.

Ab initio derived group additivity model for intramolecular hydrogen abstraction reactions

A set of group additivity values for intramolecular hydrogen abstraction reactions of alkanes, alkenes and alkynes is reported. Calculating 448 reaction rate coefficients at the CBS-QB3 level of theory for 1-2 up to 1-7 hydrogen shift reactions allowed the estimation of ΔGAV° values for 270 groups. The influence of substituents on (1) the attacking radical, (2) the attacked carbon atom, and (3) the carbon chain between the attacking and attacked reactive atom has been systematically studied. Substituents have been varied between hydrogen atoms and sp3, sp2 and sp hybridized carbon atoms. It has been assumed that substituents further away from the reactive atoms or their connecting carbon chain have negligible influences on the kinetics. This group additivity model is applicable to a wide variety of reactions in the 300-1800 K temperature range. Correlations for tunneling coefficients have been generated which are complementary to the ΔGAV°'s to obtain accurate rate coefficients without the need for imaginary frequencies or electronic energies of activation. These correlations depend on the temperature and activation energy of the exothermic step. The group additivity model has been successfully applied to a test set of reactions also calculated at the CBS-QB3 level of theory. A mean absolute deviation of 1.18 to 1.71 has been achieved showing a good overall accuracy of the model.

Ab initio study of the influence of resonance stabilization on intramolecular ring closure reactions of hydrocarbon radicals

The cyclization reactions of dieneyl radicals provide a low energy route to the formation of molecular weight growth products.

Unravelling the impact of hydrocarbon structure on the fumarate addition mechanism--a gas-phase ab initio study

The fumarate addition reaction mechanism is central to the anaerobic biodegradation pathway of various hydrocarbons, both aromatic (e.g., toluene, ethyl benzene) and aliphatic (e.g., n-hexane, dodecane). Succinate synthase enzymes, which belong to the glycyl radical enzyme family, are the main facilitators of these biochemical reactions. The overall catalytic mechanism that converts hydrocarbons to a succinate molecule involves three steps: (1) initial H-abstraction from the hydrocarbon by the radical enzyme, (2) addition of the resulting hydrocarbon radical to fumarate, and (3) hydrogen abstraction by the addition product to regenerate the radical enzyme. Since the biodegradation of hydrocarbon fuels via the fumarate addition mechanism is linked to bio-corrosion, an improved understanding of this reaction is imperative to our efforts of predicting the susceptibility of proposed alternative fuels to biodegradation. An improved understanding of the fuel biodegradation process also has the potential to benefit bioremediation. In this study, we consider model aromatic (toluene) and aliphatic (butane) compounds to evaluate the impact of hydrocarbon structure on the energetics and kinetics of the fumarate addition mechanism by means of high level ab initio gas-phase calculations. We predict that the rate of toluene degradation is ∼100 times faster than butane at 298 K, and that the first abstraction step is kinetically significant for both hydrocarbons, which is consistent with deuterium isotope effect studies on toluene degradation. The detailed computations also show that the predicted stereo-chemical preference of the succinate products for both toluene and butane are due to the differences in the radical addition rate constants for the various isomers. The computational and kinetic modeling work presented here demonstrates the importance of considering pre-reaction and product complexes in order to accurately treat gas phase systems that involve intra and inter-molecular non-covalent interactions.

Reactions of allylic radicals that impact molecular weight growth kinetics

The reactions of allylic radicals have the potential to play a critical role in molecular weight growth (MWG) kinetics during hydrocarbon oxidation and/or pyrolysis.

Investigation of free-radical copolymerization propagation kinetics of vinyl acetate and methyl methacrylate

The free-radical copolymerization propagation kinetics of vinyl acetate (VAc) and methyl methacrylate (MMA) at 50 degrees C were investigated through an experimental study combined with a computational analysis based on quantum chemistry. Copolymer composition data, obtained using pulsed laser polymerization followed by size exclusion chromatography (PLP-SEC) and proton nuclear magnetic resonance (NMR), were well represented by the terminal model using monomer reactivity ratios obtained with the computational approach (r(VAc) = 0.001 and r(MMA) = 27.9). Concerning the composition-averaged copolymerization propagation rate coefficient k(p,cop), the differences between the terminal model and the implicit penultimate unit effect (IPUE) model (s(MMA) = 0.544 and s(VAc) = 0.173) are small for VAc/MMA, with the terminal model sufficient to describe the experimental k(p,cop) data measured by PLP-SEC. Monomer and radical charge distributions determined computationally are used to explain the reactivity exhibited by the VAc/MMA system.