A study of the hydrogen abstraction reactions of C2H radical with CH3CN, C2H5CN, and C3H7CN by dual-level generalized transition state theory

The hydrogen abstraction reactions C2H + CH3CN --> products (R1), C2H + CH3CH2CN --> products (R2), and C2H + CH3CH2CH2CN --> products (R3) have been investigated by dual-level generalized transition state theory. Optimized geometries and frequencies of all the stationary points and extra points along the minimum-energy path (MEP) are performed at the BH&H-LYP and MP2 methods with the 6-311G(d, p) basis set, and the energy profiles are further refined at the MC-QCISD level of theory. The rate constants are evaluated using canonical variational transition state theory (CVT) with a small-curvature tunneling correction (SCT) over a wide temperature range 104-2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values. Our calculations show that for reaction R2, the alpha-hydrogen abstraction channel and beta-hydrogen abstraction channel are competitive over the whole temperature range. For reaction R3, the gamma-hydrogen abstraction channel is preferred at lower temperatures, while the contribution of beta-hydrogen abstraction will become more significant with a temperature increase. The branching ratio to the alpha-hydrogen abstraction channel is found negligible over the whole temperature range.     

Ab initio direct dynamics studies on the reactions of chlorine atom with CH3-nFnCH2OH (n = 1-3)

The hydrogen abstraction reactions of Cl atom with a series of fluorinated alcohols, i.e., CH(3-n)F(n)CH(2)OH + Cl (n = 1-3) (R1-R3) have been studied systematically by ab initio direct dynamics method and the canonical variational transition state theory (CVT). The potential energy surface information is calculated at the MP2/6-311G(d,p) level. Energies along the minimum energy paths are improved by a series of single-point calculations at the higher modified GAUSSIAN-2 (G2M) level of theory. Theoretical analysis shows that three kinds of hydrogen atoms can be abstracted from the reactants CH(2)FCH(2)OH and CHF(2)CH(2)OH, and for CF(3)CH(2)OH, two possible pathways are found. The rate constants for each reaction channel are evaluated by CVT with the small-curvature tunneling correction (SCT) over a wide range of temperature from 200 to 2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values for the reactions CHF(2)CH(2)OH + Cl and CF(3)CH(2)OH + Cl. However, for the reaction CH(2)FCH(2)OH + Cl, there is negative temperature dependence below 500 K, which is different from the experimental fitted. It is shown that in the low temperature ranges, the three reactions all proceed predominantly via H-abstraction from the methylene positions, and with the increase of the temperature the H-abstraction channels from the fluorinated-methyl positions should be taken into account, while the H-abstraction channels from the hydroxyl groups are negligible over the whole temperature ranges. Also, the reactivity decreases substantially with fluorine substitution at the methyl position of alcohol.     

Effects of multidimensional tunneling in the kinetics of hydrogen abstraction reactions of O (3P) with CH3OCHO

Quantum tunneling paths are important in reactions when there is a significant component of hydrogenic motion along the potential energy surface. In this study, variational transition state with multidimensional tunneling corrections are employed in the calculations of the thermal rate constants for hydrogen abstraction from the cis‐CH₃OCHO by O (³P) giving CH₃OCO + OH (R1) and CH₂OCHO + OH (R2). The structures and electronic energies are computed with the M06‐2X method. Benchmark calculations with the CBS_D–T approach give an enthalpy of reaction at 0 K for R1 (−2.8 kcal/mol) and R2 (−2.5 kcal/mol) which are in good agreement with the experiment, i.e. −2.61 and −1.81 kcal/mol. At the low and intermediate values of temperatures, small‐ and large‐curvature tunneling dominate the kinetics of R1, which is the dominant path over the range of temperature from 250 to 1200 K. This study shows the importance of multidimensional tunneling corrections for both R1 and R2, for which the total rate constant at 298 K calculated with the CVT/μOMT method is 8.2 × 10⁻¹⁵ cm³ molecule⁻¹ s⁻¹ which agrees well with experiment value of 9.3 × 10⁻¹⁵ cm³ molecule⁻¹ s⁻¹ (Mori, Bull. Inst. Chem. Res. 1981, 59, 116). © 2018 Wiley Periodicals, Inc.     

