Channel Specific Rate Constants Research Articles

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.

On the decomposition mechanism of propanal: rate constants evaluation and kinetic simulations

The reactivity of aldehydes has been the subject of considerable interest in chemical kinetics, with propanal often chosen as the representative species. Despite its relevance, the reactivity of propanal is currently estimated from analogy and fitting of experimental data measured in limited temperature and pressure ranges, while the few literature theoretical studies have focused more on the exploration the potential energy surface (PES) than on the estimation of rate constants. The purpose of this work is to reinvestigate the propanal decomposition kinetics using the ab initio transition state theory based master equation approach with the intent of: (1) Determining accurate rate constants of key reaction channels; (2) Updating and validating an existing kinetic model by simulating available experimental data on propanal pyrolysis. It is found that propanal decomposition at the initial stages of pyrolysis occurs through four unimolecular barrierless reactions to form CHO + C2H5, CH2CHO + CH3, CH3CHCHO + H, and CH3CH2CO + H, and a termolecular pathway leading to the formation of C2H4 + CO + H2. High pressure rate constants were determined for each barrierless reaction channel using Variable Reaction Coordinate Transition State Theory and used to estimate phenomenological temperature and pressure dependent channel specific rate constants integrating the 1 dimensional master equation over the whole PES. The decomposition rate constants so determined are in agreement with the few available experimental data and significantly faster than previous literature estimates. The estimated kinetic parameters were finally implemented into the CRECK kinetic mechanism, leading to an improved agreement with shock tube pyrolysis data from the literature.

Crossed-Beam and Theoretical Studies of the O(3P, 1D) + Benzene Reactions: Primary Products, Branching Fractions, and Role of Intersystem Crossing.

Reliable modeling of hydrocarbon oxidation relies critically on knowledge of the branching fractions (BFs) as a function of temperature (T) and pressure (p) for the products of the reaction of the hydrocarbon with atomic oxygen in its ground state, O(3P). During the past decade, we have performed in-depth investigations of the reactions of O(3P) with a variety of small unsaturated hydrocarbons using the crossed molecular beam (CMB) technique with universal mass spectrometric (MS) detection and time-of-flight (TOF) analysis, combined with synergistic theoretical calculations of the relevant potential energy surfaces (PESs) and statistical computations of product BFs, including intersystem crossing (ISC). This has allowed us to determine the primary products, their BFs, and extent of ISC to ultimately provide theoretical channel-specific rate constants as a function of T and p. In this work, we have extended this approach to the oxidation of one of the most important species involved in the combustion of aromatics: the benzene (C6H6) molecule. Despite extensive experimental and theoretical studies on the kinetics and dynamics of the O(3P) + C6H6 reaction, the relative importance of the C6H5O (phenoxy) + H open-shell products and of the spin-forbidden C5H6 (cyclopentadiene) + CO and phenol adduct closed-shell products are still open issues, which have hampered the development of reliable benzene combustion models. With the CMB technique, we have investigated the reaction dynamics of O(3P) + benzene at a collision energy (Ec) of 8.2 kcal/mol, focusing on the occurrence of the phenoxy + H and spin-forbidden C5H6 + CO and phenol channels in order to shed further light on the dynamics of this complex and important reaction, including the role of ISC. Concurrently, we have also investigated the reaction dynamics of O(1D) + benzene at the same Ec. Synergistic high-level electronic structure calculations of the underlying triplet/singlet PESs, including nonadiabatic couplings, have been performed to complement and assist the interpretation of the experimental results. Statistical (RRKM)/master equation (ME) computations of the product distribution and BFs on these PESs, with inclusion of ISC, have been performed and compared to experiment. In light of the reasonable agreement between the CMB experiment, literature kinetic experimental results, and theoretical predictions for the O(3P) + benzene reaction, the so-validated computational methodology has been used to predict (i) the BF between the C6H5O + H and C5H6 + CO channels as a function of collision energy and temperature (at 0.1 and 1 bar), showing that their increase progressively favors radical (phenoxy + H)-forming over molecule (C5H6 + CO and phenol stabilization)-forming channels, and (ii) channel-specific rate constants as a function of T and p, which are expected to be useful for improved combustion models.

Theoretical Study of the Extent of Intersystem Crossing in the O(3P) + C6H6 Reaction with Experimental Validation.

