The Reaction of Acetylene with Hydroxyl Radicals

Dynamics studies for the multi-well and multi-channel reaction of OH with C2H2 on a full-dimensional global potential energy surface

The C2H2 + OH reaction is an important acetylene oxidation pathway in the combustion process, as well as a typical multi-well and multi-channel reaction. Here, we report an accurate full-dimensional machine learning-based potential energy surface (PES) for the C2H2 + OH reaction at the UCCSD(T)-F12b/cc-pVTZ-F12 level, based on about 475 000 ab initio points. Extensive quasi-classical trajectory (QCT) calculations were performed on the newly developed PES to obtain detailed dynamic data and analyze reaction mechanisms. Below 1000 K, the C2H2 + OH reaction produces H + OCCH2 and CO + CH3. With increasing temperature, the product channels H2O + C2H and H + HCCOH are accessible and the former dominates above 1900 K. It is found that the formation of H2O + C2H is dominated by a direct reaction process, while other channels belong to the indirect mechanism involving long-lived intermediates along the reaction pathways. At low temperatures, the C2H2 + OH reaction behaves like an unimolecular reaction due to the unique PES topographic features, of which the dynamic features are similar to the decomposition of energy-rich complexes formed by C2H2 + OH collision. The classification of trajectories that undergo different reaction pathways to generate each product and their product energy distributions were also reported in this work. This dynamic information may provide a deep understanding of the C2H2 + OH reaction.

Propylene Oxide Addition Effect on the Chemical Speciation of a Fuel-Rich Premixed n-Heptane/Toluene Flame

1,2-Propylene oxide (PO, C₃H₆O) is considered as a promising agent for improving fuel. In this work, the effect of PO additives on the species pool in a premixed burner-stabilized fuel-rich (ϕ = 1.6) flame fueled by n-heptane/toluene mixture (7/3 by volume of liquids) at atmospheric pressure is studied by the flame-sampling molecular beam mass spectrometry and numerical modeling in order to get insight into the chemical aspects of the influence of oxygenates with an epoxy group on the formation of abundant intermediates (including PAH precursors) during combustion of fossil fuels. The flames with various loadings of PO in the fuel blend (from 0 to 16.3% in mole basis) are examined, and detailed kinetic mechanisms available in the literature are validated against the measurements of mole fraction profiles of reactants, major products, and many intermediate species. A higher reactivity of the fresh mixture and a reduction in the peak mole fractions of intermediates playing an important role in PAH formation (benzene, styrene, ethylbenzene, phenol, acetylene, diacetylene, etc.) are observed when PO is added. This was found to be due to simultaneously two factors: the partial replacement of "sooting" fuel (toluene, which is the main precursor of these species) with oxygenated additive, and the changes in the flame radical pool caused by PO addition. Propylene oxide additive was found to change the ratio between H, OH, O, and CH₃ toward an increase in the proportion of O and CH₃. The detailed kinetic mechanisms considered in the work are found to overpredict the peak mole fraction of acetylene, a key species playing a crucial role in PAH growth. Its chemistry is revisited in order to provide a better prediction of C₂H₂ and, as a result, PAHs.

Computational studies on thermo-kinetics aspects of pyrolysis of isopropyl acetate and its methyl, bromide and hydroxyl derivatives

The gas-phase decomposition kinetics of isopropyl acetate (IPA) and its methyl, bromide and hydroxyl derivatives into the corresponding acid and propene were investigated using density functional theory (DFT) with the ωB97XD and M06–2x functionals, as well as the benchmark CBS-QB3 composite method. Transition state theory (TST) and RRKM theory calculations of rate constants under atmospheric pressure and in the fall-off regime were used to supplement the measured energy profiles. The results show that the formation of propene and bromoacetic acid is the most dominant pathway at the CBS-QB3 composite method, both kinetically and thermodynamically. There was a good agreement with experimental results. Pressures greater than 0.01 bar, corresponding to larger barrier heights are insufficient to ensure saturation of the measured rate coefficient when compared to the RRKM kinetic rates.

Natural bond orbitals (NBO) charges, bond orders, bond indices, and synchronicity parameters all point to the considered pathways taking place via a homogenous, first-order concerted, as well as an asynchronous mechanism involving a non-planar cyclic six-membered transition state. The calculated data exhibit that the elongation of the C_α−O bond length and subsequent polarization of the C_α^+δ…O^−δ bond is the rate-determining step of the considered reactions in the cyclic transition state, which appears to be involved in this type of reaction.

