A quantitative kinetic analysis of CO elimination from phenoxy radicals

Investigation of Methylcyclopentadiene Reactivity: Abstraction Reactions and Methylcyclopentadienyl Radical Unimolecular Decomposition

Understanding the reactivities of methylcyclopentadiene and the methylcyclopentadienyl radical is important in order to improve our comprehension of the chemical kinetics leading to the formation, decomposition, and growth of the first aromatic ring, as it has been shown that five-membered-ring species are important intermediates in the reaction kinetics of aromatic species. In this work, the rate constants of some key H-abstraction reactions from methylcyclopentadiene to produce the methylcyclopentadienyl radical and the formation of fulvene and benzene from the latter are theoretically determined. Rate constants are evaluated using the ab initio transition state theory-based master equation approach, determining structures and Hessians of all stationary points at the ωB97X-D/aug-cc-pVTZ level, energies at the CCSD(T) level extrapolated to the complete basis set limit, RRKM rate constants using conventional and variational transition state theory, and phenomenological rate constants through the solution of the one-dimensional master equation. Variational corrections are determined in both internal and Cartesian coordinates, and it is found that the choice of the coordinate system can impact the accuracy of the calculated rate constants by up to a factor of 4 for H-abstraction reactions and 2 for the unimolecular decomposition of the methylcyclopentadienyl radical. The calculated rate constants are in good agreement with the available literature data. Prompt dissociation of methylcyclopentadienyl radicals accessed following H-abstraction from methylcyclopentadiene was also investigated, and the corresponding rate constants were determined; the results show that prompt dissociation plays a key role under combustion conditions. Finally, lumping of theoretically derived rate constants to account for methylcyclopentadiene ⇄ methylcyclopentadienyl tautomerism allowed the derivation of a simplified set of rate constants suitable to be inserted into kinetic mechanisms.

A Reaction Kinetics Study on Benzene Oxidation in the Claus Process by Sulfur Monoxide

The Claus process is used in natural gas processing plants to treat H₂S-rich acid gas to recover sulfur, but the process suffers from catalytic deactivation when aromatic contaminants such as benzene, toluene, ethylbenzene, and xylene isomers (collectively called as BTEX) are present in the acid gas feed. To safeguard the catalytic reactors, it is desired to oxidize aromatic contaminants in the furnace that are present upstream of the catalytic reactors in the process by oxidants present in it. This work develops a reaction mechanism and evaluates the reaction kinetics for the oxidation of phenyl radical by SO using CBS-QBS for reaction energetics and RRKM and transition state theory for reaction kinetics. The mechanism explores the possible products that are formed from the barrierless addition of SO on phenyl through the O atom as well as through the S atom. The exothermicity of the addition reaction is higher when the addition of SO on the aromatic structure takes place through the S atom. The major products formed from phenyl oxidation by SO are cyclopentadienyl, cyclopentadienethiol and thiopyran radicals. A remarkable similarity between the pathways for phenyl radical oxidation by O₂ and its oxidation by SO at high temperatures is observed. The proposed reactions and their rate constants are used to conduct reactor simulations to determine the important reactions that contribute to the formation of major products during phenyl-SO reactions and the temperatures suitable for benzene oxidation by SO under process conditions similar to the Claus furnace.

An automatized workflow from molecular dynamic simulation to quantum chemical methods to identify elementary reactions and compute reaction constants

We present an automatized workflow which, starting from molecular dynamics simulations, identifies reaction events, filters them, and prepares them for accurate quantum chemical calculations using, for example, Density Functional Theory (DFT) or Coupled Cluster methods. The capabilities of the automatized workflow are demonstrated by the example of simulations for the combustion of some polycyclic aromatic hydrocarbons (PAHs). It is shown how key elementary reaction candidates are filtered out of a much larger set of redundant reactions and refined further. The molecular species in question are optimized using DFT and reaction energies, barrier heights, and reaction rates are calculated. The setup is general enough to include at this stage configurational sampling, which can be exploited in the future. Using the introduced machinery, we investigate how the observed reaction types depend on the gas atmosphere used in the molecular dynamics simulation. For the re-optimization on the DFT level, we show how the additional information needed to switch from reactive force-field to electronic structure calculations can be filled in and study how well ReaxFF and DFT agree with each other and shine light on the perspective of using more accurate semi-empirical methods in the MD simulation.

Conformational Space, IR-Induced, and UV-Induced Chemistry of Carvacrol Isolated in a Low-Temperature Argon Matrix

