Dissociation of 1,1,1-trifluoroethane behind reflected shock waves: shock tube/time-of-flight mass spectrometry experiments

The dissociation of 1,1,1,-trifluoroethane, a potential non-RRKM reaction, has been studied at 600 and 1200 Torr and high temperatures (1500-1840 K) using a new shock tube/time-of-flight mass spectrometer (ST/TOF-MS). These data obtained by an independent method are in good agreement with the laser schlieren, LS, experiments of Kiefer et al. [J. Phys. Chem. A 2004, 108, 2443-2450] and extend the range of that experimental dataset. The data have been simulated by both standard RRKM calculations and the non-RRKM model reported by Kiefer et al. but with = 750 cm(-1). Both the RRKM and non-RRKM calculations provide equally good fits to the ST/TOF-MS data. Neither model simulates the combined ST/TOF-MS and LS datasets particularly well. However, the non-RRKM model predicts a pressure dependency closer to that observed in the experiments than the RRKM model.     

A shock tube and theoretical study on the pyrolysis of 1,4-dioxane

The dissociation of 1, 2 and 4% 1,4-dioxane dilute in krypton was studied in a shock tube using laser schlieren densitometry, LS, for 1550-2100 K with 56 ± 4 and 123 ± 3 Torr. Products were identified by time-of-flight mass spectrometry, TOF-MS. 1,4-dioxane was found to initially dissociate via C-O bond fission followed by nearly equal contributions from pathways involving 2,6 H-atom transfers to either the O or C atom at the scission site. The 'linear' species thus formed (ethylene glycol vinyl ether and 2-ethoxyacetaldehyde) then dissociate by central fission at rates too fast to resolve. The radicals produced in this fission break down further to generate H, CH(3) and OH, driving a chain decomposition and subsequent exothermic recombination. High-level ab initio calculations were used to develop a potential energy surface for the dissociation. These results were incorporated into an 83 reaction mechanism used to simulate the LS profiles with excellent agreement. Simulations of the TOF-MS experiments were also performed with good agreement for consumption of 1,4-dioxane. Rate coefficients for the overall initial dissociation yielded k(123Torr) = (1.58 ± 0.50) × 10(59) × T(-13.63) × exp(-43970/T) s(-1) and k(58Torr) = (3.16 ± 1.10) × 10(79) × T(-19.13) × exp(-51326/T) s(-1) for 1600 < T < 2100 K.     

High pressure pyrolysis of toluene. 1. Experiments and modeling of toluene decomposition

The pyrolysis of toluene, the simplest methyl-substituted aromatic molecule, has been studied behind reflected shock waves using a single pulse shock tube. Experiments were performed at nominal high pressures of 27 and 45 bar and spanning a wide temperature range from 1200 to 1900 K. A variety of stable species, ranging from small hydrocarbons to single ring aromatics (principal soot precursors such as phenylacetylene and indene) were sampled from the shock tube and analyzed using standard gas chromatographic techniques. A detailed chemical kinetic model with 262 reactions and 87 species was assembled to simulate the stable species profiles (specifically toluene, benzene and methane) from the current high-pressure pyrolysis data sets and shock tube-atomic resonance absorption spectrometry (ARAS) H atom profiles obtained from prior toluene pyrolysis experiments performed under similar high-temperature conditions and lower pressures from 1.5 to 8 bar. The primary steps in toluene pyrolysis represent the most sensitive and dominant reactions in the model. Consequently, in the absence of unambiguous direct experimental measurements, we have utilized recent high level theoretical estimates of the barrierless association rate coefficients for these primary reactions, C6H5 + CH3 --> C6H5CH3 (1a) and C6H5CH2 + H --> C6H5CH3 (1b) in the detailed chemical kinetic model. The available data sets can be successfully reconciled with revised values for deltaH0f(298K)(C6H5CH2) = 51.5 +/- 1.0 kcal/mol and deltaH0f(298K)(C6H5) = 78.6 +/- 1.0 kcal/mol that translate to primary dissociation rate constants, reverse of 1a and 1b, represented by k(-1a,infinity) = (4.62 x 10(25))T(-2.53)exp[-104.5 x 10(3)/RT] s(-1) and k(-1b,infinity) = (1.524 x 10(16))T(-0.04)exp[-93.5 x 10(3)/RT] s(-1) (R in units of cal/(mol K)). These high-pressure limiting rate constants suggest high-temperature branching ratios for the primary steps that vary from 0.39 to 0.52 over the temperature range 1200-1800 K.     

