Accurate and local formulation for thermodynamic properties directly from integral equation method

Free energy density functional for adsorption of fluids in nanopores

A classical free energy density functional, which is isomorphic to a usual effective hard sphere model + mean field approximation for tail contribution, is proposed for treatment of real fluids in inhomogeneous states. In the framework of the classical density functional theory (DFT), the present functional is applied to two representative model fluids, namely, a Lennard-Jones fluid and a hard core attractive Yukawa fluid, subject to influence of various external fields. A comprehensive comparison with simulation results and a detailed analysis show that the present functional holds simultaneously all of the desirable properties inherent in an excellent functional, such as high accuracy, computational simplicity, consistency with a hard wall sum rule, nonrecourse to use of adjustable parameter(s) and weighted densities, reproduction of bulk second-order direct correlation function (DCF) in bulk limit, and applicability to subcritical fluid phenomena.

A theoretical investigation on the honeycomb potential fluid

A local self-consistent Ornstein-Zernike (OZ) integral equation theory (IET) is proposed to provide a rapid route for obtaining thermodynamic and structural information for any thermodynamically stable or metastable state points in the bulk phase diagram without recourse to traditional thermodynamic integration, and extensive NVT-Monte Carlo simulations are performed on a recently proposed honeycomb potential in three dimensions to test the theory's reliability. The simulated quantities include radial distribution function (rdf) and excess internal energy, pressure, excess chemical potential, and excess Helmholtz free energy. It is demonstrated that (i) the theory reproduces the rdf very satisfactorily only if the bulk state does not enter deep into a two phases coexistence region; (ii) the excess internal energy is the only one of the four thermodynamic quantities investigated amenable to the most accurate prediction by the present theory, and the simulated pressure is somewhat overestimated by the theoretical calculations, but the deviation tends to vanish along with rising of the temperature; (iii) using the structural functions from the present local self-consistent OZ IET, a previously derived local expression, due to the present author, achieves even a higher accuracy in calculating for the excess chemical potential than the exact virial pressure formula for the pressure, and the resulting excess Helmholtz free energy is in surprisingly same with the simulation results due to offset of the errors. Based on the above observations, it is suggested that it may be a good procedure to integrate the theoretical excess internal energy along the isochors to get the excess Helmholtz free energy, which is then fitted to a polynomial to be used for calculation of all of other thermodynamic quantities in the framework of the OZ IET.

Reformulation of liquid perturbation theory for low temperatures

A perturbation strategy is proposed which allows liquid perturbation methodology to be applied to extremely low temperature situations for which the available liquid integral equation theory and traditional thermodynamic perturbation theory fail. The possibility of avoiding the low temperature problem of the thermodynamic perturbation theory not only is of relevance to the investigation of complex fluids, but also may be useful for reformulation of other liquid theories to achieve higher accuracy and avoid the respective low temperature problems.

Local structures of fluid with discrete spherical potential: Theory and grand canonical ensemble Monte Carlo simulation

Grand canonical Monte Carlo simulation and theoretical calculations based on Ornstein-Zernike (OZ) integral equation and third order+second order perturbation density functional theory (DFT) are performed to study a system of spherical particles interacting through a core-softened (CS) potential combining a repulsive square soft core and an attractive square well. Both theoretical predictions and simulation results reveal peculiar homogeneous and inhomogeneous local structures originating from the discontinuous nature of the CS potential. The bulk radial distribution function displays discontinuities at the distances coinciding with the ranges of the successive repulsive and attractive parts in the CS potential function. The density profiles of confined CS fluid show the shapes arising from the complex interplay among the steric effects and the competition between the repulsive and attractive parts of the CS potential. Satisfactory agreement between the theoretical results and simulation data leads to the following conclusions: (i) a modified hypernetted chain approximation combined with a hard sphere bridge function, which has been recently proposed by one of the authors of this study, is sufficiently reliable for the structural studies of CS fluid, and (ii) the third order+second order perturbation DFT, which has proven successful for the study of inhomogeneous structure of model fluids with continuous intermolecular potential function, posses a high adaptability to be applied for various types of interaction potentials and performs well also in the case of discontinuous CS model.