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.     

Further test of third order + second-order perturbation DFT approach: hard core repulsive Yukawa fluid subjected to diverse external fields

Grand canonical Monte Carlo simulation is used to investigate density profiles of hard-core repulsive Yukawa (HCRY) model fluid under the influence of various external fields and radial distribution function (RDF) of the bulk HCRY system. The aim of these extensive simulations is to provide exact data for purely repulsive interaction potential against which the validity of a third order + second-order perturbation DFT approach can be tested. It is found that a semiempirical parametrized bridge function due to Malijevsky and Labik performs very well for the RDF of the bulk HCRY fluid. Incorporation of a bulk second-order direct correlation function (DCF) of the HCRY fluid based on the Malijevsky-Labik bridge function into the third order + second-order perturbation DFT approach yields the resulting theoretical predictions for the density profiles of inhomogeneous HCRY fluid that are in a very good agreement with the simulation data, an exception being somewhat larger deviations appearing for the structure of the fluid around the center of a hard spherical cavity. Both theory and simulation predict layering transition and gas-liquid coexistence phenomena occurring with the HCRY model fluid under confined conditions. For the case of an inverse sixth-power repulsive potential under the influence of a flat stationary wall defined by an inverse twelfth-power repulsive potential, the present third order + second-order perturbation DFT approach is found to be superior to several existing weighted density approximations (WDA) and partitioned WDA.     

Structure of inhomogeneous attractive and repulsive hard-core yukawa fluid: grand canonical Monte Carlo simulation and density functional theory study

A density functional theory is proposed for an inhomogeneous hard-core Yukawa (HCY) fluid based on Rosenfeld's perturbative method. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for the hard-core repulsion and a quadratic functional Taylor expansion for the long-ranged attractive or repulsive interactions. To test the established theory, grand canonical ensemble Monte Carlo simulations are carried out to simulate the density profiles of attractive and repulsive HCY fluid near a wall. Comparison with the results from the Monte Carlo simulations shows that the present density functional theory gives accurate density profiles for both attractive and repulsive HCY fluid near a wall. Both the present theory and simulations suggest that there is depletion for attractive HCY fluid at low temperature, but no depletion is found for repulsive HCY fluid. The calculated results indicate that the present density functional theory is better than those of the modified version of the Lovett-Mou-Buff-Wertheim and other density functional theories. The present theory is simple in form and computationally efficient. It predicts accurate radial distribution functions of both attractive and repulsive HCY fluid except for the repulsive case at high density, where the theory overestimates the radial distribution function in the vicinity of contact.     

Density profiles and solvation forces for a Yukawa fluid in a slit pore

The effect of varying wall-particle and particle-particle interactions on the density profiles near a single wall and the solvation forces between two walls immersed in a fluid of particles is investigated by grand canonical Monte Carlo simulations. Attractive and repulsive particle-particle and particle-wall interactions are modeled by a versatile hard-core Yukawa form. These simulation results are compared to theoretical calculations using the hypernetted chain integral equation technique, as well as with fundamental measure density functional theory (DFT), where particle-particle interactions are either treated as a first order perturbation using the radial distribution function or else with a DFT based on the direct-correlation function. All three theoretical approaches reproduce the main trends fairly well, but exhibit inconsistent accuracy, particularly for attractive particle-particle interactions. We show that the wall-particle and particle-particle attractions can couple together to induce a nonlinear enhancement of the adsorption and a related "repulsion through attraction" effect for the effective wall-wall forces. We also investigate the phenomenon of bridging, where an attractive wall-particle interaction induces strongly attractive solvation forces.     

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.     

Structure and thermodynamics of hard-core Yukawa fluids: thermodynamic perturbation approaches

