Electrostatic polarization continuum

Over the last years, the basic concepts embedded within the SCRF formalism have undergone some significant improvements, and there are several commonly used variants on this idea. To exemplify the different methods and how their results differ, one recent work from this group [52] considered the sensitivity of results to the particular variant chosen. Due to its dependence upon only the dipole moment of the solute, the older approach is referred to herein as the dipole variant. The dipole method is also crude in the sense that the solute is placed in a spherical cavity within the solute medium, not a very realistic shape in most cases. The polarizable continuum method (PCM) [53,54,55] embeds the solute in a cavity that more accurately mimics the shape of the molecule, created by a series of overlapping spheres. The reaction field is represented by an apparent surface charge approach. The standard PCM approach utilizes an integral equation formulation (IEF) [56,57], A variant of this method is the conductor-polarized continuum model (CPCM) [58] wherein the apparent charges distributed on the cavity surface are such that the total electrostatic potential cancels on the surface. The self-consistent isodensity PCM procedure [59] determines the cavity self-consistently from an isodensity surface. The UAHF (United Atom model for Hartree-Fock/6-31 G ) definition [60] was used for the construction of the solute cavity. 

Generalized Bom (GB) approach. The most common implicit models used for small molecules are the Conductor-Like Screening Model (COSMO) [77,78], the DPCM [79], the Conductor-Like Modification to the Polarized Continuum Model (CPCM) [80,81], the Integral Equation Formalism Implementation of PCM (IEF-PCM) [82] PB models, and the GB SMx models of Cramer and Truhlar [23,83-86]. The newest Minnesota solvation models are the SMD universal Solvation Model based on solute electron density [26] and the SMLVE method, which combines the surface and volume polarization for electrostatic interactions model (SVPE) [87-89] with semiempirical terms that account for local electrostatics [90]. Further details on these methods can be found in Chapter 11 of Reference [23]. 

Given the diversity of different SCRF models, and the fact that solvation energies in water may range from a few kcal/mol for say ethane to perhaps 100 kcal/mol for an ion, it is difficult to evaluate just how accurately continuum methods may in principle be able to represent solvation. It seems clear, however, that molecular shaped cavities must be employed, the electrostatic polarization needs a description either in terms of atomic charges or quite high-order multipoles, and cavity and dispersion terms must be included. Properly parameterized, such models appear to be able to give absolute values with an accuracy of a few kcal/mol. Molecular properties are in many cases also sensitive to the environment, but a detailed discussion of this is outside the scope of this book. ... 

The continuum electrostatic approximation is based on the assumption that the solvent polarization density of the solvent at a position r in space is linearly related to the total local electric field at that position. The Poisson equation for macroscopic continuum media... 

The continuum model, in which solvent is regarded as a continuum dielectric, has been used to study solvent effects for a long time [2,3]. Because the electrostatic interaction in a polar system dominates over other forces such as van der Waals interactions, solvation energies can be approximated by a reaction field due to polarization of the dielectric continuum as solvent. Other contributions such as dispersion interactions, which must be explicitly considered for nonpolar solvent systems, have usually been treated with empirical quantity such as macroscopic surface tension of solvent. 

Several remedies have been suggested for improving the PB based pKa prediction methods. Most of them are based on strategies that combine conformational flexibility with the PB calculation. You and Bashford included multiple conformers by systematically scanning the side chain torsion angles [107], Alexov and Gunner used Monte-Carlo protocol to sample positions of hydroxyl and other polar protons [1], This method, referred to as the multi-conformation continuum electrostatic (MCCE), was later extended to include rotamers for residues that have strong electrostatic... 

Continuum models have a long and honorable tradition in solvation modeling they ultimately have their roots in the classical formulas of Mossotti (1850), Clausius (1879), Lorentz (1880), and Lorenz (1881), based on the polarization fields in condensed media [32, 57], Chemical thermodynamics is based on free energies [58], and the modem theory of free energies in solution is traceable to Bom s derivation (1920) of the electrostatic free energy of insertion of a monatomic ion in a continuum dielectric [59], and Kirkwood and Onsager s... 

---