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.