Energy landscapes of nucleophilic substitution reactions: a comparison of density functional theory and coupled cluster methods

We have carried out a detailed evaluation of the performance of all classes of density functional theory (DFT) for describing the potential energy surface (PES) of a wide range of nucleophilic substitution (SN2) reactions involving, amongst others, nucleophilic attack at carbon, nitrogen, silicon, and sulfur. In particular, we investigate the ability of the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA as well as hybrid DFT to reproduce high-level coupled cluster (CCSD(T)) benchmarks that are close to the basis set limit. The most accurate GGA, meta-GGA, and hybrid functionals yield mean absolute deviations of about 2 kcal/mol relative to the coupled cluster data, for reactant complexation, central barriers, overall barriers as well as reaction energies. For the three nonlocal DFT classes, the best functionals are found to be OPBE (GGA), OLAP3 (meta-GGA), and mPBE0KCIS (hybrid DFT). The popular B3LYP functional is not bad but performs significantly worse than the best GGA functionals. Furthermore, we have compared the geometries from several density functionals with the reference CCSD(T) data. The same GGA functionals that perform best for the energies (OPBE, OLYP), also perform best for the geometries with average absolute deviations in bond lengths of 0.06 A and 0.6 degrees, even better than the best meta-GGA and hybrid functionals. In view of the reduced computational effort of GGAs with respect to meta-GGAs and hybrid functionals, let alone coupled cluster, we recommend the use of accurate GGAs such as OPBE or OLYP for the study of SN2 reactions.