No abstract

We report on a computer simulation and integral equation study of a simple model of patchy spheres, each of whose surfaces is decorated with two opposite attractive caps, as a function of the fraction χ of covered attractive surface. The simple model explored -the two-patch Kern-Frenkel model -interpolates between a square-well and a hard-sphere potential on changing the coverage χ. We show that integral equation theory provides quantitative predictions in the entire explored region of temperatures and densities from the square-well limit χ = 1.0 down to χ ≈ 0.6. For smaller χ, good numerical convergence of the equations is achieved only at temperatures larger than the gasliquid critical point, where however integral equation theory provides a complete description of the angular dependence. These results are contrasted with those for the one-patch case. We investigate the remaining region of coverage via numerical simulation and show how the gas-liquid critical point moves to smaller densities and temperatures on decreasing χ. Below χ ≈ 0.3, crystallization prevents the possibility of observing the evolution of the line of critical points, providing the angular analog of the disappearance of the liquid as an equilibrium phase on decreasing the range for spherical potentials. Finally, we show that the stable ordered phase evolves on decreasing χ from a threedimensional crystal of interconnected planes to a two-dimensional independent-planes structure to a one-dimensional fluid of chains when the one-bond-per-patch limit is eventually reached.

No abstract

We study the thermodynamic and structural properties of a simple, one-patch fluid model using the reference hypernetted-chain (RHNC) integral equation and specialized Monte Carlo simulations. In this model, the interacting particles are hard spheres, each of which carries a single identical, arbitrarily-oriented, attractive circular patch on its surface; two spheres attract via a simple squarewell potential only if the two patches on the spheres face each other within a specific angular range dictated by the size of the patch. For a ratio of attractive to repulsive surface of 0.8, we construct the RHNC fluid-fluid separation curve and compare with that obtained by Gibbs ensemble and grand canonical Monte Carlo simulations. We find that RHNC provides a quick and highly reliable estimate for the position of the fluid-fluid critical line. In addition, it gives a detailed (though approximate) description of all structural properties and their dependence on patch size.

The Percus-Yevick, hypernetted-chain, and "pressure-consistent' integral equations have been solved, using numerical Hankel transforms, for a fluid of two-dimensional hard cores. The thermodynamic quantities obtained from these solutions are presented and compared among themselves and with the results of other theories; a comparison of computed pair distribution functions g with a Monte Carlo g is also presented. The Percus-Yevick equation is found to give the best over-all results.

We present an extensive integral equation study of a simple point charge model of water for a variety of thermodynamic states ranging from the vapor phase to the undercooled liquid. The calculations are carried out in the molecular reference-hypernetted chain approximation and the results are compared with extensive molecular dynamics simulations. Use of a hard sphere fluid as a reference system to provide the input reference bridge function leads to relatively good thermodynamics. However, at low temperatures the computed microscopic structure shows deficiencies that probably stem from the lack of orientational dependence in this bridge function. This is in marked contrast with results previously obtained for systems that, although similarly composed of angular triatomic molecules, do not tend to the tetrahedral coordinations that are characteristic of water.

We model a disordered planar monolayer of paramagnetic spherical particles, or ferrofluid, as a two-dimensional fluid of hard spheres with embedded three-dimensional magnetic point dipoles. This model, in which the orientational degrees of freedom are three dimensional while particle positions are confined to a plane, can be taken as a crude representation of a colloidal suspension of superparamagnetic particles confined in a water/air interface, a system that has recently been studied experimentally. In this paper, we propose an Ornstein-Zernike integral equation approach capable of describing the structure of this highly inhomogeneous fluid, including the effects of an external magnetic field. The method hinges on the use of specially tailored orthogonal polynomials whose weight function is precisely the one-particle distribution function that describes the surface- and field-induced anisotropy. The results obtained for various particle densities and external fields are compared with Monte Carlo simulations, illustrating the capability of the inhomogeneous Ornstein-Zernike equation and the proposed solution scheme to yield a detailed and accurate description of the spatial and orientational structure for this class of systems. For comparison, results from density-functional theory in the modified mean-field approximation are also presented; this latter approach turns out to yield at least qualitatively correct results.

No abstract