The extent of intersystem crossing in the O(3P) + C6H6 reaction, a prototypical system for spin-forbidden reactions in oxygenated aromatic molecules, is theoretically evaluated for the first time. Calculations are performed using nonadiabatic transition-state theory coupled with stochastic master equation simulations and Landau–Zener theory. It is found that the dominant intersystem crossing pathways connect the T2 and S0 potential energy surfaces through at least two distinct minimum-energy crossing points. The calculated channel-specific rate constants and intersystem crossing branching fractions differ from previous literature estimates and provide valuable kinetic data for the investigation of benzene and polycyclic aromatic hydrocarbons oxidation in interstellar, atmospheric, and combustion conditions. The theoretical results are supported by crossed molecular beam experiments with electron ionization mass-spectrometric detection and time-of-flight analysis at 8.2 kcal/mol collision energy. This system is a suitable benchmark for theoretical and experimental studies of intersystem crossing in aromatic species.

Combined Experimental and Theoretical Studies of the O(3P) + 1-Butene Reaction Dynamics: Primary Products, Branching Fractions, and Role of Intersystem Crossing.

Information on the detailed mechanism and dynamics (primary products, branching fractions (BFs), and channel specific rate constants as a function of temperature) for many important combustion reactions of O(3P) with unsaturated hydrocarbons is still lacking. We report synergistic experimental/theoretical studies on the mechanism and dynamics of the O(3P) + 1-C4H8 (1-butene) reaction by combining crossed molecular beam (CMB) experiments with soft electron ionization mass-spectrometric detection and time-of-flight analysis at 10.5 kcal/mol collision energy (Ec) to high-level ab initio electronic structure calculations of the underlying triplet and singlet potential energy surfaces (PESs) and statistical Rice-Ramsperger-Kassel-Marcus/Master Equation (RRKM/ME) computations of BFs including intersystem crossing (ISC). The reactive interaction of O(3P) with 1-butene is found to mainly break apart the 4-carbon atom chain, leading to the radical product channels ethyl + vinoxy (BF = 0.34 ± 0.11), methyl + C3H5O (BF = 0.28 ± 0.09), formyl + propyl (BF = 0.17 ± 0.05), ethyl + acetyl (BF = 0.014 ± 0.007), and butanal radical (ethylvinoxy) + H (BF = 0.013 ± 0.004), and molecular product channels formaldehyde + propenylidene/propene (BF = 0.15 ± 0.05) and butenone (ethyl ketene) + H2 (BF = 0.037 ± 0.015). As some of these products can only be formed via ISC from triplet to singlet PESs, from BFs an extent of ISC of 50% is inferred. This value is significantly larger than that recently observed for O(3P) + propene (22%) at similar Ec, underlying the question of how important it is to consider nonadiabatic effects for these and similar combustion reactions. Comparison of the derived BFs with those of statistical (RRKM/ME) simulations on the ab initio coupled triplet/singlet PESs shows good agreement, warranting the use of the RRKM/ME approach to provide information on the variation of BFs with temperature and to derive channel specific rate constants as a function of temperature (T) and pressure (p). Notably, ISC is predicted to decrease strongly with increasing temperature (from about 70% at 300 K to 46% at Ec = 10.5 kcal/mol, and about 1% at 2000 K). The present results lead to a detailed understanding of the complex reaction mechanism of O(3P) + 1-butene and, by providing channel specific rate constants as a function of T and p, should facilitate the improvement of current fossil-fuel (1-butene) as well as biofuel (1-butanol) combustion models.

Direct measurements of channel specific rate constants in OH + C3H8 illuminates prompt dissociations of propyl radicals

