Chemically and Thermally Activated Decomposition of Secondary Butyl 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.

Theoretical Study of the Reaction of Hydrogen Atoms with Three Pentene Isomers: 2-Methyl-1-butene, 2-Methyl-2-butene, and 3-Methyl-1-butene

This paper presents a comprehensive potential energy surface (PES) for hydrogen atom addition to and abstraction from 2-methyl-1-butene, 2-methyl-2-butene, and 3-methyl-1-butene and the subsequent ß-scission and H atom transfer reactions. Thermochemical parameters for species on the Ċ₅H₁₁ potential energy surface (PES) were calculated as a function of temperature (298-2000 K), using a series of isodesmic reactions to determine the formation enthalpies. High-pressure limiting and pressure-dependent rate constants were calculated using Rice-Ramsperger-Kassel-Marcus theory with a one-dimensional master equation. A number of studies have highlighted the fact that C₅ intermediate species play a role in polyaromatic hydrocarbon formation and that a fuel's chemical structure can be key in understanding the intermediate species formed during fuel decomposition. Rate constant recommendations for both Ḣ atom addition to, and H-atom abstraction by Ḣ atoms from, linear and branched alkenes have subsequently been proposed by incorporating our earlier work on 1- and 2-pentene, and these can be used in mechanisms of larger alkenes for which calculations do not exist. The current set of rate constants for the reactions of Ḣ atoms with both linear and branched C₅ alkenes, including their chemically activated pathways, are the first available in the literature of any reasonable fidelity for combustion modeling and are important for gasoline mechanisms. Validation of our theoretical results with pyrolysis experiments of 2-methyl-1-butene, 2-methyl-2-butene, and 3-methyl-1-butene at 2 bar in a single pulse shock tube (SPST) were carried out, with satisfactory agreement observed.

Ultraviolet photodissociation dynamics of the n-butyl, s-butyl, and t-butyl radicals

Photodissociation dynamics of the jet-cooled n-butyl radical via the 3s Rydberg state and the s-butyl radical via the 3p Rydberg states in the ultraviolet region of 233 nm-258 nm, as well as the t-butyl radical via the 3d Rydberg states at 226 nm-244 nm, are studied using the high-n Rydberg atom time-of-flight technique. The H-atom photofragment yield spectra of the n-butyl, s-butyl, and t-butyl radicals show a broad feature centered around 247 nm, 244 nm, and 234 nm, respectively. The translational energy distributions of the H + C₄H₈ products, P(E_T)'s, of the three radicals are bimodal, with a slow (low E_T) component peaking at ∼6 kcal/mol and a fast (high E_T) component peaking at ∼52 kcal/mol-57 kcal/mol, ∼43 kcal/mol, and ∼37 kcal/mol for n-butyl, s-butyl, and t-butyl, respectively. The fraction of the products' translational energy in the available energy, ⟨ f_T⟩, is 0.31, 0.30, and 0.27 for n-butyl, s-butyl, and t-butyl, respectively. The H-atom product angular distributions of the slow component are isotropic for all three radicals, while those of the fast component are anisotropic for n-butyl and s-butyl with an anisotropy parameter β ∼ 0.7 and ∼ 0.3 and that of the fast component of t-butyl is nearly isotropic. The bimodal product translational energy and angular distributions indicate two dissociation pathways to the H + C₄H₈ products in these three radicals, a direct, prompt dissociation on the repulsive potential energy surface coupling with the Rydberg excited states, and a unimolecular dissociation of the hot radical on the ground electronic state after internal conversion from the Rydberg states.

Theoretical, Experimental, and Modeling Study of the Reaction of Hydrogen Atoms with 1- and 2-Pentene

