An efficient parallel algorithm for non-equilibrium molecular dynamics simulations of very large systems in planar Couette flow

In this paper we present an efficient parallel domain decomposition algorithm for non-equilibrium molecular dynamics (NEMD) simulations of large systems under planar Couette flow. We propose a modified deforming cell method that permits NEMD simulations with negligible penalties due to the Lees-Edwards periodic boundary conditions. The algorithm was used to study large systems of the Weeks-Chandler-Andersen fluid in order to obtain better viscosity results at the low shear rate regimes where the signal-to-noise ratio is very small.