CH3CH2+O(3P)反应机理的理论研究

The reaction for CH3CH2+O(3P) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the QCISD(T)/6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2O＋CH3, CH3CHO+H and CH2CH2+OH in the reaction. For the products CH2O＋CH3 and CH3CHO+H, the major production channels are A1: (R)→IM1→TS3→(A) and B1: (R)→IM1→TS4→(B), respectively. The majority of the products CH2CH2+OH are formed via the direct abstraction channels C1 and C2: (R)→TS1(TS2)→(C). In addition, the results suggest that the barrier heights to form the CO reaction channels are very high, so the CO is not a major product in the reaction.     

Deuterium enrichment of interstellar methanol explained by atom tunneling

We calculate, down to low temperature and for different isotopes, the reaction rate constants for the hydrogen abstraction reaction H + H(3)COH → H(2) + CH(2)OH/CH(3)O. These explain the known abundances of deuterated forms of methanol in interstellar clouds, where CH(2)DOH can be almost as abundant as CH(3)OH. For abstraction from both the C- and the O-end of methanol, the barrier-crossing motion involves the movement of light hydrogen atoms. Consequently, tunneling plays a dominant role already at relatively high temperature. Our implementation of harmonic quantum transition state theory with on the fly calculation of forces and energies accounts for these tunneling effects. The results are in good agreement with previous semiclassical and quantum dynamics calculations (down to 200 K) and experimental studies (down to 295 K). Here we extend the rate calculations down to lower temperature: 30 K for abstraction from the C-end of methanol and 80 K for abstraction from the OH-group. At all temperatures, abstraction from the C-end is preferred over abstraction from the O-end, more strongly so at lower temperature. Furthermore, the tunneling behavior strongly affects the kinetic isotope effects (KIEs). D + H(3)COH → HD + CH(2)OH has a lower vibrationally adiabatic barrier than H + H(3)COH → H(2) + CH(2)OH, giving rise to an inverse KIE (k(H)/k(D) < 1) at high temperature, in accordance with previous experiments and calculations. However, since tunneling is more facile for the light H atom, abstraction by H is favored over abstraction by D below ~135 K, with a KIE k(H)/k(D) of 11.2 at 30 K. The H + D(3)COD → HD + CD(2)OD reaction is calculated to be much slower than the D + H(3)COH → HD + CH(2)OH, in agreement with low-temperature solid-state experiments, which suggests the preference for H (as opposed to D) abstraction from the C-end of methanol to be the mechanism by which interstellar methanol is deuterium-enriched.     

Theoretical studies on the unimolecular decomposition of ethylene glycol

The unimolecular decomposition processes of ethylene glycol have been investigated with the QCISD(T) method with geometries optimized at the B3LYP/6-311++G(d,p) level. Among the decomposition channels identified, the H(2)O-elimination channels have the lowest barriers, and the C-C bond dissociation is the lowest-energy dissociation channel among the barrierless reactions (the direct bond cleavage reactions). The temperature and pressure dependent rate constant calculations show that the H(2)O-elimination reactions are predominant at low temperature, whereas at high temperature, the direct C-C bond dissociation reaction is dominant. At 1 atm, in the temperature range 500-2000 K, the calculated rate constant is expressed to be 7.63 × 10(47)T(-10.38) exp(-42262/T) for the channel CH(2)OHCH(2)OH → CH(2)CHOH + H(2)O, and 2.48 × 10(51)T(-11.58) exp(-43593/T) for the channel CH(2)OHCH(2)OH → CH(3)CHO + H(2)O, whereas for the direct bond dissociation reaction CH(2)OHCH(2)OH → CH(2)OH + CH(2)OH the rate constant expression is 1.04 × 10(71)T(-16.16) exp(-52414/T).     

Theoretical study on the OH + CH3NHCOOCH3 reaction

