The performance of a variety of DFT functionals (BLYP, PBE, B3LYP, B3P86, KMLYP, B1B95, MPWPW91, MPW1B95, BB1K, MPW1K, MPWB1K, and BMK), together with the ab initio methods RHF, RMP2, and G3(MP2)-RAD, and with ONIOM methods based on combinations of these procedures, is examined for calculating the enthalpies of a range of radical reactions. The systems studied include the bond dissociation energies (BDEs) of R-X (R = CH3, CH2F, CH2OH, CH2CN, CH2Ph, CH(CH3)Ph, C(CH3)2Ph; X = H, CH3, OCH3, OH, F), RCH(Ph)-X (R = CH3, CH3CH2, CH(CH3)2, C(CH3)3, CH2F, CH2OH, CH2CN; X = H, F), R-TEMPO (R = CH3, CH2CH3, CH(CH3)2, C(CH3)3, CH2CH2CH3, CH2F, CH2OH, CH2CN, CH(CN)CH3, CH(Cl)CH3; TEMPO = 2,2,6,6,-tetramethylpiperidin-1-yloxyl) and HM1M2-X (M1, M2 = CH2CH(CH3), CH2CH(COOCH3), CH2C(CH3)(COOCH3); X = Cl, Br), the beta-scission energies of RXCH2* and RCH2CHPh* (R = CH3, CH2CH3, CH(CH3)2, C(CH3)3; X = O, S, CH2), and the enthalpies of several radical addition, ring-opening, and hydrogen- and chlorine-transfer reactions. All of the DFT methods examined failed to provide an accurate description of the energetics of the radical reactions when compared with benchmark G3(MP2)-RAD values, with all methods tested showing unpredictable deviations of up to 40 kJ mol-1 or more in some cases. RMP2 also shows large deviations from G3(MP2)-RAD in the absolute values of the enthalpies of some types of reaction and, although it fares somewhat better than the DFT methods in modeling the relative values, it fails for substituents capable of strongly interacting with the unpaired electron. However, it is possible to obtain cost-effective accurate calculations for radical reactions using ONIOM-based procedures in which a high-level method, such as G3(MP2)-RAD, is only used to model the core reaction (which should contain all substituents alpha to the reaction center), and the full system is modeled using a lower-cost procedure such as RMP2.

Conspectus This Account describes how attractive interaction of aromatic rings with other groups can influence and control the stereoselectivity of many reactions. Recent developments in theory have led to improved accuracy in the modeling of aromatic interactions. Quantum mechanical modeling can now provide insights into the roles of these interactions at a level of detail not previously accessible, both for ground-state species and for transition states of chemical reactions. In this Account, we show how transition-state modeling led to the discovery of the influence of aryl groups on the stereoselectivities of several types of organic reactions. These reaction types include asymmetric dihydroxylations, transfer hydrogenations, hetero-Diels–Alder reactions, acyl transfers, and Claisen rearrangements. Our recent studies have led to a novel mechanistic picture for two classes of (4+3) cycloadditions, both of which involve reactions of furans with oxyallyl intermediates. The first class of cycloadditions, developed by Hsung, features neutral oxyallyls containing a chiral oxazolidinone auxiliary. Originally, these cycloadditions were thought rely on differential steric crowding of the two faces of a planar intermediate. Computations reveal a different picture and show that cycloadditions with furan takes place preferentially through the more crowded transition state, with furan adding on the same side as the oxazolidinone’s Ph substituent. The crowded transition state is stabilized by a CH–π interaction between furan and Ph, worth about 2.0 kcal/mol. Stereocontrol in a second class of (4+3) cycloadditions, involving chiral alkoxy siloxyallyl cations, also is controlled by attractive interactions with aromatic rings. Alkoxy groups derived from chiral α-methylbenzyl alcohols are found to favor crowded transition states, where a CH–π interaction is again present between furan and Ar. The cationic cycloadditions are stepwise, while the Hsung cycloadditions are concerted. Our results suggest that this form of CH–π-directed stereocontrol is quite general and is likely to control the stereoselectivities of other addition reactions in which one face of a planar intermediate bears a pendant aromatic substituent.

Summary: Computational chemistry is a valuable complement to experiments in the study of polymerization processes. This article reviews the contribution of computational chemistry to understanding the kinetics and mechanism of reversible addition fragmentation chain transfer (RAFT) polymerization. Current computational techniques are appraised, showing that barriers and enthalpies can now be calculated with kcal accuracy. The utility of computational data is then demonstrated by showing how the calculated barriers and enthalpies enable appropriate kinetic models to be chosen for RAFT. Further insights are provided by a systematic analysis of structure‐reactivity trends. The development of the first computer‐designed RAFT agent illustrates the practical utility of these investigations.

CBS-QB3 enthalpies of reaction have been computed for the conjugate additions of MeSH to six α,β-unsaturated ketones. Compared with addition to methyl vinyl ketone, the reaction becomes 1-3 kcal mol(-1) less exothermic when an α-Me, β-Me, or β-Ph substituent is present on the C=C bond. The lower exothermicity for the substituted enones occurs because the substituted reactant is stabilized more by hyperconjugation or conjugation than the product is stabilized by branching. Substituent effects on the activation energies for the rate-determining step of the thiol addition (reaction of the enone with MeS(-)) were also computed. Loss of reactant stabilization, and not steric hindrance, is the main factor responsible for controlling the relative activation energies in the gas phase. The substituent effects are further magnified in solution; in water (simulated by CPCM calculations), the addition of MeS(-) to an enone is disfavored by 2-6 kcal mol(-1) when one or two methyl groups are present on the C=C bond (ΔΔG(‡)). The use of CBS-QB3 gas-phase energies in conjunction with CPCM solvation corrections provides kinetic data in good agreement with experimental substituent effects. When the energetics of the thiol additions were calculated with several popular density functional theory and ab initio methods (B3LYP, MPW1PW91, B1B95, PBE0, B2PLYP, and MP2), some substantial inaccuracies were noted. However, M06-2X (with a large basis set), B2PLYP-D, and SCS-MP2 gave results within 1 kcal mol(-1) of the CBS-QB3 benchmark values.