Structures and properties of hard sphere mixtures based on a self-consistent integral equation

Structural and thermodynamic properties of hard-sphere fluids

This Perspective article provides an overview of some of our analytical approaches to the computation of the structural and thermodynamic properties of single-component and multicomponent hard-sphere fluids. For the structural properties, they yield a thermodynamically consistent formulation, thus improving and extending the known analytical results of the Percus-Yevick theory. Approximate expressions linking the equation of state of the single-component fluid to the one of the multicomponent mixtures are also discussed.

Towards composite spheres as building blocks for structured molecules

In order to design a flexible molecular model that mimics the chemical moieties of a polyatomic molecule, we propose the 'composite-sphere' model that can assemble the essential elements to produce the structure of the target molecule. This is likened to the polymerization process where monomers assemble to form the polymer. The assemblage is built into the pair interaction potentials which can 'react' (figuratively) with selective pieces into various bonds. In addition, we preserve the spherical symmetries of the individual pair potentials so that the isotropic Ornstein-Zernike equation (OZ) for multi-component mixtures can be used as a theoretical framework. We first test our approach on generating a dumbbell molecule. An equimolar binary mixture of hard spheres and square-well spheres are allowed to react to form a dimer. As the bond length shrinks to zero, we create a site-site model of a Janus-like molecule with a repulsive moiety and an attractive moiety. We employ the zero-separation (ZSEP) closure to solve the OZ equations. The structure and thermodynamic properties are calculated at three isotherms and at several densities and the results are compared with Monte Carlo simulations. The close agreement achieved demonstrates that the ZSEP closure is a reliable theory for this composite-sphere fluid model.

Low-temperature and high-temperature approximations for penetrable-sphere fluids: comparison with Monte Carlo simulations and integral equation theories

The two-body interaction in dilute solutions of polymer chains in good solvents can be modeled by means of effective bounded potentials, the simplest of which being that of penetrable spheres (PSs). In this paper we construct two simple analytical theories for the structural properties of PS fluids: a low-temperature (LT) approximation, that can be seen as an extension to PSs of the well-known solution of the Percus-Yevick (PY) equation for hard spheres, and a high-temperature (HT) approximation based on the exact asymptotic behavior in the limit of infinite temperature. Monte Carlo simulations for a wide range of temperatures and densities are performed to assess the validity of both theories. It is found that, despite their simplicity, the HT and LT approximations exhibit a fair agreement with the simulation data within their respective domains of applicability, so that they complement each other. A comparison with numerical solutions of the PY and the hypernetted-chain approximations is also carried out, the latter showing a very good performance, except inside the core at low temperatures.

Radial distribution function of penetrable sphere fluids to the second order in density

