Benchmark ab initio energy profiles for the gas-phase Y- + CH3X → CH3Y + X- (X,Y = F,Cl,Br): validation of hybrid DFT methods

Abstract

The energetics of the gas-phase S(N)2 reactions Y- + CH3X- --> CH3Y + X- (where X,Y = F,Cl,Br), were studied using (variants on) the recent Wi and W2 ab initio computational thermochemistry methods. These calculations involve CCSD and CCSD(T) coupled cluster methods, basis sets of up to spdfgh quality, extrapolations to the one-particle basis set limit, and contributions of inner-shell correlation, scalar relativistic effects, and (where relevant) first-order spin-orbit coupling. Our computational predictions are in excellent agreement with experimental data where these have small error bars; in a number of other instances reexamination of the experimental data may be in order. Our computed benchmark data (including cases for which experimental data are unavailable altogether) are used to assess the quality of a number of compound thermochemistry schemes such as G2 theory, G3 theory, and CBS-QB3, as well as a variety of density functional theory methods. Upon applying some modifications to the level of theory used for the reference geometry (adding diffuse functions, replacing B3LYP by the very recently proposed mPW1K functional [Lynch, B. J.; Fast, P. L.; Harris, M.; Truhlar, D. G. J. Phys. Chem. A 2000, 104, 4811]), the compound methods appear to perform well. Only the "half-and-half" functionals BH&HLYP and mPWH&HPW91 and the empirical mPW1K functional, consistently find all required stationary points; the other functionals fail to find a transition state in the F/Br case. BH&HLYP and mPWH&HPW91 somewhat overcorrect for the tendency of B3LYP land, to a somewhat lesser extent, mPW1PW91) to underestimate barrier heights. The Becke97 and Becke97-1 functionals perform similarly to B3LYP for the problem under study, while the HCTH and HCTH-120 functionals both appear to underestimate central barriers. HCTH underestimates complexation energies; this problem is resolved in HCTH-120. mPW1K appears to exhibit the best performance of the functionals considered, although its energetics are still inferior to the compound thermochemistry methods. mPW1K, however, appears to be very suitable for generating reference geometries for more elaborate thermochemical methods in kinetics applications.