A Universal Approach to Solvation Modeling

Abstract

Continuum mean-field models that have been carefully designed to address the various electrostatic and nonelectrostatic interactions that develop between a molecule and a surrounding medium are particularly efficient tools for studying the effects of condensed phases on molecular structure, energetics, properties, spectra, interaction potentials, and dynamics. The SM8 model may be combined with density functional theory or Hartree-Fock theory to describe a solute's electronic structure and its self-consistent-field polarization by a solvent. A key feature is the use of class IV charge models to obtain accurate charge distributions (either in the vapor phase or in solution), even when using small basis sets that are affordable for large systems. A second key feature is that nonelectrostatic effects due to cavity formation, dispersion interactions, and changes in solvent structure are included in terms of empirical atomic surface tensions that depend on geometry but do not require atom-type assignments by the user. Use of an analytic surface area algorithm provides very stable energy gradients that allow geometry optimization in solution. The SM8 continuum model, the culmination of a series of SMx models (x = 1-8), permits the modeling of such diverse media as aqueous and organic solvents, soils, lipid bilayers, and air-water interfaces. In addition to predicting accurate transfer free energies between gaseous and condensed phases or between two different condensed phases, SMx models have been useful for predicting the significant influence of condensed phases on processes associated with a change in molecular charge, including acid/base equilibria and oxidation/reduction processes. In this Account, we provide an overview of the algorithms associated with the computation of free energies of solvation in the SM8 model. We also compare the accuracies of the SM8 model with those of other continuum solvation models. Finally, we highlight applications of the SM8 models to compute ionic solvation free energies, oxidation and reduction potentials, and pK(a) values.