Computer simulation of solid-liquid coexistence in binary hard sphere mixtures

We present the results of a computer simulation study of the solid-liquid coexistence of a binary hard sphere mixture for diameter ratios in the range 0.85 ~< ct ~< 1.00. For the solid phase we only consider substitutionally disordered FCC and HCP crystals. For 0-9425 < ct < 1"00 we find a solid-liquid coexistence curve of the 'spindle' type. For ct = 0-9425 this becomes an azeotropic and for ct = 0-875 a eutectic diagram. We compare these results with the predictions of the density functional theory of Barrat, Baus and Hansen. We observe that the density functional theory accurately predicts the point where the spindle diagram transforms into an azeotrope. However, the density functional theory differs from the simulation results on a number of counts. The greatest differences between computer simulations and theory are that the changeover from an azeotropic to a eutectic diagram is found to occur at c~ = 0-875, rather than at the predicted value of ct = 0.92, that the density difference between the solid and the liquid at liquid-solid coexistence is found to have a minimum as a function of the mole fraction of the large spheres, while density functional theory predicts a maximum, and finally that the solubility of large spheres in a solid mixture of small spheres is much larger than predicted.

[1]  M. Baus,et al.  The freezing of hard spheres: The density functional theory revisited , 1985 .

[2]  D. Frenkel,et al.  FREE ENERGY CALCULATIONS FOR SOLID SOLUTIONS BY COMPUTER SIMULATIONS , 1991 .

[3]  B. Alder,et al.  Virial Expansion for Mixtures of Rigid Spheres , 1956 .

[4]  William Hume-Rothery,et al.  The structure of metals and alloys , 1939 .

[5]  D. A. Young,et al.  Crystal Transformation for Hard Spheres , 1969 .

[6]  W. W. Wood,et al.  Equation of State of Classical Hard Spheres at High Density , 1962 .

[7]  William G. Hoover,et al.  Studies in Molecular Dynamics. V. High‐Density Equation of State and Entropy for Hard Disks and Spheres , 1968 .

[8]  W G T Kranendonk,et al.  Computer simulation of solid-liquid coexistence in binary hard-sphere mixtures , 1989 .

[9]  Z. W. Salsburg.  Erratum: Statistical‐Mechanical Properties of Small Systems in the Isothermal—Isobaric Ensemble , 1966 .

[10]  M. Baus.  Statistical mechanical theories of freezing: An overview , 1987 .

[11]  D. Frenkel,et al.  Thermodynamic properties of binary hard sphere mixtures , 1991 .

[12]  A. Haymet,et al.  Density functional theory for the freezing of 1:1 hard sphere mixtures , 1987 .

[13]  Baus,et al.  Can the thermodynamic properties of a solid be mapped onto those of a liquid? , 1990, Physical review letters.

[14]  J. S. Rowlinson,et al.  Computer simulation of mixtures of hard spheres , 1987 .

[15]  Daan Frenkel,et al.  New Monte Carlo method to compute the free energy of arbitrary solids. Application to the fcc and hcp phases of hard spheres , 1984 .

[16]  J. Hansen,et al.  Freezing of binary hard-sphere mixtures into disordered crystals: a density functional approach , 1987 .

[17]  J. Hansen,et al.  On the equilibrium structure of dense fluids , 1988 .

[18]  M. Baus,et al.  The freezing of hard spheres II: A search for structural (f.c.c.-h.c.p.) phase transitions , 1985 .

[19]  S. Rick,et al.  Density functional theory for the freezing of Lennard-Jones binary mixtures , 1989 .

[20]  Hansen,et al.  Density-functional theory of freezing of hard-sphere mixtures into substitutional solid solutions. , 1986, Physical review letters.

[21]  W. W. Wood,et al.  Molecular dynamics calculations of the hard-sphere equation of state , 1984 .

[22]  D. Oxtoby,et al.  Density functional theory for freezing of a binary hard sphere liquid , 1990 .

[23]  B. Widom,et al.  Some Topics in the Theory of Fluids , 1963 .

[24]  William G. Hoover,et al.  Melting Transition and Communal Entropy for Hard Spheres , 1968 .