Theoretical Investigations on the Hydroxyl-Initialized Oxidation of Hexafluoro-2-butyne in the Presence of Oxygen

Hexafluoro-2-butyne (C₄F₆) is a potential eco-friendly alternative gas in plasma, refrigerants, and electrical insulation applications. Mechanisms for the reactions of C₄F₆ with OH/O₂ have been revealed in detail using various theoretical methods including ROCBS-Q, RCCSD(T), multireference RS2C, and extrapolations to the complete basis-set limit with Aug-cc-pVnZ (n = T, Q, 5) basis sets. Rate coefficients and product branching ratios were predicted for a wide range of temperatures and pressures using the solution of master equations. The vibrationally adiabatic ground-state barrier for the initial C₄F₆ + OH association was best estimated to be 1.53-2.26 kcal/mol. Energetically preferable decomposition paths for the conformation-dependent C₄F₆OH adducts include six-center HF elimination, four-center proton migration, and C-C bond cleavage, but the collisional deactivation is dominant under ambient conditions. The subsequent oxidation of C₄F₆OH by O₂ bifurcates in two orientations and proceeds without any well-defined barrier followed by the successive isomerization/elimination steps, forming perfluorobiacetyl to regenerate OH radicals or trifluoroacetic acid with trifluoroacetyl radicals. The OH-recycling path accounts for a branching ratio of 70% under ambient conditions. Theoretical rate coefficients are in good agreement with the available experimental results. The effect of fluorination on the reactivity of alkynes toward OH/O₂ is discussed.

Theoretical investigations on the HOMO–LUMO gap and global reactivity descriptor studies, natural bond orbital, and nucleus-independent chemical shifts analyses of 3-phenylbenzo[d]thiazole-2(3H)-imine and its para-substituted derivatives: Solvent and substituent effects

Natural bond orbital analysis, salvation, and substituent effects of electron-releasing (–CH3, –OH) and electron-withdrawing (–Cl, –NO2, –CF3) groups at para positions on the molecular structure of synthesized 3-phenylbenzo[ d]thiazole-2(3 H)-imine and its derivatives in selected solvents (acetone, toluene, and ethanol) and in the gas phase by employing the polarizable continuum method model are studied using the M06-2x method and 6-311++G(d,p) basis set. The relative stability of the studied compounds is influenced by the possibility of intramolecular interactions between substituents and the electron donor–acceptor centers of the thiazole ring. Furthermore, atomic charges, electron density, chemical thermodynamics, energetic properties, dipole moments, and nucleus-independent chemical shifts of the studied compounds and their relative stability are considered. The dipole moment values and the highest occupied molecular orbital–lowest unoccupied molecular orbital energy gaps reveal different charge-transfer possibilities within the considered molecules. Finally, natural bond orbital analysis is carried out to picture the charge transfer between the localized bonds and lone pairs.

Investigation of Protective Properties of Organic Layers toward Reactive Oxygen Species

The antioxidant protective properties of polyaromatic organic layers were evaluated toward reactive oxygen species (ROS) using scanning electrochemical microscopy in a foot-printing strategy. The layers were prepared by electrografting of aryldiazonium salts. Where p-(methyl)phenyl films show only weak protective properties toward ROS, p-(ethynyl)phenyl films evidence efficient protection of the covered surfaces. Applied potentials and electrolytes used during oxygen reduction reaction are critical parameters to control, prevent, or reduce the influence of ROS production and hence enhance the device lifetime.

A combined photoionization time-of-flight mass spectrometry and laser absorption spectrometry flash photolysis apparatus for simultaneous determination of reaction rates and product branching

In recent years, predictions of product branching for reactions of consequence to both combustion and atmospheric chemistry have outpaced validating experiments. An apparatus is described that aims to fill this void by combining several well-known experimental techniques into one: flash photolysis for radical generation, multiple-pass laser absorption spectrometry (LAS) for overall kinetics measurements, and time-resolved photoionization time-of-flight mass spectrometry (PI TOF-MS) for product branching quantification. The sensitivity of both the LAS and PI TOF-MS detection techniques is shown to be suitable for experiments with initial photolytically generated radical concentrations of ∼1 × 10¹² molecules cm^-3. As it is fast (μs time resolution) and non-intrusive, LAS is preferred for accurate kinetics (time-dependence) measurements. By contrast, PI TOF-MS is preferred for product quantification because it provides a near-complete picture of the reactor composition in a single mass spectrum. The value of simultaneous LAS and PI TOF-MS detection is demonstrated for the chemically interesting phenyl radical + propene system.

