H-Abstraction reactions by OH, HO2, O, O2 and benzyl radical addition to O2 and their implications for kinetic modelling of toluene oxidation

Catalytic Aerobic Carbooxygenation for the Construction of Vicinal Tetrasubstituted Centers: Application to the Synthesis of Hexasubstituted γ-Lactones

Strategic design for the construction of contiguous tetrasubstituted carbon centers represents a daunting challenge in synthetic organic chemistry. Herein, we report a combined experimental and computational investigation aimed at developing catalytic aerobic carbooxygenation, involving the intramolecular addition of tertiary radicals to geminally disubstituted alkenes, followed by aerobic oxygenation. This reaction provides a straightforward route to various α,α,β,β-tetrasubstituted γ-lactones, which can be readily transformed into hexasubstituted γ-lactones through allylation/translactonization. Computational analysis reveals that the key mechanistic foundation for achieving the developed aerobic carbooxygenation involves the design of endothermic (energetically uphill) C-C bond formation followed by exothermic (energetically downhill) oxygenation. Furthermore, we highlight a unique fluorine-induced stereoelectronic effect that stabilizes the endothermic stereodetermining transition state.

Crossed molecular beam experiments and theoretical simulations on the multichannel reaction of toluene with atomic oxygen

Despite extensive experimental and theoretical studies on the kinetics of the O(3P) + C7H8 (toluene) reaction and a pioneering crossed molecular beam (CMB) investigation, the branching fractions (BFs) of the CH3C6H4O(methylphenoxy) + H, C6H5O(phenoxy) + CH3, and spin-forbidden C5H5CH3 (methylcyclopentadiene) + CO product channels remain an open question, which has hampered the proper inclusion of this important reaction in the chemical modelling of various chemical environments. We report a CMB study with universal soft electron-ionization mass-spectrometric detection of the reactions O(3P,1D) + toluene at the collision energy of 34.7 kJ mol-1. From CMB data we have inferred the reaction dynamics and quantified the BFs of the primary products and the role of intersystem crossing (ISC). The CH3-elimination channel dominates (BF = 0.69 ± 0.22) in the O(3P) reaction, while the H-displacement and CO-formation channels are minor (BF = 0.22 ± 0.07 and 0.09 ± 0.05, respectively), with ISC accounting for more than 50% of the reactive flux. Synergistic transition-state theory (TST)-based master equation simulations including nonadiabatic TST on ab initio coupled triplet/singlet potential energy surfaces were employed to compute the product BFs and assist in the interpretation of the CMB results. In the light of the good agreement between the theoretical predictions for the O(3P) + toluene reaction and the CMB results as well as the absolute rate constant as a function of temperature (T) (from literature), the so-validated computational methodology was used to predict channel-specific rate constants as a function of T at 1 atm.

Mechanism, thermochemistry, and kinetics of formation of PhC(O)H and HOPhC(O)H during unimolecular decomposition of P-PhC(O2•)HPhOH

CONTEXT

Bisphenols are one of the main components of bio-oil, produced during the pyrolysis of lignin-containing biomass. Synthetic bisphenols are used in polycarbonate plastics, epoxy resins, and thermal papers. Their mechanism of oxidation is important for the determination of the fire safety of these materials and the possibility of using them as additives in fuels for the decrease and description of ignition delays, as well as for the determination of their health risk assessment in medicine. One representative of bisphenols is p-benzylphenol (p-PhCH2PhOH), which is formed during the fast pyrolysis of lignine-containing biomass. Its thermochemistry of oxidation has been partially studied previously. It is shown that the reaction of chain oxidation of p-PhCH2PhOH is thermochemically favorable at low temperatures. However, these studies consider only two pathways of this reaction: (1) the chain oxidation of RH by RO2• and (2) the tautomerization of R'HO2• to R'O2H with following production of R'O• and OH radicals. At the same time, the reactions of intramolecular rearrangement of RO2•, produced PhC(O)H and •PhOH or HOPhC(O)H and •Ph, are not reported but can be an important part of its oxidation mechanism.

METHODS

