The freezing transition of a hard sphere fluid subject to the Percus-Yevick approximation

A classical density functional theory is applied to the calculation of the fluid-solid transition for hard spheres, using the Percus-Yevick (PY) direct correlation function. Three algebraic conditions are established for the coexistence densities and the Lindemann parameter. The terms neglected in the present analysis are small and the present theory, in our eyes, is essentially an exact solution given the PY approximation. No fluid-solid transition is found for the bcc lattice, whereas for expanded fcc lattices, the agreement with previous density functional theory-based theories is good.     

Rescaled density expansions and demixing in hard-sphere binary mixtures

The demixing transition of a binary fluid mixture of additive hard spheres is analyzed for different size asymmetries by starting from the exact low-density expansion of the pressure. Already within the second virial approximation the fluid separates into two phases of different composition with a lower consolute critical point. By successively incorporating the third, fourth, and fifth virial coefficients, the critical consolute point moves to higher values of the pressure and to lower values of the partial number fraction of the large spheres. When the exact low-density expansion of the pressure is rescaled to higher densities as in the Percus-Yevick theory, by adding more exact virial coefficients a different qualitative movement of the critical consolute point in the phase diagram is found. It is argued that the Percus-Yevick factor appearing in many empirical equations of state for the mixture has a deep influence on the location of the critical consolute point, so that the resulting phase diagram for a prescribed equation has to be taken with caution.     

Virial series for fluids of hard hyperspheres in odd dimensions

A recently derived method [R. D. Rohrmann and A. Santos, Phys. Rev. E 76, 051202 (2007)] to obtain the exact solution of the Percus-Yevick equation for a fluid of hard spheres in (odd) d dimensions is used to investigate the convergence properties of the resulting virial series. This is done both for the virial and compressibility routes, in which the virial coefficients B(j) are expressed in terms of the solution of a set of (d-1)/2 coupled algebraic equations which become nonlinear for d>/=5. Results have been derived up to d=13. A confirmation of the alternating character of the series for d>/=5, due to the existence of a branch point on the negative real axis, is found and the radius of convergence is explicitly determined for each dimension. The resulting scaled density per dimension 2eta(1/d), where eta is the packing fraction, is wholly consistent with the limiting value of 1 for d-->infinity. Finally, the values for B(j) predicted by the virial and compressibility routes in the Percus-Yevick approximation are compared with the known exact values [N. Clisby and B. M. McCoy, J. Stat. Phys. 122, 15 (2006)].     

A density functional model for the binary crystal of hard spheres with vacancies

We study the stability of a binary mixture of hard spheres in the crystalline state in which a small fraction of lattice sites in the solid structure are vacant. The optimum vacancy concentration is obtained by minimizing the free energy of the inhomogeneous solid state. We use the modified weighted density functional approximation to compute the free energy. The necessary input for the theory is the thermodynamic properties of the homogeneous state of the mixture and is obtained from the solutions of the corresponding Percus-Yevick integral equations for the binary system. We compute the free energy for the crystal having two kinds of ordered structures in which (i) both the species lie in a disordered manner on a single face-centered-cubic lattice and (ii) each of the two species lie on a separate cubic lattice. Our theoretical model obtains equilibrium vacancy fraction of O(10(-5)) near the freezing point in both cases. The vacancy concentration decreases with the increase of density in both cases.     

High temperature physical adsorption of gases on solids: a comparison of two-dimensional and three-dimensional treatments

Two theories of fluid adsorption are discussed and compared. One starts from the assumption of the existance of a strictly two-dimensional adsorbed phase on the adsorbent surface. The statistical mechanics of this model is solved in the Percus-Yevick approximation for both two dimensional and bulk phases. According to the second, three-dimensional treatment, the thermodynamic functions of adsorption are defined as surface excess functions, relative to the Gibbs dividing surface. Explicit results are given for hard spheres and a Lennard-Jones 6-12 gas in contact with a wall interacting via the Lennard-Jones 3-9 potential. The dependence of excess quantities upon the choice of the Gibbs dividing surface is also discussed.     

Particle connectedness and cluster formation in sequential depositions of particles: integral-equation theory

We applied the integral-equation theory to the connectedness problem. The method originally applied to the study of continuum percolation in various equilibrium systems was modified for our sequential quenching model, a particular limit of an irreversible adsorption. The development of the theory based on the (quenched-annealed) binary-mixture approximation includes the Ornstein-Zernike equation, the Percus-Yevick closure, and an additional term involving the three-body connectedness function. This function is simplified by introducing a Kirkwood-like superposition approximation. We studied the three-dimensional (3D) system of randomly placed spheres and 2D systems of square-well particles, both with a narrow and with a wide well. The results from our integral-equation theory are in good accordance with simulation results within a certain range of densities.     

Structure of nonuniform hard sphere fluids from shifted linear truncations of functional expansions

