Detailed kinetic study of the ring opening of cycloalkanes by CBS-QB3 calculations

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.

Flash pyrolysis vacuum ultraviolet photoionization mass spectrometry of cycloheptane: A study of the initial decomposition mechanism

Thermal decomposition of cycloheptane was studied using flash pyrolysis coupled with vacuum ultraviolet (118.2 nm) single photon ionization time-of-flight mass spectrometry at temperatures ranging from 295 K to 1380 K. C-C bond breaking of cycloheptane leading to the 1,7-heptyl diradical was considered as the initiation step. The 1,7-heptyl diradical could readily isomerize to 1-heptene and decompose into several fragments, with dissociation to •C₄H₉ and •C₃H₅ as the predominant product channel. The 1,7-heptyl diradical could undergo direct dissociation, as evidenced by the production of the C₅H₁₀ species. Quantum chemistry calculations at UCCSD(T)/cc-pVDZ//UB3LYP/cc-pVDZ level of theory on the initial reaction pathways of cycloheptane were also carried out to support the experimental observations. Other possible initiation channels, as well as some secondary reaction products, were also identified.

Genotoxic Effect of Dicyclopropanated 5-Vinyl-2-Norbornene

Dicyclopropanated 5-vinyl-2-norbornene (dcpVNB) is a strained polycyclic hydrocarbon compound with a high energy content, which makes it promising for the development of propellant components based on it. In this work, the genotoxic properties of dcpVNB were studied using whole-cell lux-biosensors based on Escherichia coli and Bacillus subtilis. It was shown that the addition of dcpVNB to bacterial cells leads to the appearance of DNA damage inducing the SOS response and Dps expression with slight activation of the OxyR-mediated response to oxidative stress. The highest toxic effect of dcpVNB is detected by the following lux-biosensors: E. coli pColD-lux, E. coli pDps, B. subtilis pNK-DinC, and B. subtilis pNK-MrgA, in which the genes of bacterial luciferases are transcriptionally fused to the corresponding promoters: P_cda, P_dps, P_dinC, and P_mrgA. It was shown that lux-biosensors based on B. subtilis, and E. coli are almost equally sensitive to dcpVNB, which indicates the same permeability to this compound of cell wall of Gram-positive and Gram-negative bacteria. The activation of P_dps after dcpVNB addition maintains even in oxyR mutant E. coli strains, which means that the P_dps induction is only partially determined by the OxyR/S regulon. Comparison of specific stress effects caused by dcpVNB and 2-ethyl(bicyclo[2.2.1]heptane) (EBH), characterized by the absence of cyclopropanated groups, shows that structural changes in hydrocarbons could significantly change the mode of toxicity.

Minimal Active Space for Diradicals Using Multistate Density Functional Theory

This work explores the electronic structure as well as the reactivity of singlet diradicals, making use of multistate density functional theory (MSDFT). In particular, we show that a minimal active space of two electrons in two orbitals is adequate to treat the relative energies of the singlet and triplet adiabatic ground state as well as the first singlet excited state in many cases. This is plausible because dynamic correlation is included in the first place in the optimization of orbitals in each determinant state via block-localized Kohn–Sham density functional theory. In addition, molecular fragment, i.e., block-localized Kohn–Sham orbitals, are optimized separately for each determinant, providing a variational diabatic representation of valence bond-like states, which are subsequently used in nonorthogonal state interactions (NOSIs). The computational procedure and its performance are illustrated on some prototypical diradical species. It is shown that NOSI calculations in MSDFT can be used to model bond dissociation and hydrogen-atom transfer reactions, employing a minimal number of configuration state functions as the basis states. For p- and s-types of diradicals, the closed-shell diradicals are found to be more reactive than the open-shell ones due to a larger diabatic coupling with the final product state. Such a diabatic representation may be useful to define reaction coordinates for electron transfer, proton transfer and coupled electron and proton transfer reactions in condensed-phase simulations.

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.

Thermal decomposition of cyclohexane by flash pyrolysis vacuum ultraviolet photoionization time-of-flight mass spectrometry: a study on the initial unimolecular decomposition mechanism

