Electrostatic energies polarizabilities

The electrostatic energy is calculated using the distributed multipolar expansion introduced by Stone [39,40], with the expansion carried out through octopoles. The expansion centers are taken to be the atom centers and the bond midpoints. So, for water, there are five expansion points (three at the atom centers and two at the O-H bond midpoints), while in benzene there are 24 expansion points. The induction or polarization term is represented by the interaction of the induced dipole on one fragment with the static multipolar field on another fragment, expressed in terms of the distributed localized molecular orbital (LMO) dipole polarizabilities. That is, the number of polarizability points is equal to the number of bonds and lone pairs in the molecule. One can opt to include inner shells as well, but this is usually not useful. The induced dipoles are iterated to self-consistency, so some many body effects are included. 

In this expression, the dipole dipole interactions are included in the electrostatic term rather than in the van der Waals interactions as in Eq. (9.43). Of the four contributions, the electrostatic energy can be derived directly from the charge distribution. As discussed in section 9.2, information on the nonelectrostatic terms can be deduced indirectly from the charge density. The polarizability a, which occurs in the expressions for the Debye and dispersion terms of Eqs. (9.41) and (9.42), can be expressed as a functional of the density (Matsuzawa and Dixon 1994), and also obtained from the quadrupole moments of the experimental charge density distribution (see section 12.3.2). However, most frequently, empirical atom-atom pair potential functions like Eqs. (9.45) and (9.46) are used in the calculation of the nonelectrostatic contributions to the intermolecular interactions. 

Dipole-Dipole Interaction. The first of the four terms in the total electrostatic energy depends on the permanent dipole moment of the solute molecule of radius a (assuming a spherical shape) immersed in a liquid solvent of static dielectric constant D. The function f(D) = 2(D - l)/(2D + 1) is known as the Onsager polarity function. The function used here is [f(D) — f(n2)] so that it is restricted to the orientational polarity of the solvent molecules to the exclusion of the induction polarity which depends on the polarizability as of the solvent molecules, related to the slightly different Debye polarity function q>(n2) according to... 

The electrostatic energy of a water molecule with the polarizability y, which acquires the polarization m in the field Ee, is... 

The electrostatic solvation energy is only a part of the total solvation energy. Cavitation, dispersion and repulsion terms must be added. We show below that the MPE method leads to similar electrostatic energies than the polarizable continuum model (PCM) of Tomasi and co-workers [10], provided the same cavities are used. Therefore, non-electrostatic terms in these methods may be computed using the same computational strategies [15]. We emphasize the fact that accurate non-electrostatic contributions are often difficult to compute since they are based on parameterized formulae that cannot be directly compared to experiment. The obtained data must therefore be used with prudence, especially if they are expected to play a major role in the process under study. Fortunately, in many circumstances, non-electrostatic terms are small and/or vary little, so that they can be neglected. Tunon et al. [80] developed a parameterized expression for the MPE method using an expression of the type... 

Electrostatic Terms. The electrostatic energy of a lattice of atoms of zero polarizability may be calculated exactly by the method of Bertaut, provided the position and charge of each atom in the structure are known, This method involves the infinite sum in reciprocal space... 

In Fig. 6, computation of the variation of [I — L W/kT)] versus the average separation of chromophores R is shown. Calculations were carried out for various values of the electrostatic energy parameter W = R W. To give a more physical meaning to these curves, it can be noted that if the chromophore polarizability and ionization potential are chosen to correspond to Disperse Red chromophores, the curves can be considered to correspond to dipole moment values of 5, 7.5, and 10... 

Fig. 1 a-e. Diagrammatic representation of some electrostatic effects. The solid lines represent electrostatic interactions, and the dashed lines represent polarizabilities, a) Electrostatic energy b) Induced moments AQ" = c) Induction energy in lowest order ) Higher-... 

---