In this work, monomers of carvacrol (5-isopropyl-2-methylphenol), a natural monoterpene exhibiting wide range bioactivity, were trapped in a cryogenic argon matrix and characterized by infrared spectroscopy, while quantum chemical calculations at the B3LYP and MP2 levels were employed to characterize the conformational landscape of the isolated molecule. Four conformers have been localized on the potential energy surface, and the factors accounting for their relative stability were analyzed. The two most stable conformers of carvacrol, differing in the relative orientation of the isopropyl group and both having the OH group pointing away from the vicinal methyl fragment, were identified in the cryomatrix for the first time. The individual spectral signatures of the two conformers were distinguished based on the change in their relative abundance induced by exposing the matrix to broadband infrared light. Matrix-isolated carvacrol was also irradiated with broadband UV light (λ > 200 nm), which resulted in the cleavage of the OH group. Recombination of the released H atom at the ortho- or para-position of the ring resulted in the formation of alkyl-substituted cyclohexadienones. These were found to undergo subsequent valence and open-ring isomerizations, leading, respectively, to the formation of a Dewar isomer and open-chain conjugated ketenes. Decarbonylation of the photoproducts was also observed for longer irradiation times. A mechanistic analysis of the observed photochemical transformations is presented.

Theoretical study of sensitive reactions in phenol decomposition

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

Energy transfer in intermolecular collisions of polycyclic aromatic hydrocarbons with bath gases He and Ar

Classical trajectory simulations of intermolecular collisions were performed for a series of polycyclic aromatic hydrocarbons interacting with the bath gases helium and argon for bath gas temperature from 300 to 2500 K. The phase-space average energy transferred per deactivating collision, ⟨∆E_down⟩, was obtained. The Buckingham pairwise intermolecular potentials were validated against high-level quantum chemistry calculations and used in the simulations. The reactive force-field was used to describe intramolecular potentials. The dependence of ⟨∆E_down⟩ on initial vibrational energy is discussed. A canonical sampling method was compared with a microcanonical sampling method for selecting initial vibrational energy at high bath gas temperatures. Uncertainties introduced by the initial angular momentum distribution were identified. The dependence of the collisional energy transfer parameters on the type of bath gas and the molecular structure of polycyclic aromatic hydrocarbons was examined.

Theoretical Study on Elementary Reaction Steps in Thermal Decomposition Processes of Syringol-Type Monolignol Compounds

This paper theoretically investigated a large number of reaction pathways and kinetics to describe the vapor-phase pyrolytic behavior of several syringol-type monolignol compounds that are derived from the primary pyrolysis of lignin: 1-(4-hydroxy-3,5-dimethoxyphenyl)prop-2-en-1-one (HDPP), sinapyl alcohol, 3-hydroxy-1-(4-hydroxy-3,5-dimethoxyphenyl)propan-1-one (HHDPP), 1-(4-hydroxy-3,5-dimethoxyphenyl)propane-1,3-diol (HDPPD), and syringol. The possible pyrolytic pathways involving unimolecular decomposition, addition, and abstraction reactions were investigated by comparing the energy barriers calculated at the B3LYP/6-311++G(d,p) level. In the proposed pathways, all syringol-type monolignols containing a side chain undergo its cleavage to form syringol through the formation of syringaldehyde or 4-vinylsyringol. Syringol is then converted into two products: (a) pyrogallol via the homolysis of the O-CH₃ bond and hydrogenation or (b) guaiacol via addition of an H atom with a carbon bearing methoxyl group in syrignol and the subsequent demethoxylation. The pyrolytic pathways of pyrogallol are classified into two processes: (a) the concerted dehydrogenation of the two hydroxyl H atoms and the unimolecular decomposition to produce acetylene (C₂H₂), ethynol (C₂HOH), and CO or (b) the displacement of an OH with H to produce catechol and resorcinol. Additionally, HDPP undergoes O-CH₃ bond cleavage to form but-1-en-3-yne. The high-pressure limit rate constants for all the proposed elementary reaction steps were evaluated on the basis of transition state theory.

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.

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.

Photodissociation dynamics of the thiophenoxy radical at 248, 193, and 157 nm

The photodissociation dynamics of the thiophenoxy radical (C6H5S) have been investigated using fast beam coincidence translational spectroscopy. Thiophenoxy radicals were produced by photodetachment of the thiophenoxide anion followed by photodissociation at 248 nm (5.0 eV), 193 nm (6.4 eV), and 157 nm (7.9 eV). Experimental results indicate two major competing dissociation channels leading to SH + C6H4 (o-benzyne) and CS + C5H5 (cyclopentadienyl) with a minor contribution of S + C6H5 (phenyl). Photofragment mass distributions and translational energy distributions were measured at each dissociation wavelength. Transition states and minima for each reaction pathway were calculated using density functional theory to facilitate experimental interpretation. The proposed dissociation mechanism involves internal conversion from the initially prepared electronic excited state to the ground electronic state followed by statistical dissociation. Calculations show that SH loss involves a single isomerization step followed by simple bond fission. For both SH and S loss, C-S bond cleavage proceeds without an exit barrier. By contrast, the CS loss pathway entails multiple transition states and minima as it undergoes five membered ring formation and presents a small barrier with respect to products. The calculated reaction pathway is consistent with the experimental translational energy distributions in which the CS loss channel has a broader distribution peaking farther away from zero than the corresponding distributions for SH loss.