Thermal decomposition of cyclohexane at temperatures up to 1310 K was performed using flash pyrolysis coupled with vacuum ultraviolet (118.2 nm) photoionization time-of-flight mass spectrometry. The experimental results revealed that the major initiation reaction of cyclohexane decomposition was C-C bond fission leading to the formation of 1,6-hexyl diradical. The 1,6-hexyl diradical could isomerize to 1-hexene and decompose into ˙C3H7 + ˙C3H5 and ˙C4H7 + ˙C2H5. The 1,6-hexyl diradical could also undergo direct dissociation; the C4H8 fragment via the 1,4-butyl diradical intermediate was observed, serving as evidence of the 1,6-hexyl diradical mechanism. Quantum chemistry calculations at UCCSD(T)/cc-pVDZ level of theory on the initial reaction pathways of cyclohexane were performed and found to be consistent with the experimental conclusions. Cyclohexyl radical was not observed as an initial intermediate in the pyrolysis. Benzene was produced from sequential H2 eliminations of cyclohexane at high temperatures.

Elucidating the differences in oxidation of high-performance α- and β- diisobutylene biofuels via Synchrotron photoionization mass spectrometry

Biofuels are a promising ecologically viable and renewable alternative to petroleum fuels, with the potential to reduce net greenhouse gas emissions. However, biomass sourced fuels are often produced as blends of hydrocarbons and their oxygenates. Such blending complicates the implementation of these fuels in combustion applications. Variations in a biofuel's composition will dictate combustion properties such as auto ignition temperature, reaction delay time, and reaction pathways. A handful of novel drop-in replacement biofuels for conventional transportation fuels have recently been down selected from a list of over 10,000 potential candidates as part of the U.S. Department of Energy's (DOE) Co-Optimization of Fuels and Engines (Co-Optima) initiative. Diisobutylene (DIB) is one such high-performing hydrocarbon which can readily be produced from the dehydration and dimerization of isobutanol, produced from the fermentation of biomass-derived sugars. The two most common isomers realized, from this process, are 2,4,4-trimethyl-1-pentene (α-DIB) and 2,4,4-trimethyl-2-pentene (β-DIB). Due to a difference in olefinic bond location, the α- and β- isomer exhibit dramatically different ignition temperatures at constant pressure and equivalence ratio. This may be attributed to different fragmentation pathways enabled by allylic versus vinylic carbons. For optimal implementation of these biofuel candidates, explicit identification of the intermediates formed during the combustion of each of the isomers is needed. To investigate the combustion pathways of these molecules, tunable vacuum ultraviolet (VUV) light (in the range 8.1-11.0 eV) available at the Lawrence Berkeley National Laboratory's Advanced Light Source (ALS) has been used in conjunction with a jet stirred reactor (JSR) and time-of-flight mass spectrometry to probe intermediates formed. Relative intensity curves for intermediate mass fragments produced during this process were obtained. Several important unique intermediates were identified at the lowest observable combustion temperature with static pressure of 93,325 Pa and for 1.5 s residence time. As this relatively short residence time is just after ignition, this study is targeted at the fuels' ignition events. Ignition characteristics for both isomers were found to be strongly dependent on the kinetics of C₄ and C₇ fragment production and decomposition, with the tert-butyl radical as a key intermediate species. However, the ignition of α-DIB exhibited larger concentrations of C₄ compounds over C₇, while the reverse was true for β-DIB. These identified species will allow for enhanced engineering modeling of fuel blending and engine design.

Genotoxic effect of 2,2'-bis(bicyclo[2.2.1] heptane) on bacterial cells

The toxic effect of strained hydrocarbon 2,2'-bis (bicyclo[2.2.1]heptane) (BBH) was studied using whole-cell bacterial lux-biosensors based on Escherichia coli cells in which luciferase genes are transcriptionally fused with stress-inducible promoters. It was shown that BBH has the genotoxic effect causing bacterial SOS response however no alkylating effect has been revealed. In addition to DNA damage, there is an oxidative effect causing the response of OxyR/S and SoxR/S regulons. The most sensitive to BBH lux-biosensor was E. coli pSoxS-lux which reacts to the appearance of superoxide anion radicals in the cell. It is assumed that the oxidation of BBH leads to the generation of reactive oxygen species, which provide the main contribution to the genotoxicity of this substance.