OH + molecules are an important class of reactions in combustion and atmospheric chemistry. Consequently, numerous studies have measured rate constants for these processes over an extended temperature range. A large majority of these experimental studies have utilized the decay of [OH] profiles (monitored either by absorption or laser-induced fluorescence) to obtain total rate constants. However, there are limited direct measurements of channel specific rate constants in this important class of reactions, particularly at combustion relevant temperatures. In the present experiments, we have directly measured site-specific rate constants for abstraction of the secondary CH bond in OH + C3H8 at high temperatures. Atomic resonance absorption spectrometry (ARAS) was used to monitor the formation of H-atoms from shock-heated mixtures of tert-butylhydroperoxide and C3H8 at high temperatures. Simulations for the experimental H-atom profiles are sensitive only to abstraction of the secondary CH bond leading to unambiguous measurements of the rate constants for this reaction. Over the T-range, 921 K < T < 1146 K, rate constants from the present experiments for OH + C3H8 → H2O + i-C3H7 can be represented by the Arrhenius expression,k=(3.935±1.387)×10−11exp(−1681±362K/T)cm3molecule−1s−1Simulations of the lower temperature data (T < 1000 K) indicate that the H-atom profiles are also influenced to a minor extent by the thermal dissociation of iso-propyl, i-C3H7 → H + C3H6, at short time-scales. Direct dynamics calculations were performed to examine in greater detail the potential role of prompt dissociations of i-C3H7 and n-C3H7 (formed from the title reaction) in interpreting the lower temperature (< 1000 K) data from the present work. These simulations suggest that prompt dissociation of propyl radicals does not influence the present experimental observations but has a minor influence on higher temperature combustion simulations.

Reaction Dynamics of O((3)P) + Propyne: II. Primary Products, Branching Ratios, and Role of Intersystem Crossing from Ab Initio Coupled Triplet/Singlet Potential Energy Surfaces and Statistical Calculations.

The mechanism of the O((3)P) + CH3CCH reaction was investigated using a combined experimental/theoretical approach. Experimentally the reaction dynamics was studied using crossed molecular beams (CMB) with mass-spectrometric detection and time-of-flight analysis at 9.2 kcal/mol collision energy. Theoretically master equation (ME) simulations were performed on a potential energy surface (PES) determined using high-level ab initio electronic structure calculations. In this paper (II) the theoretical results are described and compared with experiments, while in paper (I) are reported and discussed the results of the experimental study. The PES was investigated by determining structures and vibrational frequencies of wells and transition states at the CASPT2/aug-cc-pVTZ level using a minimal active space. Energies were then determined at the CASPT2 level increasing systematically the active space and at the CCSD(T) level extrapolating to the complete basis set limit. Two separate portions of the triplet PES were investigated, as O((3)P) can add either on the terminal or the central carbon of the unsaturated propyne bond. Minimum energy crossing points (MECPs) between the triplet and singlet PESs were searched at the CASPT2 level. The calculated spin-orbit coupling constants between the T1 and S0 electronic surfaces were ∼25 cm(-1) for both PESs. The portions of the singlet PES that can be accessed from the MECPs were investigated at the same level of theory. The system reactivity was predicted integrating stochastically the one-dimensional ME using Rice-Ramsperger-Kassel-Marcus theory to determine rate constants on the full T1/S0 PESs, accounting explicitly for intersystem crossing (ISC) using the Landau-Zener model. The computational results are compared both with the branching ratios (BRs) determined experimentally in the companion paper (I) and with those estimated in a recent kinetic study at 298 K. The ME results allow to interpret the main system reactivity: CH3CCO + H and CH3 + HCCO are the major channels active on the triplet PES and are formed from the wells accessed after O addition to the terminal and central C, respectively; (1)CH3CH + CO and C2H3 + HCO are the major channels active on the singlet PES and are formed from the methylketene and acrolein wells after ISC. However, also a large number of minor channels (∼15) are active, so that the system reactivity is quite complicated. The comparison between computational and experimental BRs is quite good for the kinetic study, while some discrepancy with the CMB estimations suggests that dynamic non-ergodic effects may influence the system reactivity. Channel specific rate constants are calculated in the 300-2250 K and 1-30 bar temperature and pressure ranges. It is found that as the temperature increases the H abstraction reaction, whose contribution is negligible in the experimental conditions, increases in relevance, and the extent of ISC decreases from ∼80% at 300 K to less than 2% at 2250 K.

Relevance of the Channel Leading to Formaldehyde + Triplet Ethylidene in the O((3)P) + Propene Reaction under Combustion Conditions.

