Adsorption of a Binary Mixture of Adhesive Fluids in Planar Pores: A Monte Carlo Study

Grand canonical Monte Carlo simulation is used to study the adsorption of a binary sticky hard-sphere fluid mixture in planar pores. The wall-component 1 and wall-component 2 contact densities are determined to calculate the pressure as a function of the composition of the mixture and the separation between the walls. From these data dependence of the solvation force between the plates on pore width is estimated. The simulation results are compared with the predictions of the Percus-Yevick approximation for planar pores. The density profiles of particular components show interesting shapes stemming from the interplay between the steric effects and the competitive adhesion among all possible species pairs. It is shown that narrowing of the pore causes selective partitioning of individual components of the mixture between the bulk phase and the interior of the pore. The agreement between the two methods is better at wider pores and for the component comprised of weakly adhesive particles.

Analysis of the validity of perturbation density functional theory: Based on extensive simulation for simple fluid at supercritical and subcritical temperature under various external potentials

Because of the scarcity of available simulation data for confined hard-core attractive Yukawa model fluid, extensive Monte Carlo (MC) simulation research for this fluid under the influence of various external potentials were carried out. The present MC simulation results were employed to test a performance of the third-order perturbation density functional theory (DFT) based on a high order direct correlation function (DCF) [S. Zhou and E. Ruckenstein, Phys. Rev. E. 61, 2704 (2000)]. It was found that the present perturbation DFT formalism is soundly structured only if the imported second-order DCF is reliable. In this case, the accuracy of the results can be satisfactory or even very high for various types of external potentials. Further, the associated adjustable parameter can be universal, i.e., independent of the particular external field responsible for the generation of a nonuniform density profile. Dependence of both the maintenance of the reliability of the formalism and holding of the universality of the adjustable parameter on the accuracy of the imported bulk second-order DCF can be strengthened by the large difference between the external field investigated and that caused by a single hard wall used for specification of the adjustable parameter. In case the gaseous density in the subcritical region is below the coexistence density, an excellent performance of the present formalism is observed even for the mean spherical approximation's second-order DCF as an input. This advantageous property, combined with the fact that the present formalism needs only the second-order DCF of fluid at the coexistence state as an input, enables the present formalism to be a very good theoretical tool for the investigations of wetting and prewetting transitions.