The multi-channel reaction of the OH radical with 5-hydroxymethylcytosine: a computational study

The hydroxyl radical may attack the new cytosine derivative 5-hydroxymethylcytosine (5-hmCyt), causing DNA oxidative damage. Two distinct mechanisms have been explored and our results provide some evidence between 5-hmCyt and tumor development.

Effects of an acid–alkaline environment on the reactivity of 5-carboxycytosine with hydroxyl radicals

The addition of ˙OH to C5C6 double bond and abstraction of H5 from 5-caCyt are more favourable in neutral, acidic and alkaline conditions.

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.

Kinetics of homoallylic/homobenzylic rearrangement reactions under combustion conditions

Homoallylic/homobenzylic radicals refer to typical radicals with the radical site located at the β position from the vinyl/phenyl group. These radicals are largely involved in combustion systems, such as the pyrolysis or oxidation of alkenes, cycloalkanes, and aromatics. The 1,2-vinyl/phenyl migration via two steps (cyclization/fission) is a peculiar reaction type for the homoallylic/homobenzylic radicals, entitled homoallylic/homobenzylic rearrangement, which has been studied by theoretical calculations including the Hirshfeld atomic charge analysis in the present work. With the help of rate constant calculations, the competition between this reaction channel and other possible pathways under combustion temperatures (500-2000 K) were evaluated. Analogous 1,3- and 1,4-vinyl/phenyl migration reactions for similar radicals with the radical sites located at the γ and δ positions from the vinyl/phenyl group were also computed. The results indicate that the 1,2-vinyl/phenyl migration is particularly important for the kinetics of unimolecular reactions of homoallylic radicals under 1500 K; nevertheless, it still has noticeable contribution at higher temperature. For those radicals with the radical site at the γ or δ positions, the respective 1,3- or 1,4-vinyl/phenyl migration channel plays an insignificant role under combustion conditions.

EXPERIMENTAL AND MODELING STUDY OF PREMIXED LAMINAR FLAMES OF ETHANOL AND METHANE

To better understand the chemistry of the combustion of ethanol, the structure of five low pressure laminar premixed flames has been investigated: a pure methane flame (φ=1), three pure ethanol flames (φ=0.7, 1.0, and 1.3), and an ethanol/methane mixture flames (φ=1). The flames have been stabilized on a burner at a pressure of 6.7 kPa using argon as dilutant, with a gas velocity at the burner of 64.3 cm/s at 333 K. The results consist of mole fraction profiles of 20 species measured as a function of the height above the burner by probe sampling followed by online gas chromatography analyses. A mechanism for the oxidation of ethanol was proposed. The reactions of ethanol and acetaldehyde were updated and include recent theoretical calculations while that of ethenol, dimethyl ether, acetone, and propanal were added in the mechanism. This mechanism was also tested against experimental results available in the literature for laminar burning velocities and laminar premixed flame where ethenol was detected. The main reaction pathways of consumption of ethanol are analyzed. The effect of the branching ratios of reaction C2H5OH+OH→Products+H2O is also discussed.

A mechanistic and kinetic study on the decomposition of morpholine

The combustion chemistry of morpholine (C(4)H(8)ONH) has been experimentally investigated recently as a representative model compound for O- and N-containing structural entities in biomass. Detailed profiles of species indicate the self-breakdown reactions prevailing over oxidative decomposition reactions. In this study, we derive thermodynamic and kinetic properties pertinent to all plausible reactions involved in the self-decomposition of morpholine and its derived morphyl radicals as a crucial task in the development of comprehensive combustion mechanism. Potential energy surfaces have been mapped out for the decomposition of morpholine and the three morphyl radicals. RRKM-based calculations predict the self-decomposition of morpholine to be dominated by 1,3-intramolecular hydrogen shift into the NH group at all temperatures and pressures. Self-decomposition of morpholine is shown to provide pathways for the formation of the experimentally detected products such as ethenol and ethenamine. Energetic requirements of all self-decomposition of morphyl radicals are predicted to be of modest values (i.e., 20-40 kcal/mol) which in turn support the occurrence of breaking-down reactions into two-heavy-atom species and the generation of doubly unsaturated four-heavy-atom segments. Calculated thermochemical parameters (in terms of standard enthalpies of formation, standard entropies, and heat capacities) and kinetic parameters (in terms of reaction rate constants at a high pressure limit) should be instrumental in building a robust kinetic model for the oxidation of morpholine.

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).

