Adaptive switching of interaction potentials in the time domain: an extended Lagrangian approach tailored to transmute force field to QM/MM simulations and back.

An extended Lagrangian formalism that allows for a smooth transition between two different descriptions of interactions during a molecular dynamics simulation is presented. This time-adaptive method is particularly useful in the context of multiscale simulation as it provides a sound recipe to switch on demand between different hierarchical levels of theory, for instance between ab initio (“QM”) and force field (“MM”) descriptions of a given (sub)system in the course of a molecular dynamics simulation. The equations of motion can be integrated straightforwardly using the usual propagators, such as the Verlet algorithm. First test cases include a bath of harmonic oscillators, of which a subset is switched to a different force constant and/or equilibrium position, as well as an all-MM to QM/MM transition in a hydrogen-bonded water dimer. The method is then applied to a smectic 8AB8 liquid crystal and is shown to be able to switch dynamically a preselected 8AB8 molecule from an all–MM to a QM/MM description ...