Alkyl radicals are prominent in combustion chemistry as they are formed by hydrocarbon decomposition or from a radical attack on hydrocarbons. Accurate determinations of the thermochemistry and kinetics of their unimolecular isomerization and decomposition reactions and related addition reactions of alkenes are therefore important in simulating the combustion chemistry of virtually all hydrocarbon fuels. In this work, a comprehensive potential energy surface (PES) for Ḣ-atom addition to and abstraction from 1- and 2-pentene, and the subsequent C-C and C-H β-scission reactions, and H-atom transfer reactions has been considered. Thermochemical values for the species on the Ċ₅H₁₁ PES were calculated as a function of temperature (298-2000 K), with enthalpies of formation determined using a network of isodesmic reactions. High-pressure limiting and pressure-dependent rate constants were calculated using the Rice-Ramsperger-Kassel-Marcus theory coupled with a one-dimensional master equation. As a validation of our theoretical results, hydrogen atomic resonance absorption spectrometry experiments were performed on the Ḣ-atom addition and abstraction reactions of 1- and 2-pentene. By incorporating our calculations into a detailed chemical kinetic model (AramcoMech 3.0), excellent agreement with these experiments is observed. The theoretical results are further validated via a comprehensive series of simulations of literature data. Our a priori model is found to reproduce important absolute species concentrations and product ratios reported therein.

100 Years of Progress in Gas-Phase Atmospheric Chemistry Research

Remarkable progress has occurred over the last 100 years in our understanding of atmospheric chemical composition, stratospheric and tropospheric chemistry, urban air pollution, acid rain, and the formation of airborne particles from gas-phase chemistry. Much of this progress was associated with the developing understanding of the formation and role of ozone and of the oxides of nitrogen, NO and NO2, in the stratosphere and troposphere. The chemistry of the stratosphere, emerging from the pioneering work of Chapman in 1931, was followed by the discovery of catalytic ozone cycles, ozone destruction by chlorofluorocarbons, and the polar ozone holes, work honored by the 1995 Nobel Prize in Chemistry awarded to Crutzen, Rowland, and Molina. Foundations for the modern understanding of tropospheric chemistry were laid in the 1950s and 1960s, stimulated by the eye-stinging smog in Los Angeles. The importance of the hydroxyl (OH) radical and its relationship to the oxides of nitrogen (NO and NO2) emerged. The chemical processes leading to acid rain were elucidated. The atmosphere contains an immense number of gas-phase organic compounds, a result of emissions from plants and animals, natural and anthropogenic combustion processes, emissions from oceans, and from the atmospheric oxidation of organics emitted into the atmosphere. Organic atmospheric particulate matter arises largely as gas-phase organic compounds undergo oxidation to yield low-volatility products that condense into the particle phase. A hundred years ago, quantitative theories of chemical reaction rates were nonexistent. Today, comprehensive computer codes are available for performing detailed calculations of chemical reaction rates and mechanisms for atmospheric reactions. Understanding the future role of atmospheric chemistry in climate change and, in turn, the impact of climate change on atmospheric chemistry, will be critical to developing effective policies to protect the planet.

Atmospheric Reaction of Cl with 4-Hydroxy-2-pentanone (4H2P): A Theoretical Study

The kinetics and the mechanism of the reaction of 4-hydroxy-2-pentanone (4H2P) with Cl atom were investigated using quantum theoretical calculations. Density functional theory, CBS-QB3, and G3B3 methods are used to explore the reaction pathways. Rice-Ramsperger-Kassel-Marcus theory is employed to obtain rate constants of the reaction at atmospheric pressure and the temperature range 278-400 K. This study provides the first theoretical and kinetic determination of Cl rate constant for reactions with 4H2P over a large temperature range. The obtained rate constant 1.47 × 10^-10 cm³ molecule^-1 s^-1 at 298 K is in reasonable agreement with those obtained for C4-C5 hydroxyketones both theoretically and experimentally. The results regarding the structure-reactivity relationship and the atmospheric implications are discussed.

Theoretical Study of the Reactions of Methane and Ethane with Electronically Excited N2(A(3)Σu(+))

Comprehensive quantum chemical analysis with the usage of density functional theory and post-Hartree-Fock approaches were carried out to study the processes in the N2(A(3)Σu(+)) + CH4 and N2(A(3)Σu(+)) + C2H6 systems. The energetically favorable reaction pathways have been revealed on the basis of the examination of potential energy surfaces. It has been shown that the reactions N2(A(3)Σu(+)) + CH4 and N2(A(3)Σu(+)) + C2H6 occur with very small or even zero activation barriers and, primarily, lead to the formation of N2H + CH3 and N2H + C2H5 products, respectively. Further, the interaction of these species can give rise the ground state N2(X(1)Σg(+)) and CH4 (or C2H6) products, i.e., quenching of N2(A(3)Σu(+)) by CH4 and C2H6 molecules is the complex two-step process. The possibility of dissociative quenching in the course of the interaction of N2(A(3)Σu(+)) with CH4 and C2H6 molecules has been analyzed on the basis of RRKM theory. It has been revealed that, for the reaction of N2(A(3)Σu(+)) with CH4, the dissociative quenching channel could occur with rather high probability, whereas in the N2(A(3)Σu(+)) + C2H6 reacting system, an analogous process was little probable. Appropriate rate constants for revealed reaction channels have been estimated by using a canonical variational theory and capture approximation. The estimations showed that the rate constant of the N2(A(3)Σu(+)) + C2H6 reaction path is considerably greater than that for the N2(A(3)Σu(+)) + CH4 one.

