A fundamental measure theory for the sticky hard sphere fluid

We construct a density functional theory (DFT) for the sticky hard sphere (SHS) fluid which, like Rosenfeld's fundamental measure theory (FMT) for the hard sphere fluid [Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989)], is based on a set of weighted densities and an exact result from scaled particle theory (SPT). It is demonstrated that the excess free energy density of the inhomogeneous SHS fluid Φ(SHS) is uniquely defined when (a) it is solely a function of the weighted densities from Kierlik and Rosinberg's version of FMT [E. Kierlik and M. L. Rosinberg, Phys. Rev. A 42, 3382 (1990)], (b) it satisfies the SPT differential equation, and (c) it yields any given direct correlation function (DCF) from the class of generalized Percus-Yevick closures introduced by Gazzillo and Giacometti [J. Chem. Phys. 120, 4742 (2004)]. The resulting DFT is shown to be in very good agreement with simulation data. In particular, this FMT yields the correct contact value of the density profiles with no adjustable parameters. Rather than requiring higher order DCFs, such as perturbative DFTs, our SHS FMT produces them. Interestingly, although equivalent to Kierlik and Rosinberg's FMT in the case of hard spheres, the set of weighted densities used for Rosenfeld's original FMT is insufficient for constructing a DFT which yields the SHS DCF.     

Scalar fundamental measure theory for hard spheres in three dimensions: application to hydrophobic solvation

Hard-sphere mixtures provide one a solvable reference system that can be used to improve the density functional theory of realistic molecular fluids. We show how the Kierlik-Rosinberg's scalar version of the fundamental measure density functional theory of hard spheres [E. Kierlik and M. L. Rosinberg, Phys. Rev. A 42, 3382 (1990)], which presents computational advantages with respect to the original Rosenfeld's vectorial formulation or its extensions, can be implemented and minimized in three dimensions to describe fluid mixtures in complex environments. This implementation is used as a basis for defining a molecular density functional theory of water around molecular hydrophobic solutes of arbitrary shape.     

Soft core thermodynamics from self-consistent hard core fluids

In an effort to generalize the self-consistent Ornstein-Zernike approximation (SCOZA)-an accurate liquid state theory that has been restricted so far to hard core systems-to arbitrary soft core systems we study a combination of SCOZA with a recently developed perturbation theory. The latter was constructed by Ben-Amotz and Stell [J. Phys. Chem. B 108, 6877 (2004)] as a reformulation of the Weeks-Chandler-Andersen [J. Chem. Phys. 54, 5237 (1971)] perturbation theory directly in terms of an arbitrary hard sphere reference system. We investigate the accuracy of the combined approach for the Lennard-Jones fluid in comparison with simulation data and pure perturbation theory predictions and determine the dependence of the thermodynamic properties and the phase behavior on the choice of the effective hard core diameter of the reference system.     

First-order mean spherical approximation for inhomogeneous fluids

The first-order mean-spherical approximation (FMSA) [Y. Tang, J. Chem. Phys., 118, 4140 (2003)] is extended to the studies of inhomogeneous fluids by combining with Rosenfeld's perturbative method [Y. Rosenfeld, J. Chem. Phys. 98, 8126 (1993)]. In the extension, the key input-direct correlation function of FMSA-is applied to constructing the free energy density functional. Preserving its high fidelity at the bulk limit, the FMSA shows satisfactory performance for Yukawa fluids near hard and attractive walls. The results are better than or comparable to several other theories reported before for the geometry. The FMSA is found, in particular, more satisfactory than the traditional mean-field theory for predicting density profiles around hard walls. The FMSA is also compared with the full MSA for inhomogeneous fluids, showing no appreciable differences. The inhomogeneous FMSA goes successfully through the self-consistency test for reproducing the radial distribution function of the bulk Yukawa fluid. As far as the computation is concerned, the FMSA can be executed much faster than any nonmean-field theories, and the speed is virtually identical to that of the mean-field theory.     

The behavior of fluids near solutes: a density functional theory and computer simulation study

The density distribution of solvent near a solute particle is studied using density functional theory and Monte Carlo simulation. The fluid atoms interact with each other via a hard sphere plus Yukawa potential, and interact with the solute via a hard sphere potential. For small solute sizes, the solvent displays liquidlike ordering near the particle. When the solute become larger, a drying transition is observed at state points near the coexistence conditions of the solvent. These predictions are similar to those of a recent theory for the hydrophobic effect by Lum, Chandler, and Weeks [J. Phys. Chem. 103, 4570 (1999)], although a comparison with simulations shows that the theory of this work is quantitatively more accurate. The connection between density functional methods and the LCW approach is also established.     