Experimental and modeling study of methyl cyclohexane pyrolysis and oxidation

Although the combustion chemistry of aliphatic hydrocarbons has been extensively documented, the oxidation of cyclic hydrocarbons has been studied to a much lesser extent. To provide a deeper understanding of the combustion chemistry of naphthenes, the oxidation of methylcyclohexane was studied in a series of high-temperature shock tube experiments. Ignition delay times for a series of mixtures, of varying methylcyclohexane/oxygen equivalence ratios (phi=0.5, 1.0, 2.0), were measured over reflected shock temperatures of 1200-2100 K and reflected shock pressures of 1.0, 2.0, and 4.0 atm. A detailed chemical kinetic mechanism has been assembled to simulate the shock tube results and flow reactor experiments, with good agreement observed.     

Thermal dissociation of ethylene glycol vinyl ether

The pyrolysis of ethylene glycol vinyl ether (EGVE), an initial product of 1,4-dioxane dissociation, was examined in a diaphragmless shock tube (DFST) using laser schlieren densitometry (LS) at 57 ± 2 and 122 ± 3 Torr over 1200-1800 K. DFST/time-of-flight mass spectrometry experiments were also performed to identify reaction products. EGVE was found to dissociate via two channels: (1) a molecular H atom transfer/C-O scission to produce C(2)H(3)OH and CH(3)CHO, and (2) a radical channel involving C-O bond fission generating ˙CH(2)CH(2)OH and ˙CH(2)CHO radicals, with the second channel being strongly dominant over the entire experimental range. A reaction mechanism was constructed for the pyrolysis of EGVE which simulates the LS profiles very well over the full experimental range. The decomposition of EGVE is clearly well into the falloff region for these conditions, and a Gorin model RRKM fit was obtained for the dominant radical channel. The results are in good agreement with the experimental data and suggest the following rate coefficient expressions: k(2,∞) = (6.71 ± 2.6) × 10(27) × T(-3.21)exp(-35512/T) s(-1); k(2)(120 Torr) = (1.23 ± 0.5) × 10(92) × T(-22.87)exp(-48?248/T) s(-1); k(2)(60 Torr) = (2.59 ± 1.0) × 10(88) × T(-21.96)exp(-46283/T) s(-1).     

Shock tube study of dissociation and relaxation in 1,1-difluoroethane and vinyl fluoride

This paper reports measurements of the thermal dissociation of 1,1-difluoroethane in the shock tube. The experiments employ laser-schlieren measurements of rate for the dominant HF elimination using 10% 1,1-difluoroethane in Kr over 1500-2000 K and 43 < P < 424 torr. The vinyl fluoride product of this process then dissociates affecting the late observations. We thus include a laser schlieren study (1717-2332 K, 75 < P < 482 torr in 10 and 4% vinyl fluoride in Kr) of this dissociation. This latter work also includes a set of experiments using shock-tube time-of-flight mass spectrometry (4% vinyl fluoride in neon, 1500-1980 K, 500 < P < 1300 torr). These time-of-flight experiments confirm the theoretical expectation that the only reaction in vinyl fluoride is HF elimination. The dissociation experiments are augmented by laser schlieren measurements of vibrational relaxation (1-20% C(2)H(3)F in Kr, 415-1975 K, 5 < P < 50 torr, and 2 and 5% C(2)H(4)F(2) in Kr, 700-1350 K, 6 < P < 22 torr). These experiments exhibit very rapid relaxation, and incubation delays should be negligible in dissociation. An RRKM model of dissociation in 1,1-difluoroethane based on a G3B3 calculation of barrier and other properties fits the experiments but requires a very large DeltaE(down) of 1600 cm(-1), similar to that found in a previous examination of 1,1,1-trifluoroethane. Dissociation of vinyl fluoride is complicated by the presence of two parallel HF eliminations, both three-center and four-center. Structure calculations find nearly equal barriers for these, and TST calculations show almost identical k(infinity). An RRKM fit to the observed falloff again requires an unusually large DeltaE(down) and the experiments actually support a slightly reduced barrier. These large energy-transfer parameters now seem routine in these large fluorinated species. It is perhaps a surprising result for which there is as yet no explanation.     