Comprehension of the detailed mechanism of O((3)P) + unsaturated hydrocarbon reactions is complicated by the existence of many possible channels and intersystem crossing (ISC) between triplet and singlet potential energy surfaces (PESs). We report synergic experimental/theoretical studies of the O((3)P) + propene reaction by combining crossed molecular beams experiments using mass spectrometric detection at 9.3 kcal/mol collision energy (Ec) with high-level ab initio electronic structure calculations of the triplet PES and RRKM/master equation computations of branching ratios (BRs) including ISC. At high Ec's and temperatures higher than 1000 K, main products are found to be formaldehyde (H2CO) and triplet ethylidene ((3)CH3CH) formed in a reaction channel that has never been identified or considered significant in previous kinetics studies at 300 K and that, as such, is not included in combustion kinetics models. Global and channel-specific rate constants were computed and are reported as a function of temperature and pressure. This study shows that BRs of multichannel reactions useful for combustion modeling cannot be extrapolated from room-temperature kinetics studies.

Fulvenallene Decomposition Kinetics

While the decomposition kinetics of the benzyl radical has been studied in depth both from the experimental and the theoretical standpoint, much less is known about the reactivity of what is likely to be its main decomposition product, fulvenallene. In this work the high temperature reactivity of fulvenallene was investigated on a Potential Energy Surface (PES) consisting of 10 wells interconnected through 11 transition states using a 1 D Master Equation (ME). Rate constants were calculated using RRKM theory and the ME was integrated using a stochastic kinetic Monte Carlo code. It was found that two main decomposition channels are possible, the first is active on the singlet PES and leads to the formation of the fulvenallenyl radical and atomic hydrogen. The second requires intersystem crossing to the triplet PES and leads to acetylene and cyclopentadienylidene. ME simulations were performed calculating the microcanonical intersystem crossing frequency using Landau-Zener theory convoluting the crossing probability with RRKM rates evaluated at the conical intersection. It was found that the reaction channel leading to the cyclopentadienylidene diradical is only slightly faster than that leading to the fulvenallenyl radical, so that it can be concluded that both reactions are likely to be active in the investigated temperature (1500-2000 K) and pressure (0.05-50 bar) ranges. However, the simulations show that intersystem crossing is rate limiting for the first reaction channel, as the removal of this barrier leads to an increase of the rate constant by a factor of 2-3. Channel specific rate constants are reported as a function of temperature and pressure.

Channel Specific Rate Constants Relevant to the Thermal Decomposition of Disilane

Rate constants for the thermal dissociation of Si2H6 are predicted with a novel transition state model. The saddle points for dissociation on the Si2H6 potential energy surface are lower in energy than the corresponding separated products, as confirmed by high level ab initio quantum mechanical calculations. Thus, the dissociations of Si2H6 to produce SiH2 + SiH4 (R1) and H3SiSiH + H2 (R2) both proceed through tight inner transition states followed by loose outer transition states. The present "dual" transition state model couples variational phase space theory treatments of the outer transition states with ab initio based fixed harmonic vibrator treatments of the inner transition states to obtain effective numbers of states for the two transition states acting in series. It is found that, at least near room temperature, such a dual transition state model is generally required for the proper description of each of the dissociations. Only at quite high temperatures, i.e., above 2000 K for (R1) and 600 K for (R2), does a single fixed inner transition state provide an adequate description. Similarly, only at quite low temperatures (below 100 and 10 K for (R1) and (R2), respectively) does a single outer transition state provide an adequate description. Pressure dependent rate constants are obtained from solutions to the multichannel master equation. These calculations confirm that dissociation channel (R2) is negligible under conditions relevant to the thermal chemical vapor deposition (CVD) processes. Rate constants for the chemical activation reactions, SiH2 + SiH4 --> Si2H6 (R-1) and SiH2 + SiH4 --> H3SiSiH + H2 (R3), are also evaluated within the dual transition state model. It is found that reaction R3 is the dominant channel for low pressures and high temperatures, i.e., below 100 Torr for temperatures above 1100 K.

Kinetics and Mechanism of the BrO Self-Reaction: Temperature- and Pressure-Dependent Studies

