Three-and many-body interactions

It must be remarked that Emt is actually a difference of energy with respect to a state of the system in whieh the internal geometry of each bonded component of the system is at equilibrium. 

The versions of MM potentials (generally called force fields) for solvent molecules may be simpler, and in fact limited libraries of MM parameters are used to describe internal geometry change effects in liquid systems. We report in Table 8.4 the definitions of the basic, and simpler, expressions used for liquids. 

The tree-body component of the interaction energy of a trimer ABC is defined as  

This function may be computed, point-by-point, over the appropriate Rabc space, either with variational and PT methods, as the dimeric interactions. The results are again affected by BSS errors, and they can be corrected with the appropriate extension of the CP procedure. A complete span of the surface is, of course, by far more demanding than for a dimer, and actually extensive scans of the decomposition of the Eabc potential energy surface have thus far been done for a very limited number of systems. 

Analogous remarks hold for the four-body component AE(ABCD Rabcd) of the cluster expansion energy, as well as for the five- and six- body components, the definition of 

It must be remarked that is actually a difference of energy with respect to a state of 

---

The second strategy we mention in this rapid survey replaces the QM description of the solvent-solvent and solute-solvent with a semiclassical description. There is a large variety of semiclassical descriptions for the interactions involving solvent molecules, but we limit ourselves to recall the (1,6,12) site formulation, the most diffuse. The interaction is composed of three terms defined in the formula by the inverse power of the corresponding interaction term (1 stays for coulombic interaction, 6 for dispersion and 12 for repulsion). Interactions are allowed for sites belonging to different molecules and are all of two-body character (in other words all the three- and many-body interactions appearing in the cluster expansion of the Hss and HMS terms of the Hamiltonian (1.1)... 

Classical molecular dynamics (MD) implementing predetermined potentials, either empirical or derived from independent electronic structure calculations, has been used extensively to investigate condensed-matter systems. An important aspect in any MD simulation is how to describe or approximate the interatomic interactions. Usually, the potentials that describe these interactions are determined a priori and the full interaction is partitioned into two-, three-, and many-body contributions, long- and short-range terms, etc., for which suitable analytical functional forms are devised. Despite the many successes with classical MD, the requirement to devise fixed potentials results in several serious problems... 

Many applications of new force fields and new QM/MM methods of necessity focus on agreement with experimental or otherwise calculated results. Also in this section we will first show that DRF indeed gives a reliable model for static and response potentials and can lead to QM/MM—or even completely MM calculations—that are as good as, e.g., SCF calculations. To that end we point at some results for simple systems like the water and benzene dimers, and the three- and four-body interactions in several systems. 

Three types of surface are in use for water simulations. The first consists of simple empirical models based on the LJ-C potential. There seems to be no purpose in continuing to develop and use such models as they give little, if any, new information. A second group attempts to improve the accuracy of the potential using semiempirical methods based on a comprehensive set of experimental data. These models allow for physical phenomena such as intramolecular relaxation, electrostatic induced terms, and many-body interactions, all of which are difficult to incorporate correctly in liquid water theories. There is room for much more work in these areas. The third group makes use of the most advanced ab initio methods to develop accurate potentials from first principles. Such calculations are now converging with parameterized surfaces based on accurate semiempirical models. Over the next few years it seems very likely that the continued application of the second and third approaches will result in a potential energy surface that achieves quantitative accuracy for water in the condensed phase. 

Much work has recently been carried out to quantify the three-body and many-body interactions in small clusters (mostly rare gases and water), which have implications on the liquid state properties, however here we consider some studies that have directly determined their influence on the bulk fluid properties. In reference87 the significant influence of three-body interactions on properties of rare gas fluids is discussed, and a recent manuscript by Szalewicz et al,107 thoroughly reviews the importance of many-body forces in general. Here we just summarise some important recent results. 

X 10 ppm amagaH. The coefficients a, G and g arise from two-, three- and four-body interactions, respectively. At low densities the shielding constant depends linearly on density, whereas at high pressures many-body collisions become important as well and cause deviation from linearity. The second virial coefficient arises from the Xe-Xe pair interactions (with the potential V(r), r is the interatomic separation) and can be presented... 

Spectral moments may be computed from expressions such as Eqs. 5.15 or 5.16. Furthermore, the theory of virial expansions of the spectral moments has shown that we may consider two- and three-body systems, without regard to the actual number of atoms contained in a sample if gas densities are not too high. Near the low-density limit, if mixtures of non-polar gases well above the liquefaction point are considered, a nearly pure binary spectrum may be expected (except near zero frequencies, where the intercollisional process generates a relatively sharp absorption dip due to many-body interactions.) In this subsection, we will sketch the computations necessary for the actual evaluation of the binary moments of low order, especially Eqs. 5.19 and 5.25, along with some higher moments. 

The extension of SCIETs to the many-body interactions is presented in Section V. Rare gases, whose constituents interact through three-body forces, are a test case to examine the validity of the SCIETs in describing real systems. Again, the problem of the thermodynamic consistency is covered in this section, since recent SANS measurements provide the structure factor S(q) at very low-q and allow us to deduce the strength of the three-body interactions. A direct comparison of the theoretical results against sharp experiments is feasible. The conclusions are given in Section VI. 

In order to account for hydrodynamic interactions among the suspended particles, Bossis and Brady (1984) used both pairwise additivity of velocities (mobilities) and forces (resistances), discussing the advantages and disadvantages of each method. While their original work did not take explicit account of three- (or more) body effects, the recent formulation of Durlofsky, Brady, and Bossis (1988) does provide a useful procedure for incorporating both the far-field, many-body interactions and near-field, lubrication forces into the grand resistance and mobility matrices. 

---