Interaction site fluids approximate theories

In this section, we review some of the important formal results in the statistical mechanics of interaction site fluids. These results provide the basis for many of the approximate theories that will be described in Section III, and the calculation of correlation functions to describe the microscopic structure of fluids. We begin with a short review of the theory of the pair correlation function based upon cluster expansions. Although this material is featured in a number of other review articles, we have chosen to include a short account here so that the present article can be reasonably self-contained. Cluster expansion techniques have played an important part in the development of theories of interaction site fluids, and in order to fully grasp the significance of these developments, it is necessary to make contact with the results derived earlier for simple fluids. We will first describe the general cluster expansion theory for fluids, which is directly applicable to rigid nonspherical molecules by a simple addition of orientational coordinates. Next we will focus on the site-site correlation functions and describe the interaction site cluster expansion. After this, we review the calculation of thermodynamic properties from the correlation functions, and then we consider the calculation of the dielectric constant and the Kirkwood orientational correlation parameters. 

The formal treatment describing the extension of perturbation theory to polar interaction site systems and to other situations where the perturbative forces are structure-determining is available and has been applied with some success to some simple models of polar diatomics. More quantitative comparisons with computer simulations need to be made. Of course, qualitatively accurate information about the structure of polar interaction site fluids has been available for some time through solutions of the SSOZ-HNC equations. However, this approach does not seem to be useful in the context of thermodynamics. Rather little attention has been paid to polarizable molecules, although these can be treated within the context of the interaction site formalism (see, for example, Chandler and, more recently, Sprik and Klein ). Although the formal treatment of the dielectric constant within the interaction site formalism is now well established, no quantitative approximations seem to emerge from any of the theories available. 

Theories based on the solution to integral equations for the pair correlation fiinctions are now well developed and widely employed in numerical and analytic studies of simple fluids [6]. Furtlier improvements for simple fluids would require better approximations for the bridge fiinctions B(r). It has been suggested that these fiinctions can be scaled to the same fiinctional fomi for different potentials. The extension of integral equation theories to molecular fluids was first accomplished by Chandler and Andersen [30] through the introduction of the site-site direct correlation fiinction c r) between atoms in each molecule and a site-site Omstein-Zemike relation called the reference interaction site... 

In this chapter we have described a theory for dynamics of polyatomic fluids based on the memory-function formalism and on the interaction-site representation of molecular liquids. Approximation schemes for memory functions appearing in the generalized Langevin equation have been developed by assuming an exponential form for memory functions and by employing the mode-coupling approach. Numerical results were presented for longitudinal current spectra of a model diatomic liquid and water, and it has been discussed how the results can be interpreted in... 

If go(r), g CrX and g (r) are known exactly, then all three routes should yield the same pressure. Since liquid state integral equation theories are approximate descriptions of pair correlation functions, and not of the effective Hamiltonian or partition function, it is well known that they are thermodynamically inconsistent [5]. This is understandable since each route is sensitive to different parts of the radial distribution function. In particular, g(r) in polymer fluids is controlled at large distance by the correlation hole which scales with the radius of gyration or /N. Thus it is perhaps surprising that the hard core equation-of-state computed from PRISM theory was recently found by Yethiraj et aL [38,39] to become more thermodynamically inconsistent as N increases from the diatomic to polyethylene. The uncertainty in the pressure is manifested in Fig. 7 where the insert shows the equation-of-state of polyethylene computed [38] from PRISM theory for hard core interactions between sites. In this calculation, the hard core diameter d was fixed at 3.90 A in order to maintain agreement with the experimental structure factor in Fig. 5. 

---