Theoretical Kinetic Study of Thermal Unimolecular Decomposition of Cyclic Alkyl Radicals

High-resolution infrared spectroscopy of jet cooled cyclobutyl in the α-CH stretch region: large-amplitude puckering dynamics in a 4-membered ring radical

Gas-phase cyclobutyl radical (c-C4H7) is generated at a rotational temperature of Trot = 26(1) K in a slit-jet discharge mixture of 70% Ne/30% He and 0.5-0.6% cyclobromobutane (c-C4H7Br). A fully rovibrationally resolved absorption spectrum of the α-CH stretch fundamental band between 3062.9 cm-1 to 3075.7 cm-1 is obtained and analyzed, yielding first precision structural and dynamical information for this novel radical species. The α-CH stretch band origin is determined to be 3068.7887(4) cm-1, which implies only a modest (≈0.8 cm-1) blue shift from rotationally unresolved infrared spectroscopic studies of cyclobutyl radicals in liquid He droplets [A. R. Brown, P. R. Franke and G. E. Douberly, J. Phys. Chem. A, 2017, 121, 7576-7587]. Of particular dynamical interest, a one-dimensional potential energy surface with respect to the ring puckering coordinate is computed at CCSD(T)/ANO2 level of theory and reveals a double minimum Cs puckered geometry, separated by an exceedingly shallow planar C2v transition state barrier (Ebarr ≈ 1 cm-1). Numerical solutions on this double minimum potential yield a zero-point energy for the ground state (Ezero-point ≈ 27 cm-1) greatly in excess of the interconversion barrier. This is indicative of highly delocalized, large amplitude motion of the four-membered ring structure, for which proper vibrationally averaging of the moment of inertia tensor reproduces the experimentally determined inertial defect remarkably well. Finally, intensity alternation in the experimental spectrum due to nuclear spin statistics upon exchange of three indistinguishable H atom pairs (IH = ½) matches Ka + Kc = even : odd = 36 : 28 predictions, implying that the unpaired electron in the radical center lies in an out-of-plane pπ orbital. Thus, high-resolution infrared spectroscopy provides first experimental confirmation of a shallow double minimum ring puckering potential with a highly delocalized ground state wave function peaked at a planar C2v transition state geometry consistent with a cyclobutyl π radical.

Photodissociation Dynamics of the Cyclohexyl Radical from the 3p Rydberg State at 248 nm

The photodissociation of jet-cooled cyclohexyl was studied by exciting the radicals to their 3p Rydberg state by using 248 nm laser light and detecting photoproducts by photofragment translational spectroscopy. Both H atom loss and dissociation to heavy fragment pairs are observed. The H atom loss channel exhibits a two-component translational energy distribution. The fast photoproduct component is attributed to impulsive cleavage directly from an excited state, likely the Rydberg 3s state, forming cyclohexene. The slow component is due to statistical decomposition of hot cyclohexyl radicals that internally convert to the ground electronic state prior to H atom loss. The fast and slow components are present in an ∼0.7:1 ratio, similar to findings in other alkyl radicals. Internal conversion to the ground state also leads to ring-opening followed by dissociation to 1-buten-4-yl + ethene in comparable yield to H-loss, with the C₄H₇ fragment containing enough internal energy to dissociate further to butadiene via H atom loss. A very minor ground-state C₅H₈ + CH₃ channel is observed, attributed predominantly to 1,3-pentadiene formation. The ground-state branching ratios agree well with RRKM calculations, which also predict C₄H₆ + C₂H₅ and C₃H₆ + C₃H₅ channels with similar yield to C₅H₈ + CH₃. If these channels were active, it was at levels too low to be observed.

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.

Ab Initio/Transition-State Theory Study of the Reactions of Ċ5H9 Species of Relevance to 1,3-Pentadiene, Part I: Potential Energy Surfaces, Thermochemistry, and High-Pressure Limiting Rate Constants

