Computer simulations of dense hard-sphere systems

We present comprehensive results of large-scale molecular dynamics and Monte Carlo simulations of systems of dense hard spheres at volume fraction φ along the disordered, metastable branch of the phase diagram from the freezing-point φf to random close packing volume φc. It is shown that many previous simulations contained deficiencies caused by crystallization and finite-size effects. We quantify the degree of local crystallization through an order parameter and study it as a function of time and initial conditions to determine the necessary conditions to obtain truly random systems. This ordering criterion is used to show that previous methods employed to ascertain the degree of randomness are inadequate. A careful study of the pressure is also carried out along the entire metastable branch. In the vicinity of the random-close packing fraction, we show that the pressure scales as (φc−φ)−γ, where γ=1 and φc=0.644±0.005. Contrary to previous studies, we find no evidence of a thermodynamic glass transition.