Papers › Optimal estimates of diffusion coefficients from molecular dynamics simulations
Optimal estimates of diffusion coefficients from molecular dynamics simulations
Jakob Tómas Bullerjahn, Sören von Bülow, Gerhard Hummer
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
Translational diffusion coefficients are routinely estimated from molecular dynamics simulations. Linear fits to mean squared displacement (MSD) curves have become the de facto standard, from simple liquids to complex biomacromolecules. Nonlinearities in MSD curves at short times are handled with a wide variety of ad hoc practices, such as partial and piece-wise fitting of the data. Here, we present a rigorous framework to obtain reliable estimates of the diffusion coefficient and its statistical uncertainty. We also assess in a quantitative manner if the observed dynamics is indeed diffusive. By accounting for correlations between MSD values at different times, we reduce the statistical uncertainty of the estimator and thereby increase its efficiency. With a Kolmogorov-Smirnov test, we check for possible anomalous diffusion. We provide an easy-to-use Python data analysis script for the estimation of diffusion coefficients. As an illustration, we apply the formalism to molecular dynamics simulation data of pure TIP4P-D water and a single ubiquitin protein. In a companion paper [J. Chem. Phys. XXX, YYYYY (2020)], we demonstrate its ability to recognize deviations from regular diffusion caused by systematic errors in a common trajectory "unwrapping" scheme that is implemented in popular simulation and visualization software.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections