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

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...

Celý popis

Uloženo v:
Podrobná bibliografie
Vydáno v:Computer physics communications Ročník 284; s. 108599
Hlavní autoři: Li, Yunguo, Ni, Huaiwei
Médium: Journal Article
Jazyk:angličtina
Vydáno: Elsevier B.V 01.03.2023
Témata:
Einstein relation
Molecular dynamics
Self-diffusion coefficient
Size dependence
Molecular dynamics
Size dependence
Einstein relation
Self-diffusion coefficient
ISSN:0010-4655, 1879-2944
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
Popis
Shrnutí: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 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.
ISSN:0010-4655
1879-2944
DOI:10.1016/j.cpc.2022.108599