PROPER STATISTICAL SAMPLING IN ISOTHERMAL DISCRETE-TIME MOLECULAR DYNAMICS

Molecular Dynamics simulations always involve a discretization of time, but the discrete-time behavior is increasingly different from that of the continuous-time physical equations as the time step is increased. This fact creates a dilemma for any simulation of a dynamical system: Use a small time step, resulting in dynamics that resemble continuous-time behavior at the expense of efficiency; or use a large time step that makes the simulation finish sooner at the expense of meaningful evolution. It is, therefore, essential to understand the features of different algorithms, such that optimal properties can be chosen for a given set of problems and objectives.

Our aim is to investigate and improve Molecular Dynamics simulation techniques for systems in thermal equilibrium. I will present our recent simple derivation [1] of a stochastic Stormer-Verlet algorithm for the evolution of Langevin equations. The method, which is as simple as conventional Verlet schemes, has been numerically tested on both low-dimensional nonlinear systems as well as more complex molecular ensembles with many degrees of freedom [2]. In light of the fundamental artifacts introduced by discrete time to dynamical simulations, I will provide a simple intuitive picture of the unique benefits of our algorithm that, unlike other algorithms, preserves proper configurational sampling (diffusion and Boltzmann distribution) in discrete time.

[1] Gronbech-Jensen & Farago, Molecular Physics Vol.111, 983 (2013).
[2] Gronbech-Jensen, Hayre, & Farago, Computer Physics Communications, Vol.185, 524 (2014).