Reflected shock tube studies of high-temperature rate constants for OH + CH4 --> CH3 + H2O and CH3 + NO2 --> CH3O + NO

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm has been used to study the reactions OH + CH(4) --> CH(3) + H(2)O and CH(3) + NO(2) --> CH(3)O + NO. Over the temperature range 840-2025 K, the rate constants for the first reaction can be represented by the Arrhenius expression k = (9.52 +/- 1.62) x 10(-11) exp[(-4134 +/- 222 K)/T] cm(3) molecule(-1) s(-1). Since this reaction is important in both combustion and atmospheric chemistry, there have been many prior investigations with a variety of techniques. The present results extend the temperature range by 500 K and have been combined with the most accurate earlier studies to derive an evaluation over the extended temperature range 195-2025 K. A three-parameter expression describes the rate behavior over this temperature range, k = (1.66 x 10(-18))T(2.182) exp[(-1231 K)/T] cm(3) molecule(-1) s(-1). Previous theoretical studies are discussed, and the present evaluation is compared to earlier theoretical estimates. Since CH(3) radicals are a product of the reaction and could cause secondary perturbations in rate constant determinations, the second reaction was studied by OH radical production from the fast reactions CH(3)O --> CH(2)O + H and H + NO(2) --> OH + NO. The measured rate constant is 2.26 x 10(-11) cm(3) molecule(-1) s(-1) and is not dependent on temperature from 233 to 1700 K within experimental error.     

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.     

Reflected shock tube studies of high-temperature rate constants for CH3 + O2, H2CO + O2, and OH + O2

The reflected shock tube technique with multipass absorption spectrometric detection of OH-radicals at 308 nm, corresponding to a total path length of approximately 2.8 m, has been used to study the reaction CH3 + O2 CH2O + OH. Experiments were performed between 1303 and 2272 K, using ppm quantities of CH3I (methyl source) and 5-10% O2, diluted with Kr as the bath gas at test pressures less than 1 atm. We have also reanalyzed our earlier ARAS measurements for the atomic channel (CH3 + O2 --> CH3O + O) and have compared both these results with other earlier studies to derive a rate expression of the Arrhenius form. The derived expressions, in units of cm3 molecule(-1) s(-1), are k = 3.11 x 10(-13) exp(-4953 K/T) over the T-range 1237-2430 K, for the OH-channel, and k = 1.253 x 10(-11) exp(-14241 K/T) over the T-range 1250-2430 K, for the O-atom channel. Since CH2O is a major product in both reactions, reliable rates for the reaction CH2O + O2 --> HCO + HO2 could be derived from [OH]t and [O]t experiments over the T-range 1587-2109 K. The combined linear least-squares fit result, k = 1.34 x 10(-8) exp(-26883 K/T) cm3 molecule(-1) s(-1), and a recent VTST calculation clearly overlap within the uncertainties in both studies. Finally, a high sensitivity for the reaction OH + O2 --> HO2 + O was noted at high temperature in the O-atom data set simulations. The values for this obtained by fitting the O-atom data sets at later times (approximately 1.2 ms) again follow the Arrhenius form, k = 2.56 x 10(-10) exp(-24145 K/T) cm3 molecule(-1) s(-1), over the T-range, 1950-2100 K.     

Reflected shock tube and theoretical studies of high-temperature rate constants for OH+CF3H<-->CF3+H2O and CF3+OH-->products

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm, using either 36 or 60 optical passes corresponding to total path lengths of 3.25 or 5.25 m, respectively, has been used to study the bimolecular reactions, OH+CF3H-->CF3+H2O (1) and CF3+H2O-->OH+CF3H (-1), between 995 and 1663 K. During the course of the study, estimates of rate constants for CF3+OH-->products (2) could also be determined. Experiments on reaction -1 were transformed through equilibrium constants to k1, giving the Arrhenius expression k1=(9.7+/-2.1)x10(-12) exp(-4398+/-275K/T) cm3 molecule(-1) s(-1). Over the temperature range, 1318-1663 K, the results for reaction 2 were constant at k2=(1.5+/-0.4)x10(-11) cm3 molecule(-1) s(-1). Reactions 1 and -1 were also studied with variational transition state theory (VTST) employing QCISD(T) properties for the transition state. These a priori VTST predictions were in good agreement with the present experimental results but were too low at the lower temperatures of earlier experiments, suggesting that either the barrier height was overestimated by about 1.3 kcal/mol or that the effect of tunneling was greatly underestimated. The present experimental results have been combined with the most accurate earlier studies to derive an evaluation over the extended temperature range of 252-1663 K. The three parameter expression k1=2.08x10(-17) T1.5513 exp(-1848 K/T) cm3 molecule(-1) s(-1) describes the rate behavior over this temperature range. Alternatively, the expression k1,th=1.78x10(-23) T3.406 exp(-837 K/T) cm3 molecule(-1) s(-1) obtained from empirically adjusted VTST calculations over the 250-2250 K range agrees with the experimental evaluation to within a factor of 1.6. Reaction 2 was also studied with direct CASPT2 variable reaction coordinate transition state theory. The resulting predictions for the capture rate are found to be in good agreement with the mean of the experimental results and can be represented by the expression k2,th=2.42x10(-11) T-0.0650 exp(134 K/T) cm3 molecule(-1) s(-1) over the 200-2500 K temperature range. The products of this reaction are predicted to be CF2O+HF.     

