Solvent models microscopic

SASA), a concept introduced by Lee and Richards [9], and the electrostatic free energy contribution on the basis of the Poisson-Boltzmann (PB) equation of macroscopic electrostatics, an idea that goes back to Born [10], Debye and Htickel [11], Kirkwood [12], and Onsager [13]. The combination of these two approximations forms the SASA/PB implicit solvent model. In the next section we analyze the microscopic significance of the nonpolar and electrostatic free energy contributions and describe the SASA/PB implicit solvent model. 

With a given set of potential functions we can evaluate various average properties of the solvent. In particular, we would like to simulate experimentally observed macroscopic properties using microscopic solvent models. To do this we have to exploit the theory of statistical mechanics... 

Development of Quantum Mechanical Solvent Effect Models—Microscopic Electrostatic Contributions. 

Other recent microscopic approaches are based on the Langevin dipoles solvent model or on the all-atom solvent model, using a standard force field with van der Waals and electrostatic terms as well as intramolecular terms, and molecular dynamics simulations of the fluctuation of the solvent and the solute, incorporating the potential from the permanent and induced solvent dipoles in the solute Hamiltonian in a self-consistent way (Luzhkov and Warshel, 1991). 

It is well-known that implicit solvent models use both discrete and continuum representations of molecular systems to reduce the number of degrees of freedom this philosophy and methodology of implicit solvent models can be extended to more general multiscale formulations. A variety of DG-based multiscale models have been introduced in an earlier paper of Wei [74]. Theory for the differential geometry of surfaces provides a natural means to separate the microscopic solute domain from the macroscopic solvent domain so that appropriate physical laws are applied to applicable domains. This portion of the chapter focuses specifically on the extension of the equilibrium electrostatics models described above to nonequilibrium transport problems that are relevant to a variety of chemical and biological S5 ems, such as molecular motors, ion channels, fuel cells, and nanofluidics, with chemically or biologically relevant behavior that occurs far from equilibrium [74-76]. 

Many additional refinements have been made, primarily to take into account more aspects of the microscopic solvent structure, within the framework of diffiision models of bimolecular chemical reactions that encompass also many-body and dynamic effects, such as, for example, treatments based on kinetic theory [35]. One should keep in mind, however, that in many cases die practical value of these advanced theoretical models for a quantitative analysis or prediction of reaction rate data in solution may be limited. 

The relation between the microscopic friction acting on a molecule during its motion in a solvent enviromnent and macroscopic bulk solvent viscosity is a key problem affecting the rates of many reactions in condensed phase. The sequence of steps leading from friction to diflfiision coefficient to viscosity is based on the general validity of the Stokes-Einstein relation and the concept of describing friction by hydrodynamic as opposed to microscopic models involving local solvent structure. In the hydrodynamic limit the effect of solvent friction on, for example, rotational relaxation times of a solute molecule is [ ]... 

It is possible to go beyond the SASA/PB approximation and develop better approximations to current implicit solvent representations with sophisticated statistical mechanical models based on distribution functions or integral equations (see Section V.A). An alternative intermediate approach consists in including a small number of explicit solvent molecules near the solute while the influence of the remain bulk solvent molecules is taken into account implicitly (see Section V.B). On the other hand, in some cases it is necessary to use a treatment that is markedly simpler than SASA/PB to carry out extensive conformational searches. In such situations, it possible to use empirical models that describe the entire solvation free energy on the basis of the SASA (see Section V.C). An even simpler class of approximations consists in using infonnation-based potentials constructed to mimic and reproduce the statistical trends observed in macromolecular structures (see Section V.D). Although the microscopic basis of these approximations is not yet formally linked to a statistical mechanical formulation of implicit solvent, full SASA models and empirical information-based potentials may be very effective for particular problems. 

The integral equation method is free of the disadvantages of the continuum model and simulation techniques mentioned in the foregoing, and it gives a microscopic picture of the solvent effect within a reasonable computational time. Since details of the RISM-SCF/ MCSCF method are discussed in the following section we here briefly sketch the reference interaction site model (RISM) theory. 