Pressure dependence and branching ratios in the decomposition of 1-pentyl radicals: shock tube experiments and master equation modeling

The decomposition and intramolecular H-transfer isomerization reactions of the 1-pentyl radical have been studied at temperatures of 880 to 1055 K and pressures of 80 to 680 kPa using the single pulse shock tube technique and additionally investigated with quantum chemical methods. The 1-pentyl radical was generated by shock heating dilute mixtures of 1-iodopentane and the stable products of its decomposition have been observed by postshock gas chromatographic analysis. Ethene and propene are the main olefin products and account for >97% of the carbon balance from 1-pentyl. Also produced are very small amounts of (E)-2-pentene, (Z)-2-pentene, and 1-butene. The ethene/propene product ratio is pressure dependent and varies from about 3 to 5 over the range of temperatures and pressures studied. Formation of ethene and propene can be related to the concentrations of 1-pentyl and 2-pentyl radicals in the system and the relative rates of five-center intramolecular H-transfer reactions and β C-C bond scissions. The 3-pentyl radical, formed via a four-center intramolecular H transfer, leads to 1-butene and plays only a very minor role in the system. The observed (E/Z)-2-pentenes can arise from a small amount of beta C-H bond scission in the 2-pentyl radical. The current experimental and computational results are considered in conjunction with relevant literature data from lower temperatures to develop a consistent kinetics model that reproduces the observed branching ratios and pressure effects. The present experimental results provide the first available data on the pressure dependence of the olefin product branching ratio for alkyl radical decomposition at high temperatures and require a value of <ΔE(down)(1000 K)> = (675 ± 100) cm(-1) for the average energy transferred in deactivating collisions in an argon bath gas when an exponential-down model is employed. High pressure rate expressions for the relevant H-transfer reactions and β bond scissions are derived and a Rice Ramsberger Kassel Marcus/Master Equation (RRKM/ME) analysis has been performed and used to extrapolate the data to temperatures between 700 and 1900 K and pressures of 10 to 1 × 10(5) kPa.

Thermal decomposition of HN(3)

The two-channel thermal decomposition of hydrogen azide, HN(3), was studied computationally. The reaction produces triplet or singlet NH and N(2). A model of the reaction was created on the basis of the theoretical study of the reaction potential-energy surface and microscopic reaction rates by Besora and Harvey (Besora, M.; Harvey, J. N. J. Chem. Phys. 2008, 129, 044303) and the experimental data on the energy-dependent rate constants reported by Foy et al. (Foy, B. R.; Casassa, M. P.; Stephenson, J. C.; King, D. S. J. Chem. Phys. 1990, 92, 2782) The properties of the model were adjusted to fit the calculated k(E) dependence to the experimental one. The experiments on thermal decomposition of HN(3) described in the literature were analyzed via kinetic modeling; the results of the analysis demonstrate that all but one of the existing studies were affected by contributions from secondary kinetics. The model of the reaction was then used in master-equation calculations of the pressure effects and the value of the critical energy transfer parameter, DeltaE(down), was adjusted based on agreement with the experimental k(T,P) data. Finally, the model was used to determine pressure- and temperature-dependent rate constants for both channels of reaction 1, which do not conform to the traditional formalism of low-pressure-limit and falloff description. Uncertainties of the model and their influence on the calculated thermal rate constant values were analyzed. Finally, parametrized expression for rate coefficients were provided for a wide range of temperatures and pressures.

Direct Measurement and Theoretical Calculation of the Rate Coefficient for Cl + CH3 in the Range from T = 202−298 K

