Alternative fundamental measure theory for additive hard sphere mixtures

The purpose of this short paper is to present an alternative fundamental measure theory (FMT) for hard sphere mixtures. Keeping the main features of the original Rosenfeld's FMT [Phys. Rev. Lett. 63, 980 (1989)] and using the dimensional and the low-density limit conditions a new functional is derived incorporating Boublik's multicomponent extension [Mol. Phys. 59, 371 (1986)] of highly accurate Kolafa's equation of state for pure hard spheres. We test the theory for pure hard spheres and hard sphere mixtures near a planar hard wall and compare the results with the original Rosenfeld's FMT and one of its modifications and with new very accurate simulation data. The test reveals an excellent agreement between the results based on the alternative FMT and simulation data for density profile near a contact and some improvement over the original Rosenfeld's FMT and its modification at the contact region.     

Structure of inhomogeneous polymer solutions: a density functional approach

The structure of polymer solutions confined between surfaces is studied using a density functional theory where the polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres. The present theory uses the concept of universality of the free energy density functional to obtain the first-order direct correlation function of the nonuniform system from that of the corresponding uniform system, calculated through the Verlet-modified type bridge function. The uniform bulk fluid direct correlation function required as input has been calculated from the reference interaction site model integral equation theory using the Percus-Yevick closure relation. The calculated results on the density profiles of the polymer as well as the solvent are shown to compare well with computer simulation results.     

Interactions between colloidal particles in polymer solutions: A density functional theory study

We present a density functional theory study of colloidal interactions in a concentrated polymer solution. The colloids are modeled as hard spheres and polymers are modeled as freely jointed tangent hard sphere chains. Our theoretical results for the polymer-mediated mean force between two dilute colloids are compared with recent simulation data for this model. Theory is shown to be in good agreement with simulation. We compute the colloid-colloid potential of mean force and the second virial coefficient, and analyze the behavior of these quantities as a function of the polymer solution density, the polymer chain length, and the colloid/polymer bead size ratio.     

On the accuracy of density functional theory for iron-sulfur clusters

A simple, yet powerful wave function manipulation method was introduced utilizing a generalized ionic fragment approach that allows for systematic mapping of the wave function space for multispin systems with antiferromagnetic coupling. The use of this method was demonstrated for developing ground state electronic wave function for [2Fe-2S] and [Mo-3Fe-4S] clusters. Using well-defined ionic wave functions for ferrous and ferric irons, sulfide, and thiolate fragments, the accuracy of various density functionals and basis sets including effective core potentials were evaluated on a [4Fe-4S] cluster by comparing the calculated geometric and electronic structures with crystallographic data and experimental atomic spin densities from X-ray absorption spectroscopy, respectively. We found that the most reasonable agreement for both geometry and atomic spin densities is obtained by a hybrid functional with 5% HF exchange and 95% density functional exchange supplemented with Perdew's 1986 correlation functional. The basis set seems to saturate only at the triple-zeta level with polarization and diffuse functions. Reasonably preoptimized structures can be obtained by employing computationally less expensive effective core potentials, such as the Stuttgart-Dresden potential with a triple-zeta valence basis set. The extension of the described calibration methodology to other biologically important and more complex iron-sulfur clusters, such as hydrogenase H-cluster and nitrogenase FeMo-co will follow.     

Perturbation theory for liquids with a hard sphere plus square well potential

Numerical data have been obtained for a series of perturbation theory in the form of Barker-Henderson discrete representation. A hard sphere liquid (864 particles) was modeled by the Monte Carlo method for 106 values of occupancy from the range η = 0.005–0.530 (step 0.005); the well width was varied from 0 to 1.5 hard sphere diameters. For each occupancy, averaging was done over 60·10⁶ ensembles. Approximating analytical expressions are given as functions of well density and width for the first three terms of the free energy decomposition. For reduced temperature T* = 1.0, the standard error in the free energy calculation (three orders of perturbation theory if the above expressions are used) is of the order of ±0.001. Analysis of four-particle correlations revealed peaks at distances of 1.35 Å and 1.65 Å (η = 0.53) in the high-density region of hard sphere liquids, which are absent in crystals and dense liquids. It is analyzed whether the results are useful for realization of WCA theory of simple liquids in second order perturbation theory.     

Structure of hard sphere packings near Bernal density

The structural features of packings of similar hard spheres in the vicinity of Bernal density corresponding to space occupancy ∼0.64 have been studied. This is maximum density for disordered packings, which give way to crystallization at higher densities. The structure was analyzed using Delaunay simplices. Aggregates of simplices approximating a regular tetrahedron (polytetrahedra) in form have been studied. They have a high local density and various morphologies, but are incompatible with translation symmetry. Percolation analysis of these clusters was carried out. The “critical” nature of structural transformations at Bernal density is related with the appearance of percolation through tetrahedra. The closest disordered packings really have no significant nucleating seeds of crystal structures, as shown by a sensitive seed detection procedure. Their appearance and fast growth are observed after passing through Bernal density. The total fraction of the crystal phase increases monotonically with the density; the fcc and hcp structures appear in different proportions in the packings.     

Thermodynamics of attractive hard rods: a test of mean field density functional theory

Mean field density functional theory (MFDFT) has been employed to calculate the free energy of a pair of attractive hard rods on a ring. The results for homogeneous and optimal inhomogeneous density profiles have been compared with the exact free energy as a test of the approach. We discuss the problems in applying MFDFT to small systems and suggest modifications which allow a reasonably accurate treatment of this particular, rather extreme, case.     

