Extending the simple weighted density approximation for a hard-sphere fluid to a Lennard-Jones fluid II. Application

A simple weighted density approximation (SWDA) was extended to nonuniform Lennard-Jones fluids by following the spirit of a partitioned density function theory [S. Zhou, Phys. Rev. E 68 (2003) 061201] and mapping the hard-core part onto an effective hard-sphere fluid whose higher order terms beyond the second order of the functional perturbation expansion are treated by the SWDA. The resultant DFT formalism performs well for Lennard-Jones fluids under the influence of diverse external fields. With the present DFT formalism, we investigate in detail the structure and adsorption properties of a low-density LJ gas in a spherical cavity with a wall consisting of hard-sphere or LJ particles. It was found that when the cavity wall exerts an attractive external potential on the LJ particles in the cavity, the excess adsorption decreases as the temperature increases, while when the cavity wall exerts a hard repulsive external potential on the LJ particles in the cavity, the excess adsorption increases as the temperature increases.     

Application of Lagrangian theorem-based density-functional approximation free of adjustable parameters to nonhard-sphere fluid

A recently proposed parameter free version of a Lagrangian theorem-based density functional approximation (LTDFA) [S. Zhou, Phys. Lett. A 319, 279 (2003)] for hard-sphere fluid is applied to hard-core attractive Yukawa model fluid by dividing bulk second-order direct correlation function (DCF) of fluid under consideration into hard-core part and tail part. The former is treated by the parameter free version of the LTDFA, while the tail part is treated by second-order functional perturbation expansion approximation as done in a recent partitioned DFA [S. Zhou, Phys. Rev. E 68, 061201 (2003)]. Two versions of mean spherical approximation (MSA) for the bulk second-order DCF are employed as input, one is the less accurate plain MSA whose tail part of the second-order DCF is strictly independent of a density argument, the other is the more accurate inverse temperature expansion version of the MSA whose tail part is not strictly independent of the density argument. Calculational results indicate that prediction based on the plain MSA is far more accurate than that based on the inverse temperature expansion version of the MSA. The reason is considered to be that the partitioned DFA requires that the tail part is highly or completely independent of the density argument, the plain MSA, by assuming that the tail part is exactly the potential itself, embodies all of the nonlinearities into the hard-core part which can be treated satisfactorily by the parameter free version of the LTDFA. The present investigation results in a universal method for constructing DFA for nonuniform any nonhard-sphere interaction potential fluids.     

Inverse density-functional theory as an interpretive tool for measuring colloid-surface interactions in dense systems

Recent advances in optical microscopy, such as total internal reflection and confocal scanning laser techniques, now permit the direct three-dimensional tracking of large numbers of colloidal particles both near and far from interfaces. A novel application of this technology, currently being developed by one of the authors under the name of diffusing colloidal probe microscopy (DCPM), is to use colloidal particles as probes of the energetic characteristics of a surface. A major theoretical challenge in implementing DCPM is to obtain the potential energy of a single particle in the external field created by the surface, from the measured particle trajectories in a dense colloidal system. In this paper we develop an approach based on an inversion of density-functional theory (DFT), where we calculate the single-particle-surface potential from the experimentally measured equilibrium density profile in a nondilute colloidal fluid. The underlying DFT formulation is based on the recent work of Zhou and Ruckenstein [Zhou and Ruckenstein, J. Chem. Phys. 112, 8079 (2000)]. For model hard-sphere and Lennard-Jones systems, using Monte Carlo simulation to provide the "experimental" density profiles, we found that the inversion procedure reproduces the true particle-surface-potential energy to an accuracy within typical DCPM experimental limitations (approximately 0.1 kT) at low to moderate colloidal densities. The choice of DFT closures also significantly affects the accuracy.     

How to extend hard sphere density functional approximation to nonuniform nonhard sphere fluids: applicable to both subcritical and supercritical temperature regions

