Coarse-grained rigid blob model for soft matter simulations

Abstract

We have developed a coarse-grained multiscale molecular simulation method for soft matter systems that directly incorporates stereochemical information. We divide the material into disjoint groups of atoms or particles that move as separate rigid bodies; we call these groups "rigid blobs," hence the name coarse-grained rigid blob model. The method is enabled by the construction of transferable interblob potentials that approximate the net intermolecular interactions, as obtained from ab initio electronic structure calculations, other all-atom empirical potentials, experimental data, or any combination of the above. We utilize a multipolar expansion to obtain the interblob potential-energy functions. The series, which contains controllable approximations that allow us to estimate the errors, approaches the original intermolecular potential as the number of terms increases. Using a novel numerical algorithm, we can calculate the interblob potentials very efficiently in terms of a few interaction moment tensors. This reduces the labor well beyond what is required in standard molecular-dynamics calculations and allows large-scale simulations for temporal scales commensurate with characteristic times of nano- and mesoscale systems. A detailed derivation of the formulas is presented, followed by illustrative applications to several systems showing that the method can effectively capture realistic microscopic details and can easily extend to large-scale simulations.