Brownian dynamics of polydisperse colloidal hard spheres: Equilibrium structures and random close packings

Recently we presented a new technique for numerical simulations of colloidal hard-sphere systems and showed its high efficiency. Here, we extend our calculations to the treatment of both 2- and 3-dimensional monodisperse and 3-dimensional polydisperse systems (with sampled finite Gaussian size distribution of particle radii), focusing on equilibrium pair distribution functions and structure factors as well as volume fractions of random close packing (RCP). The latter were determined using in principle the same technique as Woodcock or Stillinger had used. Results for the monodisperse 3-dimensional system show very good agreement compared to both pair distribution and structure factor predicted by the Percus-Yevick approximation for the fluid state (volume fractions up to 0.50). We were not able to find crystalline 3d systems at volume fractions 0.50–0.58 as shown by former simulations of Reeet al. or experiments of Pusey and van Megen, due to the fact that we used random start configurations and no constraints of particle positions as in the cell model of Hoover and Ree, and effects of the overall entropy of the system, responsible for the melting and freezing phase transitions, are neglected in our calculations. Nevertheless, we obtained reasonable results concerning concentration-dependent long-time selfdiffusion coefficients (as shown before) and equilibrium structure of samples in the fluid state, and the determination of the volume fraction of random close packing (RCP, glassy state). As expected, polydispersity increases the respective volume fraction of RCP due to the decrease in free volume by the fraction of the smaller spheres which fill gaps between the larger particles.