The five DFT (M06-2X (i = 1), B3LYP (i = 2), wB97XD (i = 3), M08HX (i = 4), MN15 (i = 5)) approaches with 6-311 +  + G(d,p) basis set are used for the determination of standard enthalpies of atomization (ΔraH°(Xi)) of considered compounds (molecules, radicals, and transition states). These values of ΔraH°(Xi) are corrected using the empirical linear calibration dependencies, reported previously. The different calibration dependencies are used for the hydrocarbons (including the aromatics and simple oxygenated derivatives) and for the peroxides. The corrected values of ΔraH°(Xi, CORR) are used according to Hess's law for the determination of ΔfH°(Xi, CORR). The most consistent values of ΔfH°(X, MEAN) are derived from the coordination of the values of ΔfH°(Xi,CORR) using the intersection of their values of standard deviations (3SDi). These values of ΔfH°(X, MEAN), as well as the B3LYP values of S°(X), which are accounting the frequencies correction and internal rotations, as well as their temperature dependencies, are used for the determination of thermochemistry of considered reactions and of the calculation, within transition state theory (TST), of the values of high pressure limits of the rate constant. The values of H°(Xi), S°(Xi), and G°(Xi) are calculated using the Gaussian 16w program. The considered mechanism is prepared using ISIS/Draw package. The temperature dependencies of thermochemical properties and the values of rate constants are determined using the ChemRate program (v.1.5). The optimized structures are visualized using the Chemcraft package.

Mechanism of formation of p-benzylenephenol peroxide radical (p-PhC(O2•)HPhOH)

CONTEXT

The reactions of radicals with O2 play the important role in the biological, medicinal, and industrial processes. The mechanism of this reaction is studied previously for the alkane, alkene, alkyne, phenols, and close-related radicals. According to these studies, the formation of intermediates in these reactions is predicted only for the aromatic radicals. Thus, the Van der Waals complexes of O2 with phenyl or benzyl radicals are predicted, as well as the π-π cluster for benzene. However, the possibility of the formation of such intermediate π-π clusters in the case of bisphenol radicals and the thermochemistry of its formation is not studied. Bisphenols are one of the main components of bio-oil, produced during pyrolysis of lignin-contained biomass. Synthetic bisphenols are used in polycarbonate plastics, epoxy resins, and thermal papers. Their mechanism of oxidation is important for the determination of fire safety of these materials, the possibility of using them as additives for fuels for the decreasing and the description of the ignition delays, as well as for the determination of its health risk assessment in medicine.

METHODS

The five DFT (M06-2X (i = 1), B3LYP (i = 2), wB97XD (i = 3), M08HX (i = 4), MN15 (i = 5)) approaches with 6-311 +  + G(d,p) basis set are used for the determination of standard enthalpies of atomization (ΔraHo(Yi)) of considered compounds (molecules, radicals, and transition states). These values of ΔraHo(Yi) are corrected using the empirical linear calibration dependencies, reported previously. The different calibration dependencies are used for the hydrocarbons (including the aromatics and simple oxygenated derivatives) and for the peroxides. The corrected values of ΔraHo(Yi, CORR) are used according to Hess's law for the determination of ΔfHo(Yi, CORR). The most consistent values of ΔfHo(Y, MEAN) are derived from the coordination of the values of ΔfHo(Yi, CORR) using the intersection of their values of standard deviations (3SDi). These values of ΔfHo(Y, MEAN), as well as the B3LYP values of So(Y), which are accounting the frequency correction and internal rotations, as well as their temperature dependencies, are used for the determination of thermochemistry of considered reactions and of the calculation, within transition state theory (TST), of the values of high-pressure limits of the rate constant. The values of Ho(Yi), So(Yi), and Go(Yi) are calculated using the Gaussian 16w program. The temperature dependencies of thermochemical properties and the values of rate constants are determined using the ChemRate program (v.1.5). The optimized structures are visualized using the Chemcraft.

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.

Experimental and numerical studies on the thermal nonequilibrium behaviors of CO with Ar, He, and H2