A methodology for the formulation of density functional approximation (DFA) for nonuniform nonhard sphere fluids is proposed by following the spirit of a partitioned density functional approximation [Zhou, Phys. Rev. E 68, 061201 (2003)] and mapping the hard core part onto an effective hard sphere whose high order part of the functional perturbation expansion is treated by existing hard sphere DFAs. The resultant density functional theory (DFT) formalism only needs a second order direct correlation function and pressure of the corresponding coexistence bulk fluid as inputs and therefore can be applicable to both supercritical and subcritical temperature cases. As an example, an adjustable parameter-free version of a recently proposed Lagrangian theorem-based DFA is imported into the present methodology; the resultant DFA is applied to Lennard-Jones fluid under the influence of external fields due to a single hard wall, two hard walls separated by a small distance, a large hard sphere, and a spherical cavity with a hard wall. By comparing theoretical predictions with previous simulation data and those recently supplied for coexistence bulk fluid situated at "dangerous" regions, it was found that the present DFA can predict subtle structure change of the density profile and therefore is the most accurate among all existing DFT approaches. A detailed discussion is given as to why so excellent DFA for nonhard sphere fluids can be drawn forth from the present methodology and how the present methodology differs from previous ones. The methodology can be universal, i.e., it can be combined with any other hard sphere DFAs to construct DFA for other nonhard sphere fluids with a repulsive core.     

Phase behavior of colloids and proteins in aqueous suspensions: Theory and computer simulations

The fluid phase behavior of colloidal suspensions with short-range attractive interactions is studied by means of Monte Carlo computer simulations and two theoretical approximations, namely, the discrete perturbation theory and the so-called self-consistent Ornstein-Zernike approximation. The suspensions are modeled as hard-core attractive Yukawa (HCAY) and Asakura-Oosawa (AO) fluids. A detailed comparison of the liquid-vapor phase diagrams obtained through different routes is presented. We confirm Noro-Frenkel's extended law of scaling according to which the properties of a short-ranged fluid at a given temperature and density are independent of the detailed form of the interaction, but just depend on the value of the second virial coefficient. By mapping the HCAY and AO fluids onto an equivalent square-well fluid of appropriate range at the critical point we show that the critical temperature as a function of the effective range is independent of the interaction potential, i.e., all curves fall in a master curve. Our findings are corroborated with recent experimental data for lysozyme proteins.     

Interfacial colloidal sedimentation equilibrium. II. Closure-based density functional theory

In Part I [R. E. Beckham and M. A. Bevan, J. Chem. Phys. 127, 164708 (2007)], results were presented for the sedimentation equilibrium of concentrated colloidal dispersions using confocal scanning laser microscopy experiments, Monte Carlo (MC) simulations, and a local density approximation perturbation theory. In this paper, we extended the modeling effort on those systems to include nonlocal density functional theory (DFT), which is capable of predicting the microstructure of the sediment at length scales comparable to the colloidal particle dimension. Specifically, we use a closure-based DFT formulation to predict interfacial colloidal sedimentation equilibrium density profiles. The colloid-colloid and colloid-surface interactions were modeled with DLVO screened electrostatic potentials using parameters taken directly from the experimental work. The DFT profiles were compared to the experimental and MC results from Part I. Good agreement was found for relatively dilute interfacial colloidal fluids, but agreement was less satisfactory as interfacial layering became more pronounced for conditions approaching the onset of interfacial crystallization. We also applied DFT in an inverse sense using the measured colloid density profile to extract the underlying colloid-surface potential; this can be thought of as a microscopic analog to the well-known procedure of using the macroscopic (coarse-grained) density profile to extract the osmotic equation of state. For the dilute interfacial fluid, the inverse DFT calculations reproduced the true colloid-surface potential to within 0.5kT at all elevations.     

Closure-based density functional theory applied to interfacial colloidal fluids

The behavior of dense colloidal fluids near surfaces can now be probed in great detail with experimental techniques like confocal microscopy. In fact, we are approaching a point where quantitative comparisons of experiment with particle-level theory, such as classical density functional theory (DFT), are appropriate. In a forward sense, we may use a known surface potential to predict a particle density distribution function from DFT; in an inverse sense, we may use an experimentally measured particle density distribution function to predict the underlying surface potential from DFT. In this paper, we tested the ability of the closure-based DFT of Zhou and Ruckenstein (J. Chem. Phys. 2000, 112, 8079-8082) to perform forward and inverse calculations on potential models commonly employed for colloidal particles and surfaces. To reduce sources of uncertainty in this initial study, Monte Carlo simulation results played the role of experimental data. The combination of Rogers-Young and modified-Verlet closures consistently performed well across the different potential models. For a reasonable range of choices of the density, temperature, and potential parameters, the inversion procedure yielded particle-surface potentials to an accuracy on the order of 0.1kT.     

