GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms
Computational Physics
2022-09-28 v2 Chemical Physics
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.
Cite
@article{arxiv.2204.02105,
title = {GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms},
author = {Stefanie Kieninger and Bettina G. Keller},
journal= {arXiv preprint arXiv:2204.02105},
year = {2022}
}