In this study, the reactions of Ċ₅H₉ radicals are theoretically investigated, with a particular emphasis on hydrogen atom addition reactions to 1,3-pentadiene (C₅H₈) to form Ċ₅H₉ radicals, although the subsequent isomerization and decomposition reactions of the Ċ₅H₉ radicals are also of direct relevance to the radicals formed from the pyrolysis and oxidation of species including pentene and cyclopentane. Moreover, H-atom abstraction reactions by hydrogen atoms from 1,3-pentadiene are also investigated. The geometries and frequencies of 63 potential energy surface (PES) minima and 88 transition states are optimized at the ωB97XD/aug-cc-pVTZ level of theory. Spin-unrestricted open-shell single-point energies for all the species are calculated at the CCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). A one-dimensional hindered rotor treatment is employed for torsional modes, with the M06-2X/6-311++G(D,P) method used to compute the potential energy as a function of the dihedral angle. The high-pressure limiting rate constants and the thermochemical properties for C₅ species are calculated using the Master Equation System Solver (MESS) with conventional transition-state theory and comparisons made with existing available literature data. A hydrogen atom can add to the terminal carbon atom of 1,3-pentadiene to form the 2,4-Ċ₅H₉ radical and/or the internal carbon atoms to form 2,5-Ċ₅H₉, 1,4-Ċ₅H₉, and 1,3-Ċ₅H₉ radicals. Among the four entrance channels for Ḣ atom addition reactions, the formation of 2,4-Ċ₅H₉ and 1,3-Ċ₅H₉ radicals is more exothermic in comparison to the other Ċ₅H₉ isomers (2,5-Ċ₅H₉, 1,4-Ċ₅H₉) because of the resonantly stabilized allylic structure. Consequently, the formation of the former is generally dominant in terms of barrier heights. Ḣ atom addition reactions to 1,3-pentadiene are compared to available C₃-C₅ alkenes and dienes, with external addition calculated to be kinetically favored over internal addition. However, the correlation between heats of formation and energy barriers for Ḣ atom addition to 1,2-dienes is different from that for 1,3- and 1,4-dienes. Hydrogen atom addition and abstraction rate constants are also compared for 1,3-pentadiene, with addition found to be dominant. The subsequent unimolecular reactions on the Ċ₅H₉ PES are found to be highly complex with reactions taking place on a multiple-well multiple-channel PES. For clarity, the reaction mechanism and kinetics of each Ċ₅H₉ radical are discussed individually in terms of the computed enthalpy of the reaction and activation, the transition-state structure/reaction class, and also in terms of the combustion species for which the reactions are of potential importance. The reactions on the Ċ₅H₉ PES are divided into three reaction classes (H-shift isomerization, cycloaddition, and β-scission reactions), and the reactivity-structure-based estimation rules for energy barriers are derived for these three reaction classes and compared to literature results for alkyl radicals.

Low Temperature Oxidation of Cyclohexane: Uncertainty of Important Thermo-Chemical Properties

The study of the standard formation enthalpy, entropy, and heat capacity for key species relevant to the low-temperature combustion of cyclohexane has been performed by applying the group additivity method of Benson. The properties of 18 Benson groups (8 of them for the first time), and 10 ring correction factors for cyclic species were estimated through different empirical and semi-empirical methods. The method validation proceeded through comparison of predicted values for certain number of newly estimated groups and available literature data derived from quantum chemistry estimations. Further validations of the estimated properties of groups have been provided by comparing estimated properties of test species with data in literature and kinetic databases. Also the standard deviation between prediction and reported values has been evaluated for each validation case. A similar approach has been applied for validation of the estimated ring correction groups. For selected well-studied cyclic molecules the predicted values and the literature data have been compared with each other, and the standard deviations have been also reported. The evaluated properties of the cyclohexane relevant species were also compared with similar ones available in other kinetic models and in databases. At the end the estimated properties have been presented in a tabulated form of NASA polynomial coefficients with extrapolation up to 3500 K.

Synthesis, Structure, and Cytotoxicity of a New Sulphanyl-Bridged Thiadiazolyl-Saccharinate Conjugate: The Relevance of S⋅⋅⋅N Interaction

Reports showing that the copper concentration is considerably higher in neoplasms than in normal tissues prompted the need to develop selective copper chelators. We disclosed recently that some N-linked tetrazole-saccharinates bind selectively to copper, forming complexes that are highly cytotoxic towards cancer cells. Because tetrazole-saccharinates are photolabile, due to the photoreactivity of tetrazoles, we proposed thiadiazolyl-saccharinates as an alternative. Herein we describe the synthesis, structure, and monomeric photochemistry of a sulphanyl-bridged thiadiazolyl-saccharinate, 3-[(5-methyl-1,3,4-thiadiazol-2-yl)sulphanyl]-1,2-benzothiazole 1,1-dioxide (MTSB). The monomeric structure, charge density analysis, and characteristic infrared spectrum of MTSB were investigated theoretically, using quantum chemical calculations, and also experimentally, using matrix-isolation infrared spectroscopy. The crystal structure was investigated by combining X-ray crystallography with infrared and Raman spectroscopies. Results show that the structure of isolated MTSB is similar to that found in the crystal, with an S⋅⋅⋅N interaction clearly contributing to the structure of the molecule and of the crystal. Matrix irradiation revealed a high photostability of MTSB, compared to parent tetrazole-saccharinates and to the 5-methyl-1,3,4-thiadiazole building block, emphasizing the photostabilizing effect of the saccharyl system. Finally, in vitro toxicity assays of MTSB showed a copper concentration-dependent toxicity against cancer cells, without affecting normal cells. In particular, MTSB was most effective towards the hepatic (HepG2), neuroblastoma (SH-SY5), and lymphoma cell lines (U937). Thus, MTSB represents a promising lead for cancer chemotherapy based on chelating agents.

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.

Kinetic and mechanistic investigations of the thermal decomposition of methyl-substituted cycloalkyl radicals

We perform systemic theoretical investigations on the thermal decomposition of 2-Me-cyclobutyl, 2-Me-cyclopentyl and 2-Me-cyclohexyl radicals at CBS-QB3 and CCSD(T) levels.

An experimental and modeling study of the low- and high-temperature oxidation of cyclohexane