The rate coefficient has been measured under pseudo-first-order conditions for the Cl+CH3 association reaction at T=202, 250, and 298 K and P=0.3-2.0 Torr helium using the technique of discharge-flow mass spectrometry with low-energy (12-eV) electron-impact ionization and collision-free sampling. Cl and CH3 were generated rapidly and simultaneously by reaction of F with HCl and CH4, respectively. Fluorine atoms were produced by microwave discharge in an approximately 1% mixture of F2 in He. The decay of CH3 was monitored under pseudo-first-order conditions with the Cl-atom concentration in large excess over the CH3 concentration ([Cl]0/[CH3]0=9-67). Small corrections were made for both axial and radial diffusion and minor secondary chemistry. The rate coefficient was found to be in the falloff regime over the range of pressures studied. For example, at T=202 K, the rate coefficient increases from 8.4x10(-12) at P=0.30 Torr He to 1.8x10(-11) at P=2.00 Torr He, both in units of cm3 molecule-1 s-1. A combination of ab initio quantum chemistry, variational transition-state theory, and master-equation simulations was employed in developing a theoretical model for the temperature and pressure dependence of the rate coefficient. Reasonable empirical representations of energy transfer and of the effect of spin-orbit interactions yield a temperature- and pressure-dependent rate coefficient that is in excellent agreement with the present experimental results. The high-pressure limiting rate coefficient from the RRKM calculations is k2=6.0x10(-11) cm3 molecule-1 s-1, independent of temperature in the range from 200 to 300 K.

Evaluation of the Kinetic Data for Intramolecular 1,x-Hydrogen Shifts in Alkyl Radicals and Structure/Reactivity Predictions from the Carbocyclic Model for the Transition State

Experimental and computational kinetic data for the intramolecular 1,x-hydrogen shift in alkyl radicals are compiled in Arrhenius format for x = 2-5. Significant experimental disparity remains, especially for x = 2 and 3. Experimental data for radicals with tert centers or bearing spectator substituents are lacking for all x, and none exist for x = 6. The common use of the strain energy of the unsubstituted (x+1)-carbocycle to coarsely model the activation energy for the 1,x-shift is extended to explore more subtle differences in progressively methyl-substituted systems by use of molecular mechanics estimates of differences in strain between radicals and carbocycles. For x = 5 and 6, a sterically driven increase in E is predicted for shifts in the tert --> tert class that apparently runs counter to the behavior of bimolecular hydrogen transfers. In contrast, a sterically driven decrease in E is predicted to result from spectator methyl groups for the prim --> prim reaction class for all x. There is no experimental basis to test these predictions; fragmentary computational evidence lends some support to the second but is ambiguous concerning the first. Possible deficiencies in the use of carbocycles as transition state models are discussed.

Experimental and theoretical investigation of rate coefficients of the reaction S(P3)+OCS in the temperature range of 298–985K

The reaction S(3P)+OCS in Ar was investigated over the pressure range of 50-710 Torr and the temperature range of 298-985 K with the laser photolysis technique. S atoms were generated by photolysis of OCS with light at 248 nm from a KrF excimer laser; their concentration was monitored via resonance fluorescence excited by atomic emission of S produced from microwave-discharged SO2. At pressures less than 250 Torr, our measurements give k(298 K)=(2.7+/-0.5)x10(-15) cm3 molecule-1 s-1, in satisfactory agreement with a previous report by Klemm and Davis [J. Phys. Chem. 78, 1137 (1974)]. New data determined for 407-985 K connect rate coefficients reported previously for T>or=860 and T<or=478 K and show a non-Arrhenius behavior. Combining our results with data reported at high temperatures, we derived an expression k(T)=(6.1+/-0.3)x10(-18) T1.97+/-0.24 exp[-(1560+/-170)/T] cm3 molecule-1 s-1 for 298<or=TK<or=1680. At 298 K and P>or=500 Torr, the reaction rate was enhanced. Theoretical calculations at the G2M(CC2) level, using geometries optimized with the B3LYP6-311+G(3df) method, yield energies of transition states and products relative to those of the reactants. Rate coefficients predicted with multichannel Rice-Ramsperger-Kassel-Marcus (RRKM) calculations agree satisfactorily with experimental observations. According to our calculations, the singlet channel involving formation of SSCO followed by direct dissociation into S2(a 1Deltag)+CO dominates below 2000 K; SSCO is formed via intersystem crossing from the triplet surface. At low temperature and under high pressure the stabilization of OCS2, formed via isomerization of SSCO, becomes important; its formation and further reaction with S atoms partially account for the observed increase in the rate coefficient under such conditions.

