Effects of skewing collision cells on transport properties in multiparticle collision dynamics simulations

Abstract

Multiparticle collision dynamics (MPCD) is a mesoscale simulation technique that uses a simplified solvent to model hydrodynamic interactions. Rather than interacting through pairwise forces, MPCD solvent particles undergo momentum-exchanging collisions within spatially localized cells according to prescribed rules. The conventional MPCD algorithm employs cubic collision cells, but this choice is not optimal for systems that are most naturally described using skewed simulation boxes. Here, we investigate the behavior of an MPCD scheme in which the collision cells are aligned with the vectors that define a triclinic (parallelepiped) simulation box. We find that skewing the collision cells has a small but statistically significant impact on the transport properties of the pure solvent. Similar but more pronounced effects are found for nearly hard spheres in solution, including a significant decrease in their nominal self-diffusion coefficient and unphysical anisotropy in their self-diffusion tensor. Thus, our analysis indicates that skewed MPCD collision cells may result in spurious behavior and should be used with caution. We posit that these artifacts may be mitigated by grid-free schemes for placing particles into collision cells.