Force fields for molecular dynamics simulations of charged dust particles with finite size in complex plasmas

Abstract

To exclude collisions leading to the overlap of finite-sized charged particles in molecular dynamics (MD) simulations in systems like complex (dusty) plasmas, we developed a scheme to generate a pair interaction force functionally depending on the pair distribution of particles. This ensures that the pair distribution function drops to zero at distances between particles corresponding to the point where the surfaces of two particles come into contact. The presented model has a wide utility for the description of various processes in dusty plasmas. As an example, using the non-equilibrium MD method, we performed calculations of the shear viscosity and thermal conductivity of charged spherical particles interacting via screened Yukawa potential. The role of finite particle size was found to be particularly important for small values of the coupling parameter, i.e., for gas-like and weakly correlated liquid-like systems. Furthermore, we have analyzed the shear viscosity and thermal conductivity in the presence of an external uniform magnetic field.