Shock tube and theory investigation of cyclohexane and 1-hexene decomposition

The decomposition of cyclohexane (c-C(6)H(12)) was studied in a shock tube using the laser-schlieren technique over the temperature range 1300-2000 K and for 25-200 Torr in mixtures of 2%, 4%, 10%, and 20% cyclohexane in Kr. Vibrational relaxation of the cyclohexane was also examined in 10 experiments covering 1100-1600 K for pressures below 20 Torr, and relaxation was found to be too fast to allow resolution of incubation times. The dissociation of 1-hexene (1- C(6)H(12)), apparently the sole initial product of cyclohexane decomposition, was also studied over 1220-1700 K for 50 and 200 Torr using 2% and 3% 1-hexene in Kr. On heating, cyclohexane simply isomerizes to 1-hexene, and this then dissociates almost entirely by a more rapid C-C scission to allyl and n-propyl radicals. This two-step reaction results in an initial small density gradient from the slight endothermicity of the isomerization. The gradient then rises strongly as the product 1-hexene dissociates. For the lower temperatures, this behavior is fully resolved here. For the higher pressures, 1-hexene decomposition generates negative gradients (exothermic reaction) as the radicals formed begin to recombine. Cyclohexane also generates such gradients, but these are now much smaller because the radical pool is depleted by abstraction from the reactant. A complete mechanism for the 1-hexene decomposition and for that of cyclohexane involving 79 reactions and 30 species is used in the final modeling of the gradients. Rate constants and RRKM fit parameters for the initial reactions are provided for the entire range of conditions. The possibility of direct reaction to allyl and n-propyl radicals, without stabilization of the intermediate 1-hexene, is examined down to pressures as low as 25 Torr, without a clear resolution of the issue. High-pressure limit rate constants from RRKM extrapolation are k(infinity)(c-C(6)H(12) --> 1-C(6)H(12)) = (8.76 x 10(17)) exp((-91.94 kcal/mol)/RT) s(-1) (T = 1300-2000 K) and k(infinity)(1-C(6)H(12) --> (*)C(3)H(7) + (*)C(3)H(5)) = (1.46 x 10(16)) exp((-69.12 kcal/mol)/RT) s(-1) (T = 1200-1700 K). This high-pressure rate for cyclohexane is entirely consistent with the notion that the isomerization involves initial C-C fission to a diradical. These extrapolated high-pressure rates are in good agreement with much of the literature.

Theoretical kinetic study of the reactions of cycloalkylperoxy radicals

Reactions of alkyl radicals with oxygen are key reactions in the low-temperature oxidation of hydrocarbons, but they have not been extensively studied yet in the case of cycloalkanes. Isomerizations of cycloalkylperoxy radicals and formation of cyclic ethers are especially important. In the present work, a theoretical study of the gas-phase reactions of cyclopentylperoxy and cyclohexylperoxy radicals has been carried out by means of quantum chemical calculations at the CBS-QB3 level. Computations on cyclopentylperoxy decomposition pathways are reported here for the first time. Thermochemical data have been obtained by means of isodesmic reactions, and the contribution of hindered rotors has been explicitly taken into account. Transition state theory has been used to calculate rate constants for all the elementary reactions. Three-parameter Arrhenius expressions have been derived in the temperature range 300-1000 K. Tunneling effects have been accounted for in the case of H-atom transfers. Our results compare well with experimental data and previous calculations available in the literature. In particular, the predicted rate constants for processes involving cyclohexylperoxy radicals, which have been introduced in a reaction mechanism scheme proposed before, exhibit excellent agreement with experiments at low and intermediate temperatures.