Fast-forward Langevin dynamics with momentum flips.

Stochastic thermostats based on the Langevin equation, in which a system is coupled to an external heat bath, are popular methods for temperature control in molecular dynamics simulations due to their ergodicity and their ease of implementation. Traditionally, these thermostats suffer from sluggish behavior in the limit of high friction, unlike thermostats of the Nosé-Hoover family whose performance degrades more gently in the strong coupling regime. We propose a simple and easy-to-implement modification to the integration scheme of the Langevin algorithm that addresses the fundamental source of the overdamped behavior of high-friction Langevin dynamics: if the action of the thermostat causes the momentum of a particle to change direction, it is flipped back. This fast-forward Langevin equation preserves the momentum distribution and so guarantees the correct equilibrium sampling. It mimics the quadratic behavior of Nosé-Hoover thermostats and displays similarly good performance in the strong coupling limit. We test the efficiency of this scheme by applying it to a 1-dimensional harmonic oscillator, as well as to water and Lennard-Jones polymers. The sampling efficiency of the fast-forward Langevin equation thermostat, measured by the correlation time of relevant system variables, is at least as good as the traditional Langevin thermostat, and in the overdamped regime, the fast-forward thermostat performs much better, improving the efficiency by an order of magnitude at the highest frictions we considered.

[1]  Generalized Langevin models and condensed-phase chemical reaction dynamics , 1982 .

[2]  P. Hünenberger.  Thermostat Algorithms for Molecular Dynamics Simulations , 2005 .

[3]  H. Trotter.  On the product of semi-groups of operators , 1959 .

[4]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[5]  Giovanni Bussi,et al.  Langevin equation with colored noise for constant-temperature molecular dynamics simulations. , 2008, Physical review letters.

[6]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[7]  Michele Parrinello,et al.  Efficient stochastic thermostatting of path integral molecular dynamics. , 2010, The Journal of chemical physics.

[8]  T. Schneider,et al.  Molecular-dynamics study of a three-dimensional one-component model for distortive phase transitions , 1978 .

[9]  M. Parrinello,et al.  Study of an F center in molten KCl , 1984 .

[10]  Thomas E. Markland,et al.  Competing quantum effects in the dynamics of a flexible water model. , 2009, The Journal of chemical physics.

[11]  Herman J. C. Berendsen,et al.  ALGORITHMS FOR BROWNIAN DYNAMICS , 1982 .

[12]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[13]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[14]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[15]  Tetsuya Morishita,et al.  Fluctuation formulas in molecular-dynamics simulations with the weak coupling heat bath , 2000 .

[16]  J. Sexton,et al.  Hamiltonian evolution for the hybrid Monte Carlo algorithm , 1992 .

[17]  Giovanni Bussi,et al.  Colored-Noise Thermostats à la Carte , 2010, 1204.0822.

[18]  M. Parrinello,et al.  Accurate sampling using Langevin dynamics. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  Artur Baumgärtner,et al.  Statics and dynamics of the freely jointed polymer chain with Lennard-Jones interaction , 1980 .

[20]  Michele Ceriotti,et al.  i-PI: A Python interface for ab initio path integral molecular dynamics simulations , 2014, Comput. Phys. Commun..

[21]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[22]  B. Leimkuhler,et al.  Robust and efficient configurational molecular sampling via Langevin dynamics. , 2013, The Journal of chemical physics.

[23]  M. Smoluchowski,et al.  MOLECULAR DYNAMICS SIMULATION OF THE LENNARD-JONES POLYMERS IN A GOOD SOLVENT , 2007 .

[24]  Michele Parrinello,et al.  Accelerating the convergence of path integral dynamics with a generalized Langevin equation. , 2011, The Journal of chemical physics.

[25]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[26]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[27]  Michele Parrinello,et al.  Stochastic thermostats: comparison of local and global schemes , 2008, Comput. Phys. Commun..