TY - JOUR
T1 - Simulation of nucleation in almost hard-sphere colloids: The discrepancy between experiment and simulation persists
AU - Filion, L.C.
AU - Ni, R.
AU - Frenkel, D.
AU - Dijkstra, M.
PY - 2011
Y1 - 2011
N2 - 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)]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)10.1063/1.3506838; S. Auer and D. Frenkel, Nature 409, 1020 (2001)
AB - 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)]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)10.1063/1.3506838; S. Auer and D. Frenkel, Nature 409, 1020 (2001)
U2 - 10.1063/1.3572059
DO - 10.1063/1.3572059
M3 - Article
SN - 0021-9606
VL - 134
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
IS - 13
M1 - 134901
ER -