Structure and adsorption of a hard-core multi-Yukawa fluid confined in a slitlike pore: grand canonical Monte Carlo simulation and density functional study

Because of the increasing interest in studying the phenomenon exhibited by charge-stabilized colloidal suspensions in confining geometry, we present a density functional theory (DFT) for a hard-core multi-Yukawa fluid. The excess Helmholtz free-energy functional is constructed by using the modified fundamental measure theory and Rosenfeld's perturbative method, in which the bulk direct correlation function is obtained from the first-order mean spherical approximation. To validate the established theory, grand canonical ensemble Monte Carlo (GCMC) simulations are carried out to determine the density profiles and surface excesses of multi-Yukawa fluid in a slitlike pore. Comparisons of the theoretical results with the GCMC data suggest that the present DFT gives very accurate density profiles and surface excesses of multi-Yukawa fluid in the slitlike pore as well as the radial distribution functions of the bulk fluid. Both the DFT and the GCMC simulations predict the depletion of the multi-Yukawa fluid near a nonattractive wall, while the mean-field theory fails to describe this depletion in some cases. Because the simple form of the direct correlation function is used, the present DFT is computationally as efficient as the mean-field theory, but reproduces the simulation data much better than the mean-field theory.     

Is perturbation DFT approach applicable to purely repulsive fluids?

A recently proposed third order + second order perturbation density functional theory (DFT) approach is tested for the validity and applicability to purely repulsive model fluids subjected to various external fields. Hard core repulsive Yukawa potential, point particle Yukawa potential, and inverse power potential are employed as sample models. Theoretical DFT results are compared with the corresponding simulation data obtained by grand canonical ensemble Monte Carlo simulation. This comparison indicates that the third order + second order perturbation DFT approach is suitable for these purely repulsive fluids only on condition of high accuracy of the imported bulk second order direct correlation function (DCF). However, in this case the origin of the successful performance somewhat differs from that observed for the mean field approximation applied to van der Waals fluids. In the present case it originates from the observation that the bulk second order DCF is strongly dependent on the density argument for the hard-core part, while for the distances exceeding the core dimension this dependence is considerably weaker.     

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.     

基于密度泛函理论研究二元排斥Yukawa流体的表面结构性质

基于自由能密度泛函理论(DFT)考察了二元排斥Yukawa (HCRY)流体在不同外场下的密度分布. 基于微扰理论, 体系的Helmholtz自由能泛函采用硬球排斥部分和长程色散部分贡献之和, 其中Kierlik和Rosinberg的加权密度近似(WDA)被用来计算硬球排斥部分, 而色散部分采用平均场理论(MFT)进行描述. 为了验证DFT计算结果的合理性, 研究中采用巨正则Monte Carlo(GCMC)模拟计算了在不同主体相密度、硬核直径和位能参数比的条件下二元HCRY混合流体的密度分布. 结果表明, 该DFT计算结果与GCMC模拟值吻合良好.     

Local solvent density augmentation around a solute in supercritical solvent bath: 1. A mechanism explanation and a new phenomenon

A recently proposed partitioned density functional (DF) approximation (Phys. Rev. E 2003, 68, 061201) and an adjustable parameter-free version of a Lagrangian theorem-based DF approximation (LTDFA: Phys. Lett. A 2003, 319, 279) are combined to propose a DF approximation for nonuniform Lennard-Jones (LJ) fluid. Predictions of the present DF approximation for local LJ solvent density inhomogeneity around a large LJ solute particle or hard core Yukawa particle are in good agreement with existing simulation data. An extensive investigation about the effect of solvent bath temperature, solvent-solute interaction range, solvent-solute interaction magnitude, and solute size on the local solvent density inhomogeneity is carried out with the present DF approximation. It is found that a plateau of solvent accumulation number as a function of solvent bath bulk density is due to a coupling between the solvent-solute interaction and solvent correlation whose mathematical expression is a convolution integral appearing in the density profile equation of the DF theory formalism. The coupling becomes stronger as the increasing of the whole solvent-solute interaction strength, solute size relative to solvent size, and the closeness to the critical density and temperature of the solvent bath. When the attractive solvent-solute interaction becomes large enough and the bulk state moves close enough to the critical temperature of the solvent bath, the maximum solvent accumulation number as a function of solvent bath bulk density appears near the solvent bath critical density; the appearance of this maximum is in contrast with a conclusion drawn by a previous investigation based on an inhomogeneous version of Ornstein-Zernike integral equation carried out only for a smaller parameter space than that in the present paper. Advantage of the DFT approach over the integral equation is discussed.     

