Electrostatic component

Fig. 1. The time evolution (top) and average cumulative difference (bottom) associated with the central dihedral angle of butane r (defined by the four carbon atoms), for trajectories differing initially in 10 , 10 , and 10 Angstoms of the Cartesian coordinates from a reference trajectory. The leap-frog/Verlet scheme at the timestep At = 1 fs is used in all cases, with an all-atom model comprised of bond-stretch, bond-angle, dihedral-angle, van der Waals, and electrostatic components, a.s specified by the AMBER force field within the INSIGHT/Discover program.

Another way of calculating the electrostatic component of solvation uses the Poisson-Boltzmann equations [22, 23]. This formalism, which is also frequently applied to biological macromolecules, treats the solvent as a high-dielectric continuum, whereas the solute is considered as an array of point charges in a constant, low-dielectric medium. Changes of the potential within a medium with the dielectric constant e can be related to the charge density p according to the Poisson equation (Eq. (41)). 

The range of systems that have been studied by force field methods is extremely varied. Some force fields liave been developed to study just one atomic or molecular sp>ecies under a wider range of conditions. For example, the chlorine model of Rodger, Stone and TUdesley [Rodger et al 1988] can be used to study the solid, liquid and gaseous phases. This is an anisotropic site model, in which the interaction between a pair of sites on two molecules dep>ends not only upon the separation between the sites (as in an isotropic model such as the Lennard-Jones model) but also upon the orientation of the site-site vector with resp>ect to the bond vectors of the two molecules. The model includes an electrostatic component which contciins dipwle-dipole, dipole-quadrupole and quadrupole-quadrupole terms, and the van der Waals contribution is modelled using a Buckingham-like function. 

Two important contributions to the study of solvation effects were made by Bom (in 192( and Onsager (in 1936). Bom derived the electrostatic component of the free energ) c solvation for placing a charge within a spherical solvent cavity [Bom 1920], and Onsagi extended this to a dipole in a spherical cavity (Figure 11.21) [Onsager 1936]. In the Bor... 

The final class of methods that we shall consider for calculating the electrostatic compone of the solvation free energy are based upon the Poisson or the Poisson-Boltzmann equatior Ihese methods have been particularly useful for investigating the electrostatic properties biological macromolecules such as proteins and DNA. The solute is treated as a body of co stant low dielectric (usually between 2 and 4), and the solvent is modelled as a continuum high dielectric. The Poisson equation relates the variation in the potential (f> within a mediu of uniform dielectric constant e to the charge density p ... 

Calculated electrostatic component" of the ionic strength effects on second-order rate constants for ions of different charges... 

Vakarelski et al. [88] also investigated the adhesive forces between a colloid particle and a flat surface in solution. In their case they investigated a sihca sphere and a mica surface in chloride solutions of monovalent cations CsCl, KCl, NaCl, and LiCl. The pH was kept at 5.6 for all the experiments. To obtain the adhesive force in the presence of an electrostatic interaction, they summed the repulsive force and the pull-off force (coined foe by the authors ) to obtain a value for the adhesive force that is independent of the electrostatic component. 

The total potential energy of the charged particles (again per mole, with charge ZjF) is the snm of a chemical and an electrostatic component ... 

It follows from Eq. (9.2) that the work function has a chemical and an electrostatic component. Its overall value can be measured, whereas an exact determination of its individual components is not possible. The chemical component depends on the interaction between the charge and the surrounding medium moreover, it includes the work performed in overcoming the image forces. 

Considering that, roughly speaking, the electrostatic component of the solvation free energy varies as the cube of the molecular dipole moment, it becomes obvious that the corrective term (13.1) should be taken into account in the determination of differential solvation properties of very polar solutes. In the computation of transfer free energies across an interface, it has been suggested that equation (13.1) be expressed as a function of the number density of one of the two media, so that the correction is zero in solvent 1 and zl,l lsl lll in solvent 2 [115]. 

Abstract This chapter reviews the theoretical background for continuum models of solvation, recent advances in their implementation, and illustrative examples of their use. Continuum models are the most efficient way to include condensed-phase effects into quantum mechanical calculations, and this is typically accomplished by the using self-consistent reaction field (SCRF) approach for the electrostatic component. This approach does not automatically include the non-electrostatic component of solvation, and we review various approaches for including that aspect. The performance of various models is compared for a number of applications, with emphasis on heterocyclic tautomeric equilibria because they have been the subject of the widest variety of studies. For nonequilibrium applications, e.g., dynamics and spectroscopy, one must consider the various time scales of the solvation process and the dynamical process under consideration, and the final section of the review discusses these issues. 

The present chapter thus provides an overview of the current status of continuum models of solvation. We review available continuum models and computational techniques implementing such models for both electrostatic and non-electrostatic components of the free energy of solvation. We then consider a number of case studies, with particular focus on the prediction of heterocyclic tautomeric equilibria. In the discussion of the latter we center attention on the subtleties of actual chemical systems and some of the dangers of applying continuum models uncritically. We hope the reader will emerge with a balanced appreciation of the power and limitations of these methods. 

CM1 charges in the calculation of AGenp, corrections for charge inadequacies appear in G " )s and it is not possible to separate the electrostatic and non-electrostatic components of the free energy of solvation. 

It is then possible, within a Noyes like charging model, to obtain the electrostatic component of the insertion energy expressed as a sum of two contribution coming from two opposite isoelectronic processes that takes... 

---