Simple liquids, density functional theory

The first way has been followed in what has become known as Car-Parrinello molecular dynamics (CPMD) (9). A solute and 60-90 solvent molecules are considered to represent the system, and the QM calculations are performed with density functionals, usually of generalised gradient approximation type (GGA), such as the Becke-Lee-Young-Parr (BLYP) (10) or the Perdew-Burke-Enzerhofer (PBE) (11,12) functionals. It is clear that the semiempirical character of concurrent density functional theory (DFT) methods and the use of these simple functionals imply a number of error sources and do not really provide a method-inherent control procedure to test the reliability of results. Recently it has been shown that these functionals even do not enable a correct description of the solvent water itself, as at ambient temperature they will describe water not as liquid but as supercooled system... 

The principal tools have been density functional theory and computer simulation, especially grand canonical Monte Carlo and molecular dynamics [17-19]. Typical phase diagrams for a simple Lennard-Jones fluid and for a binary mixture of Lennard-Jones fluids confined within cylindrical pores of various diameters are shown in Figs. 9 and 10, respectively. Also shown in Fig. 10 is the vapor-liquid phase diagram for the bulk fluid (i.e., a pore of infinite radius). In these examples, the walls are inert and exert only weak forces on the molecules, which themselves interact weakly. Nevertheless,... 

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... 

Over the last two decades the exploration of microscopic processes at interfaces has advanced at a rapid pace. With the active use of computer simulations and density functional theory the theory of liquid/vapor, liquid/liquid and vacuum/crystal interfaces has progressed from a simple phenomenological treatment to sophisticated ah initio calculations of their electronic, structural and dynamic properties [1], However, for the case of liquid/crystal interfaces progress has been achieved only in understanding the simplest density profiles, while the mechanism of formation of solid/liquid interfaces, emergence of interfacial excess stress and the anisotropy of interfacial free energy are not yet completely established [2],... 

The validity of the viscoelastic model (5.32) has been tested against experimental and molecular dynamics simulation results [26, 27, 28]. The detailed comparison has established that the viscoelastic model works remarkably well for wavenumbers k km, where km denotes the first peak position of the static structure factor S k). However, it has also been found that the situation is not so satisfactory for smaller wavenumbers, where the viscoelastic model is shown in some circumstances to yield even qualitatively incorrect results. This failure was attributed to the fact that the single relaxation time model (5.31) cannot describe both the short-time behavior of the memory function, dominated by the so-called binary collisions, and in particular the intermediate and long-time behavior where in the liquid range additional slow processes play an important role (see the next subsection). It is obvious that these conclusions demand a more rigorous consideration of the memory function, which lead to the development of the modern version of the kinetic theory. Nevertheless, the viscoelastic model provides a rather satisfactory account of the main features of microscopic collective density fluctuations in simple liquids at relatively large wavenumbers, and its value should not be undervalued. 

Curtin, W.A., and Ashcroft, N.W. Density-functional theory and freezing of simple liquids. 1986. Phys. Rev. Lett. 56 2775. 

In the case of liquid films that contain charged colloid particles (micelles), the oscillatory period. Ah, is considerably greater than the particle diameter [437,438], In this case, the theoretical prediction of Ah demands the use of density functional theory calculations and/or Monte Carlo simulations [456,512], However, the theory, simulations, and experiments showed that a simple inverse-cubic-root relation. Ah = exists between Ah and the bulk number concentration of micelles (particles), c , [437,438,498,512-514],... 

---