Quantification of Bioorthogonal Stability in Immune Phagocytes Using Flow Cytometry Reveals Rapid Degradation of Strained Alkynes

One of the areas in which bioorthogonal chemistry—chemistry performed inside a cell or organism—has become of pivotal importance is in the study of host–pathogen interactions. The incorporation of bioorthogonal groups into the cell wall or proteome of intracellular pathogens has allowed study within the endolysosomal system. However, for the approach to be successful, the incorporated bioorthogonal groups must be stable to chemical conditions found within these organelles, which are some of the harshest found in metazoans: the groups are exposed to oxidizing species, acidic conditions, and reactive thiols. Here we present an assay that allows the assessment of the stability of bioorthogonal groups within host cell phagosomes. Using a flow cytometry-based assay, we have quantified the relative label stability inside dendritic cell phagosomes of strained and unstrained alkynes. We show that groups that were shown to be stable in other systems were degraded by as much as 79% after maturation of the phagosome.

Radicals derived from acetaldehyde and vinyl alcohol

Vinyl alcohol and acetaldehyde are isoelectronic products of incomplete butanol combustion. Along with the radicals resulting from the removal of atomic hydrogen or the hydroxyl radical, these species are studied here using ab initio methods as complete as coupled cluster theory with single, double, triple, and perturbative quadruple excitations [CCSDT(Q)], with basis sets as large as cc-pV5Z. The relative energies provided herein are further refined by including corrections for relativistic effects, the frozen core approximation, and the Born-Oppenheimer approximation. The effects of anharmonic zero-point vibrational energies are also treated. The syn conformer of vinyl alcohol is predicted to be lower in energy than the anti conformer by 1.1 kcal mol^-1. The alcoholic hydrogen of syn-vinyl alcohol is found to be the easiest to remove, requiring 84.4 kcal mol^-1. Five other radicals are also carefully considered, with four conformers investigated for the 1-hydroxyvinyl radical. Beyond energetics, we have conducted an overhaul of the spectroscopic literature for these species. Our results also provide predictions for fundamental modes yet to be reported experimentally. To our knowledge, the ν₃ (3076 cm^-1) and ν₄ (2999 cm^-1) C-H stretches for syn-vinyl alcohol and all but one of the vibrational modes for anti-vinyl alcohol (ν₁-ν₁₄) are yet to be observed experimentally. For the acetyl radical, ν₆ (1035 cm^-1), ν₁₁ (944 cm^-1), ν₁₂ (97 cm^-1), and accounting for our changes to the assignment of the 1419.9 cm^-1 experimental mode, ν₁₀ (1441 cm^-1), are yet to be observed. We have predicted these unobserved fundamentals and reassigned the experimental 1419.9 cm^-1 frequency in the acetyl radical to ν₄ rather than to ν₁₀. Our work also strongly supports reassignment of the ν₁₀ and ν₁₁ fundamentals of the vinoxy radical. We suggest that the bands assigned to the overtones of these fundamentals were in fact combination bands. Our findings may be useful in constructing improved combustion models of butanol and in spectroscopically characterizing these molecules further.

An Experimental Study of the Kinetics of OH/OD(v = 1,2,3) + SO2: The Limiting High-Pressure Rate Coefficients as a Function of Temperature

The kinetics of the reaction OH/OD(v = 1,2,3) + SO₂ were studied using a photolysis/laser-induced fluorescence technique. The rate coefficients OH/OD(v = 1,2,3) + SO₂, k₁, over the temperature range of 295-810 K were used to determine the limiting high-pressure limit k₁^∞. This method is usually applicable if the reaction samples the potential well of the adduct HOSO₂ and if intramolecular vibrational relaxation is fast. In the present case, however, the rate coefficients showed an additional fast removal contribution as evidenced by the increase in k₁ with vibrational level; this behavior together with its temperature dependence is consistent with the existence of a weakly bound complex on the potential energy surface prior to adduct formation. The data were analyzed using a composite mechanism that incoporates energy-transfer mechanisms via both the adduct and the complex, and yielded a value of k₁^∞(295 K) equal to (7.2 ± 3.3) × 10^-13 cm³ molecule^-1 s^-1 (errors at 1σ), a factor of between 2 and 3 smaller than the current recommended IUPAC and JPL values of (2.0_-1.0^+2.0) and (1.6 ± 0.4) × 10^-12 cm³ molecule^-1 s^-1 at 298 K, respectively, although the error bars do overlap. k₁^∞ was observed to only depend weakly on temperature. Further evidence for a smaller k₁^∞ is presented in the companion paper.

