GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms

Abstract

Two of the most widely used Langevin integrators for molecular dynamics simulations are the GROMACS Stochastic Dynamics (GSD) integrator and the splitting method BAOAB. In this letter, we show that the GROMACS Stochastic Dynamics integrator is equal to the less frequently used splitting method BAOA. It immediately follows that GSD and BAOAB sample the same configurations and have the same high configurational accuracy. Our numerical results indicate that GSD/BAOA has higher kinetic accuracy than BAOAB.

References in corpus (6)

• Accurate sampling using Langevin dynamics
• Finite-temperature materials modeling from the quantum nuclei to the hot electrons regime
• Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
• Path probability ratios for Langevin dynamics -- exact and approximate
• The role of the dissipative and random forces in the calculation of the pressure of simple fluids with dissipative particle dynamics
• Bringing discrete-time Langevin splitting methods into agreement with thermodynamics