Variation of reaction dynamics for OH hydrogen abstraction from glycine between ab initio levels of theory

The variation in reaction dynamics of OH hydrogen abstraction from glycine between HF, MP2, CCSD(T), M05-2X, BHandHLYP, and B3LYP levels was demonstrated. The abstraction mode shows distinct patterns between these five levels and determines the barrier height, and the spin density transfer between OH radical and glycine. These differences are mainly resulted from the spin density distribution and geometry of the alpha carbon during the abstraction. The captodative effect which is commonly believed as one of the major factors to stabilize the caron-centered radical can only be observed in DFT levels but not in HF and MP2 levels. Difference in the abstraction energy were found in these calculation levels, by using the result of CCSD(T) as reference, B3LYP, BHandHLYP, and M05-2X underestimated the reaction barrier about 5.1, 0.1, and 2.4 kcal mol^-1, while HF and MP2 overestimated 19.1 kcal mol^-1 and 1.6 kcal mol^-1, respectively. These differences can be characterized by the vibration mode of imaginary frequency of transition states, which indicates the topology around transition states and determines reaction barrier height. In this model system, BHandHLYP provides the best prediction of the energy barrier among those tested methods.