Mechanism and Rate Constants for the Decomposition of 1-Pentenyl Radicals

This paper is concerned with the mechanisms and rate constants for the decomposition of 1-penten-3-yl, 1-penten-4-yl, and 1-penten-5-yl radicals. They are formed from radical attack on 1-pentene, which is an important decomposition product of normal alkyl radicals with more than 6 carbon atoms in combustion systems. This work is based on related data in the literature. These involve rate constants for the reverse radical addition process under high-pressure conditions, chemical activation experiments, and more recent direct studies. The high-pressure rate constants are based on detailed balance. The energy transfer effects and the pressure dependences of the rate constants are determined through the solution of the master equation and are projected to cover combustion conditions. The low barriers to these reactions make it necessary to treat these thermal reactions as open systems, as in chemical activation studies. The multiple reaction channels make the nature of the pressure effects different from those usually described in standard texts. The order of stability is 1-penten-3-yl approximately 1-penten-4-yl > 1-penten-5-yl and straddles those for the n-alkyl radicals. A key feature in these reactions is the effects traceable to allylic resonance. However, the 50 kJ/mol allylic resonance energy is not fully manifested. The important unsaturated products are 1,3-butadiene, the pentadienes, allyl radicals, and vinyl radicals. The results are compared with the recommendations in the literature, and significant differences are noted. Extensions to larger radicals with similar structures are discussed.

Kinetics of the unimolecular decomposition of the 2-chloroallyl radical

The thermal decomposition of the 2-chloroallyl radical, CH(2)CClCH(2) --> CH(2)CCH(2) + Cl (1), was studied using the laser photolysis/photoionization mass spectrometry technique. Rate constants were determined in time-resolved experiments as a function of temperature (720-840 K) and bath gas density ([He] = (3-12) x 10(16), [N(2)] = 6 x 10(16) molecule cm(-3)). C(3)H(4) was observed as a primary product of reaction 1. The rate constants of reaction 1 are in the falloff, close to the low-pressure limit, under the conditions of the experiments. The potential energy surface (PES) of reaction 1 was studied using a variety of quantum chemical methods. The results of the study indicate that the minimum energy path of the CH(2)CClCH(2) dissociation proceeds through a PES plateau corresponding to a weakly bound Cl-C(3)H(4) complex; a PES saddle point exists between the equilibrium CH(2)CClCH(2) structure and the Cl-C(3)H(4) complex. The results of quantum chemical calculations, the rate constant values obtained in the experimental study, and literature data on the reverse reaction of addition of Cl to allene were used to create a model of reactions 1 and -1. The experimental dependences of the rate constants on temperature and pressure were reproduced in RRKM/master equation calculations. The reaction model provides expressions for the temperature dependences of the high-pressure-limit and the low-pressure-limit rate constants and the falloff broadening factors (at T = 300-1600 K): k(infinity)(1) = 1.45 x 10(20)T(-1.75) exp(-19609 K/T) s(-1), k(infinity)(-)(1) = 8.94 x 10(-10)T(-0.40) exp(481 K/T) cm(3) molecule(-1) s(-1), k(1)(0)(He) = 5.01 x 10(-32)T(-12.02) exp(-22788 K/T) cm(3) molecule(-1) s(-1), k(1)(0)(N(2)) = 2.50 x 10(-32)T(-11.92) exp(-22756 K/T) cm(3) molecule(-1) s(-1), F(cent)(He) = 0.46 exp(-T/1001 K) + 0.54 exp(-T/996 K) + exp(-4008 K/T), and F(cent)(N(2)) = 0.37 exp(-T/2017 K) + 0.63 exp(-T/142 K) + exp(-4812 K/T). The experimental data are not sufficient to specify all the parameters of the model; consequently, some of the model parameters were obtained from quantum chemical calculations and from analogy with other reactions of radical decomposition. Thus, the parametrization is most reliable under conditions close to those used in the experiments.