Magnetophoresis, sedimentation, and diffusion of particles in concentrated magnetic fluids

A dynamic mass transfer equation for describing magnetophoresis, sedimentation, and gradient diffusion of colloidal particles in concentrated magnetic fluids has been derived. This equation takes into account steric, magnetodipole, and hydrodynamic interparticle interactions. Steric interactions have been investigated using the Carnahan-Starling approximation for a hard-sphere system. In order to study the effective interparticle attraction, the free energy of the dipolar hard-sphere system is represented as a virial expansion with accuracy to the terms quadratic in particle concentration. The virial expansion gives an interpolation formula that fits well the results of computer simulation in a wide range of particle concentrations and interparticle interaction energies. The diffusion coefficient of colloidal particles is written with regard to steric, magnetodipole and hydrodynamic interactions. We thereby laid the foundation for the formulation of boundary-value problems and for calculation of concentration and magnetic fields in the devices (for example, magnetic fluid seals and acceleration sensors), which use a concentrated magnetic fluid as a working fluid. The Monte-Carlo methods and the analytical approach are employed to study the magnetic fluid stratification generated by the gravitational field in a cylinder of finite height. The coefficient of concentration stratification of the magnetic fluid is calculated in relation to the average concentration of particles and the dipolar coupling constant. It is shown that the effective particle attraction causes a many-fold increase in the concentration inhomogeneity of the fluid if the average volume fraction of particles does not exceed 30%. At high volume concentrations steric interactions play a crucial role.     

Structures and correlation functions of multicomponent and polydisperse hard-sphere mixtures from a density functional theory

The structures of nonuniform binary hard-sphere mixtures and the correlation functions of uniform ternary hard-sphere mixtures were studied using a modified fundamental-measure theory based on the weight functions of Rosenfeld [Rosenfeld, Phys. Rev. Lett. 63, 980 (1989)] and Boublik-Mansoori-Carnahan-Starling-Leland equation of state [Boublik, J. Chem. Phys. 53, 471 (1970); Mansoori et al., J. Chem. Phys. 54, 1523 (1971)]. The theoretical predictions agreed very well with the molecular simulations for the overall density profiles, the local compositions, and the radial distribution functions of uniform as well as inhomogeneous hard-sphere mixtures. The density functional theory was further extended to represent the structure of a polydisperse hard-sphere fluid near a hard wall. Excellent agreement was also achieved between theory and Monte Carlo simulations. The density functional theory predicted oscillatory size segregations near a hard wall for a polydisperse hard-sphere fluid of a uniform size distribution.     

Structural properties of a model system with effective interparticle interaction potential applicable in modeling of complex fluids

We report grand canonical ensemble Monte Carlo (MC) simulation and theoretical studies of the structural properties of a model system described by an effective interparticle interaction potential, which incorporates basic interaction terms used in modeling of various complex fluids composed of mesoscopic particles dispersed in a solvent bath. The MC results for the bulk radial distribution function are employed to test the validity of the hard-sphere bridge function in combination with a modified hypernetted chain approximation (MHNC) in closing the Ornstein-Zernike (OZ) integral equation, while the MC data for the density profiles in different inhomogeneous environments are used to assess the validity of the third-order+second-order perturbation density functional theory (DFT). We found satisfactory agreement between the results predicted by the pure theories and simulation data, which classifies the proposed theoretical approaches as convenient tools for the investigation of complex fluids. The present investigation indicates that the bridge function approximation and density functional approximation, which are traditionally used for the study of neutral atomic fluids, also perform well for complex fluids only on condition that the underlying effective potentials include a highly repulsive core as an ingredient.     

How critical fluctuations influence adsorption properties of a van der Waals fluid onto a spherical colloidal particle

