Activation Strain Analysis of S_N2 Reactions at C, N, O, and F Centers

Fundamental principles that determine chemical reactivity and reaction mechanisms are the very foundation of chemistry and many related fields of science. Bimolecular nucleophilic substitutions (S_N2) are among the most common and therefore most important reaction types. In this report, we examine the trends in the S_N2 reactions with respect to increasing electronegativity of the reaction center by comparing the well-studied backside S_N2 Cl^- + CH₃Cl with similar Cl^- substitutions on the isoelectronic series with the second period elements N, O, and F in place of C. Relativistic (ZORA) DFT calculations are used to construct the gas phase reaction potential energy surfaces (PES), and activation strain analysis, which allows decomposition of the PES into the geometrical strain and interaction energy, is employed to analyze the observed trends. We find that S_N2@N and S_N2@O have similar PES to the prototypical S_N2@C, with the well-defined reaction complex (RC) local minima and a central barrier, but all stationary points are, respectively, increasingly stable in energy. The S_N2@F, by contrast, exhibits only a single-well PES with no barrier. Using the activation strain model, we show that the trends are due to the interaction energy and originate mainly from the decreasing energy of the empty acceptor orbital (σ*_A-Cl) on the reaction center A in the order of C, N, O, and F. The decreasing steric congestion around the central atom is also a likely contributor to this trend. Additional decomposition of the interaction energy using Kohn-Sham molecular orbital (KS-MO) theory provides further support for this explanation, as well as suggesting electrostatic energy as the primary reason for the distinct single-well PES profile for the FCl reaction.