Shape factors in equations of state. Part II. Repulsion phenomena in multicomponent chain fluids

An equation of state for the multicomponent fluid phase of nonattracting rigid particles of arbitrary shape is presented. The equation is a generalization of a previously presented equation of state for pure fluids of rigid particles; the approach describes the volumetric properties of a pure fluid in terms of a shape factor, zeta, which can be back calculated by scaling the volumetric properties of pure fluids to that of a hard sphere. The performance of the proposed equation is tested against mixtures of chain fluids immersed in a "monomeric" solvent of hard spheres of equal and different sizes. Extensive new Monte Carlo simulation data are presented for 19 binary mixtures of hard homonuclear tangent freely-jointed hard sphere chains (pearl-necklace) of various lengths (three to five segments), with spheres of several size ratios and at various compositions. The performance of the proposed equation is compared to the hard-sphere SAFT approach and found to be of comparable accuracy. The equation proposed is further tested for mixtures of spheres with spherocylinders. In all cases, the equation proved to be accurate and simple to use.     

A new generalization of the Carnahan-Starling equation of state to additive mixtures of hard spheres

We introduce an expansion of the equation of state for additive hard-sphere mixtures in powers of the total packing fraction with coefficients which depend on a set of weighted densities used in scaled particle theory and fundamental measure theory. We demand that the mixture equation of state recovers the quasiexact Carnahan-Starling [J. Chem. Phys. 51, 635 (1969)] result in the case of a one-component fluid and show from thermodynamic considerations and consistency with an exact scaled particle relation that the first and second orders of the expansion lead unambiguously to the Boublik-Mansoori-Carnahan-Starling-Leland [J. Chem. Phys. 53, 471 (1970); J. Chem. Phys. 54, 1523 (1971)] equation and the extended Carnahan-Starling equation introduced by Santos et al. [Mol. Phys. 96, 1 (1999)]. In the third order of the expansion, our approach allows us to define a new equation of state for hard-sphere mixtures which we find to be more accurate than the former equations when compared to available computer simulation data for binary and ternary mixtures. Using the new mixture equation of state, we calculate expressions for the surface tension and excess adsorption of the one-component fluid at a planar hard wall and compare its predictions to available simulation data.     

Density functional theory of liquid crystals and surface anchoring: hard Gaussian overlap-sphere and hard Gaussian overlap-surface potentials

This article applies the density functional theory to confined liquid crystals, comprised of ellipsoidal shaped particles interacting through the hard Gaussian overlap (HGO) potential. The extended restricted orientation model proposed by Moradi and co-workers [J. Phys.: Condens. Matter 17, 5625 (2005)] is used to study the surface anchoring. The excess free energy is calculated as a functional expansion of density around a reference homogeneous fluid. The pair direct correlation function (DCF) of a homogeneous HGO fluid is approximated, based on the optimized sum of Percus-Yevick and Roth DCF for hard spheres; the anisotropy introduced by means of the closest approach parameter, the expression proposed by Marko [Physica B 392, 242 (2007)] for DCF of HGO, and hard ellipsoids were used. In this study we extend an our previous work [Phys. Rev. E 72, 061706 (2005)] on the anchoring behavior of hard particle liquid crystal model, by studying the effect of changing the particle-substrate contact function instead of hard needle-wall potentials. We use the two particle-surface potentials: the HGO-sphere and the HGO-surface potentials. The average number density and order parameter profiles of a confined HGO fluid are obtained using the two particle-wall potentials. For bulk isotropic liquid, the results are in agreement with the Monte Carlo simulation of Barmes and Cleaver [Phys. Rev. E 71, 021705 (2005)]. Also, for the bulk nematic phase, the theory gives the correct density profile and order parameter between the walls.     

Scaled particle theory for hard sphere pairs. II. Numerical analysis

We use the extension of scaled particle theory presented in the accompanying paper [F. H. Stillinger et al., J. Chem. Phys. 125, 204504 (2006)] to calculate numerically the pair correlation function of the hard sphere fluid over the density range 0< or =rhosigma(3)< or =0.96. Comparison with computer simulation results reveals that the new theory is able to capture accurately the fluid's structure across the entire density range examined. The pressure predicted via the virial route is systematically lower than simulation results, while that obtained using the compressibility route is lower than simulation predictions for rhosigma(3)< or =0.67 and higher than simulation predictions for rhosigma(3)> or =0.67. Numerical predictions are also presented for the surface tension and Tolman length of the hard sphere fluid.     