The multiple-channel reactions OH + CH3NHC(O)OCH3 --> products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6-311+G(d,p) level, and energetic information is further refined by the BMC-CCSD (single-point) method. The rate constants for every reaction channels, R1, R2, R3, and R4, are calculated by canonical variational transition state theory with small-curvature tunneling correction over the temperature range 200-1000 K. The total rate constants are in good agreement with the available experimental data and the two-parameter expression k(T) = 3.95 x 10(-12) exp(15.41/T) cm3 molecule(-1) s(-1) over the temperature range 200-1000 K is given. Our calculations indicate that hydrogen abstraction channels R1 and R2 are the major channels due to the smaller barrier height among four channels considered, and the other two channels to yield CH3NC(O)OCH3 + H2O and CH3NHC(O)(OH)OCH3 + H2O are minor channels over the whole temperature range.     

Branching ratios in the hydroxyl reaction with propene

The branching ratios for the reactions of attachment of hydroxyl radical to propene and hydrogen-atom abstraction were measured at 298 K over the buffer gas pressure range 60-400 Torr (N(2)) using a subatmospheric pressure turbulent flow reactor coupled with a chemical ionization quadrupole mass spectrometer. Isotopically enriched water H(2)(18)O was used to produce (18)O-labeled hydroxyl radicals in reaction with fluorine atoms. The β-hydroxypropyl radicals formed in the attachment reactions 1a and 1b , OH + C(3)H(6) → CH(2)(OH)C(?)HCH(3) (eq 1a ) and OH + C(3)H(6) → C(?)H(2)CH(OH)CH(3) (eq 1b ), were converted to formaldehyde and acetaldehyde in a sequence of secondary reactions in O(2)- and NO-containing environment. The (18)O-labeling propagates to the final products, allowing determination of the branching ratio for the attachment channels of reaction 1. The measured branching ratio for attachment is β(1b) = k(1b)/(k(1a) + k(1b)) = 0.51 ± 0.03, independent of pressure over the 60-400 Torr pressure range. An upper limit on the hydrogen-abstraction channel, OH + C(3)H(6) → H(2)O + C(3)H(5) (eq 1c ), was determined by measuring the water yield in reactions of OH and OD radicals (produced via H(D) + NO(2) → OH(OD) + NO reactions) with C(3)H(6) as k(1c)/(k(1a) + k(1b) + k(1c)) < 0.05 (at 298 K, 200 Torr N(2)).     

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.     

Thermal decomposition mechanisms of the methoxyphenols: formation of phenol, cyclopentadienone, vinylacetylene, and acetylene

The pyrolyses of the guaiacols or methoxyphenols (o-, m-, and p-HOC(6)H(4)OCH(3)) have been studied using a heated SiC microtubular (μ-tubular) reactor. The decomposition products are detected by both photoionization time-of-flight mass spectroscopy (PIMS) and matrix isolation infrared spectroscopy (IR). Gas exiting the heated SiC μ-tubular reactor is subject to a free expansion after a residence time of approximately 50-100 μs. The PIMS reveals that, for all three guaiacols, the initial decomposition step is loss of methyl radical: HOC(6)H(4)OCH(3) → HOC(6)H(4)O + CH(3). Decarbonylation of the HOC(6)H(4)O radical produces the hydroxycyclopentadienyl radical, C(5)H(4)OH. As the temperature of the μ-tubular reactor is raised to 1275 K, the C(5)H(4)OH radical loses a H atom to produce cyclopentadienone, C(5)H(4)═O. Loss of CO from cyclopentadienone leads to the final products, acetylene and vinylacetylene: C(5)H(4)═O → [CO + 2 HC≡CH] or [CO + HC≡C-CH═CH(2)]. The formation of C(5)H(4)═O, HCCH, and CH(2)CHCCH is confirmed with IR spectroscopy. In separate studies of the (1 + 1) resonance-enhanced multiphoton ionization (REMPI) spectra, we observe the presence of C(6)H(5)OH in the molecular beam: C(6)H(5)OH + λ(275.1?nm) → [C(6)H(5)OH ?] + λ(275.1nm) → C(6)H(5)OH(+). From the REMPI and PIMS signals and previous work on methoxybenzene, we suggest that phenol results from a radical/radical reaction: CH(3) + C(5)H(4)OH → [CH(3)-C(5)H(4)OH]* → C(6)H(5)OH + 2H.     