Percus showed that approximate theories for the structure of nonuniform hard sphere fluids can be generated by linear truncations of functional expansions of the nonuniform density rho(r) about that of an appropriately chosen uniform system. We consider the most general such truncation, which we refer to as the shifted linear response (SLR) equation, where the density response rho(r) to an external field phi(r) is expanded to linear order at each r about a different uniform system with a locally shifted chemical potential. Special cases include the Percus-Yevick (PY) approximation for nonuniform fluids, with no shift of the chemical potential, and the hydrostatic linear response (HLR) equation, where the chemical potential is shifted by the local value of phi(r). The HLR equation gives exact results for very slowly varying phi(r) and reduces to the PY approximation for hard core phi(r), where generally accurate results are found. We show that a truncated expansion about the bulk density (the PY approximation) also gives exact results for localized fields that are nonzero only in a "tiny" region whose volume V(phi) can accommodate at most one particle. The SLR equation can also exactly describe a limit where the fluid is confined by hard walls to a very narrow slit. This limit can be related to the localized field limit by a simple shift of the chemical potential, leading to an expansion about the ideal gas. We then try to develop a systematic way of choosing an optimal local shift in the SLR equation for general phi(r) by requiring that the predicted rho(r) is insensitive to small variations about the appropriate local shift, a property that the exact expansion to all orders would obey. The resulting insensitivity criterion (IC) gives a theory that reduces to the HLR equation for slowly varying phi(r) and is much more accurate than HLR both for very narrow slits, where the IC agrees with exact results, and for fields confined to "tiny" regions, where the IC gives very accurate (but not exact) results. However, the IC is significantly less accurate than the PY and HLR equations for single hard core fields. Only a small change in the predicted reference density is needed to correct this remaining limit.     

Poisson's ratio of the fcc hard sphere crystal at high densities

Elastic constants and the Poisson ratio of the fcc hard-sphere crystalline phases, free of defects and with vacancies, are determined by two Monte Carlo methods: (i) the analysis of the box fluctuations in the constant pressure ensemble with variable box shape (N-P-T) and (ii) by the free-energy differentiation with respect to deformation in the fixed box ensemble (N-V-T). Very good agreement is observed for the extrapolated to the infinitely large system limit results of both the methods. The coefficients of the leading singularities of the elastic constants near close packing are estimated; they are well described by the free volume approximation. Two mechanisms influencing the Poisson ratio are studied. (i) It is shown that at high densities particle motions decrease the Poisson ratio with respect to the static case which corresponds to zero temperature. Simulations performed for systems of soft spheres, interacting through n-inverse-power potentials, r-n, show that the elastic constants of the hard spheres can be obtained in the limit n-->infinity. When T-->0 the elastic constants of the soft spheres tend to those of the static model. (ii) It is also shown that vacancies decrease C11 and C44 and increase C12 and, hence, increase the Poisson ratio with respect to the defect-free state of the system.     

Depletion potential in the infinite dilution limit

The depletion force and depletion potential between two in principle unequal "big" hard spheres embedded in a multicomponent mixture of "small" hard spheres are computed using the rational function approximation method for the structural properties of hard-sphere mixtures [S. B. Yuste, A. Santos, and M. Lopez de Haro, J. Chem. Phys. 108, 3683 (1998)]. The cases of equal solute particles and of one big particle and a hard planar wall in a background monodisperse hard-sphere fluid are explicitly analyzed. An improvement over the performance of the Percus-Yevick theory and good agreement with available simulation results are found.     

Thermodynamically consistent closure approximation for hard spheres systems

We present a new closure relation that is an extension of the Percus-Yevick approximation. In the proposed closure, we introduce an additional term and a mixing coefficient that can be determined by imposing a condition of thermodynamic self-consistency. Moreover, the mixing coefficient is calculated analytically within a linear approximation. In the case of a monodisperse system of hard spheres, we compare the results of our model to well-established thermodynamic expressions and also to the structural properties of fairly known closure approximations. In the second case, and using an equivalent scheme, the new closure relation is extended to the depletion potential between two large hard spheres immersed in a liquid of small hard spheres. In both cases, the results of our model are in good agreement with numerical simulations performed at intermediate concentrations.     

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.     

Polydisperse telechelic polymers at interfaces: analytic results and density functional theory

We use a recently developed continuum theory to expand on an exact treatment of the interfacial properties of telechelic polymers displaying Schulz-Flory polydispersity. Our results are remarkably compact and can be derived from the properties of equilibrium, ideal polymers at interfaces. A new surface adsorption transition is identified for ideal telechelic chains, wherein the central block is an equilibrium polymer. This transition occurs in the limit of strong end adsorption. Additionally, closed expressions are derived for the ideal continuum telechelic chain in contact with two large spheres, using the Derjaguin approximation. We analyze the interactions between colloids as a function of polydispersity and molecular weight, and the results are compared with polymer density functional theory in the dilute limit. Significant variations in polymer mediated forces are observed as a function of polydispersity, molecuar weight, and chain stiffness.     

