Effective particle size from molecular dynamics simulations in fluids

We report molecular dynamics simulations designed to investigate the effective size of colloidal particles suspended in a fluid in the vicinity of a rigid wall where all interactions are defined by smooth atomic potential functions. These simulations are used to assess how the behavior of this system at the atomistic length scale compares to continuum mechanics models. In order to determine the effective size of the particles, we calculate the solvent forces on spherical particles of different radii as a function of different positions near and overlapping with the atomistically defined wall and compare them to continuum models. This procedure also then determines the effective position of the wall. Our analysis is based solely on forces that the particles sense, ensuring self-consistency of the method. The simulations were carried out using both Weeks–Chandler–Andersen and modified Lennard-Jones (LJ) potentials to identify the different contributions of simple repulsion and van der Waals attractive forces. Upon correction for behavior arising the discreteness of the atomic system, the underlying continuum physics analysis appeared to be correct down to much less than the particle radius. For both particle types, the effective radius was found to be \(\sim 0.75\sigma \), where \(\sigma \) defines the length scale of the force interaction (the LJ diameter). The effective “hydrodynamic” radii determined by this means are distinct from commonly assumed values of \(0.5\sigma \) and \(1.0\sigma \), but agree with a value developed from the atomistic analysis of the viscosity of such systems.