The time-dependent rotational and vibrational temperatures were measured to study the shock-heated thermal nonequilibrium behaviors of CO with Ar, He, and H2 as collision partners. Three interference-free transition lines in the fundamental vibrational band of CO were applied to the fast, in situ, and state-specific measurements. Vibrational relaxation times of CO were summarized over a temperature range of 1110-2820 K behind reflected shocks. The measured rotational temperature instantaneously reached an equilibrium state behind shock waves. The measured vibrational temperature experienced a relaxation process before reaching the equilibrium state. The measured vibrational temperature time histories were compared with predictions based on the Landau-Teller model and the state-to-state approach. The state-to-state approach treats the vibrational energy levels of CO as pseudo-species and accurately describes the detailed thermal nonequilibrium processes behind shock waves. The datasets of state-specific inelastic rate coefficients of CO-Ar, CO-He, CO-CO, and CO-H2 collisions were calculated in this study using the mixed quantum-classical method and the semiclassical forced harmonic oscillator model. The predictions based on the state-to-state approach agreed well with the measured data and nonequilibrium (non-Boltzmann) vibrational distributions were found in the post-shock regions, while the Landau-Teller model predicted slower vibrational temperature time histories than the measured data. Modifications were applied to the Millikan-White vibrational relaxation data of the CO-Ar and CO-H2 systems to improve the performance of the Landau-Teller model. In addition, the thermal nonequilibrium processes behind incident shocks, the acceleration effects of H2O on the relaxation process of CO, and the characterization of vibrational temperature were highlighted.

Multiwavelength Speciation in Pyrolysis of n-Pentane and Experimental Determination of the Rate Coefficient of nC5H12 = nC3H7 + C2H5 in a Shock Tube

We report the application of a multiwavelength speciation strategy to the study of n-pentane (nC₅H₁₂) pyrolysis behind reflected shock waves in a shock tube. Experiments were conducted with 2% nC₅H₁₂/0.8%CO₂/Ar (by mole) between 1150 and 1520 K in the pressure range of 1-2 atm. Utilization of laser absorption spectroscopy at eight wavelengths allowed time-resolved measurements of n-pentane, ethylene, methane, heavy alkenes, and temperature. The measured time histories were compared against the predictions of four recently developed chemical kinetic models for heavy hydrocarbons. It was found that none of the models reconciled the measured species time histories simultaneously. Sensitivity analysis was conducted to identify key reactions influencing the evolution of ethylene and other major pyrolysis products. The analysis revealed that the unimolecular decomposition of n-pentane into n-propyl and ethyl radicals has a dominating influence over the evolution of ethylene in the temperature range of 1150-1450 K. The rate coefficient of this reaction was then adjusted to match the measured ethylene time histories for each experiment. The rate coefficients thus determined, were fit against temperature using an Arrhenius expression given by k₁(T) = 3.5 × 10¹⁴ exp(-67.2 kcal/RT) s^-1. The average overall 2σ uncertainty of the measured rate coefficient was found to be ±35%, resulting primarily from uncertainties in the rate coefficients of secondary reactions. The measured rate coefficient, when used with the models, leads to a significant improvement in the prediction of species time histories. Further improvements in the model are possible if the rate coefficients of relevant reactions pertaining to small hydrocarbon chemistry are determined with an improved accuracy, and less uncertainty. To the best knowledge of the authors, this is the first experimental determination of the rate coefficient of C₅H₁₂ → nC₃H₇ + C₂H_5.

A Theoretical Study of H-Abstractions of Benzaldehyde by H, O3(P), 3O2, OH, HO2, and CH3 Radicals: Ab Initio Rate Coefficients and Their Uncertainty Quantification

