The dispersion term

The concept of dispersion was introduced by London (1930), using by far simpler arguments based on the application of the perturbation theory, as will be shown in the following subsection. A different but related interpretation puts the emphasis on the correlation in the motions of electrons. 

It is worth spending some words on electron correlation. 

We have already considered similar eoneepts in discussing die Pauh exclusion principle and die antisymmetry of the electronie wave funetions. Actually, die Pauli principle holds for particles bearing die same set of values for the characterizing quantum numbers, including spin. It says nodiing about two electrons with different spin. 

This fact has important consequences on the structure of the Hartree-Foek (HF) description of electrons in a molecule or in a dimer. The UF wave function and the eorrespond-ing electron distribution function take into aeeount the correlation of motions of electrons widi the same spin (there is a description of a hole in die probability, ealled a Fermi hole), but do not correlate motions of electrons of different spin (there is no the seeond component of the electron probability hole, called a Coulomb hole). 

This remark is important because almost all the calculations thus far performed to get molecular interaction energies have been based on the UF proeedure, which still remains the basic starting approach for all the ab initio ealeulations. The HF procedure gives the best definition of the molecular wave function in terms of a single antisymmetrized product of molecular orbitals (MO). To improve the HF description, one has to introduce in the caleu-lations other antisymmetrized products obtained from the basic one by replacing one or more MOs with others (replacement of occupied MOs with virtual MOs). This is a proce- 

---

Ihrig and Smith extended their study by running a regression analysis including reaction field terms, dispersion terms and various combinations of the solvent refractive index and dielectric constant. The best least squares fit between VF F and solvent parameters was found with a linear function of the reaction field term and the dispersion term. The reaction field term was found to be approximately three times as important as the dispersion term and the coefficients of the terms were opposite in sign. 

Here, is the distance between atoms i andj, C(/ is a dispersion coefficient for atoms i andj, which can be calculated directly from tabulated properties of the individual atoms, and /dampF y) is a damping function to avoid unphysical behavior of the dispersion term for small distances. The only empirical parameter in this expression is S, a scaling factor that is applied uniformly to all pairs of atoms. In applications of DFT-D, this scaling factor has been estimated separately for each functional of interest by optimizing its value with respect to collections of molecular complexes in which dispersion interactions are important. There are no fundamental barriers to applying the ideas of DFT-D within plane-wave DFT calculations. In the work by Neumann and Perrin mentioned above, they showed that adding dispersion corrections to forces... 

The expression of the three-body exchange energy was proposed in . For the four-body exchange energy we used a similar term. The analytical form of the dispersion terms was taken from perturbation theory up to fourth order. For the two-body dispersion energy, the dipole-dipole (r ), dipole-quadrupole... 

Throughout this summary we have neglected the effect of dispersion on the overall transport of mass and heat. This is due to the fact that if dispersion is included, dispersion tensors must be determined before the equation can be solved. This can be done by solving the appropriate transport equation within a unit cell. Because a unit cell cannot be defined in most reinforcements used in polymer matrix composites, however, dispersion tensors cannot be accurately determined, so we have left dispersion effects out of our equations. In general, we anticipate dispersion to play a minor role in the IP, AP, and RTM processes. This assumption can be checked, however, by evaluating the dispersion terms using an approach similar to [16] where experiments and correlations are used to determine the importance of dispersion. 

ST. D. Sokolov (Moscow) By which method was the dispersion term of the energy of the hydrogen bond calculated If it was calculated by London s formula then, first, it should be remembered that it is not applicable near the equilibrium state and, second, that in reality the dispersion interaction (the correlation effects) does not enter additively the energy of the electronic shell... 

Bandyopadhyay and Yashonath (31), in an extension of their work on MD studies of noble gas diffusion, presented MD results for methane diffusion in NaY and NaCaA zeolites. The zeolite models were the same as those used in the noble gas simulations (13, 15, 17, 18, 20, 28, 29) and the zeolite lattice was held rigid. The methane molecule was approximated as a single interaction center and the guest-host potential parameters were calculated from data of Bezus et al. (49) (for the dispersive term) and by setting the force on a pair of atoms equal to zero at the sum of their van der Waals radii (for the repulsive term). Simulations were run for 600 ps with a time step of 10 fs. 

For line shape calculations, purely exponential dipole models, and exponential models combined with a dispersion term, Dn/R 1, have been used, see Chapter 4 for details. The variation of the line shape with these dipole models has been studied. In all cases, the purely exponential dipole model gives inferior results when compared with the exp-7 models the influence of the dispersion term, while small, is nevertheless significant in line shape analyses moment analyses, in contrast, have reportedly not been able to demonstrate the significance of the dispersion term. It would seem that the accuracy of quantum line shape calculations of the absorption by rare gas pairs has reached the point where further progress must await more accurate experiments at low gas densities and over a wider range of frequencies and temperatures. 

For a complete description of the solvatochromic shift the dispersion term should be added, but in many cases all terms in f(n2) can be neglected if the... 

In deriving the material balance equations, the dispersed plug flow model will first be used to obtain the general form but, in the numerical calculations, the dispersion term will be omitted for simplicity. As used previously throughout, the basis for the material balances will be unit volume of the whole reactor space, i.e. gas plus liquid plus solids. Thus in the equations below, for the transfer of reactant A kLa is the volumetric mass transfer coefficient for gas-liquid transfer, and k,as is the volumetric mass transfer coefficient for liquid-solid transfer. 

Table 4.4.1 lists typical values of dipole moments and polarizabilities of some simple molecules and the three coefficients of the r 6 term at 300 K. Except for H2O which is small and highly polar, the dispersion term dominates the long-range energy. The induction term is always the least significant. 

