Comparison of structure and transport properties of concentrated hard and soft sphere fluids

Abstract

Using Newtonian and Brownian dynamics simulations, the structural and transport properties of hard and soft spheres have been studied. The soft spheres were modeled using inverse power potentials (V～r−n, with 1/n the potential softness) Although, at constant density, the pressure, diffusion coefficient, and viscosity depend on the particle softness up to extremely high values of n, we show that scaling the density with the freezing point for every system effectively collapses these parameters for n≥ 18 (including hard spheres) or large densities. At the freezing points, the long range structure of all systems is identical, when length is measured in units of the interparticle distance, but differences appear at short distances (due to the different shapes of the interaction potential). This translates into differences at short times in the velocity and stress autocorrelation functions, although they concur to give the same value of the corresponding transport coefficient (for the same density to freezing ratio); the microscopic dynamics also affects the short time behavior of the correlation functions and absolute values of the transport coefficients, but the same scaling with the freezing density works for Newtonian or Brownian dynamics. For hard spheres, the short time behavior of the stress autocorrelation function has been studied in detail, confirming quantitatively the theoretical forms derived for it.