Benzaldehyde is a vital intermediate during the oxidation of toluene and oxygenated aromatics, but benzaldehyde's combustion chemistry currently gains little attention. The H-atom abstraction of benzaldehyde is one critical reaction class, yet its rate coefficients rely on the analogy of molecular similarity. For this reason, we have employed the ab initio calculations at the CCSD(T)/cc-pVTZ//M06-2X/6-311+g(d,p) level of theory, along with the RRKM master equation to predict the rate coefficients of H-abstraction reactions of benzaldehyde by H, O³(P), ³O₂, OH, HO₂, and CH₃ radicals. The calculated rate coefficients of benzaldehyde + OH generally agree with literature measurements. From the branching ratio analysis, all the H-abstractions from ring sites can be neglected in kinetic model construction except by OH radical at temperature above 1700 K, where the ring sites contribute a relatively large branching ratio to consume benzaldehyde. A random sampling method has been used to estimate the global uncertainty in the calculated rate coefficients. The logarithm of uncertainty is proportional to the reciprocal temperature, and the global uncertainty is primarily derived from the energy errors and is almost independent of the attacking species. Comparison between benzaldehyde and acetaldehyde reveals that their rate coefficients are consistent only in H-abstractions by H and HO₂ but are conflictive in other reactions, especially in the case of OH. By incorporating the new calculations, existing models show a faster prediction in autoignition delay times of benzaldehyde. This study reports the first high-level ab initio calculations for the H-atom abstraction reaction class of benzaldehyde. It is necessary for the future comprehensive kinetic modeling of aromatic aldehyde.

Radical-Mediated C-C Coupling of Alcohols Induced by Plasmonic Hot Carriers

The C-C coupling reactions of aliphatic alcohols to aromatics and larger-mass compounds have large endothermicities and activation energies, calling for catalysts operating at high temperatures. Here, we demonstrate that plasmon-excited nanoparticles catalyze the C-C coupling of aliphatic alcohols at room temperature to produce polyaromatic hydrocarbons and graphene oxide. The conversion is quenched by radical and electron scavengers and by the surface passivation of metals, suggesting that the reaction proceeds through alkoxy, peroxyl, hydroxyalkyl, and alkyl radical intermediates created by the metal to molecule transfer of plasmonic hot carriers. Besides being the first realization of C-C coupling of aliphatic alcohols at room temperature, the result constitutes a rare example of an endothermic plasmon-induced reaction producing new bonds and a new method for photogenerating graphene derivatives. More importantly, the result demonstrates the facile generation of organic radicals directly from alcohols, which may be used as precursors for radical-based organic reactions.

Investigating the aging mechanism of asphaltene and its dependence on environmental factors through AIMD simulations and DFT calculations

To understand the complex aging mechanism of asphalt and its dependence on environmental factors, the chemical reactivity of asphaltene during aging under different environmental conditions was studied through first-principles molecular simulations and density functional theory calculations. The aging of asphaltene was demonstrated to involve a series of subreactions along different pathways on the asphaltene molecules, including hydrogen abstraction from carbon, formation of polar groups, aromatization of cycloalkanes, and homolysis of side chains. These subreactions occurred with different free-energy barriers and, therefore, had different kinetic rates. Asphaltene aging was found to be slightly accelerated in the presence of water owing to the improved electron transfer ability of the asphaltene molecule in an aqueous solvent. Under ultraviolet radiation, the asphaltene molecule transitioned to an excited state with an excitation energy of 348.7 kJ/mol, significantly increasing its aging rate. This work bridges the gap between electronic-scale modeling and diversified experimental observations related to asphalt aging and is expected to provide theoretical guidance for strategies to prevent or delay the aging-induced failure of asphalt pavements.

An Experimental Kinetics Study of Isopropanol Pyrolysis and Oxidation behind Reflected Shock Waves

