Kinetics of the hydrogen abstraction OH + alkane --> H2O + alkyl reaction class: an application of the reaction class transition state theory

This paper presents an application of the reaction class transition state theory (RC-TST) to predict thermal rate constants for hydrogen abstraction reactions of the type OH + alkane --> HOH + alkyl. We have derived all parameters for the RC-TST method for this reaction class from rate constants of 19 representative reactions, coupling with linear energy relationships (LERs), so that rate constants for any reaction in this class can be predicted from its reaction energy calculated at either the AM1 semiempirical or BH&HLYP/cc-pVDZ level of theory. The RC-TST/LER thermal rate constants for selected reactions are in good agreement with those available in the literature. Detailed analyses of the results show that the RC-TST/LER method is an efficient method for accurately estimating rate constants for a large number of reactions in this class. Analysis of the LERs leads to the discovery of the beta-carbon radical stabilization effect that stabilizes the transition state of any reaction in this class that yields products having one or more beta-carbons, and thus leads to the lower barrier for such a reaction.     

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.     

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.     

Proton transfer dynamics of the reaction H3O(+)(NH3,H2O)NH4+ studied using the crossed molecular beam technique

The proton transfer reaction of H3O+ and NH3 was studied using the crossed molecular beam technique at relative energies of 0.41, 0.81, and 1.27 eV. At all three energies, the center-of-mass flux distribution of the product ion NH4+ exhibits sharply asymmetry, and the maximum is close to the velocity and direction of the precursor ammonia beam. The reaction transforms almost all of the 1.69 eV exothermicity into internal excitation of the products at all three collision energies. At the lowest collision energy of 0.41 eV, nearly 77% of the total energy appears in NH4+ internal excitation. However, almost 100% of the incremental translational energy in the two higher-energy experiments appears in the product translational energy. Such an observation provides a classic example of the "induced repulsive energy release" mechanism that is expected to be operative on the highly skewed potential energy surfaces characteristic of the heavy+light-heavy mass combination. These results indicate that the proton transfer proceeds through a direct reaction mechanism; a Rice-Ramsperger-Kassel-Marcus theory calculation shows that the lifetime of the intermediate complex [NH3-H-H2O]+ is about 100 fs. Proton transfer occurs early on the reaction coordinate, when the incipient N-H bond is extended, and results in highly vibrationally excited NH4+ products, with excitation primarily in N-H stretching modes.     

Quantum calculations of the O(3P)+H2-->OH+H reaction

