On The Importance of Accurate Algorithms for Reliable Molecular Dynamics Simulations

31 December 2019, Version 1
This content is a preprint and has not undergone peer review at the time of posting.


Molecular dynamics is expected to produce accurate results over a wide range of conditions and timescales. However, this is not always the case since the field has been too reluctant to abandon historically popular techniques known to introduce artefacts. Two recent papers have suggested there are reliability issues in the GROMACS code since it no longer uses a legacy twin-range algorithm. Here, we show there are order-of-magnitude differences in accuracy favoring the modern Trotter decomposition, and that a force field relying on the old algorithm will have errors parametrized into the force field. Similarly, the suggestions about incorrect virial calculations turn out to be explained by insufficient accuracy in the default SHAKE settings used for GROMOS, while the GROMACS default choices are accurate. This highlights the importance of being more critical to error cancellation in simulations in order for algorithms and parameters to both gradually converge to more perfect ones.


molecular dynamics
force field
Multiple Time Step


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.