Shock tube and theoretical studies on the thermal decomposition of propane: evidence for a roaming radical channel

The thermal decomposition of propane has been studied using both shock tube experiments and ab initio transition state theory-based master equation calculations. Dissociation rate constants for propane have been measured at high temperatures behind reflected shock waves using high-sensitivity H-ARAS detection and CH(3) optical absorption. The two major dissociation channels at high temperature are C(3)H(8) → CH(3) + C(2)H(5) (eq 1a) and C(3)H(8) → CH(4) + C(2)H(4) (eq 1b). Ultra high-sensitivity ARAS detection of H-atoms produced from the decomposition of the product, C(2)H(5), in (1a), allowed measurements of both the total decomposition rate constants, k(total), and the branching to radical products, k(1a)/k(total). Theoretical analyses indicate that the molecular products are formed exclusively through the roaming radical mechanism and that radical products are formed exclusively through channel 1a. The experiments were performed over the temperature range 1417-1819 K and gave a minor contribution of (10 ± 8%) due to roaming. A multipass CH(3) absorption diagnostic using a Zn resonance lamp was also developed and characterized in this work using the thermal decomposition of CH(3)I as a reference reaction. The measured rate constants for CH(3)I decomposition agreed with earlier determinations from this laboratory that were based on I-atom ARAS measurements. This CH(3) diagnostic was then used to detect radicals from channel 1a allowing lower temperature (1202-1543 K) measurements of k(1a) to be determined. Variable reaction coordinate-transition state theory was used to predict the high pressure limits for channel (1a) and other bond fission reactions in C(3)H(8). Conventional transition state theory calculations were also used to estimate rate constants for other tight transition state processes. These calculations predict a negligible contribution (<1%) from all other bond fission and tight transition state processes, indicating that the bond fission channel (1a) and the roaming channel (1b) are indeed the only active channels at the temperature and pressure ranges of the present experiments. The predicted reaction exo- and endothermicities are in excellent agreement with the current version of the Active Thermochemical Tables. Master equation calculations incorporating these transition state theory results yield predictions for the temperature and pressure dependence of the dissociation rate constants for channel 1a. The final theoretical results reliably reproduce the measured dissociation rate constants that are reported here and in the literature. The experimental data are well reproduced over the 500-2500 K and 1 × 10(-4) to 100 bar range (errors of ～15% or less) by the following Troe parameters for Ar as the bath gas: k(∞) = 1.55 × 10(24)T(-2.034) exp(-45?490/T) s(-1), k(0) = 7.92 × 10(53)T(-16.67) exp(-50?380/T) cm(3) s(-1), and F(c) = 0.190 exp(-T/3091) + 0.810 exp(-T/128) + exp(-8829/T).     

High-temperature shock tube measurements of methyl radical decomposition

