Integral equations hard-sphere fluid models

In Sec. 3 our presentation is focused on the most important results obtained by different authors in the framework of the rephca Ornstein-Zernike (ROZ) integral equations and by simulations of simple fluids in microporous matrices. For illustrative purposes, we discuss some original results obtained recently in our laboratory. Those allow us to show the application of the ROZ equations to the structure and thermodynamics of fluids adsorbed in disordered porous media. In particular, we present a solution of the ROZ equations for a hard sphere mixture that is highly asymmetric by size, adsorbed in a matrix of hard spheres. This example is relevant in describing the structure of colloidal dispersions in a disordered microporous medium. On the other hand, we present some of the results for the adsorption of a hard sphere fluid in a disordered medium of spherical permeable membranes. The theory developed for the description of this model agrees well with computer simulation data. Finally, in this section we demonstrate the applications of the ROZ theory and present simulation data for adsorption of a hard sphere fluid in a matrix of short chain molecules. This example serves to show the relevance of the theory of Wertheim to chemical association for a set of problems focused on adsorption of fluids and mixtures in disordered microporous matrices prepared by polymerization of species. 

The theories developed for calculating the oscillatory force are based on modeling by means of the integral equations of statistical mechanics [449-453] or numerical simulations [454-457]. As a rule, these approaches are related to complicated theoretical expressions or numerical procedures, in contrast with the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, one of its main advantages being its simplicity [36], To overcome this difficulty, some relatively simple semiempirical expressions have been proposed [458,459] on the basis of fits of theoretical results for hard-sphere fluids. 

The present chapter is organized as follows. We focus first on a simple model of a nonuniform associating fluid with spherically symmetric associative forces between species. This model serves us to demonstrate the application of so-called first-order (singlet) and second-order (pair) integral equations for the density profile. Some examples of the solution of these equations for associating fluids in contact with structureless and crystalline solid surfaces are presented. Then we discuss one version of the density functional theory for a model of associating hard spheres. All aforementioned issues are discussed in Sec. II. 

Integral equation methods provide another approach, but their use is limited to potential models that are usually too simple for engineering use and are moreover numerically difficult to solve. They are useful in providing equations of state for certain simple reference fluids (e.g., hard spheres, dipolar hard spheres, charged hard spheres) that can then be used in the perturbation theories or density functional theories. 

The so-called product reactant Ornstein-Zernike approach (PROZA) for these systems was developed by Kalyuzhnyi, Stell, Blum, and others [46-54], The theory is based on Wertheim s multidensity Ornstein-Zernike (WOZ) integral equation formalism [55] and yields the monomer-monomer pair correlation functions, from which the thermodynamic properties of the model fluid can be obtained. Based on the MSA closure an analytical theory has been developed which yields good agreement with computer simulations for short polyelectrolyte chains [44, 56], The theory has been recently compared with experimental data for the osmotic pressure by Zhang and coworkers [57], In the present paper we also show some preliminary results for an extension of this model in which the solvent is now treated explicitly as a separate species. In this first calculation the solvent molecules are modelled as two fused charged hard spheres of unequal radii as shown in Fig. 1 [45],... 

Application of the GDI method to the coexistence lines requires establishment of a coexistence datum on each. A point on the vapor-liquid line can be determined by a GE simulation. At high temperature the model behaves as a system of hard spheres, and the liquid-solid coexistence line approaches the fluid-solid transition for hard spheres, which is known [76,77]. Integration of liquid-solid coexistence from the hard-sphere transition proceeds much as described in Section III.C.l for the Lennard-Jones example. The limiting behavior (fi - 0) finds that /IP is well behaved and smoothly approaches the hard-sphere value [76,77] of 11.686 at f = 0 (unlike the LJ case, we need not work with j81/2). Thus the appropriate governing equation for the GDI procedure is... 

---