MD2D: a python module for accurate determination of diffusion coefficient from molecular dynamics

Yunguo Li,Huaiwei Ni
DOI: https://doi.org/10.1016/j.cpc.2022.108599
IF: 4.717
2022-11-25
Computer Physics Communications
Abstract:Self-diffusion coefficient can be derived from molecular dynamics (MD) simulations by fitting the mean squared displacement (MSD) into the Einstein relation. However, the finite system size, nonfulfillment of the Brownian motion, and finite simulation time may bring in significant uncertainties that need to be estimated. We present a python module to facilitate the accurate determination of self-diffusion coefficient from the Einstein relation. We show that the ballistic stage can be clearly recognized and excluded to improve the accuracy and efficiency of self-diffusion coefficient calculation. The correct self-diffusion coefficient and its uncertainty can be conveniently obtained by taking the ensemble average of diffusion coefficients calculated at different time intervals. At the meantime, the module calculates viscosity that can correct the MD-derived self-diffusion coefficient to the thermodynamic limit. Program summary Program Title: MD2D VERSION 1.2.0 CPC Library link to program files: https://doi.org/10.17632/d2x8rw83jb.1 Code Ocean capsule: https://codeocean.com/capsule/2259048 Licensing provisions: GNU General Public License, version 3 Programming language: Python Nature of problem: Accurate determination of self-diffusion coefficient from molecular dynamics by using the Einstein relation is hampered by the finite system size. Solution method: Self-diffusion coefficient can be accurately determined from molecular dynamics simulations by fitting the Einstein relation with correct choice of parameters and correcting to the thermodynamic limit.
physics, mathematical,computer science, interdisciplinary applications
What problem does this paper attempt to address?