13C, 18O, and D fractionation effects in the reactions of CH3OH isotopologues with Cl and OH radicals

A relative rate experiment is carried out for six isotopologues of methanol and their reactions with OH and Cl radicals. The reaction rates of CH2DOH, CHD2OH, CD3OH, (13)CH3OH, and CH3(18)OH with Cl and OH radicals are measured by long-path FTIR spectroscopy relative to CH3OH at 298 +/- 2 K and 1013 +/- 10 mbar. The OH source in the reaction chamber is photolysis of ozone to produce O((1)D) in the presence of a large excess of molecular hydrogen: O((1)D) + H2 --> OH + H. Cl is produced by the photolysis of Cl2. The FTIR spectra are fitted using a nonlinear least-squares spectral fitting method with measured high-resolution infrared spectra as references. The relative reaction rates defined as alpha = k(light)/k(heavy) are determined to be: k(OH + CH3OH)/k(OH + (13)CH3OH) = 1.031 +/- 0.020, k(OH + CH3OH)/k(OH + CH3(18)OH) = 1.017 +/- 0.012, k(OH + CH3OH)/k(OH + CH2DOH) = 1.119 +/- 0.045, k(OH + CH3OH)/k(OH + CHD2OH) = 1.326 +/- 0.021 and k(OH + CH3OH)/k(OH + CD3OH) = 2.566 +/- 0.042, k(Cl + CH3OH)/k(Cl + (13)CH3OH) = 1.055 +/- 0.016, k(Cl + CH3OH)/k(Cl + CH3(18)OH) = 1.025 +/- 0.022, k(Cl + CH3OH)/k(Cl + CH2DOH) = 1.162 +/- 0.022 and k(Cl + CH3OH)/k(Cl + CHD2OH) = 1.536 +/- 0.060, and k(Cl + CH3OH)/k(Cl + CD3OH) = 3.011 +/- 0.059. The errors represent 2sigma from the statistical analyses and do not include possible systematic errors. Ground-state potential energy hypersurfaces of the reactions were investigated in quantum chemistry calculations at the CCSD(T) level of theory with an extrapolated basis set. The (2)H, (13)C, and (18)O kinetic isotope effects of the OH and Cl reactions with CH3OH were further investigated using canonical variational transition state theory with small curvature tunneling and compared to experimental measurements as well as to those observed in CH4 and several other substituted methane species.     

Theoretical Study on Mechanism of Reaction of OH with HO2NO2

The reaction of HO2NO2 (peroxynitric acid, PNA) with OH was studied by the hybrid density functional B3LYP and CBS-QB3 methods. Based on the calculated potential energy surface, five reaction channels, H2O+NO2+O2, HOOH+NO3, NO2+HO3H, HO2+HONO2 and HO2+HOONO, were examined in detail. The major reaction channel is PNA+OH→M1→TS1→H2O+NO2+O2. Taking a pre-equilibrium approximation and using the CBS-QB3 energies, the theoretical rate constant of this channel was calculated as 1.13×10-12 cm3/(molecule s) at 300 K, in agreement with the experimental result. Comparison between reactions of HOONO2+OH and HONO2+OH was carried out. For HOR+OH reactions, the total rate constants increase from R=NO2 to R=ONO2, which is consistent with experimental measurements.     

Low-energy paths for the unimolecular decomposition of CH3OH: a G2M/statistical theory study