Quantum mechanical rate constants for H + O2 <--> O + OH and H + O2 --> HO2 reactions

Canonical rate constants for both the forward and reverse H + O(2) <--> O + OH reactions were calculated using a quantum wave packet-based statistical model on the DMBE IV potential energy surface of Varandas and co-workers. For these bimolecular reactions, the results show reasonably good agreement with available experimental and theoretical data up to 1500 K. In addition, the capture rate for the H + O(2) --> HO(2) addition reaction at the high-pressure limit was obtained on the same potential using a time-independent quantum capture method. Excellent agreement with experimental and quasi-classical trajectory results was obtained except for at very low temperatures, where a reaction threshold was found and attributed to the centrifugal barrier of the orbital motion.     

A new ab initio potential-energy surface of HO2(X2A") and quantum studies of HO2 vibrational spectrum and rate constants for the H + O2 <--> O + OH reactions

A new global potential-energy surface for the ground electronic state of HO(2)(X(2)A(")) has been developed by three-dimensional cubic spline interpolation of more than 15 000 ab initio points, which were calculated at the multireference configuration-interaction level with Davidson correction using the augmented correlation-consistent polarized valence quadruple zeta basis set. Low-lying vibrational states were obtained in this new potential using the Lanczos method and assigned. The calculated vibrational frequencies are in much better agreement with the available experimental band origins than those obtained from a previous potential. In addition, rate constants for the H+O(2) <--> O + OH reactions were obtained using a wave-packet-based statistical model. Reasonably good agreement with experimental data was obtained. These results demonstrate the accuracy of the potential.     

High-temperature rate constants for CH3OH + Kr --> products, OH + CH3OH --> products, OH + (CH3)(2)CO --> CH2COCH3 + H2O, and OH + CH3 --> CH) + H2O

The reflected shock tube technique with multipass 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 dissociation of methanol between 1591 and 2865 K. Rate constants for two product channels [CH3OH + Kr --> CH3 + OH + Kr (1) and CH3OH + Kr --> 1CH2 + H2O + Kr (2)] were determined. During the course of the study, it was necessary to determine several other rate constants that contributed to the profile fits. These include OH + CH3OH --> products, OH + (CH3)2CO --> CH2COCH3 + H2O, and OH + CH3 --> 1,3CH2 + H2O. The derived expressions, in units of cm(3) molecule(-1) s(-1), are k(1) = 9.33 x 10(-9) exp(-30857 K/T) for 1591-2287 K, k(2) = 3.27 x 10(-10) exp(-25946 K/T) for 1734-2287 K, kOH+CH3OH = 2.96 x 10-16T1.4434 exp(-57 K/T) for 210-1710 K, k(OH+(CH3)(2)CO) = (7.3 +/- 0.7) x 10(-12) for 1178-1299 K and k(OH+CH3) = (1.3 +/- 0.2) x 10(-11) for 1000-1200 K. With these values along with other well-established rate constants, a mechanism was used to obtain profile fits that agreed with experiment to within <+/-10%. The values obtained for reactions 1 and 2 are compared with earlier determinations and also with new theoretical calculations that are presented in the preceding article in this issue. These new calculations are in good agreement with the present data for both (1) and (2) and also for OH + CH3 --> products.     

Ab initio investigation on the reaction path and rate for the gas-phase reaction of HO + H2O <--> H2O + OH

This article describes an ab initio investigation on the potential surfaces for one of the simplest hydrogen atom abstraction reactions, that is, HO + H2O <--> H2O + OH. In accord with the findings in the previously reported theoretical investigations, two types of the hydrogen-bonding complexes [HOH--OH] and [H2O--HO] were located on the potential energy surface. The water molecule acts as a hydrogen donor in the [HOH--OH] complex, while the OH radical acts as a hydrogen donor in the [H2O--HO] complex. The energy evaluations at the MP2(FC) basis set limit, as well as those through the CBS-APNO procedure, have provided estimates for enthalpies of association for these complexes at 298 K as -2.1 approximately -2.3 and -4.1 approximately -4.3 kcal/mol, respectively. The IRC calculations have suggested that the [H2O--HO] complex should be located along the reaction coordinate for the hydrogen abstraction. Our best estimate for the classical barrier height for the hydrogen abstraction is 7.8 kcal/mol, which was obtained from the CBS-APNO energy evaluations. After fitting the CBS-APNO potential energy curve to a symmetrical Eckart function, the rate constants were calculated by using the transition state theory including the tunneling correction. Our estimates for the Arrhenius parameters in the temperature region from 300 to 420 K show quite reasonable agreement with the experimentally derived values.     

Including anharmonicity in the calculation of rate constants. II. The OH+H2-->H2O+H reaction

A recently developed method for calculating anharmonic vibrational energy levels at nonstationary points along a reaction path that is based on second-order perturbation theory in curvilinear coordinates is combined with variational transition state theory with semiclassical multidimensional tunneling approximations to calculate thermal rate constants for the title reaction. Two different potential energy surfaces were employed for these calculations, an improved version of the author's surface 5 and the WSLFH surface of Wu et al. [J. Chem. Phys. 113, 3150 (2000)]. We present detailed comparisons of rate constants computed for the two surfaces with and without anharmonicity and with various approximations for incorporating tunneling along the reaction path. The results for this system are quite sensitive to the surface employed, the choice of coordinates (curvilinear versus rectilinear), and the inclusion of anharmonicity. A comparison with experiment provides information on the accuracy of these surfaces.     

Quasiclassical trajectory simulations of OH(v) + NO2 --> HONO2* --> OH(v') + NO2: capture and vibrational deactivation rate constants

Quasiclassical trajectory calculations are used to investigate the dynamics of the OH(v) + NO(2) --> HONO(2) --> OH(v') + NO(2) recombination/dissociation reaction on an analytic potential energy surface (PES) that gives good agreement with the known structure and vibrational frequencies of nitric acid. The calculated recombination rate constants depend only weakly on temperature and on the initial vibrational energy level of OH(v). The magnitude of the recombination rate constant is sensitive to the potential function describing the newly formed bond and to the switching functions in the PES that attenuate inter-mode interactions at long range. The lifetime of the nascent excited HONO(2) depends strongly not only on its internal energy but also on the identity of the initial state, in disagreement with statistical theory. This disagreement is probably due to the effects of slow intramolecular vibrational energy redistribution (IVR) from the initially excited OH stretching mode. The vibrational energy distribution of product OH(v') radicals is different from statistical distributions, a result consistent with the effects of slow IVR. Nonetheless, the trajectory results predict that vibrational deactivation of OH(v) via the HONO(2) transient complex is approximately 90% efficient, almost independent of initial OH(v) vibrational level, in qualitative agreement with recent experiments. Tests are also carried out using the HONO(2) PES, but assuming the weaker O-O bond strength found in HOONO (peroxynitrous acid). In this case, the predicted vibrational deactivation efficiencies are significantly lower and depend strongly on the initial vibrational state of OH(v), in disagreement with experiments. This disagreement suggests that the actual HOONO PES may contain more inter-mode coupling than found in the present model PES, which is based on HONO(2). For nitric acid, the measured vibrational deactivation rate constant is a useful proxy for the recombination rate, but IVR randomization of energy is not complete, suggesting that the efficacy of the proxy method must be evaluated on a case-by-case basis.     

Kinetics and rate constants of the reaction CH2OH+O2-->CH2O+HO2 in the temperature range of 236-600 K

The kinetics and absolute rate constants of the gas-phase reaction of the hydroxymethyl radical (CH2OH) with molecular oxygen have been studied using laser photolysis/near-IR absorption spectroscopy. The reaction was tracked by monitoring the time-dependent changes in the production of the hydroperoxy radical (HO2) concentration. For sensitive detection of HO2, two-tone frequency modulation absorption spectroscopy was used in combination with a Herriott-type optical multipass absorption cell. Rate constants were determined as a function of temperature (236 K相似文献   

A shock tube study of the OH + OH → H2O + O reaction

The rate coefficient for the reaction has been determined in mixtures of nitric acid (HNO₃) and argon in incident shock wave experiments. Quantitative OH time-histories were obtained by cw narrow-linewidth uv laser absorption of the R₁(5) line of the A² σ⁺ ← X² Π_i (0,0) transition at 32606.56 cm^?1 (vacuum). The experiments were conducted over the temperature range 1050–2380 K and the pressure range 0.18–0.60 atm. The second-order rate coefficient was determined to be with overall uncertainties of +11%, ?16% at high temperatures and +25%, ?22% at low temperatures. By incorporating data from previous investigations in the temperature range 298–578 K, the following expression is determined for the temperature range 298–2380 K © 1994 John Wiley & Sons, Inc.     

Planar study of H2O+Cl<-->HO+HCl reactions

The first four dimensional (4D) quantum scattering calculations on the tetra-atomic H2O+Cl<-->HO+HCl reactions are reported. With respect to a full (6D) treatment, only the planar constraint and a fixed length for the HO spectator bond are imposed. This work explicitly accounts for the bending and local HO stretching vibrations in H2O, for the vibration of HCl and for the in-plane rotation of the H2O, HO and HCl molecules. The calculations are performed with the potential energy surface of Clary et al. and use a Born-Oppenheimer type separation between the motions of the light and the heavy nuclei. State-to-state cross sections are reported for a collision energy range 0-1.8 eV measured with respect to H2O+Cl. For the H2O+Cl reaction, present results agree with previous (3D) non planar calculations and confirm that excitation of the H2O stretching promotes more reactivity than excitation of the bending. New results are related to the rotation of the H2O molecule: the cross sections are maximal for planar rotational states corresponding to 10相似文献   

Bayesian uncertainty quantification of recent shock tube determinations of the rate coefficient of reaction H + O2→ OH + O

We analyze the ignition delay in hydrogen–oxygen combustion and the important chain ‐branching reaction H + O₂→ OH + O that occurs behind the shock waves in shock tube experiments. We apply a stochastic Bayesian approach to quantify uncertainties in the theoretical model and experimental data. The approach involves a statistical inverse problem, which has four “components” as input information: (a) model, (b) prior joint probability density function (PDF) of the uncertain parameters, (c) experimental data, and (d) uncertainties in the scenario parameters. The solution of this statistical inverse problem is a posterior joint PDF of the uncertain parameters from which we can easily extract statistical information. We first perform a parametric study to investigate how the level of the total uncertainty (which we define as the sum of model uncertainty and experimental uncertainty) affects the uncertainty in the rate coefficient k₁ of the reaction H + O₂→ OH + O, which is “most likely” expressed by k₁=1.73×10²³T^?2.5exp(?11550/T) cm³ mol^?1 s^?1 over the experimental temperature range 1100–1472 K. We also introduce the idea of “irreducible” uncertainty when considering other parameters in the system. After statistically calibrating the parameters modeling the rate coefficient k₁, we predict its 95% confidence interval (CI) for different temperature regimes and compare the CI against the values of k₁ obtained deterministically. Our results show that a small uncertainty in gas temperature (±5 K) introduces appreciable uncertainty in k₁. © 2012 Wiley Periodicals, Inc. Int J Chem Kinet 44: 586–597, 2012     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

Differential and integral cross sections for the H + O2 --> OH + O combustion reaction

We present accurate differential and integral cross sections for the H + O2 --> OH + O reaction obtained on a newly developed ab initio potential energy surface using time-independent and time-dependent quantum mechanical methods. The product angular distributions near the reaction threshold show pronounced forward and backward peaks, reflecting the complex-forming mechanism. However, the asymmetry of these peaks suggests certain nonstatistical behaviors, presumably due to some relatively short-lived resonances. The integral cross section increases monotonically with the collision energy above a reaction threshold.     

Direct dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2

We present a direct ab initio dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2, which is predicted to have four possible reaction channels caused by different attacking orientations of HO2 radical to CH2O. The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of the four reaction channels are calculated at the B3LYP/cc-pVTZ level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of some single-point multilevel energy calculations (HL). The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that, in the whole temperature range, the more favorable reaction channels are Channels 1 and 3. The total ICVT/SCT rate constants of the four channels at the HL//B3LYP/cc-pVTZ level of theory are in good agreement with the available experiment data over the measured temperature ranges, and the corresponding three-parameter expression is k(ICVT/SCT) = 3.13 x 10(-20) T(2.70) exp(-11.52/RT) cm3 mole(-1) s(-1) in the temperature range of 250-3000 K. Additionally, the flexibility of the dihedral angle of H2O2 is also discussed to explain the different experimental values.