Pairwise

To first order, the dispersion (a-a) interaction is independent of the structure in a condensed medium and should be approximately pairwise additive. Qualitatively, this is because the dispersion interaction results from a small perturbation of electronic motions so that many such perturbations can add without serious mutual interaction. Because of this simplification and its ubiquity in colloid and surface science, dispersion forces have received the most significant attention in the past half-century. The way dispersion forces lead to long-range interactions is discussed in Section VI-3 below. Before we present this discussion, it is useful to recast the key equations in cgs/esu units and SI units in Tables VI-2 and VI-3. 

In this section we consider electromagnetic dispersion forces between macroscopic objects. There are two approaches to this problem in the first, microscopic model, one assumes pairwise additivity of the dispersion attraction between molecules from Eq. VI-15. This is best for surfaces that are near one another. The macroscopic approach considers the objects as continuous media having a dielectric response to electromagnetic radiation that can be measured through spectroscopic evaluation of the material. In this analysis, the retardation of the electromagnetic response from surfaces that are not in close proximity can be addressed. A more detailed derivation of these expressions is given in references such as the treatise by Russel et al. [3] here we limit ourselves to a brief physical description of the phenomenon. 

The total interaction between two slabs of infinite extent and depth can be obtained by a summation over all atom-atom interactions if pairwise additivity of forces can be assumed. While definitely not exact for a condensed phase, this conventional approach is quite useful for many purposes [1,3]. This summation, expressed as an integral, has been done by de Boer [8] using the simple dispersion formula, Eq. VI-15, and following the nomenclature in Eq. VI-19 ... 

There are tliree important varieties of long-range forces electrostatic, induction and dispersion. Electrostatic forces are due to classical Coulombic interactions between the static charge distributions of the two molecules. They are strictly pairwise additive, highly anisotropic, and can be either repulsive or attractive. 

A direct and transparent derivation of the second virial coefficient follows from the canonical ensemble. To make the notation and argument simpler, we first assume pairwise additivity of the total potential with no angular contribution. The extension to angularly-mdependent non-pairwise additive potentials is straightforward. The total potential... 

The correlation functions provide an alternate route to the equilibrium properties of classical fluids. In particular, the two-particle correlation fimction of a system with a pairwise additive potential detemrines all of its themiodynamic properties. It also detemrines the compressibility of systems witir even more complex tliree-body and higher-order interactions. The pair correlation fiinctions are easier to approximate than the PFs to which they are related they can also be obtained, in principle, from x-ray or neutron diffraction experiments. This provides a useful perspective of fluid stmcture, and enables Hamiltonian models and approximations for the equilibrium stmcture of fluids and solutions to be tested by direct comparison with the experimentally detennined correlation fiinctions. We discuss the basic relations for the correlation fiinctions in the canonical and grand canonical ensembles before considering applications to model systems. 

This leads to the third virial coefficient for hard spheres. In general, the nth virial coefficient of pairwise additive potentials is related to the coefficient7) in the expansion of g(r), except for Coulombic systems for which the virial coefficients diverge and special teclmiques are necessary to resiim the series. 

The thennodynamic properties of a fluid can be calculated from the two-, tln-ee- and higher-order correlation fiinctions. Fortunately, only the two-body correlation fiinctions are required for systems with pairwise additive potentials, which means that for such systems we need only a theory at the level of the two-particle correlations. The average value of the total energy... 

For a pairwise additive potential, each temi in the sum of pair potentials gives the same result in the above expression and there are N(N - l)/2 such temis. It follows that... 

The integral is easily simplified for a pairwise additive system, and one finds... 

Successive n and n + 1 particle density fiinctions of fluids with pairwise additive potentials are related by the Yvon-Bom-Green (YBG) hierarchy [6]... 

Kirkwood derived an analogous equation that also relates two- and tlnee-particle correlation fiinctions but an approximation is necessary to uncouple them. The superposition approximation mentioned earlier is one such approximation, but unfortunately it is not very accurate. It is equivalent to the assumption that the potential of average force of tlnee or more particles is pairwise additive, which is not the case even if the total potential is pair decomposable. The YBG equation for n = 1, however, is a convenient starting point for perturbation theories of inliomogeneous fluids in an external field. 