Isopropanol has potential as a future bio-derived fuel and is a promising substitute for ethanol in gasoline blends. Even so, little has been done in terms of high-temperature chemical kinetic speciation studies of this molecule. To this end, experiments were conducted in a shock tube using simultaneous CO and H2O laser absorption measurements. Water and CO formation during isopropanol pyrolysis was also examined at temperatures between 1127 and 2162 K at an average pressure of 1.42 atm. Species profiles were collected at temperatures between 1332 and 1728 K and at an average pressure of 1.26 atm for equivalence ratios of 0.5, 1.0, and 2.0 in highly diluted mixtures of 20% helium and 79.5% argon. Species profiles were also compared to four modern C3 alcohol mechanisms, including the impact of recent rate constant measurements. The Li et al. (2019) and Saggese et al. (2021) models both best predict CO and water production under pyrolysis conditions, while the AramcoMech 3.0 and Capriolo and Konnov models better predict the oxidation experimental profiles. Additionally, previous studies have collected ignition delay time (τign) data for isopropanol but are limited to low pressures in highly dilute mixtures. Therefore, real fuel–air experiments were conducted in a heated shock tube with isopropanol for stoichiometric and lean conditions at 10 and 25 atm between 942 and 1428 K. Comparisons to previous experimental results highlight the need for real fuel–air experiments and proper interpretation of shock-tube data. The AramcoMech 3.0 model over predicts τign values, while the Li et al. model severely under predicts τign. The models by Capriolo and Konnov and Saggese et al. show good agreement with experimental τign values. A sensitivity analysis using these two models highlights the underlying chemistry for isopropanol combustion at 25 atm. Additionally, modifying the Li et al. model with a recently measured reaction rate shows improvement in the model’s ability to predict CO and water profiles during dilute oxidation. Finally, a regression analysis was performed to quantify τign results from this study.

π-Diradical Aromatic Soot Precursors in Flames

Soot emitted from incomplete combustion of hydrocarbon fuels contributes to global warming and causes human disease. The mechanism by which soot nanoparticles form within hydrocarbon flames is still an unsolved problem in combustion science. Mechanisms proposed to date involving purely chemical growth are limited by slow reaction rates, whereas mechanisms relying on solely physical interactions between molecules are limited by weak intermolecular interactions that are unstable at flame temperatures. Here, we show evidence for a reactive π-diradical aromatic soot precursor imaged using non-contact atomic force microscopy. Localization of π-electrons on non-hexagonal rings was found to allow for Kekulé aromatic soot precursors to possess a triplet diradical ground state. Barrierless chain reactions are shown between these reactive sites, which provide thermally stable aromatic rim-linked hydrocarbons under flame conditions. Quantum molecular dynamics simulations demonstrate physical condensation of aromatics that survive for tens of picoseconds. Bound internal rotors then enable the reactive sites to find each other and become chemically cross-linked before dissociation. These species provide a rapid, thermally stable chain reaction toward soot nanoparticle formation and could provide molecular targets for limiting the emission of these toxic combustion products.

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.

Rational Design and Testing of Anti-Knock Additives

An innovative and informed methodology for the rational design and testing of anti-knock additives is reported. Interaction of the additives with OH● and HO2● is identified as the key reaction pathway by which non-metallic anti-knock additives are proposed to operate. Based on this mechanism, a set of generic design criteria for anti-knock additives is outlined. It is suggested that these additives should contain a weak X-H bond and form stable radical species after hydrogen atom abstraction. A set of molecular structural, thermodynamic, and kinetic quantities that pertain to the propensity of the additive to inhibit knock by this mechanism are identified and determined for a set of 12 phenolic model compounds. The series of structural analogues was carefully selected such that the physical thermodynamic and kinetic quantities could be systematically varied. The efficacy of these molecules as anti-knock additives was demonstrated through the determination of the research octane number (RON) and the derived cetane number(DCN), measured using an ignition quality tester (IQT), of a RON 95 gasoline treated with 1 mole % of the additive. The use of the IQT allows the anti-knock properties of potential additives to be studied on one tenth of the scale, compared to the analogous RON measurement. Using multiple linear regression, the relationship between DCN/RON and the theoretically determined quantities is studied. The overall methodology reported is proposed as an informed alternative to the non-directed experimental screening approach typically adopted in the development of fuel additives.

Reinvestigation of the Deceptively Simple Reaction of Toluene with OH and the Fate of the Benzyl Radical: The "Hidden" Routes to Cresols and Benzaldehyde

