Field integrated electrostatic approach

In the above formula, Q is the nuclear coordinate, p, and I/r are the ground state and excited electronic terms. Here Kv is provided through the traditional Rayleigh-Schrodinger perturbation formula and K0 have an electrostatic meaning. This expression will be called traditional approach, which has, in principle, quantum correctness, but requires some amendments when different particular approaches of electronic structure calculation are employed (see the Bersuker s work in this volume). In the traditional formalism the vibronic constants P0 dH/dQ Pr) can be tackled with the electric field integrals at nuclei, while the K0 is ultimately related with electric field gradients. Computationally, these are easy to evaluate but the literally use of equations (1) and (2) definitions does not recover the total curvature computed by the ab initio method at hand. 

Theoretical considerations leading to a density functional theory (DFT) formulation of the reaction field (RF) approach to solvent effects are discussed. The first model is based upon isolelectronic processes that take place at the nucleus of the host system. The energy variations are derived from the nuclear transition state (ZTS) model. The solvation energy is expressed in terms of the electrostatic potential at the nucleus of a pseudo atom having a fractional nuclear charge. This procedure avoids the introduction of arbitrary ionic radii in the calculation of insertion energy, since all integrations involved are performed over [O.ooJ The quality of the approximations made are discussed within the frame of the Kohn-Sham formulation of density functional theory. 

The bottleneck of a calculation in solution is the evaluation of the polarization which, in the case of PCM, corresponds to the evaluation of the apparent surface charges. In particular, the bottleneck is represented by the evaluation of the products between the integral matrices of the electrostatic potential (matrix S in Equation (1.8.6)) or of the normal component of the electric field (matrix D in Equation (1.92)) and the apparent charges vector q. Thus the criterion we use to compare the standard and the simultaneous approach is based on the number of matrix products (Sq or D q) necessary in the whole optimization process. We also remind the reader that the dimension of the matrices is equal to the square of the number of the surface elements. 

An alternative to integral equation theories of the nonprimitive inhomogeneous electric double layer is a mean electrostatic field analysis of an ion-solvent dipole mixture against a charged wall [83-90]. Although this approach has been successful with the primitive model and avoids the difficult problem with the bridge function, it is still in the early stages of development with the nonprimitive electric double layer model. 

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. 

For QM solutes, volume polarization is treated approximately (but accurately [89]) by Eq. (11.17), and Chipman has called this approach surface and simulation of volume polarization for electrostatics [SS(V)PE] [15]. An equivalent form of Eq. (11.17) was actually derived prior to Chipman s work, where it was called the integral equation formalism (lEF) [10, 58]. The equivalence is not obvious, as the original lEF requires the solute s electric field as an input in addition to its electrostatic potential, but it was later shown that the former could be eliminated in order to obtain Eq. (11.17) [9]. The operator K can similarly be manipulated into different forms, by means of the identity [15]... 

Implicit solvation models developed for condensed phases represent the solvent by a continuous electric field, and are based on the Poisson equation, which is valid when a surrounding dielectric medium responds linearly to the charge distribution of the solute. The Poisson equation is actually a special case of the Poisson-Boltzmann (PB) equation PB electrostatics applies when electrolytes are present in solution, while the Poisson equation applies when no ions are present. Solving the Poisson equation for an arbitrary equation requires numerical methods, and many researchers have developed an alternative way to approximate the Poisson equation that can be solved analytically, known as the Generalized Born (GB) approach. The most common implicit models used for small molecules are the Conductor-like Screening Model (COSMO) [96,97], the Dielectric Polarized Continuum Model (DPCM) [98], the Conductor-like modification to the Polarized Continuum Model (CPCM) [99], the Integral Equation Formalism implementation of PCM (lEF-PCM) [100] PB models and the GB SMx models of Cramer and Truhlar [52,57,101,102]. The newest Miimesota solvation models are the SMD (universal Solvation Model based on solute electron Density [57]) and the SMLVE method, which combines the surface and volume polarization for electrostatic interactions model (SVPE) [103-105] with semiempirical terms that account for local electrostatics [106]. Further details on these methods can be found in Chapter 11 of reference 52. 

The simple virtual charge model discussed by Constanciel and Tapia [6] has been developed into an extended generalized Born (EGB) approach. Different approximations have been proposed. Constanciel [40] has analyzed the theoretical basis used as foundations for empirical reaction field approximations through the continuum model to the surrounding medium. Artifacts in the EGB scheme have been clearly identified. The new approximate formulation proposed derives from an exact integral equation of classical electrostatics following a well defined procedure. It is shown there how the wavefunction of solvated species imbedded in cavities formed by interlocking sphere in a polarizable continuum can be computed. 

---