We have studied the two-channel thermal decomposition of methyl radicals in argon, involving the reactions CH3 + Ar --> CH + H2 + Ar (1a) and CH3 + Ar --> CH2 + H + Ar (1b), in shock tube experiments over the 2253-3527 K temperature range, at pressures between 0.7 and 4.2 atm. CH was monitored by continuous-wave, narrow-line-width laser absorption at 431.1311 nm. The collision-broadening coefficient for CH in argon, 2gamma(CH-Ar), was measured via repeated single-frequency experiments in the ethane pyrolysis system behind reflected shock waves. The measured 2gamma(CH-Ar) value and updated spectroscopic and molecular parameters were used to calculate the CH absorption coefficient at 431.1311 nm (23194.80 cm(-1)), which was then used to convert raw traces of fractional transmission to quantitative CH concentration time histories in the methyl decomposition experiments. The rate coefficient of reaction 1a was measured by monitoring CH radicals generated upon shock-heating highly dilute mixtures of ethane, C2H6, or methyl iodide, CH3I, in an argon bath. A detailed chemical kinetic mechanism was used to model the measured CH time histories. Within experimental uncertainty and scatter, no pressure dependence could be discerned in the rate coefficient of reaction 1a in the 0.7-4.2 atm pressure range. A least-squares, two-parameter fit of the current measurements, applicable between 2706 and 3527 K, gives k(1a) (cm(3) mol(-1) s(-1)) = 3.09 x 1015 exp[-40700/T (K)]. The rate coefficient of reaction 1b was determined by shock-heating dilute mixtures of C2H6 or CH3I and excess O2 in argon. During the course of reaction, OH radicals were monitored using the well-characterized R(1)(5) line of the OH A-X (0,0) band at 306.6871 nm (32606.52 cm(-1)). H atoms generated via reaction 1b rapidly react with O2, which is present in excess, forming OH. The OH traces are primarily sensitive to reaction 1b, reaction 9 (H + O2 --> OH + O) and reaction 10 (CH3 + O2 --> products), where the rate coefficients of reactions 9 and 10 are relatively well-established. No pressure dependence could be discerned for reaction 1b between 1.1 and 3.9 atm. A two-parameter, least-squares fit of the current data, valid over the 2253-2975 K temperature range, yields the rate expression k(1b) (cm(3) mol(-1) s(-1)) = 2.24 x 10(15) exp[-41600/T (K)]. Theoretical calculations carried out using a master equation/RRKM analysis fit the measurements reasonably well.     

A shock tube study of the reaction of methyl radicals with hydroxyl radicals

The reaction of CH₃ with OH has been studied near 1200 K and 1 atmosphere pressure in shock tube experiments in which UV absorption was used to monitor [OH]. A rate coefficient of (1.1 ± 0.3) × 10¹³ cm³/mol-s was measured for removal of OH by CH₃. This measured value is compared with previous experimental data and calculations. Several possible reaction channels are discussed, and although products were not monitored, it seems probable, on the basis of other work and theoretical estimates, that the primary mechanism (?75%) for the removal of OH by CH₃ at these conditions is their combination to form CH₃OH. Rate coefficients of (5.3 ± 0.8) × 10¹² and (9.0 ± 1.4) × 10¹² cm³/mol-s were measured for the reactions of OH with acetone and ethane, respectively, at the same temperature and pressure.     

Reaction of hydrogen atoms with propyne at high temperatures: an experimental and theoretical study

The kinetics of the reaction of hydrogen atoms with propyne (pC3H4) was experimentally studied in a shock tube at temperatures ranging from 1200 to 1400 K and pressures between 1.3 and 4.0 bar with Ar as the bath gas. The hydrogen atoms (initial mole fraction 0.5-2.0 ppm) were produced by pyrolysis of C2H5I and monitored by atomic resonance absorption spectrometry under pseudo-first-order conditions with respect to propyne (initial mole fraction 5-20 ppm). From the hydrogen atom time profiles, overall rate coefficients k(ov) identical with -([pC3H4][H])(-1) x d[H]/dt for the reaction H + pC3H4 --> products ( not equal H) were deduced; the following temperature dependence was obtained: kov = 1.2 x 10(-10) exp(-2270 K/T) cm(3) s(-1) with an estimated uncertainty of +/-20%. A pressure dependence was not observed. The results are analyzed in terms of statistical rate theory with molecular and transition state data from quantum chemical calculations. Geometries were optimized using density functional theory at the B3LYP/6-31G(d) level, and single-point energies were computed at the QCISD(T)/cc-pVTZ level of theory. It is confirmed that the reaction proceeds via an addition-elimination mechanism to yield C2H2 + CH3 and via a parallel direct abstraction to give C3H3 + H2. Furthermore, it is shown that a hydrogen atom catalyzed isomerization channel to allene (aC3H4), H + pC3H4 --> aC3H4 + H, is also important. Kinetic parameters to describe the channel branching of these reactions are deduced.     

