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. [ABSTRACT FROM AUTHOR]
Copyright of Journal of Chemical Physics is the property of American Institute of Physics and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