The potential energy surface (PES) of the CH3OH system has been characterized by ab initio molecular orbital theory calculations at the G2M level of theory. The mechanisms for the decomposition of CH3OH and the related bimolecular reactions, CH3 + OH and 1CH2 + H2O, have been elucidated. The rate constants for these processes have been calculated using variational RRKM theory and compared with available experimental data. The total decomposition rate constants of CH3OH at the high- and low-pressure limits can be represented by k infinity = 1.56 x 10(16) exp(-44,310/T) s-1 and kAr0 = 1.60 x 10(36) T-12.2 exp(-48,140/T) cm3 molecule-1 s-1, respectively, covering the temperature range 1000-3000 K, in reasonable agreement with the experimental values. Our results indicate that the product branching ratios are strongly pressure dependent, with the production of CH3 + OH and 1CH2 + H2O dominant under high (P > 10(3) Torr) and low (P < 1 atm) pressures, respectively. For the bimolecular reaction of CH3 and OH, the total rate constant and the yields of 1CH2 + H2O and H2 + HCOH at lower pressures (P < 5 Torr) could be reasonably accounted for by the theory. For the reaction of 1CH2 with H2O, both the yield of CH3 + OH and the total rate constant could also be satisfactorily predicted theoretically. The production of 3CH2 + H2O by the singlet to triplet surface crossing, predicted to occur at 4.3 kcal mol-1 above the H2C...OH2 van der Waals complex (which lies 82.7 kcal mol-1 above CH3OH), was neglected in our calculations.     

Reaction of carbon monoxide and hydrogen on neutral Nb8 clusters in the gas phase

Reactions of neutral V(n), Nb(n), and Ta(n) metal clusters (n< or =11) with CO+H(2) mixed gases and CH(3)OH in a flow tube reactor (1-50 Torr) are studied by time of flight mass spectroscopy and density functional theory calculations. Metal clusters are generated by laser ablation, and reactants and products are ionized by low fluence (approximately 200 microJ/cm(2)) 193 nm excimer laser light. Nb(n) clusters exhibit strong size dependent reactivity in reactions both with CO+H(2) and CH(3)OH compared with V(n) and Ta(n) clusters. A "magic number" (relatively intense) mass peak at Nb(8)COH(4) is observed in the reaction of Nb(n) clusters with CO+H(2), and CH(3)OH is suggested to be formed. This feature at Nb(8)COH(4) remains the most intense peak independent of the relative concentrations of CO and H(2) in the flow tube reactor. No other Nb(n), Ta(n), or V(n) feature behaves in this manner. In reactions of CH(3)OH with metal clusters M(n) (M=V, Nb, and Ta, n=3-11), nondehydrogenated products M(n)COH(4)/M(n)CH(3)OH are only observed on Nb(8) and Nb(10), whereas dehydrogenated products M(n)CO/CM(n)O are observed for all other clusters. These observations support the suggestion that CH(3)OH can be formed on Nb(8) in the reaction of Nb(n) with CO+H(2). A reaction mechanism is suggested based on the experimental results and theoretical calculations of this work and of those in the literature. Methanol formation from CO+H(2) on Nb(8) is overall barrierless and thermodynamically and kinetically favorable.     

Dual-level direct dynamics studies for the reactions of dimethyl ether with hydrogen atom and methyl radical

The dual-level direct dynamics approach is employed to study the dynamics of the CH(3)OCH(3) + H (R1) and CH(3)OCH(3) + CH(3) (R2) reactions. Low-level calculations of the potential energy surface are carried out at the MP2/6-311+G(d,p) level of theory. High-level energetic information is obtained at the QCISD(T) level of theory with the 6-311+G(3df,3pd) basis set. The dynamics calculations are performed using variational transition state theory (VTST) with the interpolated single-point energies (ISPE) method, and small-curvature tunneling (SCT) is included. It is shown that the reaction of CH(3)OCH(3) with H (R1) may proceed much easier and with a lower barrier height than the reaction with CH(3) radical (R2). The calculated rate constants and activation energies are in good agreement with the experimental values. The calculated rate constants are fitted to k(R1) = 1.16 x 10(-19) T(3) exp(-1922/T) and k(R2) = 1.66 x 10(-28) T(5) exp(-3086/T) cm(3) mol(-1) s(-1) over a temperature range 207-2100 K. Furthermore, a small variational effect and large tunneling effect in the lower temperature range are found for the two reactions.     

Theoretical Investigation of CH3CF2O2+HOO Reaction