In a previous work, we have investigated the initial steps of the reaction of toluene with the hydroxyl radical using several quantum chemical approaches including density functional and composite post-Hartree-Fock models. Comparison of H-abstraction from the methyl group and additions at different positions of the phenyl ring showed that the former reaction channel is favored at room temperature. This conclusion appears at first sight incompatible with the experimental observation of a lower abundance of the product obtained from abstraction (benzaldehyde) with respect to those originating from addition (cresols). Further reactions of the intermediate radicals with oxygen, water, and additional OH radicals are explored in this paper through theoretical calculations on more than 120 species on the corresponding potential energy surface. The study of the addition reactions, to obtain the cresols through hydroxy methylcyclodienyl intermediate radicals, showed that only in the case of o-cresol the reaction proceeds by addition of O₂ to the ring, internal H-transfer, and hydroperoxyl abstraction and not through direct H-abstraction. For both p- and m-cresol, instead, the reaction occurs through a higher-energy direct H-abstraction, thus explaining in part the observed larger concentration of the ortho isomer in the final products. It was also found that the benzyl radical, formed by H-abstraction from the methyl group, is able to react further if additional OH is present. Two reaction paths leading to o-cresol, two leading to p-cresol, and one leading to m-cresol were determined. Moreover, in this situation, the benzyl radical is predicted to produce benzyl alcohol, as was found in some experiments. The commonly accepted route to benzaldehyde was found to be not the energetically favored one. Instead, a route leading to the benzoyl radical (and ultimately to benzoic acid) with the participation of one water molecule was clearly more favorable, both thermodynamically and kinetically.

OH-Initiated Reactions of para-Coumaryl Alcohol Relevant to the Lignin Pyrolysis. Part III. Kinetics of H-Abstraction by H, OH, and CH3 Radicals

Lignin is the most complex component of biomass, and development of a detailed chemical kinetic model for biomass pyrolysis mainly relies on the understanding of the lignin decomposition kinetics. para-Coumaryl alcohol (p-CMA, HOPh-CH═CH-CH₂OH), the focus of our analysis, is the simplest of the lignin monomers (monolignols) containing a typical side-chain double bond and both alkyl- and phenolic-type OH-groups. In parts I and II of our work (Asatryan, R. J. Phys. Chem. A 2019, 123, 2570-2585; Hudzik, J. M. J. Phys. Chem. A 2020, current issue), we created a detailed potential energy surface (PES) and performed a kinetic analysis of chemically activated, unimolecular, and bimolecular reactions pathways for p-CMA + OH. Reaction pathways analyzed include dissociation, intramolecular abstraction, group transfer, and elimination processes. The α- and β-carbon addition reactions generate 1,3- (RA1) and 1,2-diol (RB1) adduct radicals, respectively. Well depths are approximately 29 and 41 kcal/mol below the p-CMA + OH entrance level. Kinetic analysis aides in determining the major pathways for our conventional and fractional pyrolysis experiments. The current paper focuses on the H-abstraction reactions via H, OH, and CH₃ light ("pool") radicals from p-CMA. The thermochemical properties of all stable, radical, and transition-state species were determined using the ωB97XD density functional theory (DFT) and higher-level CBS-QB3 composite methods. Barrier heights from the prereaction complexes, for OH-radical abstractions, to the transition states for the propanoid side chain are compared to the model H-abstraction reactions of allyl alcohol (AA) with OH and p-CMA with H and CH₃ radicals. The lowest-energy, most stable, p-CMA radical formed is at the C9 allylic position (p-CMA-C9j) with exothermicity of 26.63, 41.32, and 27.34 kcal/mol for H, OH, and CH₃, respectively. For OH-radical abstraction at this position, our findings are consistent with corresponding data on AA + OH at 37.44 kcal/mol and similar to that of RB1. A similar stable radical with an exothermicity of 34.95 kcal/mol occurs for the phenol hydroxyl group, generating the p-CMA-O4j radical. H-abstraction pathways are considered in relation to other major pathways previously considered for p-CMA + OH reactions including H-atom shifts, dehydration, and β-scission reactions. Derived rate coefficients for substituted phenols can be utilized in detailed kinetic models for lignin/biomass pyrolysis.

Role of ring-enlargement reactions in the formation of aromatic hydrocarbons

Ring-enlargement reactions can provide a fast route towards the formation of six-membered single-ring or polycyclic aromatic hydrocarbons (PAHs).

Theoretical study of sensitive reactions in phenol decomposition

