Self-diffusion coefficient of the square-well fluid from molecular dynamics simulations within the constant force approach

Abstract

We present a systematic study of the self-diffusion coefficient for a fluid of particles interacting via the square-well pair potential by means of molecular dynamics simulations in the canonical (N, V, T) ensemble. The discrete nature of the interaction potential is modeled by the constant force approximation, and the self-diffusion coefficient is determined for several fluid densities at supercritical thermodynamic states. The dependence of the self-diffusion coefficient on the potential range λ is analyzed in the range of 1.1 ≤ λ ≤ 1.5. The obtained simulation results are in agreement with the self-diffusion coefficient predicted by the Enskog method. Additionally, we show that the diffusion coefficient is very sensitive to the potential range λ. Our results for the self-diffusion coefficient times density extrapolate well to the values in the zero-density limit obtained from the Chapman-Enskog theory for dilute gases. The constant force approximation used in this work to model the discrete pair potentials has shown to be an excellent scheme to compute the transport properties of square-well fluids using molecular dynamics simulations. Finally, the simulation results presented here are useful for improving theoretical approaches, such as the Enskog method.