These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
hess2019.pdf (240.44 kB)

On The Importance of Accurate Algorithms for Reliable Molecular Dynamics Simulations

submitted on 30.12.2019, 10:44 and posted on 31.12.2019, 09:17 by Berk Hess, David van der Spoel, Mark J. Abraham, Erik Lindahl

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.



BioExcel Centre of Excellence for ComputationalBiomolecular Research

European Commission

Find out more...

SNIC 2.0: Swedish National Infrastructure for Computing

Swedish Research Council

Find out more...


Email Address of Submitting Author


KTH Royal Institute of Technology



ORCID For Submitting Author


Declaration of Conflict of Interest