Thermal decomposition of monomethylhydrazine: Shock tube experiments and kinetic modeling

The thermal decomposition of gaseous monomethylhydrazine (MMH) was studied by recording MMH absorption at 220 nm of the reacting gas behind a reflected shock wave at temperatures of 900–1370 K, pressures of 140–450 kPa, and in mixtures containing 97.5–99 mol% argon. Based on previous work (Sun and Law; J Phys Chem A 2007, 111(19), 3748–3760), a kinetic mechanism was developed over extended temperature and pressure ranges to model these experimental data. Specifically, the temperature and pressure dependence of the unimolecular rate coefficients on the dissociation of MMH and the associated radicals were calculated by the QRRK/Master equation analysis at temperatures of 300–2000 K and pressures of 1–100 atm based on published thermochemical and kinetic parameters. They were then fitted using the Troe formalism and incorporated in the kinetic model. This unadjusted model was then used to predict the MMH decomposition profiles at different temperatures and pressures for seven groups of MMH/Ar mixtures and the half‐life decomposition times from shock tube experiments. Good agreement was observed below 940 K and above 1150 K for the diluted MMH/Ar mixtures. The model predictions further show that the overall MMH decomposition rate follows first‐order kinetics, and that the N–N bond scission is the most sensitive reaction path for the modeling of the homogeneous decomposition of MMH at elevated pressures. However, the model predictions deviate from the experimental data with the incubation period of ca. 100 μs observed in the 1030–1090 K temperature range, and it also predicts longer ignition delays for highly concentrated MMH/Ar mixtures. The discrepancy between the model predictions and experimental data at these special conditions of MMH decomposition was analyzed. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 41: 176–186, 2009     

Ion-molecule rate constants and branching ratios for the reaction of N(3) (+)+O(2) from 120 to 1400 K

The kinetics of the reaction of N(3) (+) with O(2) has been studied from 120 to 1400 K using both a selected ion flow tube and high-temperature flowing afterglow. The rate constant decreases from 120 K to approximately 1200 K and then increases slightly up to the maximum temperature studied, 1400 K. The rate constant compares well to most of the previous measurements in the overlapping temperature range. Comparing the results to drift tube data shows that there is not a large difference between increasing the translational energy available for reaction and increasing the internal energy of the reactants over much of the range, i.e., all types of energies drive the reactivity equally. The reaction produces both NO(+) and NO(2) (+), the latter of which is shown to be the higher energy NOO(+) linear isomer. The ratio of NOO(+) to NO(+) decreases from a value of over 2 at 120 K to less than 0.01 at 1400 K because of dissociation of NOO(+) at the higher temperatures. This ratio decreases exponentially with increasing temperature. High-level theoretical calculations have also been performed to compliment the data. Calculations using multi-reference configuration interaction theory at the MRCISD(Q)/cc-pVTZ level of theory show that singlet NOO(+) is linear and is 4.5 eV higher in energy than ONO(+). A barrier of 0.9 eV prevents dissociation into NO(+) and O((1)D); however, a crossing to a triplet surface connects to NO(+) and O((3)P) products. A singlet and a triplet potential energy surface leading to products have been determined using coupled cluster theory at the CCSD(T)/aug-cc-pVQZ level on structures optimized at the Becke3-Lee, Yang, and Parr (B3LYP)/aug-cc-pVTZ level of theory. The experimental results and reaction mechanism are evaluated using these surfaces.     

Kinetic modeling of propane oxidation and pyrolysis