Assuming the perturbing potential is pairwise additive, an argument virtnally identical to the calcnlation of = shows that... 

The relationship between g(r) and the interparticle potential energy is most easily seen if we assume that the interparticle energy can be factorized into pairwise additive potentials as... 

Here we assumed pairwise additivity V= and defined w(r) = r(dv(r)/dr). Also easily derived in the... 

The intennolecular forces between water molecules are strongly non-additive. It is not realistic to expect any pair potential to reproduce the properties of both the water dimer and the larger clusters, let alone liquid water. There has therefore been a great deal of work on developing potential models with explicit pairwise-additive and nonadditive parts [44, 50, 51]. It appears that, when this is done, the energy of the larger clusters and ice has a nonadditive contribution of about 30%. 

Similarly, van der Waals forces operate between any two colloidal particles in suspension. In the 1930s, predictions for these interactions were obtained from the pairwise addition of molecular interactions between two particles [38]. The interaction between two identical spheres is given by... 

Tables of single substance surface tensions y can be built up tlirough tire measurement of contact angles [134] a few examples are collected in table C2.14.1. These can be combined pairwise according to ...

More recently Andrews and Juzeliunas [6, 7] developed a unified tlieory that embraces botli radiationless (Forster) and long-range radiative energy transfer. In otlier words tliis tlieory is valid over tire whole span of distances ranging from tliose which characterize molecular stmcture (nanometres) up to cosmic distances. It also addresses tire intennediate range where neitlier tire radiative nor tire Forster mechanism is fully valid. Below is tlieir expression for tire rate of pairwise energy transfer w from donor to acceptor, applicable to transfer in systems where tire donor and acceptor are embedded in a transparent medium of refractive index ... 

Knowledge of tire pairwise energy trairsfer rates fonrrs a basis for finding tire average rate of energy trairsfer in air ensemble of molecules. To tlris end, a system of master equations is commonly employed [15,16 aird 17]. Then, tire probability, to find excitation on molecule cair be calculated as ... 

Here t. is the intrinsic lifetime of tire excitation residing on molecule (i.e. tire fluorescence lifetime one would observe for tire isolated molecule), is tire pairwise energy transfer rate and F. is tire rate of excitation of tire molecule by the external source (tire photon flux multiplied by tire absorjDtion cross section). The master equation system (C3.4.4) allows one to calculate tire complete dynamics of energy migration between all molecules in an ensemble, but tire computation can become quite complicated if tire number of molecules is large. Moreover, it is commonly tire case that tire ensemble contains molecules of two, tliree or more spectral types, and experimentally it is practically impossible to distinguish tire contributions of individual molecules from each spectral pool. 

Parallel molecular dynamics codes are distinguished by their methods of dividing the force evaluation workload among the processors (or nodes). The force evaluation is naturally divided into bonded terms, approximating the effects of covalent bonds and involving up to four nearby atoms, and pairwise nonbonded terms, which account for the electrostatic, dispersive, and electronic repulsion interactions between atoms that are not covalently bonded. The nonbonded forces involve interactions between all pairs of particles in the system and hence require time proportional to the square of the number of atoms. Even when neglected outside of a cutoff, nonbonded force evaluations represent the vast majority of work involved in a molecular dynamics simulation. 

In order to improve parallelism and load balancing, a hybrid force-spatial decomposition scheme was adopted in NAMD 2. Rather than decomposing the nonbonded computation into regions of space or pairwise atomic interactions, the basic unit of work was chosen to be interactions between atoms... 

Consideration of stereochem-iitry. The parity or handedness - R/S or chjirans - of a stcreoccnter can be obtained by considering the sequence of the Morgan numbers of die atoms, similarly to CIP. Then the number of pairwise interclianges is counted until the numbers arc in ascending order (see Section 2,8,5). 

---