Quantum scattering calculations are reported for the O(3P)+H2(v=0,1) reaction using chemically accurate potential energy surfaces of 3A' and 3A" symmetry. We present state-to-state reaction cross sections and rate coefficients as well as thermal rate coefficients for the title reaction using accurate quantum calculations. Our calculations yield reaction cross sections that are in quantitative accord with results of recent crossed molecular beam experiments. Comparisons with results obtained using the J-shifting calculations show that the J-shifting approximation is quite reliable for this system. Thermal rate coefficients from the exact calculations and the J-shifting approximation agree remarkably well with experimental results. Our calculations also reproduce the markedly different OH(v'=0)/OH(v'=1) branching in O(3P)+H2(v=1) reaction, observed in experiments that use different O(3P) atom sources. In particular, we show that the branching ratio is a strong function of the kinetic energy of the O(3P) atom.     

Quantum dynamics of the H + O2 --> O + OH reaction on an accurate ab initio potential energy surface

We report exact time-dependent and time-independent quantum mechanical studies of the title reaction on an accurate ab initio potential energy surface of Xu et al. (J. Chem. Phys. 2005, 122, 24305). The J = 0 reaction probabilities for several reactant states show sharp resonance structures superimposed on relatively low backgrounds, and they are remarkably different from existing quantum results on an earlier potential energy surface (DMBE-IV). The new findings reported here suggest that our current understanding of this important reaction might require significant revision.     

Product spin-orbit state resolved dynamics of the H+H2O and H+D2O abstraction reactions

The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.     

Angular momentum conservation in the O + OH ↔ O2 + H reaction

We have studied the O + OH ↔ O₂ + H reaction on Varandas's DMBE IV potential using a variety of statistical methods, all involving the RRKM assumption for the HO₂* complex. Comparing our results using microcanonical variational transition‐state theory (μVT) with those using microcanonical/fixed‐J variational transition‐state theory (μVT‐J), we find that the effect of angular momentum conservation on the rate coefficient is imperceptible up to a temperature of about 700 K. Above 700 K angular momentum conservation increasingly reduces the rate coefficient, but only by approximately 21% even at 5000 K. Comparing our μVT‐J calculations with the quasi‐classical trajectory (QCT) results of Miller and Garrett [ 1 ], we confirm their conclusion that non‐RRKM dynamics of the HO₂* complex reduces the rate coefficient by about a factor of 2 independent of temperature. Our calculations of k^(c), the rate coefficient for HO₂* formation from O + OH, are in excellent agreement with the QCT results of Miller and Garrett. Although the differences are not large, we find k_CVT^(c) > k_μVT^(c) > k_μVT‐J^(c) > k_QCT^(c), where CVT stands for canonical variational transition‐state theory. © 1999 John Wiley & Sons, Inc. Int J Chem Kinet 31: 753–756, 1999     

A quantum mechanical approach to the kinetics of the hydrogen abstraction reaction H2O2 + •OH → HO2 + H2O

The kinetics of the hydrogen abstraction from H₂O₂ by ^?OH has been modeled with MP2/6‐31G*//MP2/6‐31G*, MP2‐SAC//MP2/6‐31G*, MP2/6‐31+G**//MP2/6‐31+G**, MP2‐SAC// MP2/6‐31+G**, MP4(SDTQ)/6‐311G**//MP2/6‐31G*, CCSD(T)/6‐31G*//CCSD(T)/6‐31G*, CCSD(T)/6‐31G**//CCSD(T)/6‐31G**, CCSD(T)/6‐311++G**//MP2/6‐31G* in the gas phase. MD simulations have been used to generate initial geometries for the stationary points along the potential energy surface for hydrogen abstraction from H₂O₂. The effective fragment potential (EFP) has been used to optimize the relevant structures in solution. Furthermore, the IEFPCM model has been used for the supermolecules generated via MD calculations. IEFPCM/MP2/6‐31G* and IEFPCM/CCSD(T)/6‐31G* calculations have also been performed for structures without explicit water molecules. Experimentally, the rate constant for hydrogen abstraction by ^?OH drops from 1.75 × 10^?12 cm³ molecule^?1 s^?1 in the gas phase to 4.48 × 10^?14 cm³ molecule^?1 s^?1 in solution. The same trend has been reproduced best with MP4 (SDTQ)/6‐311G**//MP2/6‐31G* in the gas phase (0.415 × 10^?12 cm³ molecule^?1 s^?1) and with EFP (UHF/6‐31G*) in solution (3.23 × 10^?14 cm³ molecule^?1 s^?1). © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 502–514, 2005     

水催化HBrO与OH自由基反应的微观机理

采用密度泛函理论,分别在B3LYP/6-311++g(d,p)和B3LYP/aug-cc-PVTZ理论水平下,系统研究了无水和水催化的OH自由基与HBrO反应,即HBrO+OH和HBrO+OH+H_2O 2个反应的微观反应机理,给出了所有可能发生的反应路径,并指出能量最低的反应通道.对于没有水参与的反应,由于OH自由基进攻HBrO方式不同,存在顺式方向和反式方向2种进攻方式的反应路径;当有一分子水参与反应时,考虑HBrO H_2O复合物与OH自由基的反应和HBrO与H_2O OH复合物2种反应情况,共发现4条不同的反应路径.这2种反应的所有路径均是在OH自由基提取氢之前以氢键复合物形式存在,反应过程均为无势垒加合过程,总反应为放热反应.水对目标反应起催化作用,有效地降低了反应的势垒,可以加快OH自由基和HBrO的消耗速度.     

Full-dimensional quantum dynamics study of the H + H2O and H + HOD exchange reactions

The initial state-selected time-dependent wave packet approach is employed to study the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H, HOH' + D exchange reactions with both OH bonds in the H(2)O reactant and OH(D) bond in the HOD reactant treated as reactive bonds. The total reaction probabilities for different partial waves, as well as the integral cross sections, which are the exact CC (coupled-channel) results, are first obtained in this study for the H(2)O(HOD) reactant initially in the ground rovibrational state. Because of the shallow C(3v) minimum along the reaction path, the reaction probabilities for the three reactions present several resonance peaks, with one dominant resonance peak just above the threshold. The cross sections for the H' + HOD → HOH' + D reaction are substantially smaller than those for the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H reactions, indicating that the H'/H exchange reactions are much more favored. In the CC calculations, the resonance peaks in the reaction probabilities diminish quickly with the increase in total angular momenta J, resulting in the existence of a clear step-like feature just above the threshold in the cross sections for the title reactions, which manifests the signature of shape resonances in these reactions. In the CS calculations, the resonance peaks on reaction probabilities persist in many partial waves, and thus the resonance structures can no longer survive the partial-wave summation and are washed out completely in the CS cross sections for the title reactions.     

State-to-state dynamics of the Cl + CH3OH --> HCl + CH2OH reaction

Molecular chlorine, methanol, and helium are co-expanded into a vacuum chamber using a custom designed "late-mixing" nozzle. The title reaction is initiated by photolysis of Cl2 at 355 nm, which generates monoenergetic Cl atoms that react with CH3OH at a collision energy of 1960 +/- 170 cm(-1) (0.24 +/- 0.02 eV). Rovibrational state distributions of the nascent HCl products are obtained via 2 + 1 resonance enhanced multiphoton ionization, center-of-mass scattering distributions are measured by the core-extraction technique, and the average internal energy of the CH3OH co-products is deduced by measuring the spatial anisotropy of the HCl products. The majority (84 +/- 7%) of the HCl reaction products are formed in HCl(v = 0) with an average rotational energy of [Erot] = 390 +/- 70 cm(-1). The remaining 16 +/- 7% are formed in HCl(v = 1) and have an average rotational energy of [Erot] = 190 +/- 30 cm(-1). The HCl(v = 1) products are primarily forward scattered, and they are formed in coincidence with CH2OH products that have little internal energy. In contrast, the HCl(v = 0) products are formed in coincidence with CH2OH products that have significant internal energy. These results indicate that two or more different mechanisms are responsible for the dynamics in the Cl + CH3OH reaction. We suggest that (1) the HCl(v = 1) products are formed primarily from collisions at high impact parameter via a stripping mechanism in which the CH2OH co-products act as spectators, and (2) the HCl(v = 0) products are formed from collisions over a wide range of impact parameters, resulting in both a stripping mechanism and a rebound mechanism in which the CH2OH co-products are active participants. In all cases, the reaction of fast Cl atoms with CH3OH is with the hydrogen atoms on the methyl group, not the hydrogen on the hydroxyl group.     

Thermal decomposition of CF3 and the reaction of CF2 + OH --> CF2O + H

The reflected shock tube technique with multipass absorption spectrometric detection (at a total path length of approximately 1.75 m) of OH-radicals at 308 nm has been used to study the dissociation of CF3-radicals [CF3 + Kr --> CF2 + F + Kr (a)] between 1,803 and 2,204 K at three pressures between approximately 230 and 680 Torr. The OH-radical concentration buildup resulted from the fast reaction F + H2O --> OH + HF (b). Hence, OH is a marker for F-atoms. To extract rate constants for reaction (a), the [OH] profiles were modeled with a chemical mechanism. The initial rise in [OH] was mostly sensitive to reactions (a) and (b), but the long time values were additionally affected by CF2 + OH --> CF2O + H (c). Over the experimental temperature range, rate constants for (a) and (c) were determined from the mechanistic fits to be kCF3+Kr = 4.61 x 10-9 exp(-30,020 K/T) and kCF2+OH = (1.6 +/- 0.6) x 10-10, both in units of cm3 molecule-1 s-1. Reaction (a), its reverse recombination reaction reaction (-a), and reaction (c) are also studied theoretically. Reactions (c) and (-a) are studied with direct CASPT2 variable reaction coordinate transition state theory. A master equation analysis for reaction (a) incorporating the ab initio determined reactive flux for reaction (-a) suggests that this reaction is close to but not quite in the low-pressure limit for the pressures studied experimentally. In contrast, reaction (c) is predicted to be in the high-pressure limit due to the high exothermicity of the products. A comparison with past and present experimental results demonstrates good agreement between the theoretical predictions and the present data for both (a) and (c).     

Variational transition‐state theory study of the atmospheric reaction OH + O3 → HO2 + O2

We report variational transition‐state theory calculations for the OH + O₃→ HO₂ + O₂ reaction based on the recently reported double many‐body expansion potential energy surface for ground‐state HO₄ [Chem Phys Lett 2000, 331, 474]. The barrier height of 1.884 kcal mol^?1 is comparable to the value of 1.77–2.0 kcal mol^?1 suggested by experimental measurements, both much smaller than the value of 2.16–5.11 kcal mol^?1 predicted by previous ab initio calculations. The calculated rate constant shows good agreement with available experimental results and a previous theoretical dynamics prediction, thus implying that the previous ab initio calculations will significantly underestimate the rate constant. Variational and tunneling effects are found to be negligible over the temperature range 100–2000 K. The O₁? O₂ bond is shown to be spectator like during the reactive process, which confirms a previous theoretical dynamics prediction. © 2007 Wiley Periodicals, Inc. 39: 148–153, 2007     

Quantum dynamics study of the dissociative photodetachment of HOCO-

Six-dimensional wave packet calculations are carried out to study the behavior of HOCO subsequent to the photodetachment of an electron from the negative anion, HOCO-. It is possible to form stable and/or long-lived HOCO complexes, as well as the dissociative products OH+CO and H+CO2. A variety of observables are determined: the electron kinetic energy (eKE) distributions associated with the OH+CO and H+CO2 channels, the correlated eKE and product translational energy distribution for the OH+CO channel, and product branching ratios. Most of our results are in good accord with the experimental results of Clements, Continetti, and Francisco [J. Chem. Phys. 117, 6478 (2002)], except that the calculated eKE distribution for the H+CO2 channel is noticeably colder than experiment. Reasons for this discrepancy are suggested.     

A flash photolysis resonance fluorescence investigation of the reaction OH + H2O2 → HO2 + H2O

The flash photolysis resonance fluorescence technique has been used to measure the rate constant for the reaction over the temperature range of 250–370 K. The present results are in excellent agreement with three very recent studies, and the combined data set can been used to derive the expression similar to that currently used in atmospheric modeling applications. A summary of our computer simulation of this reaction system is presented. The results of the computations indicate the absence of secondary reaction complications in the present work while revealing significant problems in the earlier (pre-1980) studies of the title reaction.     

Effects of a single water molecule on the OH + H2O2 reaction

The effect of a single water molecule on the reaction between H(2)O(2) and HO has been investigated by employing MP2 and CCSD(T) theoretical approaches in connection with the aug-cc-PVDZ, aug-cc-PVTZ, and aug-cc-PVQZ basis sets and extrapolation to an ∞ basis set. The reaction without water has two elementary reaction paths that differ from each other in the orientation of the hydrogen atom of the hydroxyl radical moiety. Our computed rate constant, at 298 K, is 1.56 × 10(-12) cm(3) molecule(-1) s(-1), in excellent agreement with the suggested value by the NASA/JPL evaluation. The influence of water vapor has been investigated by considering either that H(2)O(2) first forms a complex with water that reacts with hydroxyl radical or that H(2)O(2) reacts with a previously formed H(2)O·OH complex. With the addition of water, the reaction mechanism becomes much more complex, yielding four different reaction paths. Two pathways do not undergo the oxidation reaction but an exchange reaction where there is an interchange between H(2)O(2)·H(2)O and H(2)O·OH complexes. The other two pathways oxidize H(2)O(2), with a computed total rate constant of 4.09 × 10(-12) cm(3) molecule(-1) s(-1) at 298 K, 2.6 times the value of the rate constant of the unassisted reaction. However, the true effect of water vapor requires taking into account the concentration of the prereactive bimolecular complex, namely, H(2)O(2)·H(2)O. With this consideration, water can actually slow down the oxidation of H(2)O(2) by OH between 1840 and 20.5 times in the 240-425 K temperature range. This is an example that demonstrates how water could be a catalyst in an atmospheric reaction in the laboratory but is slow under atmospheric conditions.     

Kinetics study of the OH + alkene --> H2O + alkenyl reaction class

In this paper we report on the kinetics of hydrogen abstraction for the OH + alkene reaction class, using the reaction class transition state theory (RC-TST) combined with the linear energy relationship (LER) and the barrier height grouping (BHG) approaches. Parameters for the RC-TST were derived from theoretical calculations using a set of 15 reactions representing the hydrogen abstractions from the terminal and nonterminal carbon sites of the double bond of alkene compounds. Both the RC-TST/LER, where only reaction energy is needed at either density functional theory BH&HLYP or semiempirical AM1 levels, and RC-TST/BHG, where no additional information is required, are found to be promising methods for predicting rate constants for a large number of reactions in this reaction class. Detailed error analyses show that, when compared to explicit theoretical calculations, the averaged systematic errors in the calculated rate constants using both the RC-TST/LER and RC-TST/BHG methods are less than 25% in the temperature range 300-3000 K. The estimated rate constants using these approaches are in good agreement with available data in the literature.     

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.