The thermodynamic perturbation theories, which are based on the power series of a coupling constant (λ-expansion), have been proposed for studying the structural and thermodynamic properties of a hard-core Yukawa (HCY) fluid: one (A1-approximation) is the perturbation theory based on the hard-sphere repulsion as a reference system. The other (A2-approximation) is the perturbation theory based on the reference system which incorporates both the repulsive and short-range attractive interactions. The first-order mean-spherical approximation (FMSA) provided by Tang and Lu [J. Chem. Phys. 99, 9828 (1993)] has been employed for investigating the thermodynamic properties of a HCY fluid using the alternative method via the direct correlation function. The calculated results show that (i) the A1 and A2 approximations are in excellent agreements with previous computer simulation results in the literature and compare with the semi-empirical works of Shukla including the higher-order free energy terms, (ii) the A1 and A2 approximations are better than the FMSA and the mean-spherical approximation, (iii) the A2-approximation compares with the A1-approximation, even though the perturbation effect of an A2-approximation is much smaller than that of an A1-approximation, and that (iv) the FMSA study is particularly of advantage in providing the structure and thermodynamics in a simple and analytic manner.     

Non-hard sphere thermodynamic perturbation theory

A non-hard sphere (HS) perturbation scheme, recently advanced by the present author, is elaborated for several technical matters, which are key mathematical details for implementation of the non-HS perturbation scheme in a coupling parameter expansion (CPE) thermodynamic perturbation framework. NVT-Monte Carlo simulation is carried out for a generalized Lennard-Jones (LJ) 2n-n potential to obtain routine thermodynamic quantities such as excess internal energy, pressure, excess chemical potential, excess Helmholtz free energy, and excess constant volume heat capacity. Then, these new simulation data, and available simulation data in literatures about a hard core attractive Yukawa fluid and a Sutherland fluid, are used to test the non-HS CPE 3rd-order thermodynamic perturbation theory (TPT) and give a comparison between the non-HS CPE 3rd-order TPT and other theoretical approaches. It is indicated that the non-HS CPE 3rd-order TPT is superior to other traditional TPT such as van der Waals/HS (vdW/HS), perturbation theory 2 (PT2)/HS, and vdW/Yukawa (vdW/Y) theory or analytical equation of state such as mean spherical approximation (MSA)-equation of state and is at least comparable to several currently the most accurate Ornstein-Zernike integral equation theories. It is discovered that three technical issues, i.e., opening up new bridge function approximation for the reference potential, choosing proper reference potential, and/or using proper thermodynamic route for calculation of f(ex-ref), chiefly decide the quality of the non-HS CPE TPT. Considering that the non-HS perturbation scheme applies for a wide variety of model fluids, and its implementation in the CPE thermodynamic perturbation framework is amenable to high-order truncation, the non-HS CPE 3rd-order or higher order TPT will be more promising once the above-mentioned three technological advances are established.     

Use of the Yukawa fluid as the reference system in perturbation theory

The Weeks-Chandler-Anderson (WCA) perturbation theory is studied utilising recent results for the Yukawa model fluid. Replacing the attractive tail of the Lennard-Jones potential with a Yukawa tail, where the Yukawa parameters are chosen using a least squares fit, it is shown that accurates field dstribution functions can be generated via the EXP approximation of the WCA optimized cluster theory. The comparative case and accuracy with which the correlation functions for the Yukawa fluid can be compared render this a very useful method for studying the equilibrium properties of simple liquids.     

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

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

The structural properties of a two-Yukawa fluid: Simulation and analytical results

Standard Monte Carlo simulations are carried out to assess the accuracy of theoretical predictions for the structural properties of a model fluid interacting through a hard-core two-Yukawa potential composed of a short-range attractive well next to a hard repulsive core, followed by a smooth, long-range repulsive tail. Theoretical calculations are performed in the framework provided by the Ornstein-Zernike equation, solved either analytically with the mean spherical approximation (MSA) or iteratively with the hypernetted-chain (HNC) closure. Our analysis shows that both theories are generally accurate in a thermodynamic region corresponding to a dense vapor phase around the critical point. For a suitable choice of potential parameters, namely, when the attractive well is deep and/or large enough, the static structure factor displays a secondary low-Q peak. In this case HNC predictions closely follow the simulation results, whereas MSA results progressively worsen the more pronounced this low-Q peak is. We discuss the appearance of such a peak, also experimentally observed in colloidal suspensions and protein solutions, in terms of the formation of equilibrium clusters in the homogeneous fluid.     