The reactivity of phenol is of utmost importance in combustion systems.

New theoretical insights into the reaction kinetics of toluene and hydroxyl radicals

This work provides theoretical insights into the kinetics of toluene + OH, focusing on the anharmonic effect and the accuracy of barrier heights.

Radical-Pair Formation in Hydrocarbon (Aut)Oxidation

The reaction profiles for the uni- and bimolecular decomposition of benzyl hydroperoxide have been studied in the context of initiation reactions for the (aut)oxidation of hydrocarbons. The unimolecular dissociation of benzyl hydroperoxide was found to proceed through the formation of a hydrogen-bonded radical-pair minimum located +181 kJ mol^-1 above the hydroperoxide substrate and around 15 kJ mol^-1 below the separated radical products. The reaction of toluene with benzyl hydroperoxide proceeds such that O-O bond homolysis is coupled with a C-H bond abstraction event in a single kinetic step. The enthalpic barrier of this molecule-induced radical formation (MIRF) process is significantly lower than that of the unimolecular O-O bond cleavage. The same type of reaction is also possible in the self-reaction between two benzyl hydroperoxide molecules forming benzyloxyl and hydroxyl radical pairs along with benzaldehyde and water as co-products. In the product complexes formed in these MIRF reactions, both radicals connect to a centrally placed water molecule through hydrogen-bonding interactions.

Combined Experimental and Computational Study on the Reaction Dynamics of the 1-Propynyl (CH3CC)-1,3-Butadiene (CH2CHCHCH2) System and the Formation of Toluene under Single Collision Conditions