Theoretical direct correlation function for two-dimensional fluids of monodisperse hard spheres

The direct correlation function plays an important role in describing the effects of the structure of particle systems with respect to light diffraction, x-ray diffraction as well as transmission and transmission fluctuations of radiation through a dense suspension. In this paper, the direct correlation function for a monolayer of monodisperse hard spheres or disks is derived theoretically. Based on the approximation of Baus and Colot [Phys. Rev. A 36, 3912 (1987)] and the equation of state for a fluid of hard disks by Santos et al. [J. Chem. Phys. 103, 4622 (1995)], we propose a new direct correlation function, which compares well to the approximate analytical expressions and gives a good prediction of the structure factor in a wide range of monolayer density or suspension concentration. The resulting radial distribution function also agrees well with Monte Carlo computer simulation data. The corresponding contact values of the radial distribution function compare well with the results of analytic approximations, numerical solutions, and computer simulations. Our proposed direct correlation function is applied to the transmission fluctuation spectrometric study. Experimental results show good agreement with the theory.     

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.     

A consistent integral equation theory for hard spheres

The standard integral equation approach is used to extract the bridge function and other correlation functions of hard spheres fluid. To achieve this, we first use a recent consistent closure relation proposed by Bomont et al. [J. Chem. Phys. 119, 2188 (2003)] that has already proven to be accurate to describe the Lennard-Jones fluid properties. Second, we take advantage of the coherent scheme derived by Bomont [J. Chem. Phys. 119, 11484 (2003)] to calculate the excess chemical potential, the entropy and some relative transport properties. Very good agreement is obtained for structural quantities and thermodynamic properties as compared to exact data at densities ranging from 0.1 to 0.9.     

Structures and adsorption of binary hard-core Yukawa mixtures in a slitlike pore: grand canonical Monte Carlo simulation and density-functional study

The grand canonical ensemble Monte Carlo simulation and density-functional theory are applied to calculate the structures, local mole fractions, and adsorption isotherms of binary hard-core Yukawa mixtures in a slitlike pore as well as the radial distribution functions of bulk mixtures. The excess Helmholtz energy functional is a combination of the modified fundamental measure theory of Yu and Wu [J. Chem. Phys. 117, 10156 (2002)] for the hard-core contribution and a corrected mean-field theory for the attractive contribution. A comparison of the theoretical results with the results from the Monte Carlo simulations shows that the corrected theory improves the density profiles of binary hard-core Yukawa mixtures in the vicinity of contact over the original mean-field theory. Both the present corrected theory and the simulations suggest that depletion and desorption occur at low temperature, and the local segregation can be observed in most cases. For binary mixtures in the hard slitlike pore, the present corrected theory predicts more accurate surface excesses than the original one does, while in the case of the attractive pore, no improvement is found in the prediction of a surface excess of the smaller molecule.     

一个新的桥泛函及其在非均一流体密度泛函理论中的应用

基于对OZ 方程的渐近行为与Taylor级数展开的分析，提出了一个新的桥泛函，桥泛函被表达为间接相关函数的函数，Taylor级数展开的重整化导致了一个可调参数，通过将所提出的桥泛函与一个最近提出的密度泛函理论方法学，以及单个硬墙的sum 规则结合，可以确定可调参数.所提出的桥泛函能预言如下非均一流体的密度分布：硬球流体接近一个硬墙与在球形空隙内，Lennard-Jones 流体与缔合硬球流体在两个硬墙之内.理论预言与文献所报导的模拟数据符合很好.     

Crystal nucleation of hard spheres using molecular dynamics, umbrella sampling, and forward flux sampling: a comparison of simulation techniques

Over the last number of years several simulation methods have been introduced to study rare events such as nucleation. In this paper we examine the crystal nucleation rate of hard spheres using three such numerical techniques: molecular dynamics, forward flux sampling, and a Bennett-Chandler-type theory where the nucleation barrier is determined using umbrella sampling simulations. The resulting nucleation rates are compared with the experimental rates of Harland and van Megen [Phys. Rev. E 55, 3054 (1997)], Sinn et al. [Prog. Colloid Polym. Sci. 118, 266 (2001)], Sch?tzel and Ackerson [Phys. Rev. E 48, 3766 (1993)], and the predicted rates for monodisperse and 5% polydisperse hard spheres of Auer and Frenkel [Nature 409, 1020 (2001)]. When the rates are examined in units of the long-time diffusion coefficient, we find agreement between all the theoretically predicted nucleation rates, however, the experimental results display a markedly different behavior for low supersaturation. Additionally, we examined the precritical nuclei arising in the molecular dynamics, forward flux sampling, and umbrella sampling simulations. The structure of the nuclei appears independent of the simulation method, and in all cases, the nuclei contains on average significantly more face-centered-cubic ordered particles than hexagonal-close-packed ordered particles.     