We recently proposed a new method referred to as RISM-SCF/MCSCF based on the ab initio electronic structure theory and the integral equation theory of molecular liquids (RISM). Ten-no et al. [12,13] proposed the original RISM-SCF method in 1993. The basic idea of the method is to replace the reaction field in the continuum models with a microscopic expression in terms of the site-site radial distribution functions between solute and solvent, which can be calculated from the RISM theory. Exploiting the microscopic reaction field, the Fock operator of a molecule in solution can be expressed by... 

Solvent effects on chemical equilibria and reactions have been an important issue in physical organic chemistry. Several empirical relationships have been proposed to characterize systematically the various types of properties in protic and aprotic solvents. One of the simplest models is the continuum reaction field characterized by the dielectric constant, e, of the solvent, which is still widely used. Taft and coworkers [30] presented more sophisticated solvent parameters that can take solute-solvent hydrogen bonding and polarity into account. Although this parameter has been successfully applied to rationalize experimentally observed solvent effects, it seems still far from satisfactory to interpret solvent effects on the basis of microscopic infomation of the solute-solvent interaction and solvation free energy. 

Since the interface behaves like a capacitor, Helmholtz described it as two rigid charged planes of opposite sign [2]. For a more quantitative description Gouy and Chapman introduced a model for the electrolyte at a microscopic level [2]. In the Gouy-Chapman approach the interfacial properties are related to ionic distributions at the interface, the solvent is a dielectric medium of dielectric constant e filling the solution half-space up to the perfect charged plane—the wall. The ionic solution is considered as formed... 

Continuum models of solvation treat the solute microscopically, and the surrounding solvent macroscopically, according to the above principles. The simplest treatment is the Onsager (1936) model, where aspirin in solution would be modelled according to Figure 15.4. The solute is embedded in a spherical cavity, whose radius can be estimated by calculating the molecular volume. A dipole in the solute molecule induces polarization in the solvent continuum, which in turn interacts with the solute dipole, leading to stabilization. 

In Eq. (6) Ecav represents the energy necessary to create a cavity in the solvent continuum. Eel and Eydw depict the electrostatic and van-der-Waals interactions between solute and the solvent after the solute is brought into the cavity, respectively. The van-der-Waals interactions divide themselves into dispersion and repulsion interactions (Ed sp, Erep). Specific interactions between solute and solvent such as H-bridges and association can only be considered by additional assumptions because the solvent is characterized as a structureless and polarizable medium by macroscopic constants such as dielectric constant, surface tension and volume extension coefficient. The use of macroscopic physical constants in microscopic processes in progress is an approximation. Additional approximations are inherent to the continuum models since the choice of shape and size of the cavity is arbitrary. Entropic effects are considered neither in the continuum models nor in the supermolecule approximation. Despite these numerous approximations, continuum models were developed which produce suitabel estimations of solvation energies and effects (see Refs. 10-30 in 68)). 

The elasticity can be related to very different contributions to the energy of the interface. It includes classical and nonclassical (exchange, correlation) electrostatic interactions in ion-electron systems, entropic effects, Lennard-Jones and van der Waals-type interactions between solvent molecules and electrode, etc. Therefore, use of the macroscopic term should not hide its relation to microscopic reality. On the other hand, microscopic behavior could be much richer than the predictions of such simplified electroelastic models. Some of these differences will be discussed below. 

The structure of the interface between two immiscible electrolyte solutions (ITIES) has been the matter of considerable interest since the beginning of the last century [1], Typically, such a system consists of water (w) and an organic solvent (o) immiscible with it, each containing an electrolyte. Much information about the ITIES has been gained by application of techniques that involve measurements of the macroscopic properties, such as surface tension or differential capacity. The analysis of these properties in terms of various microscopic models has allowed us to draw some conclusions about the distribution and orientation of ions and neutral molecules at the ITIES. The purpose of the present chapter is to summarize the key results in this field. 

---