Three-body term

The computational effort is significantly increased if three-body terms are included in the model. Even with a simple pairwise model, the non-bonded interactions usually require by far the greatest amount of computational effort. The number of bond, angle and torsional terms increases approximately with the number of atoms (N) in the system, but the number of non-bonded interactions increases with N. There are N(N —l)/2 distinct pairs of... 

The analytic PES function is usually a summation of two- and three-body terms. Spline functions have also been used. Three-body terms are often polynomials. Some of the two-body terms used are Morse functions, Rydberg... 

The term Uss is the solvent-solvent interaction term [the Unb and Uqq terms of eq. (3.1)] and t/ind is the induced dipoles three-body term which includes now the field both from the solute and the solvent. With a potential surface for a solvated solute we can address the important issue of evaluating solvation energies. In principle, one can try to evaluate the average poten-... 

Table IV. Comparison of stability and structure of Ain clusters between ab initio and parameterized interaction results with two- and three-body terms (2+3-b) as well as using only the two-body (2-b) interaction. Binding energies (Dc in eV) per atom, and bond distances (rg in ao) are given...

One additional important reason why nonbonded parameters from quantum chemistry cannot be used directly, even if they could be calculated accurately, is that they have to implicitly account for everything that has been neglected three-body terms, polarization, etc. (One should add that this applies to experimental parameters as well A set of parameters describing a water dimer in vacuum will, in general, not give the correct properties of bulk liquid water.) Hence, in practice, it is much more useful to tune these parameters to reproduce thermodynamic or dynamical properties of bulk systems (fluids, polymers, etc.) [51-53], Recently, it has been shown, how the cumbersome trial-and-error procedure can be automated [54-56A],... 

Extension to a molecule with more than four atoms or to a solid is straightforward. Usually the two-body terms are much larger than the three-body terms, which in turn are greater than the four-body. For ionic solids, for example, the three-body and four-body terms are often neglected. In contrast, for metals and semiconductors including only two-body terms leads to very poor results (see Sutton (Further reading)). 

Sorbie and Murrell have developed a function for triatomic systems based on a many body expansion of PE. It has been applied to repulsive potentials as well and has become an attractive option for fitting ah initio PES. It has sufficient flexibility and extendable to higher polyatomic systems and to two-valued surfaces. The potential for triatomic system is decomposed into one-, two- and three-body terms. The potential is given as... 

The monoatomic terms are constant and can be taken as zero in their ground states. For diatomic terms extended Rydberg (ER) can be used. The three-body term is expressed as a polynomial multiplied by a product of switching functions... 

It should also be noted that ternary and higher order polymer-polymer interactions persist in the theta condition. In fact, the three-parameter theoretical treatment of flexible chains in the theta state shows that in real polymers with finite units, the theta point corresponds to the cancellation of effective binary interactions which include both two body and fundamentally repulsive three body terms [26]. This causes a shift of the theta point and an increase of the chain mean size, with respect to Eq. (2). However, the power-law dependence, Eq. (3), is still valid. The RG calculations in the theta (tricritical) state [26] show that size effect deviations from this law are only manifested in linear chains through logarithmic corrections, in agreement with the previous arguments sketched by de Gennes [16]. The presence of these corrections in the macroscopic properties of experimental samples of linear chains is very difficult to detect. 

Let us make clear now the correspondence between our treatment here and Erdahl s 1978 treatment [4, Sec. 8]. Erdahl works in general Fock space and his operators conserve only the parity of the number of nuclei. He exhibits two families of operators that are polynomials in the annihilation and creation operators containing a three-body and a one-body term. Generic instances of these operators are denoted y and w. The coefficients are real, and Erdahl stresses that this is essential for his treatment. The one-body term is otherwise unrestricted, but the three-body term must satisfy conditions to guarantee that y+y or H +w does not contain a six-body term. For the first family the conditions amount to the three-body term being even under taking the adjoint, and for... 

The A scp term is calculated using the standard CP-method. At the correlated MP2 level, we have shown for several systems [7-10], that the AE terms are usually and systematically smaller than the dominant ( )+ Ecj) terms. The sum of these two terms provides a good approximation to the total interaction energy at the correlated level. It is important to emphasize that the AE values were obtained by making the difference with the values of the CP-corrected subsystems i.e. taking into consideration the "benefit effect" of the superposition of the basis set [3, 6]. As the charge-transfer components are of importance in the two-body interaction, (see a discussion in ref. 10), we will also investigate them separately for the three-body terms in the studied systems. 

Let us now turn to the discussion of the three-body terms at the correlated level. The calculation were performed at the MP3 and MP4 levels, using the SMO-LMBPT formalism. Certain problems arise as to the definition. Namely, there is no general agreement in the literature how to define the three-body... 

A well-known prototype of three-body interaction is the helium trimer [25-28]. The three-body interaction energy is discussed in detail in the above-mentioned works. It is interesting that the results, certainly due to the different approaches used, are not uniform. As demonstrated this by one example, the three-body term is found to be negative at the correlated level for a linear He-trimer (in ref.28.) while other authors (ref 27.) obtained a positive value for the same term. 

For water molecules, there is a three-body term of the form... 

For atoms with sp3 hybridization (e.g., Si in Si02), the parameter cosdg = —1/3 corresponds to the angle in the tetrahedron. The three-body term with an angular dependence is also used in the Stillinger-Weber potential [63] for systems with the diamond lattice... 

In case of complexes consisting of more than two constituents it has been shown that three-body terms are of significant magnitudes 120,121) i.e. ... 

In principle, the most general representation of the matrix elements of equation (68) includes a many-body expansion with one-, two-, and three-body terms. However, only two conditions [(a) and (b) in equations (66), (67)] may be used in the determination of the one- and two-body energy terms. If more conditions were to exist, the dimension of the matrix necessary to represent the adiabatic potential would be larger.16 In what follows we shall argue that the definition of V12 must depend on the type of conical intersection, namely on whether its locus is finite or infinite in extent (Section II.C). For the cases of HzO and 03 [equations (66), (67)] the conical intersection occurs along the C projection line, but only for finite values of the molecular perimeter (for H20, there is also a simple intersection at the O + H2 asymptotic channel, which is avoided for finite values of the OH distances117). For example, one gets for the linear dissociation of 03... 

The most general many-body expansion of the potential for the ground doublet state of an s3 system may, in principle, be obtained by adding a three-body term to the diagonal and nondiagonal elements of the hamiltonian matrix (75), say K(31), and V ... 

Simulations of the liquid water properties have been the subject of many papers, see Ref. (374) for a review. Recently a two-body potential for the water dimer was computed by SAPT(DFT)375. Its accuracy was checked375 by comparison with the experimental second virial coefficients at various temperatures. As shown on Figure 1-16, the agreement between the theory and experiment is excellent. Given an accurate pair potential, and three-body terms computed by SAPT376, simulations of the radial 0-0, 0-H, and H-H distribution functions could be... 

---