Adaptive partitioning molecular dynamics using an extended Hamiltonian approach Available to Purchase

A recently proposed extended Hamiltonian approach to switching interaction potentials is generalized to enable adaptive partitioning molecular dynamics simulations. Switching is performed along a fictitious classical degree of freedom whose value determines the mixing ratio of the two potentials on a time scale determined by its associated mass. We propose to choose this associated fictitious mass adaptively so as to ensure a constant time scale for all switching processes. For different model systems, including a harmonic oscillator and a Lennard-Jones fluid, we investigate the window of switching time scales that guarantees the conservation of the extended Hamiltonian for a large number of switching events. The methodology is first applied in the microcanonical ensemble and then generalized to the canonical ensemble using a Nosé–Hoover chain thermostat. It is shown that the method is stable for thousands of consecutive switching events during a single simulation, with constant temperature and a conserved extended Hamiltonian. A slight modification of the original Hamiltonian is introduced to avoid accumulation of small numerical errors incurred after each switching process.