Langevin thermostat for robust configurational and kinetic sampling

Research output: Contribution to journalArticlepeer-review

7 Scopus citations

Abstract

We reformulate the algorithm of Grønbech-Jensen and Farago (GJF) for Langevin dynamics simulations at constant temperature. The GJF algorithm has become increasingly popular in molecular dynamics simulations because it provides robust (i.e., insensitive to variations in the time step) and accurate configurational sampling of the phase space with larger time steps than other Langevin thermostats. In the original derivation (Grønbech-Jensen and Farago, 2013), the algorithm was formulated as a velocity-Verlet type integrator with an in-site velocity variable. Here, we reformulate it as a leap frog scheme with a half-step velocity variable. In contrast to the original form, the reforumlated one also provides robust and accurate estimations of kinetic measures such as the average kinetic energy. We analytically prove that the newly presented algorithm gives the exact configurational and kinetic temperatures of a harmonic oscillator for any time step smaller than the Verlet stability limit, and use computer simulations to demonstrate the configurational and kinetic robustness of the algorithm in strongly non-linear systems. This property of the new formulation of the GJF thermostat makes it very attractive for implementation in computer simulations.

Original languageEnglish
Article number122210
JournalPhysica A: Statistical Mechanics and its Applications
Volume534
DOIs
StatePublished - 15 Nov 2019

Keywords

  • Discrete-time integration
  • Langevin thermostat
  • Molecular dynamics
  • Verlet algorithm

ASJC Scopus subject areas

  • Statistics and Probability
  • Condensed Matter Physics

Fingerprint

Dive into the research topics of 'Langevin thermostat for robust configurational and kinetic sampling'. Together they form a unique fingerprint.

Cite this