The reactions of CH3CF2O2 with HOO are important chemical cyclic processes of photochemical contamination. In this paper, the reaction pathways and reaction mechanism of CH3CF2O2+HOO are investigated extensively with the Gaussian 98 package at the B3LYP/6-311++G** basis sets. The use of vibrational mode analysis and electron population analysis to reveal the reaction mechanism is firstly reported. The study shows that CH3CF2CO2+HOO→IM1→TS1→CH3CF2O2H+O2 channel is the energetically most favorable, CH3CF2CO2H and O2 are the principal products, and the formation of CH3OH and CF2O is also possible.     

Theoretical study on the reaction CX3 + SiH(CH3)3 (X = H, F)

Theoretical investigations are carried out on the multiple-channel reactions, CH(3) + SiH(CH(3))(3) → products and CF(3) + SiH(CH(3))(3) → products. The minimum energy paths (MEP) are calculated at the MP2/6-311 + G(d,p) level, and energetic information is further refined by the MC-QCISD (single point) method. The rate constants for major reaction channels are calculated by the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) correction over the temperature range 200-1500 K. The theoretical rate constants are in good agreement with the available experimental data and are found to be k(1a)(T) = 1.93 × 10(-24) T(3.15) exp(-1214.59/T) and k(2a)(T) = 1.33 × 10(-25) T(4.13) exp(-397.94/T) (in unit of cm(3) molecule(-1) s(-1)). Our calculations indicate that hydrogen abstraction channel from SiH group is the major channel due to the smaller barrier height among five channels considered.     

Implementation of a microcanonical variational transition state theory for direct dynamics calculations of rate constants

Calculation of microcanonical rate constants has been an important field in chemical dy-namic studies for many years because it can be used not only to give good prediction of rate con-stants in microcanonical assembly, but also to calculate rate constants with certain conserved quantum numbers such as the total angular momentum, and in turn, can be easily converted into thermal rate constants[1—3]. The widely used method for calculating microcanonical rate constants of unimolecular reac-tions…     

Full-dimensional quantum dynamics study of exchange processes for the D + H2O and D + HOD reactions

The exchange processes of D + H(2)O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H(2)O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H(2)O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D(') + H(2)O → D(')OH + H and D(') + HOD → D(')OD + H reactions are substantially larger than those for the D(') + HOD → HOD(') + D reaction, indicating that the D(')/H exchange reactions are much more favored than the D(')/D exchange.     

Reactions of hydrogen atom with hydrogen peroxide

Rate coefficients are calculated using canonical variational transition state theory with multidimensional tunneling (CVT/SCT) for the reactions H + H2O2 --> H2O + OH (1a) and H + H2O2 --> HO2 + H2 (1b). Reaction barrier heights are determined using two theoretical approaches: (i) comparison of parametrized rate coefficient calculations employing CVT/SCT to experiment and (ii) high-level ab initio methods. The evaluated experimental data reveal considerable variations of the barrier height for the first reaction: although the zero-point-exclusive barrier for (1a) derived from the data by Klemm et al. (First Int. Chem. Kinet. Symposium 1975, 61) is 4.6 kcal/mol, other available measurements result in a higher barrier of 6.2 kcal/mol. The empirically derived zero-point-exclusive barrier for (1b) is 10.4 kcal/mol. The electronic structure of the system at transition state geometries in both reactions was found to have "multireference" character; therefore special care was taken when analyzing electronic structure calculations. Transition state geometries are optimized by multireference perturbation theory (MRMP2) with a variety of one-electron basis sets, and by a multireference coupled cluster (MR-AQCCSD) method. A variety of single-reference benchmark-level calculations have also been carried out; included among them are BMC-CCSD, G3SX(MP3), G3SX, G3, G2, MCG3, CBS-APNO, CBS-Q, CBS-QB3, and CCSD(T). Our data obtained at the MRMP2 level are the most complete; the barrier height for (1a) using MRMP2 at the infinite basis set limit is 4.8 kcal/mol. Results are also obtained with midlevel single-reference multicoefficient correlation methods, such as MC3BB, MC3MPW, MC-QCISD/3, and MC-QCISD-MPWB, and with a variety of hybrid density functional methods, which are compared with high-level theory. On the basis of the evaluated experimental values and the benchmark calculations, two possible recommended values are given for the rate coefficients.