Molecular dynamics simulations were carried out to investigate the diffusion behavior of penetrable-sphere model fluids characterized by a finite energy barrier ε. The self-diffusion coefficient was evaluated from the time-dependent velocity autocorrelation function and mean-square displacement. Detailed insights into the cluster formation for penetrable spheres were gained from the Enskog factor, the effective particle volume fraction, the mean free path, and the collision frequency for both the soft-type penetrable and the hard-type reflective collisions. The simulation data were compared to theoretical predictions from the Boltzmann kinetic equation and from a simple extension to finite ε of the Enskog prediction for impenetrable hard spheres (ε→∞). A reasonable agreement between theoretical and simulation results was found in the cases of ε - ε/kT = 0.2, 0.5 and 1.0. However, for dense systems (packing fraction >0.6) with a highly repulsive energy barrier (ε - = 3.0), a poorer agreement was observed due to metastable static effects of clustering formation and dynamic effects of correlated collision processes among these cluster-forming particles.

Molecular Dynamics Simulation Study of Self-Diffusion for Penetrable-Sphere Model Fluids. Suh, S.H., Kim, C.H., Kim, S.C., Santos, A.: Physical Review E, 2010, 82[5], 051202