Solid-fluid equilibrium density functional theories

An alternative approach is by the application of an approximate theory. At present, the most useful theoretical treatment for the estimation of the equilibrium properties is generally considered to be the density functional theory (DFT). This involves the derivation of the density profile, p(r), of the inhomogeneous fluid at a solid surface or within a given set of pores. Once p(r) is known, the adsorption isotherm and other thermodynamic properties, such as the energy of adsorption, can be calculated. The advantage of DFT is its speed and relative ease of calculation, but there is a risk of oversimplification through the introduction of approximate forms of the required functionals (Gubbins, 1997). 

Over the years, vapour adsorption and condensation in porous materials continue to attract a great deal of attention because of (i) the fundamental physics of low-dimension systems due to confinement and (ii) the practical applications in the field of porous solids characterisation. Particularly, the specific surface area, as in the well-known BET model [I], is obtained from an adsorbed amount of fluid that is assumed to cover uniformly the pore wall of the porous material. From a more fundamental viewpoint, the interest in studying the thickness of the adsorbed film as a function of the pressure (i.e. t = f (P/Po) the so-called t-plot) is linked to the effort in describing the capillary condensation phenomenon i.e. the gas-Fadsorbed film to liquid transition of the confined fluid. Indeed, microscopic and mesoscopic approaches underline the importance of the stability of such a film on the thermodynamical equilibrium of the confined fluid [2-3], In simple pore geometry (slit or cylinder), numerous simulation works and theoretical studies (mainly Density Functional Theory) have shown that the (equilibrium) pressure for the gas/liquid phase transition in pores greater than 8 nm is correctly predicted by the Kelvin equation provided the pore radius Ro is replaced by the core radius of the gas phase i.e. (Ro -1) [4]. Thirty year ago, Saam and Cole [5] proposed that the capillary condensation transition is driven by the instability of the adsorbed film at the surface of an infinite... 

The formalism of density functional theory (DFT) has received considerable attention as a way to describe the adsorption process at the fluid-solid interfece. The older approach was to treat the adsorbate as a separate, two-dimensional phase existing in equilibrium with the bulk gas phase. This model works well... 

All three areas will be addressed here. The application of classical density functional theory has led to some of the most important recent theoretical advances in SFE and these have been the subject of several authoritative review articles [10-16]. On the other hand, we know of no recent comprehensive review addressing theoretical approaches other than density functional theories (DFT) and the other two subject areas, particularly the last one, and it was this that motivated us to write this chapter. We hope that the somewhat broader coverage of molecular modeling research in SFE given in this chapter will be of benefit to researchers new to the field. We should mention that this Chapter is written from a perspective that is more strongly influenced by liquid-state statistical mechanics than by solid-state theory. The interests of the authors in the problem at hand are an outgrowth of their previous work on phase equilibrium in fluids and fluid mixtures. 

The classical models of adsorption processes like Langmuir, BET, DR or Kelvin treatments and their numerous variations and extensions, contain several uncontrolled approximations. However, the classical theories are convenient and their usage is very widespread. On the other hand, the aforementioned classical theories do not start from a well - defined molecular model, and the result is that the link between the molecular behaviour and the macroscopic properties of the systems studied are blurred. The more developed and notable descriptions of the condensed systems include lattice models [408] which are solved by means of the mean - field or other non-classical techniques [409]. The virial formalism of low -pressure adsorption discussed above, integral equation method and perturbation theory are also useful approaches. However, the state of the art technique is the density functional theory (DFT) introduced by Evans [410] and Tarazona [411]. The DFT method enables calculating the equilibrium density profile, p (r), of the fluid which is in contact with the solid phase. The main idea of the DFT approach is that the free energy of inhomogeneous fluid which is a function of p (r), can be... 

Figure 3. The equilibrium vapor and liquid densities of an associating fluid with one square-well bonding site. The circles are data from RCMC-Gibbs ensemble simulations, and the lines are calculations from three different implementations of a theory for associating fluids. The solid line uses exact values of the reference fluid radial distribution function the dashed and long dashed-short dashed lines use the WCA and modified WCA approximations to the radial distribution function, respectively. (Reprinted with permission from Muller et al. [43]. Copyright 1995 American Institute of Physics.)...

Over the p t several years we and our collaborators have pursued a continuous space liquid state approach to developing a computationally convenient microscopic theory of the equilibrium properties of polymeric systems. Integral equations method [5-7], now widely employed to understand structure, thermodynamics and phase transitions in atomic, colloidal, and small molecule fluids, have been generalized to treat macromolecular materials. The purpose of this paper is to provide the first comprehensive review of this work referred to collectively as Polymer Reference Interaction Site Model (PRISM) theory. A few new results on polymer alloys are also presented. Besides providing a unified description of the equilibrium properties of the polymer liquid phase, the integral equation approach can be combined with density functional and/or other methods to treat a variety of inhomogeneous fluid and solid problems. 

The simultaneous calculation of the amount adsorbed and the distribution of adsorbate within the pores of a reconstructed solid can be impl ented afier the i lication of a Draisity Functional Theory (DFT) mean field model, which is particularly suited for on-lattice simulations on di tised structures [31]. In a lattice model the spatial distribution of adsorbate can be described at each site by the local density fimction. The equilibrium density profile for a given matrix realization and chemical potential ft (xin then be determined by minimising the grand potential 0 p with respect to the fluid density on each lattice site, / (x), leading to ... 

---