Propane oxidation in jet-stirred reactor was modeled using a comprehensive kinetic reaction mechanism including the most recent findings concerning the kinetics of the reactions involved in the oxidation of C₁? C₄ hydrocarbons. The present detailed mechanism is able to reproduce experimental species concentration profiles obtained in our high-pressure jet-stirred reactor (900 ? T/K ? 1200; 1 ? P/atm ? 10; 0.15 ? ? ? 4) and in a turbulent flow reactor at 1 atm; ignition delay times measured in shock tube (1200 ? T/K ? 1700; 2 ? P/atm ? 15; 0.125 ? ? ? 2); H-atoms concentrations measured in shock tube during the pyrolysis of propane and burning velocities of freely propagating premixed propane-air laminar flames. The computed results are discussed in terms of pressure and equivalence ratio (?) effects on propane oxidation. The same detailed kinetic reaction mechanism can also be used to model the oxidation of methane, ethylene, ethane, and propene in similar conditions. © John Wiley & Sons, Inc.     

Kinetics, mechanism, and thermochemistry of the gas-phase reaction of atomic chlorine with pyridine

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the reaction of atomic chlorine with pyridine (C(5)H(5)N) as a function of temperature (215-435 K) and pressure (25-250 Torr) in nitrogen bath gas. At T> or = 299 K, measured rate coefficients are pressure independent and a significant H/D kinetic isotope effect is observed, suggesting that hydrogen abstraction is the dominant reaction pathway. The following Arrhenius expression adequately describes all kinetic data at 299-435 K for C(5)H(5)N: k(1a) = (2.08 +/- 0.47) x 10(-11) exp[-(1410 +/- 80)/T] cm(3) molecule(-1) s(-1) (uncertainties are 2sigma, precision only). At 216 K < or =T< or = 270 K, measured rate coefficients are pressure dependent and are much faster than computed from the above Arrhenius expression for the H-abstraction pathway, suggesting that the dominant reaction pathway at low temperature is formation of a stable adduct. Over the ranges of temperature, pressure, and pyridine concentration investigated, the adduct undergoes dissociation on the time scale of our experiments (10(-5)-10(-2) s) and establishes an equilibrium with Cl and pyridine. Equilibrium constants for adduct formation and dissociation are determined from the forward and reverse rate coefficients. Second- and third-law analyses of the equilibrium data lead to the following thermochemical parameters for the addition reaction: Delta(r)H = -47.2 +/- 2.8 kJ mol(-1), Delta(r)H = -46.7 +/- 3.2 kJ mol(-1), and Delta(r)S = -98.7 +/- 6.5 J mol(-1) K(-1). The enthalpy changes derived from our data are in good agreement with ab initio calculations reported in the literature (which suggest that the adduct structure is planar and involves formation of an N-Cl sigma-bond). In conjunction with the well-known heats of formation of atomic chlorine and pyridine, the above Delta(r)H values lead to the following heats of formation for C(5)H(5)N-Cl at 298 K and 0 K: Delta(f)H = 216.0 +/- 4.1 kJ mol(-1), Delta(f)H = 233.4 +/- 4.6 kJ mol(-1). Addition of Cl to pyridine could be an important atmospheric loss process for pyridine if the C(5)H(5)N-Cl product is chemically degraded by processes that do not regenerate pyridine with high yield.     

Shock tube investigation of CH3 + CH3OCH3

The title reaction has been investigated in a diaphragmless shock tube by laser schlieren densitometry over the temperature range 1163-1629 K and pressures of 60, 120, and 240 Torr. Methyl radicals were produced by dissociation of 2,3-butanedione in the presence of an excess of dimethyl ether. Rate coefficients for CH(3) + CH(3)OCH(3) were obtained from simulations of the experimental data yielding the following expression which is valid over the range 1100-1700 K: k = (10.19 ± 3.0)T(3.78)?exp((-4878/T)) cm(3) mol(-1)s(-1). The experimental results are in good agreement with estimates by Curran and co-workers [Fischer, S. L.; Dryer, F. L.; Curran, H. J. Int. J. Chem. Kinet.2000, 32 (12), 713-740. Curran, H. J.; Fischer, S. L.; Dryer, F. L. Int. J. Chem. Kinet.2000, 32 (12), 741-759] but about a factor of 2.6 lower than those of Zhao et al. [Zhao, Z.; Chaos, M.; Kazakov, A.; Dryer, F. L. Int. J. Chem. Kinet.2008, 40 (1), 1-18].     