Effect of quenched size polydispersity on the fluid-solid transition in charged colloidal suspensions

We study the effect of quenched size polydispersity on the phase behavior of charged colloidal suspensions using free-energy calculations in Monte Carlo simulations. The colloids are assumed to interact with a hard-core repulsive Yukawa (screened-Coulomb) interaction with constant surface potential, so that the particles are polydisperse both in size and charge. In addition, we take the size distribution to be fixed in both the fluid and crystal phase (no size fractionation is allowed). We study the fluid-solid transition for various screening lengths and surface potentials, finding that upon increasing the size polydispersity the freezing transition shifts toward higher packing fractions and the density discontinuity between the two coexisting phases diminishes. Our results provide support for a terminal polydispersity above which the freezing transition disappears.     

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.     

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.     

Yukawa流体的相平衡和界面张力研究

应用二阶微扰理论, Duh-Mier-Y-Teran状态方程和在平均球近似(mean spherical approximation, MSA)的基础上获得的直接相关函数, 建立了适用于均匀流体和非均匀流体的状态方程. 结合此状态方程, 重整化群理论(renormalization group theory, RG)和密度泛函理论(density functional theory, DFT), 分别研究了Yukawa流体的相平衡和界面张力. 结果与分子模拟数据吻合良好.     

Sedimentation equilibrium of colloidal suspensions in a planar pore based on density functional theory and the hard-core attractive Yukawa model

The sedimentation equilibrium of colloidal suspensions modeled by hard-core attractive Yukawa (HCAY) fluids in a planar pore is studied. The density profile of the HCAY fluid in a gravitational field and its distribution between the pore and uniform phases are investigated by a density functional theory (DFT) approach, which results from employing a recently proposed parameter-free version of the Lagrangian theorem-based density functional approximation (Zhou, S. Phys. Lett. A 2003, 319, 279) for hard-sphere fluids to the hard-core part of the HCAY fluid, and the second-order functional perturbation expansion approximation to the tail part as was done in a recent partitioned density functional approximation (Zhou, S. Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 2003, 68, 061201). The resultant DFT approach is, thus, the first adjustable parameter-free DFT for HCAY fluids. The validity of the present DFT for HCAY fluids of reduced range parameter z(red) = 1.8 under various external potentials is established in the first of the papers cited previously. The present DFT for HCAY fluids can predict the radial distribution function for the bulk HCAY fluid accurately in the colloidal limit (large value of z(red)), and in the hard-sphere limit, its prediction for the density profile of the hard-sphere fluid in a gravitational field is in very good agreement with the existing simulation data. The dependence of the density profile and distribution coefficient on the magnitude of the interparticle attraction, gravitational field, and degree of confinement is investigated in detail by the present DFT approach. Intuitive and qualitative analyses are also compared with the quantitative DFT calculational results.     

Thermodynamic and structural properties of mixed colloids represented by a hard-core two-Yukawa mixture model fluid: Monte Carlo simulations and an analytical theory

The interaction between colloidal particles is well represented by a hard-core two-Yukawa potential. In order to assess the accuracy of theoretical predictions for the thermodynamic and structural properties of mixed colloids, standard Monte Carlo simulations are carried out for the hard-core two-Yukawa mixtures. In the simulations, one range parameter in the two-Yukawa potential is taken as 1.8 or 2.8647, and another is taken as 4, 8, or 13.5485. Both attractive and repulsive dominant cases of the potential outside the hard core are considered. The effects of temperature, density, composition, size and energy parameter ratios on internal energy, compressibility factor, and radial distribution function are investigated extensively. Theoretical calculations are performed in the framework of analytical solution for the Ornstein-Zernike equation with the first-order mean spherical approximation (FMSA). Our analysis shows that the FMSA is very accurate for the prediction of the compressibility factor of the hard-core two-Yukawa mixtures at all conditions studied. The FMSA generally predicts accurate internal energy, but overestimates the internal energy of the systems at lower temperatures. Furthermore, we found that a simplified exponential version of the FMSA predicts fairly good radial distribution function at contact for the mixed two-Yukawa fluids. The comparison of the theoretical compressibility factor with that from the Monte Carlo simulations suggests that the FMSA can be used to investigate the fluid-fluid equilibria of hard-core two-Yukawa mixtures.     