Theoretical Investigation of the Radical-Radical Reaction of O((3)P) + C2H3 and Comparison with Gas-Phase Crossed-Beam Experiments

Herein, we present an ab initio study of the prototypal radical-radical reactions of ground-state atomic oxygen [O((3)P)] with the vinyl (C2H3) radical using density functional theory and a complete basis set model. Two distinctive pathways on the lowest doublet potential energy surfaces (PESs) were predicted to be in competition: addition and abstraction. The barrierless addition of O((3)P) to the hydrocarbon radicals leads to energy-rich intermediate formation followed by subsequent isomerization and decomposition to yield various products: CH2CO (ketene) + H, CO + CH3, C2HOH (acetylenol) + H, (3,1)CCHOH + H, H2O + C2H, (3,1)CH2 + HCO, H2CO (formaldehyde) + CH, C2H2 (acetylene) + OH, and (3,1)CCH2 + OH. The competing but minor H-atom abstraction mechanisms produce C2H2 + OH and (1,3)CCH2 + OH. The optimized structures of the reactants, products, intermediates, and transition states and the reaction mechanisms were obtained on the lowest doublet PESs. The major pathway was predicted to be the formation of CH2CO + H through the low-barrier, single-step cleavages of the addition intermediates. The Levine-Bernstein prior method, statistical surprisal approach, and microcanonical Rice-Ramsperger-Kassel-Marcus theory were applied to deduce the energy distributions of H atoms and OH products and quantitative rate constants. On the basis of the statistical theory and the population analysis, the predicted energy distributions were compared to the kinetic energy release of H and the preferential population of the Π(A') component of OH products reported in recent gas-phase crossed-beam investigations (Park, M. J.; Jang, S. C.; Choi, J. H. J. Chem. Phys. 2012, 137, 204311), and their kinetic and dynamic characteristics were discussed.

Probing the kinetic energy-release dynamics of H-atom products from the gas-phase reaction of O(3P) with vinyl radical C2H3

The gas-phase radical-radical reaction dynamics of ground-state atomic oxygen O((3)P) with vinyl radicals C2H3 has been studied by combining the results of vacuum-ultraviolet laser-induced fluorescence spectroscopy in a crossed beam configuration with ab initio calculations. The two radical reactants O((3)P) and C2H3 were produced by photolysis of NO2 and supersonic flash pyrolysis of C2H3I, respectively. Doppler profile analysis of the kinetic energy release of the nascent H-atom products from the title reaction O((3)P) + C2H3→ H((2)S) + CH2CO (ketene) revealed that the average translational energy of the products and the average fraction of the total available energy were 7.03 ± 0.30 kcal mol(-1) and 7.2%. The empirical data combined with CBS-QB3 level ab initio theory and statistical calculations demonstrated that the title oxygen-hydrogen exchange reaction is a major reaction channel, through an addition-elimination mechanism involving the formation of a short-lived, dynamical complex on the doublet potential energy surface. On the basis of systematic comparison with several exchange reactions of hydrocarbon radicals, the observed kinetic energy release can be explained in terms of the weak impulse at the moment of decomposition in the loose transition state with a product-like geometry and a small reverse barrier along the exit channel.

Theoretical study on the gas phase reaction of propargyl alcohol with hydroxyl radical

The reaction of propargyl alcohol with hydroxyl radical has been studied extensively at CCSD(T)/aug-cc-pVTZ//MP2/cc-pVTZ level. This is the first time to gain a conclusive insight into the reaction mechanism and kinetics for this important reaction in detail. Two reaction mechanisms were revealed, namely addition/elimination and hydrogen abstraction mechanism. The reaction mechanism confirms that OH addition to C≡C triple bond forms the chemically activated adducts, IM1 (·CHCOHCH2OH) and IM2 (CHOH·CCH2OH), and the hydrogen abstraction pathways (-CH2OH bonded to the carbon atom and alcohol hydrogen) may occur via low barriers. Harmonic model of Rice-Ramsperger-Kassel-Marcus theory and variational transition state theory are used to calculate the overall and individual rate constants over a wide range of temperatures and pressures. The calculated rate constants are in good agreement with the experimental data. At atmospheric pressure with Ar as bath gas, IM1 (·CHCOHCH2OH) and IM2 (CHOH·CCH2OH) formed by collisional stabilization are dominant in the low temperature range. The production of CHCCHOH + H2O via hydrogen abstraction becomes dominate at higher temperature. The fraction of IM3 (CH2COHCH2·O) is very significant over the moderate temperature range.

A combined crossed-beam and theoretical study of the reaction dynamics of O(3P) + C2H3 → C2H2 + OH: analysis of the nascent OH products with the preferential population of the Π(A') component

The gas-phase reaction dynamics of ground-state atomic oxygen [O((3)P) from the photo-dissociation of NO(2)] with vinyl radicals [C(2)H(3) from the supersonic flash pyrolysis of vinyl iodide, C(2)H(3)I] has been investigated using a combination of high-resolution laser-induced fluorescence spectroscopy in a crossed-beam configuration and ab initio calculations. Unlike the previous gas-phase bulk kinetic experiments by Baulch et al. [J. Phys. Chem. Ref. Data 34, 757 (2005)], a new exothermic channel of O((3)P) + C(2)H(3) → C(2)H(2) + OH (X (2)Π: υ" = 0) has been identified for the first time, and the population analysis shows bimodal nascent rotational distributions of OH products with low- and high-N" components with a ratio of 2.4:1. No spin-orbit propensities were observed, and the averaged ratios of Π(A('))∕Π(A") were determined to be 1.66 ± 0.27. On the basis of computations at the CBS-QB3 theory level and comparison with prior theory, the microscopic mechanisms responsible for the nascent populations can be understood in terms of two competing dynamical pathways: a direct abstraction process in the low-N" regime as the major pathway and an addition-complex forming process in the high-N" regime as the minor pathway. Particularly, during the bond cleavage process of the weakly bound van der Waals complex C(2)H(2)-OH, the characteristic pathway from the low dihedral-angle geometry was consistent with the observed preferential population of the Π(A') component in the nascent OH products. A molecular-level discussion of the reactivity, mechanism, and dynamical features of the title reaction are presented together with a comparison to gas-phase oxidation reactions of a series of prototypical hydrocarbon radicals.

Interception of excited vibrational quantum states by O2 in atmospheric association reactions

Bimolecular reactions in Earth's atmosphere are generally assumed to proceed between reactants whose internal quantum states are fully thermally relaxed. Here, we highlight a dramatic role for vibrationally excited bimolecular reactants in the oxidation of acetylene. The reaction proceeds by preliminary adduct formation between the alkyne and OH radical, with subsequent O(2) addition. Using a detailed theoretical model, we show that the product-branching ratio is determined by the excited vibrational quantum-state distribution of the adduct at the moment it reacts with O(2). Experimentally, we found that under the simulated atmospheric conditions O(2) intercepts ~25% of the excited adducts before their vibrational quantum states have fully relaxed. Analogous interception of excited-state radicals by O(2) is likely common to a range of atmospheric reactions that proceed through peroxy complexes.

Interactions between theory and experiment in the investigation of elementary reactions of importance in combustion

Elementary reactions are a central component of models of combustion processes. Rate constants and channel yields are needed for those models. Both experimental and theoretical methods used to determine such rate data are discussed in this tutorial review, which is of interest to reaction kinetics and combustion engineering communities. Applications to combustion present particular problems because the conditions required can be well outside the ranges of temperature and pressure accessible to experiment, and the rate data can show a complex dependence on conditions. Under these conditions, the interplay between theory and experiment becomes important.

Master equation methods in gas phase chemical kinetics

In this article, we discuss the application of master equation methods to problems in gas phase chemical kinetics. The focus is on reactions that take place over multiple, interconnected potential wells and on the dissociation of weakly bound free radicals. These problems are of paramount importance in combustion chemistry. To illustrate specific points, we draw on our experience with reactions we have studied previously.

Reaction of ethylene with hydroxyl radicals: a theoretical study

Ab initio calculations of portions of the C2H5O potential energy surface critical to the title reaction are presented. These calculations are based on QCISD geometries and frequencies and RQCISD(T) energies extrapolated to the complete-basis-set limit. Rate coefficients for the reaction of C2H4 with OH are calculated using this surface and the two transition-state model of Greenwald and co-workers [J. Phys. Chem. A 2005, 109, 6031] for the association of OH with C2H4. The present calculations reproduce most of the experimental data, including the temperature and pressure dependence of the rate coefficients, with only a small (0.4 kcal/mol) adjustment to the energy barrier for direct hydrogen abstraction. We confirm the importance of this channel above 800 K and find that a significant fraction of the total rate coefficient (approximately 10%) is due to the formation of vinyl alcohol above this temperature. Calculations of the vinyl alcohol channel are consistent with the recent observation of this molecule in low-pressure flames [Taatjes, C. A.; Hansen, N.; McIlroy, A.; Miller, J. A.; Senosiain, J. P.; Klippenstein, S. J.; Qi, F.; Sheng, L.; Zhang, Y.; Cool, T. A.; Wang, J.; Westmoreland, P. R.; Law, M. E.; Kasper, T.; Kohse-Höinghaus, K. Science 2005, 308, 1887] and suggest that this reaction should be included in hydrocarbon oxidation mechanisms.

H Atom Yields from the Reactions of CN Radicals with C2H2, C2H4, C3H6, trans-2-C4H8, and iso-C4H8

The kinetics and H atom channel yield at both 298 and 195 K have been determined for reactions of CN radicals with C2H2 (1.00+/-0.21, 0.97+/-0.20), C2H4 (0.96+/-0.032, 1.04+/-0.042), C3H6 (pressure dependent), iso-C4H8 (pressure dependent), and trans-2-C4H8 (0.039+/-0.019, 0.029+/-0.047) where the first figure in each bracket is the H atom yield at 298 K and the second is that at 195 K. The kinetics of all reactions were studied by monitoring both CN decay and H atom growth by laser-induced fluorescence at 357.7 and 121.6 nm, respectively. The results are in good agreement with previous studies where available. The rate coefficients for the reaction of CN with trans-2-butene and iso-butene have been measured at 298 and 195 K for the first time, and the rate coefficients are as follows: k298K=(2.93+/-0.23)x10(-10) cm3 molecule(-1) s(-1), k195K=(3.58+/-0.43)x10(-10) cm3 molecule(-1) s(-1) and k298K=(3.17+/-0.10)x10(-10) cm3 molecule(-1) s(-1), k195K=(4.32+/-0.35)x10(-10) cm3 molecule(-1) s(-1), respectively, where the errors represent a combination of statistical uncertainty (2sigma) and an estimate of possible systematic errors. A potential energy surface for the CN+C3H6 reaction has been constructed using G3X//UB3LYP electronic structure calculations identifying a number of reaction channels leading to either H, CH3, or HCN elimination following the formation of initial addition complexes. Results from the potential energy surface calculations have been used to run master equation calculations with the ratio of primary:secondary addition, the average amount of downward energy transferred in a collision DeltaEd, and the difference in barrier heights between H atom elimination and an H atom 1, 2 migration as variable parameters. Excellent agreement is obtained with the experimental 298 K H atom yields with the following parameter values: secondary addition complex formation equal to 80%, DeltaEd=145 cm(-1), and the barrier height for H atom elimination set 5 kJ mol(-1) lower than the barrier for migration. Finally, very low temperature master equation simulations using the best fit parameters have been carried out in an increased precision environment utilizing quad-double and double-double arithmetic to predict H and CH3 yields for the CN+C3H6 reaction at temperatures and pressures relevant to Titan. The H and CH3 yields predicted by the master equation have been parametrized in a simple equation for use in modeling.

The Reaction of n- and i-C4H5 Radicals with Acetylene

In this article, we discuss the reactions of i-C4H5 and n-C4H5 with acetylene. Both have been proposed as possible cyclization steps, forming benzene or fulvene, in rich flames burning aliphatic fuels. The relevant parts of the potential energy surface were determined from rQCISD(T) calculations extrapolated to the infinite-basis-set limit. Using this information in a Rice-Ramsperger-Kassel-Marcus-based master equation, we have calculated thermal rate coefficients and product distributions for both reactions as a function of temperature and pressure. The results are cast in forms that can be used in modeling, and the implications of the results for flame chemistry are discussed.

Experimental and Master Equation Study of the Kinetics of OH + C2H2: Temperature Dependence of the Limiting High Pressure and Pressure Dependent Rate Coefficients

The kinetics of the reaction OH + C2H2 have been studied using laser flash photolysis at 248 nm to generate OH radicals and laser-induced fluorescence to monitor OH removal. An attempt was made to use the rate coefficients OH (v = 1,2) + C2H2 to determine the limiting high-pressure rate coefficient, k(1a)(infinity), over the temperature range of 195-823 K. This method is usually applicable if the reaction samples the potential energy well of the adduct, HOC2H2, and if intramolecular vibrational relaxation is fast. In the present case, however, the rate coefficients for loss of the vibrationally excited states by reaction with C2H2 also contain a substantial contribution from nonreactive vibrational relaxation, which occurs via a mechanism that does not sample the adduct potential energy well but involves, at least at low temperatures, collisions that access a shallower, longer range van der Waals well. The data were analyzed using a composite mechanism that incorporates both reactive and nonreactive energy transfer mechanisms, which allows the determination of k(1a)(infinity)(T) for OH + C2H2 with satisfactory accuracy over the temperature range 195-823 K. The kinetics of the reaction OH (v = 0) + C2H2 were also studied in He over the range of conditions: 210-373 K and 5-760 Torr. A one-dimensional master equation (ME) analysis of the experimental data provided a further determination of k(1a)(infinity)(T) and also <DeltaE>(down) for He. Combining the two sets of results gives a consistent dataset for k(1a)(infinity) and the Arrhenius parameters A1ainfinity = 7.3 x 10(-12) cm(3) molecule(-1) s(-1) and E(1a)(infinity) = 5.3 kJ mol(-1), with <DeltaE>(down) = 150(T/300 K) cm(-1). Additional experiments were conducted at room temperature in N(2) and SF(6) by laser flash photolysis with cavity ring down spectroscopy, and ME calculations were then optimized for the pressure falloff in N(2) by varying the average downward energy transfer parameter (<DeltaE>(down)). The output from the best fit ME was parametrized using a modified Troe expression to provide rate data for use in atmospheric modeling.

Kinetics of the Reaction of Methyl Radical with Hydroxyl Radical and Methanol Decomposition

The CH3 + OH bimolecular reaction and the dissociation of methanol are studied theoretically at conditions relevant to combustion chemistry. Kinetics for the CH3 + OH barrierless association reaction and for the H + CH2OH and H + CH3O product channels are determined in the high-pressure limit using variable reaction coordinate transition state theory and multireference electronic structure calculations to evaluate the fragment interaction energies. The CH3 + OH --> 3CH2 + H2O abstraction reaction and the H2 + HCOH and H2 + H2CO product channels feature localized dynamical bottlenecks and are treated using variational transition state theory and QCISD(T) energies extrapolated to the complete basis set limit. The 1CH2 + H2O product channel has two dynamical regimes, featuring both an inner saddle point and an outer barrierless region, and it is shown that a microcanonical two-state model is necessary to properly describe the association rate for this reaction over a broad temperature range. Experimental channel energies for the methanol system are reevaluated using the Active Thermochemical Tables (ATcT) approach. Pressure dependent, phenomenological rate coefficients for the CH3 + OH bimolecular reaction and for methanol decomposition are determined via master equation simulations. The predicted results agree well with experimental results, including those from a companion high-temperature shock tube determination for the decomposition of methanol.

CH and C2 Measurements Imply a Radical Pool within a Pool in Acetylene Flames

Measured CH and C2 profiles show a striking resemblance as a function of time in a series of seven well-characterized fuel-rich (phi=1.2-2.0) non-sooting acetylene flames. This implied commonality and interrelationship are unexpected as these radicals have dissimilar chemical kinetic natures. As a result, a rigorous examination was undertaken of the behavior of each of the hydrocarbon species known to be present, C, CH, CH2, CH3, CH4, CHO, CHOH, CH2O, CH2OH, CH3O, CH3OH, C2, C2H, C2H2, CHCO, CH2CO, and C2O. This emphasized the main region where CH and C2 are observed (50-600 micros) and reduced the kinetic reactions to only those that operate efficiently and are dominant. It was immediately apparent that this region of the flame reflects the nature of a hydrogen flame heavily doped with CO and CO2 and containing traces of hydrocarbons. The radical species, H, OH, O, along with H2, H2O, and O2, form an important controlling radical pool that is in partial equilibrium, and the concentrations of each of the hydrocarbon radicals are minor to this, playing secondary roles. As a result, the dominant fast reactions are those between the hydrocarbons and the basic hydrogen/oxygen radicals. Hydrocarbon-hydrocarbon reactions are unimportant here at these equivalence ratios. CH and C2 are formed and destroyed on a sub-microsecond time scale so that their flame profiles are the reflection of a complex kinetically dynamic system. This is found to be the case for all of the hydrocarbon species examined. As might be expected, these rapidly form steady-state distributions. However, with the exceptions of C, CHO, CHOH, and CH2O, which are irreversibly being oxidized, the others all form an interconnected hydrocarbon pool that is under the control of the larger hydrogen radical pool. The hydrocarbon pool can rapidly adjust, and the CH and C2 decay together as the pool is drained. This is either by continuing oxidation in less rich mixtures, or in richer flames where this is negligible by the onset of hydrocarbon-hydrocarbon reactions. The implications of such a hydrocarbon pool are significant. It introduces a buffering effect on their distribution and provides the indirect connection between CH and C2. Moreover, because they are members of this radical pool, flame studies alone cannot answer questions concerning their specific importance in combustion other than their contributing role to this pool. The presence of such a pool modifies the exactness that is needed for kinetic mechanisms, and knowledge of every species in the system no longer is necessary. Furthermore, as rate constants become refined, it will allow for the calculation of the relative concentrations of the hydrocarbon species and facilitate reduced kinetic mechanisms. It provides an explanation for previous isotopically labeled experiments and illustrates the difficulty of exactly identifying in flames the role of individual species. It resolves the fact that differing kinetic models can show similar levels of accuracy and has implications for sensitivity analyses. It finally unveils the mechanism of the flame ionization detector and has implications for the differing interpretations of diamond formation mechanisms.

Methyl Vinyl Ketone+OH and Methacrolein+OH Oxidation Reactions: A Master Equation Analysis of the Pressure- and Temperature-Dependent Rate Constants

High-level electronic structure calculations and master equation analyses were carried out to obtain the pressure- and temperature-dependent rate constants of the methyl vinyl ketone+OH and methacrolein+OH reactions. The balance between the OH addition reactions at the high-pressure limit, the OH addition reactions in the fall-off region, and the pressure-independent hydrogen abstractions involved in these multiwell and multichannel systems, has been shown to be crucial to understand the pressure and temperature dependence of each global reaction. In particular, the fall-off region of the OH addition reactions contributes to the inverse temperature dependence of the rate constants in the Arrhenius plots, leading to pressure-dependent negative activation energies. The pressure dependence of the methyl vinyl ketone+OH reaction is clearly more important than in the case of the methacrolein+OH reaction owing to the weight of the hydrogen abstraction process in this second system. Comparison of the theoretical rate constants and the experimental measurements shows quite good agreement.

Reflected shock tube studies of high-temperature rate constants for OH + C2H2 and OH + C2H4

The reflected shock tube technique with multi-pass absorption spectrometric detection of OH-radicals at 308 nm (corresponding to a total path length of approximately 4.9 m) has been used to study the reactions, OH + C(2)H(2)--> products (1) and OH + C(2)H(4)--> C(2)H(3) + H(2)O (2). The present optical configuration gives a S/N ratio of approximately 1 at approximately 0.5-1.0 x 10(12) radicals cm(-3). Hence, kinetics experiments could be performed at [OH](0) = approximately 4-20 ppm thereby minimizing secondary reactions. OH was produced rapidly from the dissociations of either CH(3)OH or NH(2)OH (hydroxylamine). A mechanism was then used to obtain profile fits that agreed with the experiment to within <+/-5%. The derived Arrhenius expressions, in units of cm(3) molecule(-1) s(-1) are: k(1) = (1.03 +/- 0.24) x 10(-10) exp(-7212 +/- 417 K/T) for 1509-2362 K and k(2) = (10.2 +/- 5.8) x 10(-10) exp(-7411 +/- 871 K/T) for 1463-1931 K. The present study is the first ever direct measurement for reaction (1) at temperatures >1275 K while the present results extend the temperature range for (2) by approximately 700 K. These values are compared with earlier determinations and with recent theoretical calculations. The calculations agree with the present data for both reactions to within +/-10% over the entire T-range.