GSHMC: An efficient method for molecular simulation. The hybrid Monte Carlo (HMC) method is a popular and rigorous method for sampling from a canonical ensemble. The HMC method is based on classical molecular dynamics simulations combined with a Metropolis acceptance criterion and a momentum resampling step. While the HMC method completely resamples the momentum after each Monte Carlo step, the generalized hybrid Monte Carlo (GHMC) method can be implemented with a partial momentum refreshment step. This property seems desirable for keeping some of the dynamic information throughout the sampling process similar to stochastic Langevin and Brownian dynamics simulations. It is, however, ultimate to the success of the GHMC method that the rejection rate in the molecular dynamics part is kept at a minimum. Otherwise an undesirable Zitterbewegung in the Monte Carlo samples is observed. In this paper, we describe a method to achieve very low rejection rates by using a modified energy, which is preserved to high-order along molecular dynamics trajectories. The modified energy is based on backward error results for symplectic time-stepping methods. The proposed generalized shadow hybrid Monte Carlo (GSHMC) method is applicable to NVT as well as NPT ensemble simulations.
Keywords for this software
References in zbMATH (referenced in 9 articles , 1 standard article )
Showing results 1 to 9 of 9.
- Escribano, Bruno; Akhmatskaya, Elena; Mujika, Jon I.: Combining stochastic and deterministic approaches within high efficiency molecular simulations (2013)
- Arizumi, Nana; Bond, Stephen D.: On the estimation and correction of discretization error in molecular dynamics averages (2012)
- Bou-Rabee, Nawaf; Vanden-Eijnden, Eric: A patch that imparts unconditional stability to explicit integrators for Langevin-like equations (2012)
- Lelièvre, Tony; Rousset, Mathias; Stoltz, Gabriel: Langevin dynamics with constraints and computation of free energy differences (2012)
- Bou-Rabee, Nawaf; vanden-Eijnden, Eric: Pathwise accuracy and ergodicity of metropolized integrators for SDEs (2010)
- Akhmatskaya, Elena; Bou-Rabee, Nawaf; Reich, Sebastian: A comparison of generalized hybrid Monte Carlo methods with and without momentum flip (2009)
- Akhmatskaya, Elena; Bou-Rabee, Nawaf; Reich, Sebastian: Erratum to “A comparison of generalized hybrid Monte Carlo methods with and without momentum flip” (2009)
- Leimkuhler, Benedict; Reich, Sebastian: A Metropolis adjusted Nosé-Hoover thermostat (2009)
- Akhmatskaya, Elena; Reich, Sebastian: GSHMC: An efficient method for molecular simulation (2008)