Simple bond length dependence: a correspondence between reactive fluid theories

Two elementary models of reactive fluids are examined, the first being a standard construction assuming molecular dissociation at infinite separation; the second is an open mixture of nondissociative molecules and free atoms in which the densities of free atoms and molecules are coupled. An approximation to the density of molecules, to low order in site density, is derived in terms of the classical associating fluid theory variously described by Wertheim [J. Chem. Phys. 87, 7323 (1987)] and Stell [Physica A 231, 1 (1996)]. The results are derived for a fluid of dimerizing hard spheres, and predict dependence of the molecular density on the total site density, the hard sphere diameter, and the bond length of the dimer. The results for the two reactive models are shown to be qualitatively similar, and lead to equivalent predictions of the molecular density for the infinitely short and infinitely long bond lengths.     

Molecular dynamics and theory for the contact values of the radial distribution functions of hard-disk fluid mixtures

We report molecular dynamics results for the contact values of the radial distribution functions of binary additive mixtures of hard disks. The simulation data are compared with theoretical predictions from expressions proposed by Jenkins and Mancini [J. Appl. Mech. 54, 27 (1987)] and Santos et al. [J. Chem. Phys. 117, 5785 (2002)]. Both theories agree quantitatively within a very small margin, which renders the former still a very useful and simple tool to work with. The latter (higher-order and self-consistent) theory provides a small qualitative correction for low densities and is superior especially in the high-density domain.     

Density functional theory of adsorption of mixtures of charged chain particles and spherical counterions

We propose a microscopic density functional theory to describe nonuniform ionic fluids composed of chain molecules with charged "heads" and spherical counterions. The chain molecules are modeled as freely jointed chains of hard spheres, the counterions are oppositely charged spheres of the same diameter as all segments of chain molecules. The theory is based on the approach of Yu and Wu [J. Chem. Phys. 117, 2368 (2002)] of adsorption of chain molecules and on theory of adsorption of electrolytes [O. Pizio, A. Patrykiejew, and S. Sokolowski, J. Chem. Phys. 121, 11957 (2004)]. As an application of the proposed formalism we investigate the structure and adsorption of fluids containing segments of different length in a slitlike pore.     

The equation of state of isotropic fluids of hard convex bodies from a high-level virial expansion

We have calculated virial coefficients up to seventh order for the isotropic phases of a variety of fluids composed of hard aspherical particles. The models studied were hard spheroids, hard spherocylinders, and truncated hard spheres, and results are obtained for a variety of length-to-width ratios. We compare the predicted virial equations of state with those determined by simulation. We also use our data to calculate the coefficients of the y expansion [B. Barboy and W. M. Gelbart, J. Chem. Phys. 71, 3053 (1979)] and to study its convergence properties. Finally, we use our data to estimate the radius of convergence of the virial series for these aspherical particles. For fairly spherical particles, we estimate the radius of convergence to be similar to that of the density of closest packing. For more anisotropic particles, however, the radius of convergence decreases with increased anisotropy and is considerably less than the close-packed density.     

Theoretical aspects and computer simulations of flexible charged oligomers in salt-free solutions

The structural and thermodynamic properties of a model solution containing flexible charged oligomers and an equivalent number of counterions were studied by means of the canonical Monte Carlo simulation and integral equation theory. The oligomers were represented as freely jointed chains of charged hard spheres. In accordance with the primitive model of electrolyte solutions, the counterions were modeled as charged hard spheres and the solvent as a dielectric continuum. Simulations were performed for a set of model parameters, independently varying the chain length and concentration of the oligomers. Structural properties in the form of pair distribution functions were calculated as functions of model parameters. In addition, thermodynamic properties such as the excess energy of solution and the excess chemical potential of counterions were obtained. These properties were correlated with the conformational averages of oligomers as reflected in the end-to-end distances and radii of gyration obtained from the simulations. The relation with the experimental data for heats of dilution and for the activity coefficient is discussed. Finally, theories based on Wertheim's integral equation approach (product reactant Ornstein-Zernike approach) [J. Stat. Phys. 42, 477 (1986)] in the so-called polymer mean spherical and polymer hypernetted chain approximations were tested against the new and existing computer simulations. For the values of parameters examined in this study, the integral equation theory yields semiquantitative agreement with computer simulations.