Ab initio study of the two competitive channels HO2 + H → H2 + O2 and HO2 + H → H2O + O

Saddle point geometries and barrier heights have been calculated for the H abstraction reaction HO₂(²A″)+H(²S) → H₂(¹Σ⁺_g)+O₂(³Σ⁻_g) and the concerted H approach-O removing reaction HO₂ (²A″)+H(²S) → H₂O(¹A₁)+O(³P) by using SDCI wavefunctions with a valence double-zeta plus polarization basis set. The saddle points are found to be of C_s symmetry and the barrier heights are respectively 5.3 and 19.8 kcal by including size consistent correction. Moreoever kinetic parameters have been evaluated within the framework of the TST theory. So activation energies and the rate constants are estimated to be respectively 2.3 kcal and 0.4×10⁹ ℓ mol⁻¹ s⁻¹ for the first reaction, 20.0 kcal and 5.4.10⁻⁵ ℓ mol⁻¹ s⁻¹ for the second. Comparison of these results with experimental determinations shows that hydrogen abstraction on HO₂ is an efficient mechanism for the formation of H₂ + O₂, while the concerted mechanism envisaged for the formation of H₂O + O is highly unlikely.     

O(3P) + CO2 collisions at hyperthermal energies: dynamics of nonreactive scattering, oxygen isotope exchange, and oxygen-atom abstraction