硬核Yukawa流体中胶体粒子间的耗尽势

基于Roth、Evans和Dietrich有关耗尽势的密度泛函理论研究了硬核Yukawa(HCY)流体中胶体粒子间的耗尽势.在极稀溶液条件下,通过计算两个胶体粒子在不同条件下的耗尽势,分析了HCY流体的相关因素对耗尽势的影响.结果表明,胶体粒子与溶剂分子的尺寸比率、HCY流体分子间的相互作用、HCY流体分子的体相密度以及胶体粒子与流体分子之间的相互作用等因素均可对胶体粒子间耗尽势产生显著影响.研究结果可为实验上调控胶体粒子间的相互作用提供可能的线索.     

Local structure and thermodynamics of a core-softened potential fluid: theory and simulation.

Phase behavior and structural properties of homogeneous and inhomogeneous core-softened (CS) fluid consisting of particles interacting via the potential, which combines the hard-core repulsion and double attractive well interaction, are investigated. The vapour-liquid coexistence curves and critical points for various interaction ranges of the potential are determined by discrete molecular dynamics simulations to provide guidance for the choice of the bulk density and potential parameters for the study of homogeneous and inhomogeneous structures. Spatial correlations in the homogeneous CS system are studied by the Ornstein-Zernike integral equation in combination with the modified hypernetted chain (MHNC) approximation. The local structure of CS fluid subjected to diverse external fields maintaining the equilibrium with the bulk CS fluid are studied on the basis of a recently proposed third order+second order perturbation density functional approximation (DFA). The accuracy of DFA predictions is tested against the results of a grand canonical ensemble Monte Carlo simulation. Reasonable agreement between the results of both methods proves that the DFA theory applied in this work is a convenient theoretical tool for the investigation of the CS fluid, which is practically applicable for modeling numerous real systems.     

Pair force distributions in simple fluids

Analytic expressions are derived for the frequency distribution, P(f), of pair forces, f, and those of their α-Cartesian component, f(α), or P(f(α)), for some typical model simple fluids, expressed in terms of the radial distribution function and known constants. For strongly repulsive inverse power (IP), exponential and Yukawa purely repulsive potentials, P(f) diverges at the origin approximately as ～f(-1), but with different limiting analytic forms. P(f(α)) is also shown to diverge as ～f(-1) as f → 0 for the IP fluid. For the Lennard-Jones potential fluid, P(f) is finite for all f ≥ 0 but has two singularities for negative f, corresponding to the zero force limit (i.e., f → 0(-)) and the point of inflection in the potential. The corresponding component force distribution is singular as f(α) → 0 from both positive and negative force sides. The large force limit of P(f), which originates from the close neighbor interactions, is nearly exponential for the IP and LJ fluids, as is also found for granular materials. A more complete picture of force distributions in off-lattice particulate systems as a function of force law and state point (particularly the extent of "thermalization" of the particles) is provided.     

Self-assembly of patchy particles into polymer chains: a parameter-free comparison between Wertheim theory and Monte Carlo simulation

The authors numerically study a simple fluid composed of particles having a hard-core repulsion, complemented by two short-ranged attractive (sticky) spots at the particle poles, which provides a simple model for equilibrium polymerization of linear chains. The simplicity of the model allows for a close comparison, with no fitting parameters, between simulations and theoretical predictions based on the Wertheim perturbation theory. This comparison offers a unique framework for the analytic prediction of the properties of self-assembling particle systems in terms of molecular parameters and liquid state correlation functions. The Wertheim theory has not been previously subjected to stringent tests against simulation data for ordering across the polymerization transition. The authors numerically determine many of the thermodynamic properties governing this basic form of self-assembly (energy per particle, order parameter or average fraction of particles in the associated state, average chain length, chain length distribution, average end-to-end distance of the chains, and the static structure factor) and find that predictions of the Wertheim theory accord remarkably well with the simulation results.