A shock tube study of the reactions of the hydroxyl radical with several combustion species

Reactions of the hydroxyl radical, OH, with several organic species of interest in combustion chemistry have been studied near 1200 K and 1 atm in shock tube experiments in which UV absorption was used to monitor the OH concentration. Rate coefficients were measured for the reactions of OH with 2,3-dimethylbutane, isooctane, neooctane, ethylene, propylene, acetylene, formaldehyde, methanol, and ethanol. The values were found to be (in units of 10¹² cm³/mol-s): 21, 22, 18, 2.6, 9.6, 0.28, 12, 5.2, and 5.3. These measured values are compared with previous experimental results and, where appropriate, transition-state theory calculations.     

Reactions of 1-naphthyl radicals with ethylene. Single pulse shock tube experiments, quantum chemical, transition state theory, and multiwell calculations

The reactions of 1-naphthyl radicals with ethylene were studied behind reflected shock waves in a single pulse shock tube, covering the temperature range 950-1200 K at overall densities behind the reflected shocks of approximately 2.5 x 10(-5) mol/cm3. 1-Iodonaphthalene served as the source for 1-naphthyl radicals as its C-I bond dissociation energy is relatively small. It is only approximately 65 kcal/mol as compared to the C-H bond strength in naphthalene which is approximately 112 kcal/mol and can thus produce naphthyl radicals at rather low reflected shock temperatures. The [ethylene]/[1-iodo-naphthalene] ratio in all of the experiments was approximately 100 in order to channel the free radicals into reactions with ethylene rather than iodonaphthalene. Four products resulting from the reactions of 1-naphthyl radicals with ethylene were found in the post shock samples. They were vinyl naphthalene, acenaphthene, acenaphthylene, and naphthalene. Some low molecular weight aliphatic products at rather low concentrations, resulting from the attack of various free radicals on ethylene were also found in the shocked samples. In view of the relatively low temperatures employed in the present experiments, the unimolecular decomposition rate of ethylene is negligible. Three potential energy surfaces describing the production of vinyl naphthalene, acenaphthene, and acenaphthylene were calculated using quantum chemical methods and rate constants for the elementary steps on the surfaces were calculated using transition state theory. Naphthalene is not part of the reactions on the surfaces. Acenaphthylene is obtained only from acenaphthene. A kinetics scheme containing 27 elementary steps most of which were obtained from the potential energy surfaces was constructed and computer modeling was performed. An excellent agreement between the experimental yields of the four major products and the calculated yields was obtained.     

Experimental investigation of toluene + H --> benzyl + H2 at high temperatures

The reaction of toluene with hydrogen atoms yielding benzyl and molecular hydrogen, C(6)H(5)CH(3) + H --> C(6)H(5)CH(2) + H(2), was investigated using UV laser absorption of benzyl radicals at 266 nm in shock tube experiments. Test gas mixtures of toluene and ethyl iodide, an H-atom source, diluted in argon were heated in reflected shock waves to temperatures ranging from 1256 to 1667 K at total pressures around 1.7 bar. Measurement of laser absorption at 266 nm due to benzyl radicals allowed determination of the rate coefficient of the title reaction, reaction 1. A two-parameter best-fit Arrhenius expression for the rate determinations over the temperature range of these experiments is given by k(1)(T) = 1.33 x 10(15) exp(-14880 [cal/mol]/RT) [cm(3) mol(-1) s(-1)]. With the use of both the high-temperature shock tube measurements reported here and the rate coefficient determination of Ellis et al. (Ellis, C.; Scott, M. S.; Walker, R. W. Combust. Flame 2003, 132, 291) at 773 K the best-fit rate coefficient for reaction 1 can be described using a three-parameter Arrhenius expression by k(1)(T) = 6.47T (3.98) exp(-3384 [cal/mol]/RT) [cm(3) mol(-1) s(-1)].     

High-temperature measurements of the reactions of OH with toluene and acetone

