Ab Initio Study of Hydrogen Abstraction Reactions

Abstract

Ab initio electronic structure calculations have been carried out on six hydrogen abstraction reactions of the form, X−H + Y• → X• + H−Y, where X, Y = CH3, NH2, and OH. Geometric structures for the reactants, reactant complexes, transition states, product complexes, and products of each reaction have been gradient optimized at both the UMP2 and DFT(B3LYP) theory levels using the 6-311++G(2d,p) basis set. The character of each stationary state as a minimum or saddle point was determined by a harmonic force field calculation. PMP2 and CCSD(T) energies were also calculated at the UMP2 optimized geometries. CCSD(T) geometry optimizations were carried out for selected cases to resolve differences in results between UMP2 and DFT. The calculated reaction energies, barrier heights, and weak complex stabilities are compared to experiment, where possible, and among the different theory levels. The geometric structures and wave function properties are compared between the UMP2, CCSD(T), and DFT(B3LYP) methods. In general, all the methods predict reaction energies to within several kcal/mol of experiment. Calculated DFT(B3LYP) activation barriers are usually equal to or smaller than fitted Arrehenius equation activation energies (E a), with the gap increasing as X and Y are more electronegative. The DFT(B3LYP) activation barrier for the CH4 + CH3 • reaction is, therefore, in the best agreement with experiment. The UMP2, PMP2//UMP2, and CCSD(T)//UMP2 methods give activation energies that are somewhat greater than E a, as expected. The geometric structures of the three exchange reaction (X ≠ Y) transition states (TS) generally agree reasonably well between the UMP2 and DFT(B3LYP) methods, with active site bond length differences reflecting variations in the calculated exothermicities of the reactions. The UMP2 and DFT(B3LYP) methods give different conformations for the NH3 + OH• transition state, and CCSD(T) endorses the DFT results. For the H2O + OH• reaction UMP2 gives the symmetric structure as an energy minimum. DFT(B3LYP) predicts the symmetric structure to be the TS, in agreement with the optimized CCSD(T) result. For the NH3 + NH2 • reaction the UMP2 transition state is symmetric but with an unusually low imaginary reaction path frequency, while the DFT(B3LYP) frequency is reasonable. CCSD(T)//UMP2 and DFT(B3LYP) hydrogen-bonding energies for the reaction and product complexes are very similar and generally smaller than the UMP2 and PMP2//UMP2 results. The equilibrium geometric structures at the UMP2 and DFT(B3LYP) levels generally agree reasonably well. However, for the NH2···H−OH product complex DFT(B3LYP) gives a planar conformation, while UMP2 predicts a C s nonplanar configuration.

Communication: An accurate full 15 dimensional permutationally invariant potential energy surface for the OH + CH4 → H2O + CH3 reaction

“…This is consistent with the recent theoretical work. 16,38 The P-vdW complex, on the other hand, was found in a H 3 C ... HOH configuration, stabilized by 1.84 kcal/mol with respect to the product asymptote, consistent with the previous theoretical reports. 21,31 It is interesting that some small imaginary frequencies persist and are hard to eliminate at the ab initio level, as shown in the supplementary material, 37 due apparently to the floppy nature of the system.…”

Assessment of theoretical methods for the study of hydrogen abstraction kinetics of global warming gas species during their degradation and byproduct formation (IUPAC Technical Report)

“…The contradiction was explained based on the possibility of hydrogen bonding in condensed phases, which blocks the reaction site for the hydrogen abstraction from the hydroxyl group. Hydrogen abstraction from 2-propanol by hydroxyl radical was studied using MP2/6-31G(d) with scaling of the correlation energy [35]. The calculated rate constants for the combined hydrogen atom abstractions were in good agreement with the available experimental data.…”

Density functional, single and multireference perturbation theory study of the reaction Σg3O2+HOCH2CH2⋅→HOO⋅+HOCH=CH2, modeling an important step in tropospheric benzene oxidation

“…All DFT barriers appear to be rather low ͑3 to 5 kcal mol Ϫ1 ͒, and it is possible that this method is prone to underestimate the energy of hydrogen abstraction barriers. Indeed, this feature had been observed in a recent DFT study on H abstraction reactions, 30 for which suitable experimental and high-level computational references were available. The underestimation appeared to be larger when more electronegative atoms were involved in hydrogen exchange.…”