Multipole contributions potential model calculations

Another method is to calculate the molecular electronic electrostatic potential by replacing p(r ) in Eq. 19 by its multipole formulation (Eq. 8). The quantity obtained represents the electrostatic potential of a molecule removed from the crystal lattice. First calculations have been performed by the Pittsburgh group (Stewart, Craven, He, and co-workers) [43] electrostatic potential calculations were also derived from the Hansen Coppens [lib] electron density model [41,44], The atomic total electrostatic potential including nuclear contribution may be calculated as ... 

Solvent continuum models are now routinely used in quantum mechanical (QM) studies to calculate solvation effects on molecular properties and reactivity. In these models, the solvent is represented by a dielectric continuum that in the presence of electronic and nuclear charges of the solute polarizes, creating an electrostatic potential, the so-called reaction field . The concept goes back to classical electrostatic schemes by Martin [1], Bell [2] and Onsager [3] who made fundamental contributions to the theory of solutions. Scholte [4] and Kirkwood [5] introduced the use of multipole moment distributions. The first implementation in QM calculations was reported in a pioneer work by Rivail and Rinaldi [6,7], Other fundamental investigations were carried out by Tapia and Goscinski [8], Hilton-McCreery et al. [9] and Miertus et al. [10], Many improvements have been made since then (for a review,... 

Errors (/>(calc)-/ (expt)) in the relaxed crystal structures using a fixed repulsion-dispersion potential (as in Table 11.2) and various electrostatic models derived from a DMA of a 6-31G SCF wave function DMA, full multipoles up to hexadecapole SDMA, all multipoles scaled by 0.9 CHAR, just the charge component of the DMA. The r.m.s. % errors were calculated over the three cell lengths. The electrostatic contribution, (7estat, to the total lattice energy, Us, is given at both the experimental (for the SDMA model) and relaxed crystal structures. This can be compared with the experimental heat of sublimation Af/sutl (Chickos 1987), where available. 

A useful alternative approach is to isolate the components of the perturbation expansion, namely the repulsion, electrostatic interaction, induction, and dispersion terms, and to calculate each of them independently by the most appropriate technique. Thus the electrostatic interaction can be calculated accurately from distributed multipole descriptions of the individual molecules, while the induction and dispersion contributions may be derived from molecular polarizabilities. This approach has the advantage that the properties of the monomers have to be calculated only once, after which the interactions may be evaluated easily and efficiently at as many dimer geometries as required. The repulsion is not so amenable, but it can be fitted by suitable analytic functions much more satisfactorily than the complete potential. The result is a model of the intermolecular potential that is capable of describing properties to a high level of accuracy. 

If the charge distribution is described by a set of distributed multipoles, as described in Section 4.2.3, the coulombic contributions to the intermolecular potential energy are calculated as multipole-multipole terms. The main disadvantage of even a rigorous distributed multipole model is that such a representation is still very localized, so that coulombic energies miss a large part of the penetration contribution. For use in a complete representation of intermolecular interactions, the dispersion, polarization, and repulsion terms must be evaluated separately by some semi-empirical or fiilly empirical method, for example the approximate atom-atom formulations of equations 4.38. 39. This approach has been extensively exploited by S. L. Price and coworkers over the years, in applications to molecular crystals [48]. 

---