Multibody Molecular Dynamics I: Theoretical Development
This is the first paper in a series of two papers on using multibody dynamics algorithms and methods for coarse-grained molecular dynamics simulations. This paper presents the underlying framework for multi-scale modelling of biomolecules and polymers. In this framework, the system to be simulated is sub-structured into a hierarchy of multi-resolution models that are simulated using efficient multibody dynamics algorithms. The algorithms work in a unified framework, enabling efficient multi-scale (or multi-resolution) simulations. A discussion of the hierarchy of models with different resolutions along with the salient features of the appropriate multibody dynamics algorithms used for simulating them is presented. The unified scheme and the qualitative advantages of the method are discussed. Important implementation details such as boundary conditions, temporal integration schemes, interaction force field calculations and solvent models are also presented. In the next paper applications and results are discussed.