The third component to the electrostatic interaction is caused by the motion of the electron cloud, which creates an oscillating field. It couples to the oscillating field of the neighboring molecules, which gives an attractive contribution to the total energy. This should be obvious, from the following example. Consider the frozen electron distribution of a nonpolar molecule (e.g., a noble gas). The instantaneous distribution possesses a dipole moment, which for the same reason as described above, induces a dipole in neighboring molecules, which in turn act on the first molecule, etc. This contribution is denoted the dispersion term or the London term [9]. Note that this contribution is only approximately pairwise additive. 

The density susceptibilities of the monomers can be expanded in terms of a single set of atomic orbitals, making of aA(r, r ) a two index quantity and thus greatly simplifying calculations of the dispersion term for large systems. See Section 7 for a more detailed discussion of this technique applied to the calculations of the dispersion energy. 

Thole, B.T. and Duijnen P.Th. van, The direct reaction field hamiltonian analysis of the dispersion term and application to the water dimer. Chem.Phys. (1982) 71 211-220. 

Figure 5a shows CD spectra of tartaric acid, which has an absorption i the short wavelength region and thus is prone to suffer from dispersion effect as compared with transition metal complexes. Two solution spectra in solvent of different polarity, water and dioxane, are similar to each other, but the CD C a nujol mull is quite different from that in solution. A KBr disc prepared t avoid dispersion effects gave a solid-state tartaric acid spectrum similar to thi in solution (Fig. 5b). Thus the difference between the nujol mull CD and solutid CD is not due to the different molecular conformation or intermolecular intera tion in the two phases. Most likely, it is due to the dispersion effect in the cas of the nujol mull form. Many nujol mull CD spectra of organic compound have been reported recently, but most of them appear to suffer from substanth dispersion effects. It is to be noted that the dispersion terms for molecules C...

Figure 3. Plots of the dispersion function for the equilibrium He dimer from Ref. 36. The electron coordinates, xj and X2 are defined in fig. (a). The accurate dispersion function (calculated with Gaussian geminals) is shown in plot (b). If the dispersion term is reduced to the dipole-dipole component, one obtains plot (c). If every electron is allowed to use the complete dimer basis set as well as the bond functions one obtains plot (d).

The multi-center form of the dispersion term is more in the spirit of the other terms, and despite its simplicity proved to perform better in some applications, e.g. in calculations of tunneling splittings between nonsuperimposable forms of the water dimer [56,57]. 

In the previous section, stability criteria were obtained for gas-hquid bubble columns, gas-solid fluidized beds, liquid-sohd fluidized beds, and three-phase fluidized beds. Before we begin the review of previous work, let us summarize the parameters that are important for the fluid mechanical description of multiphase systems. The first and foremost is the dispersion coefficient. During the derivation of equations of continuity and motion for multiphase turbulent dispersions, correlation terms such as esv appeared [Eqs. (3) and (10)]. These terms were modeled according to the Boussinesq hypothesis [Eq. (4)], and thus the dispersion coefficients for the sohd phase and hquid phase appear in the final forms of equation of continuity and motion [Eqs. (5), (6), (14), and (15)]. However, for the creeping flow regime, the dispersion term is obviously not important. 

The dispersion in equation of motion (12) appears in the second-order derivative, while Batchelor (1988) has considered mobility to derive the dispersion term, which is a first-order derivative. 

Batchelor (1988) has derived the particle phase governing equation directly without any need for writing a hquid phase momentum equation. This enabled elimination of the pressure term. The dispersion term was not included in the equation of continuity. However, it was included in the equation of motion. A revisit to the derivation of governing equations is expected to be useful in future development. 

This criterion involves the assumption that the gas phase stress terms are negligible. This assumption may not be valid in case of solid-liquid fluidized beds or liquid-liquid dispersions. In this case, the criterion is of the same form as Eq. (172), with different definitions of the parameters Ml, M2, and M3, which are given in Table VII. Table VII also gives the parameters of the criterion when the dispersion terms are not included in the continuity equations of both the phases. 

The interaction between water molecules and silica substrate is described in the framework of the PN-TrAZ model [15] which has proven to model successfiilly the adsorption of simple adsorbates on various zeolites [20]. In this model, the pair potential decomposes in two parts a repulsion term Ae" due to electronic clouds and the attractive dispersion terms. The repulsive parameters (A,b) for silica atoms (Si, O, H) are those obtained from studies of adsorption of simple gazes on various zeolites [20] and mesoporous glass [21]. Those for water oxygen are chosen to fit the repulsive part of Lennard-Jones from SPC model in the range around equilibriiun distance, and those for water hydrogen are taken equal to the parameters for surface hydrogen of vycor. The cross repulsive parameters A and b are obtained by Bohm and Ahlrichs [22] combination mles. The dispersion terms are calculated from polarizabilities and effective niunber of electron Neff according to the PN-TrAZ model up to order r °. Values are listed in table 1. 

Note that is called the axial-dispersion coejficient, and that the dispersion term of the molar flow rate is formulated by analogy to molecular diffusion. Pick s First Law states that the flux of species A (moles/area/time) can be formulated as ... 

For these, empirical values are usually not available. One calculates first the wave-function of the excited state. A, using SCF-CI techniques. The charge distribution of this state is obtained, and expressed as a sum of multipolar, multicentric terms, up to quadrupole. From this expression, the electrostatic interactions are calculated as in the ground state. The polarization term is harder to calculate, it is often taken to be equal to the ground-state term. The dispersion term may be calculated approximately by considering the interaction between two dipolar distributions due to the excited states of A and D. This approximation is valid when the first excited state... 

---