Conservative Algorithm for an Adaptive Change of Resolution in Mixed Atomistic / Coarse-Grained Multiscale Simulations

Understanding complex materials often requires investigating multiple, tightly coupled time and length scales. Neither atomistic nor coarse-grained simulations are often able to adequately capture all the relevant scales. To combine the efficiency of coarse-grained models with the accuracy of atomistic models for systems that require atomistic resolution only locally, for example at an interface, mixed-resolution models have been developed. These models use a coarse-grained description for the part of the system distant from an active site and atomistic description for the active site and its direct environment. Since the active zone may diffuse during a simulation, the simulation algorithm needs to permit an on-the-fly reclassication of atoms as they transition between the high- and low-resolution regimes. In this paper, we derive a conservative Hamiltonian and present an explicit symplectic integrator for mixed-resolution systems that allows for such a change in resolution of selected groups of atoms during a MD simulation.