Macromolecule-Induced Clustering of Hard Spheres

The connectivity Ornstein-Zernike formalism, together with the polymer reference interaction site model (PRISM), is employed to describe connectivity and network formation in mixtures of spheres and polymers. Results are presented for the percolation of spheres induced by both flexible coil-like and rigid rod-like linear polymers; the Percus-Yevick (PY) approximation is used throughout. Our results are compared with predictions based on the adhesive hard sphere (AHS) model, and correlations with the polymer-mediated second virial coefficient between spheres are discussed. Copyright 2001 Academic Press.     

Hard-sphere radial distribution function again

A theoretically based closed-form analytical equation for the radial distribution function, g(r), of a fluid of hard spheres is presented and used to obtain an accurate analytic representation. The method makes use of an analytic expression for the short- and long-range behaviors of g(r), both obtained from the Percus-Yevick equation, in combination with the thermodynamic consistency constraint. Physical arguments then leave only three parameters in the equation of g(r) that are to be solved numerically, whereas all remaining ones are taken from the analytical solution of the Percus-Yevick equation.     

Structure and ordering kinetics of micelles in triblock copolymer solutions in selective solvents

We have used small-angle x-ray scattering (SAXS), and small-angle neutron scattering (SANS) to study the micelle structure of a polystyrene-block-poly(ethene-co-butene)-block-polystyrene triblock copolymer in dilute - semidilute solutions in solvents selective for either the outer styrene block (dioxane) or for the middle block (heptane or tetradecane). Measurements of equilibrium structure factors showed that micelles were formed in both types of selective solvents. In the case of dioxane, the micelles are isolated whereas in the case of heptane or tetradecane, a bridged micellar structure may be formed at higher copolymer concentrations. In both cases we observed an ordered cubic structure of insoluble domains (micellar cores) at high concentrations (> 8 %). The micellar scattering function was fit to the Percus-Yevick interacting hard-sphere model. The temperature dependence of the core radius, the hard-sphere interaction radius and the volume fraction of hard spheres were obtained. We also used synchrotron-based time-resolved SAXS to examine the kinetics of ordering of the micelles on a cubic lattice for many different temperature jumps into the ordered cubic phase starting from the disordered micellar fluid phase in different solvents at different concentrations. The time evolution of the structure changes was determined by fitting the data with Gaussians to describe the structure factor of the ordered Bragg peaks and the Percus-Yevick structure factor was used to describe the micellar fluid. The time dependence of the peak intensities and widths as well as of the micellar parameters will be presented. The results showing the kinetics of the transformation from the fluid to the ordered phase were analyzed using the Mehl-Johnson-Avrami theory of nucleation.     

Isothermal compressibility and correlation radius near the critical point: Numerical analysis of the Ornstein-Zernike equation

Solution of the Ornstein-Zernike equation is analyzed numerically in the Percus-Yevick and hyperchain approximations for a system of Lennard-Jones particles in a critical region. The temperature dependences of correlation functions, isothermal compressibility η, and correlation radius of density fluctuations ζ are investigated at a critical density; the corresponding critical indices are determined. It is shown that the Percus-Yevick approximation yields satisfactory results when the correlation functions are calculated within a range corresponding to approximately 50 atomic (molecular) diameters. In this case, with ≈5% deviations from the critical temperature, the calculated and experimental values of η and critical indices are in good agreement. Tver State University. Translated fromZhurnal Strukturnoi Khimii, Vol. 36, No. 5, pp. 799–807, September–October, 1995. Translated by I. Izvekova     

Solvation forces in ionic and neutral liquids: A density functional approach

The solvation forces between two planar charged surfaces in ionic solutions, corresponding to charged and neutral hard spheres representing the ions and the solvent, respectively, are studied here using a weighted density functional theory for inhomogeneous Coulomb systems developed by us recently. The hard sphere contributions to the one-particle correlation function are evaluated nonperturbatively using a position-dependent effective density, while the electrical contributions are obtained through a perturbative expansion around this weighted density. The calculated results on the solvation forces between two charged hard walls compare well with available simulation results for ionic systems. For a neutral system, the present results show good agreement with the experimentally observed oscillating forces for two mica surfaces in octamethylcyclotetrasiloxane. The present approach thus provides a direct route to the calculation of interaction energies between colloidal particles.     

Percus-Yevick theory of adsorption

We study adsorption of hard sphere particles on to a plane surface with a delta function adsorption potential. The calculation takes account of exclusion via the Percus-Yevick approximation. At low and intermediate bulk adsorbate densities, both type II and type III BET adsorption isotherms can be found for the surface excess density and for the monolayer surface density. The surface excess isotherm agrees with an expansion of the exact surface excess isotherm to second order in the density. We mention some biochemical ramifications of the results.