No abstract

We have carried out ab initio calculations for the reaction OH+CH4→H2O+CH3 using second-order Mo/ller–Plesset perturbation theory, employing a very large basis set and scaling all correlation energy for the final energy calculation, but optimizing the equilibrium and transition state structures without scaling (MP-SAC2//MP2). We found that inclusion of correlation energy has an important effect on the geometry, barrier height, and vibrational frequencies of the transition state. The final calculated values for the forward and reverse classical barrier heights are 7.4 and 20.6 kcal/mol, respectively. We have used these with interpolated canonical variational transition state theory and the centrifugal-dominant small-curvature tunneling approximation, including information at the reactants, products, transition state, and two other points along the minimum energy path, to predict the rate constants for the above reaction in the temperature range from 223 to 2400 K. The calculated rate constants agree well with experiment over a wide temperature range.

Ab initio and density functional theory molecular orbital calculations were carried out at both the HF/6-31 +G(d) and B3LYP/6-31+G(d) levels for the four antioxidants, p-hydroxycinnamic acid derivatives, namely, the p-coumaric, caffeic, ferulic, and sinapinic acid and the corresponding radicals, in an attempt to explain the structural dependency of the antioxidant activity of these compounds. Optimized resulting geometries, vibrational frequencies, absolute infrared intensities, and electron-donating ability are discussed. Both the high degree of conjugation and the extended spin delocalization in the phenoxyl radicals offer explanation for the scavenging activity of the four acids. In structurally related compounds, the calculated heat of formation value in radical formation appears as a meaningful molecular descriptor of antioxidant activity in accordance with experimental data. This becomes more clear at the B3LYP level.

In this paper we optimize several algorithms for the computation of reaction rates based on information calculated along minimum energy reaction paths and we evaluate the efficiencies of the optimized algorithms. The investigations are based on the calculation of chemical reaction rate constants using variational transition state theory and multidimensional semiclassical transmission coefficients including reaction path curvature. Several methods are evaluated and compared by a systematic set of applications to test cases involving the hydrogen-atom transfer reactions CH3+H2→CH4+H and OH+H2→H2O+H. For each method we present general recommendations for all algorithmic choices other than gradient step size so that future calculations may be carried out reasonably efficiently by varying only one parameter. In the process of these optimizations we have found that the accuracy of the Euler stabilization method can be significantly increased by choosing the auxiliary parameters differently than in previous work; the optimized algorithm is called ES1*. Our final recommendations for future work are (i) when the Hessian/gradient computational cost ratio is low (≲3): the Page–McIver algorithm with the Hessian recalculated at every step, with a cubic starting step, and with curvature calculated from the derivative of the gradient, and (ii) when the Hessian/gradient computational cost ratio is moderate or large: the ES1* algorithm with a Hessian step size three times larger than the gradient step size, with a quadratic starting step, and with curvature calculated from the derivative of the gradient.