The reaction of hydroxyl [OH] radicals with toluene [C6H5CH3] was studied at temperatures between 911 and 1389 K behind reflected shock waves at pressures of approximately 2.25 atm. OH radicals were generated by rapid thermal decomposition of shock-heated tert-butyl hydroperoxide [(CH3)3-CO-OH], and monitored by narrow-line width ring dye laser absorption of the well-characterized R1(5) line of the OH A-X (0,0) band near 306.7 nm. OH time histories were modeled by using a comprehensive toluene oxidation mechanism. Rate constants for the reaction of C6H5CH3 with OH were extracted by matching modeled and measured OH concentration time histories in the reflected shock region. Detailed error analyses yielded an uncertainty estimate of +/-30% at 1115 K for the rate coefficient of this reaction. The current high-temperature data were fit with the lower temperature measurements of Tully et al. [J. Phys. Chem. 1981, 85, 2262-2269] to the following two-parameter form, applicable over 570-1389 K: k3 = (1.62 x 10(13)) exp(-1394/T [K]) [cm3 mol(-1) s(-1)]. The reaction between OH radicals and acetone [CH3COCH3] was one of the secondary reactions encountered in the toluene + OH experiments. Direct high-temperature measurements of this reaction were carried out at temperatures ranging from 982 to 1300 K in reflected shock wave experiments at an average total pressure of 1.65 atm. Uncertainty limits were estimated to be +/-25% at 1159 K. A two-parameter fit of the current data yields the following rate expression: k6 = (2.95 x 10(13)) exp(-2297/T [K]) [cm3 mol(-1) s(-1)].     

Experimental and theoretical studies of the reaction of the OH radical with alkyl sulfides: 1. Direct observations of the formation of the OH-DMS adduct-pressure dependence of the forward rate of addition and development of a predictive expression at low temperature

A pulsed laser photolysis-pulsed laser-induced fluorescence (PLP-PLIF) system was employed to study the kinetics and mechanisms of reactions (1) OH + h6-DMS --> products and (2) OH + d6-DMS --> products. We report direct observations of the rate coefficients for the formation and dissociation of the h6-OHDMS and d6-OHDMS adducts over the pressure range 50-650 Torr and between 240 and 245 K, together with measurements of the oxygen dependence of the effective rate coefficients for reactions 1 and 2 under similar conditions. The effective rate coefficients increased as a function of O2 concentration reaching their limiting values in each case. The values of the adduct formation rate, obtained from the O2 dependencies, were in excellent agreement with values determined from direct observation of adduct equilibration in N2. OH regeneration is insignificant. The rate coefficients for the formation of the adduct isotopomers showed slight differences in their falloff behavior and do not approach the high-pressure limit in either case. The equilibrium constants obtained show no dependence on isotopomer and are in good agreement with previous work. A "second-law" analysis of the temperature dependence of the equilibrium constant gives an adduct bond strength (DeltaH degrees =-10.9 +/- 1.0 kcal mol(-1)), also in good agreement with previously reported values. Using the entropy calculated from the ab initio vibrational frequencies, we obtain a "third-law" value for the reaction enthalpy at 240 K, DeltaH(240K) degrees = -10.5 kcal mol(-1) in good agreement with the other approach. The rate coefficient for the reactions of the adducts with O2 was obtained from an analysis of the O2 dependence and was determined to be 6.3 +/- 1.2 x 10(-13) cm3 molecule(-1) s(-1), with no dependence on pressure or isotopomer. The pressure and temperature dependence for all of the elementary processes in the initial steps of the dimethylsulfide (DMS) oxidation mechanism have been characterized in the range 238-245 K, allowing the formulation of an expression which can be used to calculate the effective rate coefficient for reaction 1 at any pressure and oxygen concentration. The expression can calculate the effective rate coefficient for reaction 1 to +/- 40% over the range 220-260 K, with the largest errors at the extremes of this range. Gaussian 03 has been used to calculate the structure of the OH-DMS adduct and its deuterated isotopomer. We find similar bound structures for both isotopomers. The calculated enthalpies of formation of the adducts are lower than the experimentally determined values.