The dynamics of O((3)P) + CO(2) collisions at hyperthermal energies were investigated experimentally and theoretically. Crossed-molecular-beams experiments at = 98.8 kcal mol(-1) were performed with isotopically labeled (12)C(18)O(2) to distinguish products of nonreactive scattering from those of reactive scattering. The following product channels were observed: elastic and inelastic scattering ((16)O((3)P) + (12)C(18)O(2)), isotope exchange ((18)O + (16)O(12)C(18)O), and oxygen-atom abstraction ((18)O(16)O + (12)C(18)O). Stationary points on the two lowest triplet potential energy surfaces of the O((3)P) + CO(2) system were characterized at the CCSD(T)/aug-cc-pVTZ level of theory and by means of W4 theory, which represents an approximation to the relativistic basis set limit, full-configuration-interaction (FCI) energy. The calculations predict a planar CO(3)(C(2v), (3)A') intermediate that lies 16.3 kcal mol(-1) (W4 FCI excluding zero point energy) above reactants and is approached by a C(2v) transition state with energy 24.08 kcal mol(-1). Quasi-classical trajectory (QCT) calculations with collision energies in the range 23-150 kcal mol(-1) were performed at the B3LYP/6-311G(d) and BMK/6-311G(d) levels. Both reactive channels observed in the experiment were predicted by these calculations. In the isotope exchange reaction, the experimental center-of-mass (c.m.) angular distribution, T(θ(c.m.)), of the (16)O(12)C(18)O products peaked along the initial CO(2) direction (backward relative to the direction of the reagent O atoms), with a smaller isotropic component. The product translational energy distribution, P(E(T)), had a relatively low average of = 35 kcal mol(-1), indicating that the (16)O(12)C(18)O products were formed with substantial internal energy. The QCT calculations give c.m. P(E(T)) and T(θ(c.m.)) distributions and a relative product yield that agree qualitatively with the experimental results, and the trajectories indicate that exchange occurs through a short-lived CO(3)* intermediate. A low yield for the abstraction reaction was seen in both the experiment and the theory. Experimentally, a fast and weak (16)O(18)O product signal from an abstraction reaction was observed, which could only be detected in the forward direction. A small number of QCT trajectories leading to abstraction were observed to occur primarily via a transient CO(3) intermediate, albeit only at high collision energies (149 kcal mol(-1)). The oxygen isotope exchange mechanism for CO(2) in collisions with ground state O atoms is a newly discovered pathway through which oxygen isotopes may be cycled in the upper atmosphere, where O((3)P) atoms with hyperthermal translational energies can be generated by photodissociation of O(3) and O(2).     

Kinetics and mechanism of the beta-alanine + OH gas phase reaction: a quantum mechanical approach

The OH hydrogen abstraction reaction from beta-alanine has been studied using the BHandHLYP hybrid HF-density functional and 6-311G(d,p) basis sets. The energies have been improved by single point calculations at the CCSD(T)/6-311G(d,p) level of theory. The structures of the different stationary points are discussed. Reaction profiles are modeled including the formation of pre-reactive and product complexes. Negative net activation energy is obtained for the overall reaction. A complex mechanism is proposed, and the rate coefficients are calculated using transition state theory over the temperature range of 250-400 K. The rate coefficients are proposed for the first time and it was found that in the gas phase the hydrogen abstraction occurs mainly from the CH(2) group next to the amino end. The following expressions, in cm(3) mol(-1) s(-1), are obtained for the overall rate constants, at 250-400 and 290-310 K, respectively: k(250-400)= 2.36 x 10(-12) exp(340/T), and k(290-310)= 1.296 x 10(-12) exp(743/T). The three parameter expression that best describes the studied reaction is k(250-400)= 1.01 x 10(-21)T(3.09) exp(1374/T). The beta-alanine + OH reaction was found to be 1.5 times faster than the alpha-alanine + OH reaction.     

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.     

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.     

Trajectory surface hopping study of the O(3P) + ethylene reaction dynamics

Spin-orbit coupling (SOC) induced intersystem crossing (ISC) has long been believed to play a crucial role in determining the product distributions in the O(3P) + C2H4 reaction. In this paper, we present the first nonadiabatic dynamics study of the title reaction at two center-of-mass collision energies: 0.56 eV, which is barely above the H-atom abstraction barrier on the triplet surface, and 3.0 eV, which is in the hyperthermal regime. The calculations were performed using a quasiclassical trajectory surface hopping (TSH) method with the potential energy surface generated on the fly at the unrestricted B3LYP/6-31G(d,p) level of theory. To simplify our calculations, nonadiabatic transitions were only considered when the singlet surface intersects the triplet surface. At the crossing points, Landau-Zener transition probabilities were computed assuming a fixed spin-orbit coupling parameter, which was taken to be 70 cm-1 in most calculations. Comparison with a recent crossed molecular beam experiment at 0.56 eV collision energy shows qualitative agreement as to the primary product branching ratios, with the CH3 + CHO and H + CH2CHO channels accounting for over 70% of total product formation. However, our direct dynamics TSH calculations overestimate ISC so that the total triplet/singlet ratio is 25:75, compared to the observed 43:57. Smaller values of SOC reduce ISC, resulting in better agreement with the experimental product relative yields; we demonstrate that these smaller SOC values are close to being consistent with estimates based on CASSCF calculations. As the collision energy increases, ISC becomes much less important and at 3.0 eV, the triplet to singlet branching ratio is 71:29. As a result, the triplet products CH2 + CH2O, H + CH2CHO and OH + C2H3 dominate over the singlet products CH3 + CHO, H2 + CH2CO, etc.     

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.     

Reactions of Nb2 and Nb3 with CO, D2, N2, and O2: reconciling experimental kinetics with density functional theory-calculated reaction profiles

Calculated energy profiles for the reactions of neutral Nb(2) and Nb(3) metal clusters with CO, D(2), N(2), and O(2) are presented. In each reaction path, both a physisorption energy minimum, where the reactant remains intact, and a chemisorption energy minimum, where the reactant has dissociated, are calculated and linked by saddle points. We calculate branching ratios for the forward (dissociative) and reverse reactions which we compare with the experimental kinetic data. It is found that a combination of average thermal energies and barrier heights leads to wide variation in branching ratios which compares favourably to previously determined experimental reaction rates.     

Full dimensional calculations of vibrational energies of H3O+ and D3O+

We report full dimensional calculations of vibrational energies of H3O+ and D3O+ using two implementations of the code MULTIMODE. In one implementation, the reference geometry is the minimum of the potential (the standard choice for MULTIMODE). This implementation is not able to readily describe splittings in the vibrational energies due to motion through the inversion barrier. A second implementation, in which the reference geometry is the inversion saddle point, is able to describe the splittings. These full dimensional calculations are done using the realistic, though not spectroscopically accurate, potential of Ojamae, Singer and Shavitt, and the results are compared with experiment.     

Direct dynamics simulations of O(3P) + HCl at hyperthermal collision energies

The dynamics of the O(3P) + HCl reaction at hyperthermal collision energies were investigated using the quasiclassical trajectory method. Stationary points on the OClH 3A" and 3A' potential energy surfaces (PESs) were also examined. The lowest transition state leading to OCl + H on the 3A" surface is 2.26 eV above the reagents at the CCSD(T)/cc-pVTZ level of theory. This saddle point is bent and product-like. Direct dynamics calculations at the MP2/cc-pVTZ level of theory were used to investigate the excitation functions for OH + Cl, OCl + H, and O + H + Cl formation. OCl is formed mainly from small-impact-parameter collisions, and the OCl + H excitation function peaks around 5 eV, where it is similar in magnitude to the OH + Cl excitation function. The shape of the OCl + H excitation function is discussed, and features are identified that should be general to hyperthermal collision dynamics.     

Quasi-classical trajectory study of Si+O2-->SiO+O reaction

Quasi-classical trajectory calculations for the Si(3P)+O2(X 3Sigmag-)-->SiO(X 1Sigma+)+O(1D) reaction have been carried out using the analytical ground 1A' potential energy surface (PES) recently reported by Dayou and Spielfiedel [J. Chem. Phys. 119, 4237 (2003)]. The reaction has been studied for a wide range of collision energies (0.005-0.6 eV) with O2 in its ground rovibrational state. The barrierless PES leads to a decrease of the total reaction cross section with increasing collision energy. It has been brought to evidence that the reaction proceeds through different reaction mechanisms whose contributions to reactivity are highly dependent on the collision energy range. At low collision energy an abstraction mechanism occurs involving the collinear SiOO potential well. The associated short-lived intermediate complex leads to an inverted vibrational distribution peaked at v'=3 and low rotational excitation of SiO(v',j') with a preferentially backward scattering. At higher energies the reaction proceeds mainly through an insertion mechanism involving the bent and linear OSiO deep potential wells and associated long-lived intermediate complexes, giving rise to nearly statistical energy disposals into the product modes and a forward-backward symmetry of the differential cross section.     

Mechanistic and kinetic study of the CH3CO + O2 reaction

Potential-energy surface of the CH3CO + O2 reaction has been calculated by ab initio quantum chemistry methods. The geometries were optimized using the second-order Moller-Plesset theory (MP2) with the 6-311G(d,p) basis set and the coupled-cluster theory with single and double excitations (CCSD) with the correlation consistent polarized valence double zeta (cc-pVDZ) basis set. The relative energies were calculated using the Gaussian-3 second-order Moller-Plesset theory with the CCSD/cc-pVDZ geometries. Multireference self-consistent-field and MP2 methods were also employed using the 6-311G(d,p) and 6-311++G(3df,2p) basis sets. Both addition/elimination and direct abstraction mechanisms have been investigated. It was revealed that acetylperoxy radical [CH3C(O)OO] is the initial adduct and the formation of OH and alpha-lactone [CH2CO2(1A')] is the only energetically accessible decomposition channel. The other channels, e.g., abstraction, HO2 + CH2CO, O + CH3CO2, CO + CH3O2, and CO2 + CH3O, are negligible. Multichannel Rice-Ramsperger-Kassel-Marcus theory and transition state theory (E-resolved) were employed to calculate the overall and individual rate coefficients and the temperature and pressure dependences. Fairly good agreement between theory and experiments has been obtained without any adjustable parameters. It was concluded that at pressures below 3 Torr, OH and CH2CO2(1A') are the major nascent products of the oxidation of acetyl radicals, although CH2CO2(1A') might either undergo unimolecular decomposition to form the final products of CH2O + CO or react with OH and Cl to generate H2O and HCl. The acetylperoxy radicals formed by collisional stabilization are the major products at the elevated pressures. In atmosphere, the yield of acetylperoxy is nearly unity and the contribution of OH is only marginal.     

H2O3 as a reactive oxygen species: formation of 8-oxoguanine from its reaction with guanine

Reaction of guanine with H2O3 in the absence and presence of a water molecule leading to the formation of 8-oxoguanine (8-oxoG) was investigated. Initial calculations were performed using imidazole (Im) as a model for the five-membered ring of guanine. The reactant, intermediate, and product complexes as well as transition states were obtained in gas phase at the B3LYP/6-31+G* and B3LYP/AUG-cc-pVDZ levels of theory. In all the cases, except for the reactions involving imidazole, single-point energy calculations were performed in gas phase at the MP2/AUG-cc-pVDZ level of theory. Solvation calculations in aqueous media were carried out using the polarizable continuum model (PCM) of the self-consistent reaction field (SCRF) theory. Vibrational frequency analysis and intrinsic reaction coordinate (IRC) calculations were performed to ensure that the transition states connected the reactant and product complexes properly. Zero-point energy (ZPE)-corrected total energies and Gibbs free energies at 298.15 K in gas phase and aqueous media were obtained. When a reaction of H2O3 in place of H2O2 with guanine is considered, the major barrier energy which is encountered at the first step is almost halved showing that H2O3 would be much more reactive than H2O2. Considering the reaction schemes investigated here and the observed fact that H2O3 is dissociated easily under ambient conditions, it appears that H2O3 would serve as an effective reactive oxygen species.     

Ab initio potential energy surfaces of the ion-molecule reaction: C2H2 + O+

High level ab initio calculations using complete active space self-consistent field and multi reference single and double excitation configuration interaction methods with cc-pVDZ (correlation consistent polarized valence double zeta) and cc-pVTZ (triple zeta) basis sets have been performed to elucidate the reaction mechanism of the ion-molecule reaction, C2H2(1Sigmag+) + O+(4S), for which collision experiment has been performed by Chiu et al. [J. Chem. Phys. 109, 5300 (1998)]. The minor low-energy process leading to the weak spin-forbidden product C2H2+ (2Piu) + O(1D) has been studied previously and will not be discussed here. The major pathways to form charge-transfer (CT) products, C2H2+ (2Piu) + O(3P) (CT1) and C2H2+ (4A2) + O(3P) (CT2), and the covalently bound intermediates are investigated. The approach of the oxygen atom cation to acetylene goes over an energy barrier TS1 of 29 kcal/mol (relative to the reactant) and adiabatically leads the CT2 product or a weakly bound intermediate Int1 between CT2 products. This transition state TS1 is caused by the avoided crossing between the reactant and CT2 electronic states. As the C-O distance becomes shorter beyond the above intermediate, the C1 reaction pathway is energetically more favorable than the Cs pathway and goes over the second transition state TS2 of a relative energy of 39 kcal/mol. Although this TS connects diabatically to the covalent intermediate Int2, there are many states that interact adiabatically with this diabatic state and these lead to the other charge-transfer product CT1 via either of several nonadiabatic transitions. These findings are consistent with the experiment, in which charge transfer and chemical reaction products are detected above 35 and 39 kcal/mol collision energies, respectively.     

A systematic computational study of the reactions of HO2 with RO2: the HO2 + C2H5O2 reaction

The reaction of HO2 with C2H5O2 has been studied using the density functional theory (B3LYP) and the coupled-cluster theory [CCSD(T)]. The reaction proceeds on the triplet potential energy surface via hydrogen abstraction to form ethyl hydroperoxide and oxygen. On the singlet potential energy surface, the addition-elimination mechanism is revealed. Variational transition state theory is used to calculate the temperature-dependent rate constants in the range 200-1000 K. At low temperatures (e.g., below 300 K), the reaction takes place predominantly on the triplet surface. The calculated low-temperature rate constants are in good agreement with the experimental data. As the temperature increases, the singlet reaction mechanism plays more and more important role, with the formation of OH radical predominantly. The isotope effect of the reaction (DO2 + C2D5O2 vs HO2 + C2H5O2) is negligible. In addition, the triplet abstraction energetic routes for the reactions of HO2 with 11 alkylperoxy radicals (CnHmO2) are studied. It is shown that the room-temperature rate constants have good linear correlation with the activation energies for the hydrogen abstraction.     

Accurate ab initio potential energy surfaces for the 3A' and 3A' electronic states of the O(3P)+HBr system

In this work, we report the construction of potential energy surfaces for the (3)A(') and (3)A(') states of the system O((3)P) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O((3)P) + HBr → OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A(') electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A(') surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A(') and 4.16 kcal/mol for the (3)A(') state.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

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.