The crossed beams reactions of the 1-propynyl radical (CH₃CC; X²A₁) with 1,3-butadiene (CH₂CHCHCH₂; X¹A_g), 1,3-butadiene- d₆ (CD₂CDCDCD₂; X¹A_g), 1,3-butadiene- d₄ (CD₂CHCHCD₂; X¹A_g), and 1,3-butadiene- d₂ (CH₂CDCDCH₂; X¹A_g) were performed under single collision conditions at collision energies of about 40 kJ mol^-1. The underlying reaction mechanisms were unraveled through the combination of the experimental data with electronic structure calculations at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) + ZPE(B3LYP/6-311G(d,p) level of theory along with statistical Rice-Ramsperger-Kassel-Marcus (RRKM) calculations. Together, these data suggest the formation of the thermodynamically most stable C₇H₈ isomer-toluene (C₆H₅CH₃)-via the barrierless addition of 1-propynyl to the 1,3-butadiene terminal carbon atom, forming a low-lying C₇H₉ intermediate that undergoes multiple isomerization steps resulting in cyclization and ultimately aromatization following hydrogen atom elimination. RRKM calculations predict that the thermodynamically less stable isomers 1,3-heptadien-5-yne, 5-methylene-1,3-cyclohexadiene, and 3-methylene-1-hexen-4-yne are also synthesized. Since the 1-propynyl radical may be present in cold molecular clouds such as TMC-1, this pathway could potentially serve as a carrier of the methyl group incorporating itself into methyl-substituted (poly)acetylenes or aromatic systems such as toluene via overall exoergic reaction mechanisms that are uninhibited by an entrance barrier. Such pathways are a necessary alternative to existing high energy reactions leading to toluene that are formally closed in the cold regions of space and are an important step toward understanding the synthesis of polycyclic aromatic hydrocarbons (PAHs) in space's harsh extremes.

Kinetics of the Hydrogen Abstraction PAH + •OH → PAH Radical + H2O Reaction Class: An Application of the Reaction Class Transition State Theory (RC-TST) and Structure-Activity Relationship (SAR)

A reaction class transition state theory (RC-TST) augmented with structure-activity relationship (SAR) methodology is applied to predict high-pressure limit thermal rate constants for hydrogen abstraction by ^•OH radical from polycyclic aromatic hydrocarbons (PAHs) reaction class in the temperature range of 300-3000 K. The rate constants for the reference reaction of C₆H₆ + ^•OH → C₆H₅ + H₂O is calculated by the canonical variational transition state theory (CVT) with small curvature tunneling (SCT). Only the reaction energy is needed to predict RC-TST rates for other processes within the family, the parameters needed were obtained from M06-2X/cc-pVTZ data for a training set of 34 reactions. The systematic error of the resulting RC-TST rates is smaller than 50% in comparison with explicit rate calculations, which facilitates application of the proposed methodology to the automated reaction mechanism generators (ARMGs) schemes.

A new era for combustion research

Current topics in combustion chemistry include aspects of a changing fuel spectrum with a focus on reducing emissions and increasing efficiency. This article is intended to provide an overview of selected recent work in combustion chemistry, especially addressing reaction pathways from fuel decomposition to emissions. The role of the molecular fuel structure will be emphasized for the formation of certain regulated and unregulated species from individual fuels and their mixtures, exemplarily including fuel compounds such as alkanes, alkenes, ethers, alcohols, ketones, esters, and furan derivatives. Depending on the combustion conditions, different temperature regimes are important and can lead to different reaction classes. Laboratory reactors and flames are prime sources and targets from which such detailed chemical information can be obtained and verified with a number of advanced diagnostic techniques, often supported by theoretical work and simulation with combustion models developed to transfer relevant details of chemical mechanisms into practical applications. Regarding the need for cleaner combustion processes, some related background and perspectives will be provided regarding the context for future chemistry research in combustion energy science.

Accurate entropy calculation for large flexible hydrocarbons using a multi-structural 2-dimensional torsion method

Entropy is one of the key thermodynamic parameters in combustion kinetic modeling. Accurate entropy prediction needs to account for the conformational torsional anharmonicity, which could be solved by the state-of-the-art multi-structural torsion (MS-T) method. However, this method is computationally expensive or even not feasible for large flexible molecules. To address this issue, we proposed a multi-structural 2-dimensional torsion (MS-2DT) method that adopts minimally coupled torsions to reduce the computational cost. In this method, a series of 2-dimensional coupled torsion combinations were used to generate an initial conformer space with a size of CN2·9 (N is the number of torsions). The standard entropy (and the heat capacity) values of 18 C6-C8 alkanes with 5-7 torsions were computed at 200-2000 K. The MS-2DT calculation is in good agreement with the benchmark MS-T method: only a small deviation of -0.19 ± 0.15 cal mol-1 K-1 in standard entropy and -0.10 ± 0.21 cal mol-1 K-1 in heat capacity. Additionally, a further application of MS-2DT to n-decane with 9 torsions implies an improved accuracy in entropy (and heat capacity) prediction compared to other conventional simplified treatments. This method provides an affordable and accurate solution to treat the conformational torsional anharmonicity of large flexible alkanes.

Kinetics of the Toluene Reaction with OH Radical

We calculated the kinetics of chemical activation reactions of toluene with hydroxyl radical in the temperature range from 213 K to 2500 K and the pressure range from 10 Torr to the high-pressure limit by using multistructural variational transition state theory with the small-curvature tunneling approximation (MS-CVT/SCT) and using the system-specific quantum Rice-Ramsperger-Kassel method. The reactions of OH with toluene are important elementary steps in both combustion and atmospheric chemistry, and thus it is valuable to understand the rate constants both in the high-pressure, high-temperature regime and in the low-pressure, low-temperature regime. Under the experimental pressure conditions, the theoretically calculated total reaction rate constants agree well with the limited experimental data, including the negative temperature dependence at low temperature. We find that the effect of multistructural anharmonicity on the partition functions usually increases with temperature, and it can change the calculated reaction rates by factors as small as 0.2 and as large as 4.2. We also find a large effect of anharmonicity on the zero-point energies of the transition states for the abstraction reactions. We report that abstraction of H from methyl should not be neglected in atmospheric chemistry, even though the low-temperature results are dominated by addition. We calculated the product distribution, which is usually not accessible to experiments, as a function of temperature and pressure.

Detailed kinetics of substituted phenolic species in pyrolysis bio-oils

A comprehensive kinetic model for the pyrolysis and combustion of substituted phenolic species, key components of fast pyrolysis bio-oils.