Scaled particle theory of mixtures of hard spheres with negatively non-additive diameters

The scaled particle theory is solved for mixtures of hard spheres with arbitrary negative non-additive diameters in three dimensions. The results obtained with the Gibbs—Duhem relation show excellent agreement with the Monte-Carlo calculations. The results from the virial relation for high densities and high non-additivity are less satisfactory.     

Exact density functional theory for ideal polymer fluids with nearest neighbor bonding constraints

We present a new density functional theory of ideal polymer fluids, assuming nearest-neighbor bonding constraints. The free energy functional is expressed in terms of end site densities of chain segments and thus has a simpler mathematical structure than previously used expressions using multipoint distributions. This work is based on a formalism proposed by Tripathi and Chapman [Phys. Rev. Lett. 94, 087801 (2005)]. Those authors obtain an approximate free energy functional for ideal polymers in terms of monomer site densities. Calculations on both repulsive and attractive surfaces show that their theory is reasonably accurate in some cases, but does differ significantly from the exact result for longer polymers with attractive surfaces. We suggest that segment end site densities, rather than monomer site densities, are the preferred choice of "site functions" for expressing the free energy functional of polymer fluids. We illustrate the application of our theory to derive an expression for the free energy of an ideal fluid of infinitely long polymers.     

Modeling of the hydration sphere of gadolinium(III) ion using density functional theory

The hydrated gadolinium(III) ion cluster Gd³⁺ (H₂O)_x, with x = 8,9, was studied using density functional theory. The different conformations of the first hydration shell were calculated. For x = 8, the results for the cubic conformation correspond to previously published Hartree-Fock and MP2 results, whereas much lower energies were found for the square antiprismatic and dodecahedral conformations. For x = 9, the energy of the tricapped trigonal prism is nearly identical to the one of the antiprism plus an extra free water molecule. © 1997 John Wiley & Sons, Inc.     

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.     

Dielectric properties of a dipolar hard sphere liquid: Association equilibria theory

The hindered rotation model is used to generalize the association equilibria theory for the case of a dense dipolar hard sphere (DHS) liquid. The derived expressions for the orientation constant of association and the Kirkwood factor of a DHS liquid are well consistent with machine experiments in the high temperature region.     

The density profile of hard sphere liquid system under gravity

The density profile of hard sphere liquid under gravity is calculated by using density functional theory and Monte Carlo simulation method. The two methods give consistent results for a wide range of parameters. Meanwhile, the validity range of the density functional theory is also established. The results are quite different from the barometric height distribution rho(z)=rho(0) exp(-zL(G)) in almost all cases studied, which indicates that the interaction between particles plays an important role in the density distribution under external fields. Moreover, the crystallizing phenomenon is also predicted at the bottom part of the system under strong gravitation.     

Surface forces in polymer fluids: a comparison between simulations and density functional theory

A polymer density functional theory is evaluated in terms of its ability to predict interactions between large surfaces in a polymer fluid. Comparisons are made with results from simulations in an expanded isotension ensemble. The variation of the net surface-surface interaction with adsorption strength is examined. Cases where the monomers interact via a pure hard sphere potential are investigated, but we have also studied the effect of attractions between the monomers. In all cases, we obtain an almost quantitative agreement between the simulated results and the predictions from the polymer density functional theory.     

Effect of attractions on the structure of polymer solutions confined between surfaces: a density functional approach

A density functional theory is presented to study the effect of attractions on the structure of polymer solutions confined between surfaces. The polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres, both having Yukawa-type attractions and the mixture being confined between attractive Yukawa-type surfaces. The present theory treats the ideal gas free energy functional exactly and uses weighted density approximation for the hard chain and hard sphere contributions to the excess free energy functional. The attractive interactions are calculated using the direct correlation function obtained from the polymer reference interaction site model theory along with the mean spherical approximation closure. The theoretical predictions on the density profiles of the polymer and the solvent molecules are found to agree quite well with the Monte Carlo simulation results for varying densities, chain lengths, wall separations, and different sets of interaction potentials.     

Interactions between polymer brushes in solvents of variable quality: a density functional theory study

We present a density functional theory study of interactions between sterically stabilized colloidal particles in solvents of variable quality. Both flat and spherical polymer brushes are considered, as well as both monatomic and polymeric solvents. It is shown that the interaction between sterically stabilized particles can be tuned from repulsive to attractive by varying the solvent quality, the relative length of free and grafted chains, and by employing a mixed brush consisting of both well and poorly solvated chains.     

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.     

Mesophase formation in a system of top-shaped hard molecules: density functional theory and Monte Carlo simulation

We present the phase diagram of a system of mesogenic top-shaped molecules based on the Parsons-Lee density functional theory and Monte Carlo simulation. The molecules are modeled as a hard spherocylinder with a hard sphere embedded in its center. The stability of five different phases is studied, namely, isotropic, nematic, smectic A, smectic C, and columnar phases. The positionally ordered phases are investigated only for the case of parallel alignment. It is found that the central spherical unit destabilizes the nematic with respect to the isotropic phase, while increasing the length of the cylinder has the opposite effect. Also, the central hard sphere has a strong destabilizing effect on the smectic A phase, due the inefficient packing of the molecules into layers. For large hard sphere units the smectic A phase is completely replaced by a smectic C structure. The columnar phase is first stabilized with increasing diameter of the central unit, but for very large hard sphere units it becomes less stable again. The density functional results are in good agreement with the simulations.