The experimental study of the oxidation of cyclohexane has been performed in a jet-stirred reactor at temperatures ranging from 500 to 1100 K (low- and intermediate temperature zones including the negative temperature-coefficient area), at a residence time of 2 s and for dilute mixtures with equivalence ratios of 0.5, 1, and 2. Experiments were carried out at quasi-atmospheric pressure (1.07 bar). The fuel and reaction product mole fractions were measured using online gas chromatography. A total of 34 reaction products have been detected and quantified in this study. Typical reaction products formed in the low-temperature oxidation of cyclohexane include cyclic ethers (1,2-epoxycyclohexane and 1,4-epoxycyclohexane), 5-hexenal (formed from the rapid decomposition of 1,3-epoxycyclohexane), cyclohexanone, and cyclohexene, as well as benzene and phenol. Cyclohexane displays high low-temperature reactivity with well-marked negative temperature-coefficient (NTC) behavior at equivalence ratios 0.5 and 1. The fuel-rich system (ϕ = 2) is much less reactive in the same region and exhibits no NTC. To the best of our knowledge, this is the first jet-stirred reactor study to report NTC in cyclohexane oxidation. Laminar burning velocities were also measured by the heated burner method at initial gas temperatures of 298, 358, and 398 K and at 1 atm. The laminar burning velocity values peak at ϕ = 1.1 and are measured as 40 and 63.1 cm/s for Ti = 298 and 398 K, respectively. An updated detailed chemical kinetic model including low-temperature pathways was used to simulate the present (jet-stirred reactor and laminar burning velocity) and literature experimental (laminar burning velocity, rapid compression machine, and shock tube ignition delay times) data. Reasonable agreement is observed with most of the products observed in our reactor, as well as the literature experimental data considered in this paper.

Shock tube and chemical kinetic modeling study of the oxidation of 2,5-dimethylfuran

A detailed kinetic model describing the oxidation of 2,5-dimethylfuran (DMF), a potential second-generation biofuel, is proposed. The kinetic model is based upon quantum chemical calculations for the initial DMF consumption reactions and important reactions of intermediates. The model is validated by comparison to new DMF shock tube ignition delay time measurements (over the temperature range 1300-1831 K and at nominal pressures of 1 and 4 bar) and the DMF pyrolysis speciation measurements of Lifshitz et al. [ J. Phys. Chem. A 1998 , 102 ( 52 ), 10655 - 10670 ]. Globally, modeling predictions are in good agreement with the considered experimental targets. In particular, ignition delay times are predicted well by the new model, with model-experiment deviations of at most a factor of 2, and DMF pyrolysis conversion is predicted well, to within experimental scatter of the Lifshitz et al. data. Additionally, comparisons of measured and model predicted pyrolysis speciation provides validation of theoretically calculated channels for the oxidation of DMF. Sensitivity and reaction flux analyses highlight important reactions as well as the primary reaction pathways responsible for the decomposition of DMF and formation and destruction of key intermediate and product species.

PROGRESS IN DETAILED KINETIC MODELING OF THE COMBUSTION OF OXYGENATED COMPONENTS OF BIOFUELS

Due to growing environmental concerns and diminishing petroleum reserves, a wide range of oxygenated species has been proposed as possible substitutes to fossil fuels: alcohols, methyl esters, acyclic and cyclic ethers. After a short review the major detailed kinetic models already proposed in the literature for the combustion of these molecules, the specific classes of reactions considered for modeling the oxidation of acyclic and cyclic oxygenated molecules respectively, are detailed.

Towards cleaner combustion engines through groundbreaking detailed chemical kinetic models

In the context of limiting the environmental impact of transportation, this critical review discusses new directions which are being followed in the development of more predictive and more accurate detailed chemical kinetic models for the combustion of fuels. In the first part, the performance of current models, especially in terms of the prediction of pollutant formation, is evaluated. In the next parts, recent methods and ways to improve these models are described. An emphasis is given on the development of detailed models based on elementary reactions, on the production of the related thermochemical and kinetic parameters, and on the experimental techniques available to produce the data necessary to evaluate model predictions under well defined conditions (212 references).

Mass spectra of cyclic ethers formed in the low-temperature oxidation of a series of n-alkanes

Cyclic ethers are important intermediate species formed during the low-temperature oxidation of hydrocarbons. Along with ketones and aldehydes, they could consequently represent a significant part of the heavy oxygenated pollutants observed in the exhaust gas of engines. Apart a few of them such as ethylene oxide and tetrahydrofuran, cyclic ethers have not been much studied and very few of them are available for calibration and identification. Electron impact mass spectra are available for very few of them, making their detection in the exhaust emissions of combustion processes very difficult. The main goal of this study was to complete the existing set of mass spectra for this class of molecules. Thus cyclic ethers have been analyzed in the exhaust gases of a jet-stirred reactor in which the low-temperature oxidation of a series of n-alkanes was taking place. Analyzes were performed by gas chromatography coupled to mass spectrometry and to MS/MS. The second goal of this study was to derive some rules for the fragmentation of cyclic ethers in electron impact mass spectrometry and allow the identification of these species when no mass spectrum is available.