The simplest bounded potential is that of penetrable spheres, which takes a positive finite value epsilon if the two spheres are overlapped, being zero otherwise. In this paper we derive the cavity function to second order in density and the fourth virial coefficient as functions of T* identical with k(B)T/epsilon (where k(B is the Boltzmann constant and T is the temperature) for penetrable sphere fluids. The expressions are exact, except for the function represented by an elementary diagram inside the core, which is approximated by a polynomial form in excellent agreement with accurate results obtained by Monte Carlo integration. Comparison with the hypernetted-chain (HNC) and Percus-Yevick (PY) theories shows that the latter is better than the former for T* < or similar to 1 only. However, even at zero temperature (hard sphere limit), the PY solution is not accurate inside the overlapping region, where no practical cancellation of the neglected diagrams takes place. The exact fourth virial coefficient is positive for T* < or similar to 0.73, reaches a minimum negative value at T* approximately 1.1, and then goes to zero from below as 1/T(*4) for high temperatures. These features are captured qualitatively, but not quantitatively, by the HNC and PY predictions. In addition, in both theories the compressibility route is the best one for T* < or similar to 0.7, while the virial route is preferable if T* > or similar to 0.7.

Thermodynamic and dynamic anomalies for a three-dimensional isotropic core-softened potential

Using molecular-dynamics simulations and integral equations (Rogers-Young, Percus-Yevick, and hypernetted chain closures) we investigate the thermodynamics of particles interacting with continuous core-softened intermolecular potential. Dynamic properties are also analyzed by the simulations. We show that, for a chosen shape of the potential, the density, at constant pressure, has a maximum for a certain temperature. The line of temperatures of maximum density (TMD) was determined in the pressure-temperature phase diagram. Similarly the diffusion constant at a constant temperature, D, has a maximum at a density rho(max) and a minimum at a density rho(min) < rho(max). In the pressure-temperature phase diagram the line of extrema in diffusivity is outside of the TMD line. Although this interparticle potential lacks directionality, this is the same behavior observed in simple point charge/extended water.

Structure of penetrable-rod fluids: Exact properties and comparison between Monte Carlo simulations and two analytic theories

Bounded potentials are good models to represent the effective two-body interaction in some colloidal systems, such as the dilute solutions of polymer chains in good solvents. The simplest bounded potential is that of penetrable spheres, which takes a positive finite value if the two spheres are overlapped, being 0 otherwise. Even in the one-dimensional case, the penetrable-rod model is far from trivial, since interactions are not restricted to nearest neighbors and so its exact solution is not known. In this paper the structural properties of one-dimensional penetrable rods are studied. We first derive the exact correlation functions of the penetrable-rod fluids to second order in density at any temperature, as well as in the high-temperature and zero-temperature limits at any density. It is seen that, in contrast to what is generally believed, the Percus-Yevick equation does not yield the exact cavity function in the hard-rod limit. Next, two simple analytic theories are constructed: a high-temperature approximation based on the exact asymptotic behavior in the limit T--> infinity and a low-temperature approximation inspired by the exact result in the opposite limit T--> 0. Finally, we perform Monte Carlo simulations for a wide range of temperatures and densities to assess the validity of both theories. It is found that they complement each other quite well, exhibiting a good agreement with the simulation data within their respective domains of applicability and becoming practically equivalent on the borderline of those domains. A comparison with numerical solutions of the Percus-Yevick and the hypernetted-chain approximations is also carried out. Finally, a perspective on the extension of our two heuristic theories to the more realistic three-dimensional case is provided.

Structure of highly asymmetric hard-sphere mixtures: An efficient closure of the Ornstein-Zernike equations

A simple modification of the reference hypernetted chain (RHNC) closure of the multicomponent Ornstein-Zernike equations with bridge functions taken from Rosenfeld's hard-sphere bridge functional is proposed. Its main effect is to remedy the major limitation of the RHNC closure in the case of highly asymmetric mixtures--the wide domain of packing fractions in which it has no solution. The modified closure is also much faster, while being of similar complexity. This is achieved with a limited loss of accuracy, mainly for the contact value of the big sphere correlation functions. Comparison with simulation shows that inside the RHNC no-solution domain, it provides a good description of the structure, while being clearly superior to all the other closures used so far to study highly asymmetric mixtures. The generic nature of this closure and its good accuracy combined with a reduced no-solution domain open up the possibility to study the phase diagram of complex fluids beyond the hard-sphere model.

Potential of mean force in confined colloids: Integral equations with fundamental measure bridge functions

The potential of mean force for uncharged macroparticles suspended in a fluid confined by a wall or a narrow pore is computed for solvent-wall and solvent-macroparticle interactions with attractive forces. Bridge functions taken from Rosenfeld's density-functional theory are used in the reference hypernetted chain closure of the Ornstein-Zernike integral equations. The quality of this closure is assessed by comparison with simulation. As an illustration, the role of solvation forces is investigated. When the "residual" attractive tails are given a range appropriate to "hard sphere-like" colloids, the unexpected role of solvation forces previously observed in bulk colloids is confirmed in the confinement situation.