A recently proposed 3rd-order thermodynamic perturbation theory (TPT) is extended to its 5th-order version and non-uniform counterpart by supplementing with density functional theory (DFT) and a number of ansatzs for a bulk 2nd-order direct correlation function (DCF). Employment of the ansatzs DCF enables the resultant non-uniform formalism devoid of any adjustable parameter and free from numerically solving an Ornstein–Zernike integral equation theory. Density profiles calculated by the present non-uniform formalism for a hard core attractive Yukawa (HCAY) fluid near a spherical geometry are favorably compared with corresponding simulation data available in literature, and are more accurate than those based on a previous 3rd + 2nd-order perturbation DFT. The non-uniform 5th-order TPT is employed to investigate adsorption of the HCAY fluid onto a colloidal particle; it is disclosed that a depletion adsorption can be induced when the coexistence bulk fluid is situated in neighborhood of a critical point or near a bulk vapor–liquid coexistence gaseous phase or liquid phase density. A physical interpretation is given for such depletion adsorption and for its connection with parameters of the potential under consideration, which is ascribed to critical density fluctuations existing within a wide region of the bulk diagram. For a large spherical external potential inducing wetting transition, it is found that only round wetting transition is found instead of 1st-order pre-wetting transition in the case of a planar wall external potential, and the wetting transition temperature increases relative to that for the planar wall external potential. The present theoretical results for wetting transitions are supported by previous investigation based on thermodynamic considerations and a phenomenological Landau mean field theory, and are also in conformity with the present qualitative physical interpretation.     

Theory of repulsive charged colloids in slit-pores

Using classical density functional theory (DFT) we analyze the structure of the density profiles and solvation pressures of negatively charged colloids confined in slit pores. The considered model, which was already successfully employed to study a real colloidal (silica) suspension [S. H. L. Klapp et al., Phys. Rev. Lett. 100, 118303 (2008)], involves only the macroions which interact via the effective Derjaguin-Landau-Verwey-Overbeek (DLVO) potential supplemented by a hard core interaction. The solvent enters implicitly via the screening length of the DLVO interaction. The free energy functional describing the colloidal suspension consists of a hard sphere contribution obtained from fundamental measure theory and a long range contribution which is treated using two types of approximations. One of them is the mean field approximation (MFA) and the remaining is based on Rosenfeld's perturbative method for constructing the Helmholtz energy functional. These theoretical calculations are carried out at different bulk densities and wall separations to compare finally to grand canonical Monte Carlo simulations. We also consider the impact of charged walls. Our results show that the perturbative DFT method yields generally qualitatively consistent and, for some systems, also quantitatively reliable results. In MFA, on the other hand, the neglect of charge-induced correlations leads to a breakdown of this approach in a broad range of densities.     

Density functional theory of fluids in nanopores: analysis of the fundamental measures theory in extreme dimensional-crossover situations

Two density functional theories, the fundamental measures theory of Rosenfeld [Phys. Rev. Lett. 63, 980 (1989)] and a subsequent approximation by Tarazona [Phys. Rev. Lett. 84, 694 (2000)] are applied to the study of the hard-sphere fluid in two situations: the cylindrical pore and the spherical cavity. The results are compared with those obtained with grand canonical ensemble Monte Carlo simulations. The differences between both theories are evaluated and interpreted in the terms of the dimensional crossover from three to one and zero dimensions.     

Hard sphere fluids at a soft repulsive wall: a comparative study using Monte Carlo and density functional methods

Hard-sphere fluids confined between parallel plates at a distance D apart are studied for a wide range of packing fractions including also the onset of crystallization, applying Monte Carlo simulation techniques and density functional theory. The walls repel the hard spheres (of diameter σ) with a Weeks-Chandler-Andersen (WCA) potential V(WCA)(z) = 4ε[(σ(w)/z)(12) - (σ(w)/z)(6) + 1/4], with range σ(w) = σ/2. We vary the strength ε over a wide range and the case of simple hard walls is also treated for comparison. By the variation of ε one can change both the surface excess packing fraction and the wall-fluid (γ(wf)) and wall-crystal (γ(wc)) surface free energies. Several different methods to extract γ(wf) and γ(wc) from Monte Carlo (MC) simulations are implemented, and their accuracy and efficiency is comparatively discussed. The density functional theory (DFT) using fundamental measure functionals is found to be quantitatively accurate over a wide range of packing fractions; small deviations between DFT and MC near the fluid to crystal transition need to be studied further. Our results on density profiles near soft walls could be useful to interpret corresponding experiments with suitable colloidal dispersions.