Simulation of nucleation in almost hard-sphere colloids: The discrepancy between experiment and simulation persists

In this paper we examine the phase behavior of the Weeks–Chandler–Andersen (WCA) potential with βε = 40. Crystal nucleation in this model system was recently studied by Kawasaki and Tanaka [Proc. Natl. Acad. Sci. U.S.A. 107, 14036 (2010)] https://doi.org/10.1021/pr100656u, who argued that the computed nucleation rates agree well with experiment, a finding that contradicted earlier simulation results. Here we report an extensive numerical study of crystallization in the WCA model, using three totally different techniques (Brownian dynamics, umbrella sampling, and forward flux sampling). We find that all simulations yield essentially the same nucleation rates. However, these rates differ significantly from the values reported by Kawasaki and Tanaka and hence we argue that the huge discrepancy in nucleation rates between simulation and experiment persists. When we map the WCA model onto a hard-sphere system, we find good agreement between the present simulation results and those that had been obtained for hard spheres [L. Filion, M. Hermes, R. Ni, and M. Dijkstra, J. Chem. Phys. 133, 244115 (2010) https://doi.org/10.1063/1.3506838; S. Auer and D. Frenkel, Nature 409, 1020 (2001) https://doi.org/10.1038/35059035].