DOE OSTI2016
We recently showed for a large dataset of pK a s and reduction potentials that free energies calculated directly within the SMD continuum model compares very well with corresponding thermodynamic cycle calculations in both aqueous and organic solvents (Phys. Chem. Chem. Phys. 2015, 17, 2859). In this paper, we significantly expand the scope of our study to examine the suitability of this approach for the calculation of general solution phase kinetics and thermodynamics, in conjunction with several commonly used solvation models (SMDM062X, SMD-HF, CPCM-UAKS, and CPCM-UAHF) for a broad range of systems and reaction types. This includes cluster-continuum schemes for pK a calculations, as well as various neutral, radical and ionic reactions such as enolization, cycloaddition, hydrogen and chlorine atom transfer, and bimolecular SN2 and E2 reactions. On the basis of this benchmarking study, we conclude that the accuracies of both approaches are generally very similar – the mean errors for Gibbs free energy changes of neutral and ionic reactions are approximately 5 kJ mol -1 and 25 kJ mol -1 respectively. In systems where there are significant structural changes due to solvation, as is the case for certain ionic transition states and amino acids, the direct approach generally afford free energy changes that are in better agreement with experiment. The results indicate that when appropriate combinations of electronic structure methods are employed, the direct approach provides a reliable alternative to the thermodynamic cycle calculations of solution phase kinetics and thermodynamics across a broad range of organic reactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