The flash photolysis/UV absorption technique has been used to study the self-reaction of BrO radicals over the temperature range 222−298 K and the pressure range 100−760 Torr of N2 or O2. Two chemical sources of BrO radicals were used: photolysis of Br2 in the presence of excess ozone and photolysis of O2 in the presence of excess Br2. The overall rate constant, k1, for the BrO self-reaction (defined by −d[BrO]/dt = 2k1[BrO]2) was found to be temperature and pressure independent at T ≥ 250 K, with k1 = (2.88 ± 0.20) × 10-12 cm3 molecule-1 s-1. At temperatures below 250 K, k1 was found to be pressure dependent, due to the emergence of a new termolecular channel of the BrO self-reaction 1c, −1c forming the BrO dimer, Br2O2 (BrO + BrO + M ⇌ Br2O2 + M). Channel-specific rate constants were determined for the two bimolecular channels of the BrO self-reaction above 250 K, giving for (1a) (BrO + BrO → 2Br + O2) k1a = (5.31 ± 1.17) × 10-12 exp{(−211 ± 59)/T} cm3 molecule-1 s-1 and for (1b) (BrO + BrO → Br2 + O2) k1b = (1.13 ± 0.47) × 10-14 exp{(983 ± 111)/T} cm3 molecule-1 s-1. Below 250 K, the overall rate coefficient of the two bimolecular channels is reduced as the dimer forming channel emerges. At 235 and 222 K, rate constants for the formation (k1c) and decomposition (k-1c) of Br2O2 were recorded. Using the values for K1c, ΔHr for reaction 1c was estimated as −58.6 ± 0.1 kJ mol-1. A UV absorption spectrum attributed to Br2O2 was also recorded over the wavelength range 300−390 nm. The cross section of the smooth Br2O2 spectrum was found to be 1.2 × 10-17cm2 molecule-1 at 320 nm. These results are rationalized in terms of a mechanism of the BrO self-reaction that shows competition, at low temperatures, between collisional quenching and unimolecular dissociation of an excited BrOOBr* intermediate. The rate constant for the reaction of oxygen atoms with molecular bromine was also determined in the course of these experiments [O + Br2 → BrO + Br (5)], giving k5 = (5.12 ± 1.86) × 10-13 exp{(989 ± 91)/T} cm3 molecule-1 s-1. All errors are 1σ. The atmospheric implications of these results are discussed.

Collisions of O(1D) with HCl, Cl2, and COCl2: Total quenching, channel specific rate constants, and yields of Cl(P1/22)

The method of time-resolved laser magnetic resonance was used to monitor Cl atoms produced by the reactions of O(D1) with HCl, Cl2, and COCl2 at room temperature. For these reactions, the measured relative quantum yields of spin-orbit excited atoms are [Cl(P1/22)]/([Cl(P1/22)]+[Cl(P3/22)])=(0.10±0.04), &amp;lt;0.05, &amp;lt;0.10, respectively. The rate constants for total removal of O(D1) by HCl, Cl2, and COCl2 have been determined to be (1.5±0.3), (2.5±0.5) and (2.6±0.5), respectively. The reactive rate constants have been determined to be k2=(0.94±0.2), k4=(1.9±0.3), and k5a+2k5b=(2.1±0.27); here k2, k4, k5a, and k5b are rate constants of O(D1)+HCl → OH+Cl (2), O(D1)+Cl2 → ClO+Cl (4), and O(D1)+ COCl2 → CO+ClO+Cl (5a), CO2 +2Cl (5b), respectively. The rate constants are in units of 10−10 cm3/molecule⋅s, quoted errors are 2σ.

Channel specific rate constants for reactions of O(1D) with HCl and HBr

Time resolved resonance fluorescence detection of O(3P) and H(2S) has been employed in conjunction with pulsed laser photolysis of O3/HX/He mixtures to determine absolute rate coefficients and product yields for reactions of O(1D) with HCl(1) and HBr(2) at 297 K. Total rate coefficients for O(1D) removal are found to be, in units of 10−10 cm3 molecule−1 s−1, k1=1.50±0.18 and k2=1.48±0.16, where the quoted errors are 2σ and represent precision only. The absolute accuracy of these rate coefficients is estimated to be ±20%. Formation of O(3P)+HX is found to account for 9%±5% of the overall rate of reaction (1) and 20%±7% of the overall rate of reaction (2). Formation of H(2S)+XO is found to account for 24%±5% of the overall rate of reaction (1) but &amp;lt;4.5% of the overall rate of reaction (2). By difference, the OH+X(2P) channel accounts for 67%±10% of the overall rate of reaction (1) and 80%±12% of the overall rate of reaction (2).