Metadata-Version: 2.1
Name: md2d
Version: 1.2.0
Summary: A python module for accurate determination of diffusion coefficient from molecular dynamics
Home-page: https://github.com/lyg-ucl/md2d
Author: Yunguo Li
Author-email: liyungo@ustc.edu.cn
License: GNU General Public License v3 (GPLv3)
Keywords: diffusion coefficient,molecular dynamics
Platform: UNKNOWN
Classifier: Development Status :: 3 - Alpha
Requires-Python: >=3.6, <4
Description-Content-Type: text/markdown
License-File: LICENSE.txt

Installing
============

MD2D needs the dependencies:


    numpy
    scipy
    matplotlib
    mdanaysis

It can be installed locally. Under the package directory


    pip install .

It can also be installed from pypi


    pip install md2d

Usage
=====

MD2D can be used as a standalone code


    md2d

or imported as a module


    >>> import md2d
    >>> md2d.diffusion(xdatcar="XDATACR",outcar="OUTCAR")

Files
=====

Input:
------
MD2D reads the output of molecular dynamics simulations in the VASP format. Only two input files will be read: XDATCAR and OUTCAR. MD2D uses MDAnalysis to transform the outputs from other molecular dynamics codes to the VASP format.

..

    POTIM  = 1.0000    time-step for ionic-motion
    NBLOCK =      1;   KBLOCK =      1    inner block; outer block
    Mass of Ions in am
       POMASS =  55.85 
    volume of cell :     711.71
    kin. lattice  EKIN_LAT=         93.098925  (temperature 7275.16 K)
    Total+kin.     3416.08  3605.40  3585.13   -42.72  -109.39   -38.68
    
..

The keys words “POTIM”, “NBLOCK”, “Mass of Ions in am”, “number of ions“, “volume of cell”, “EKIN_LAT=”, “Total+kin.” are used for searching.

Output:
-------
Calculated diffusion coefficients, viscosities, and correction to diffusion coefficient are streamed in standard output. Besides, D-tau, mean squared displacement (MSD) and stress auto-correlation function (SACF) will be written into file as X_D_tau.dat (X is the element name), msd.dat and sacf_average.dat, and their plotsare saved as D_tau.pdf, msd.pdf and SACF.pdf. In each X_D_tau.dat, the first column is the time interval, the left three columns are the Ds for x, y and z directions. The columns in msd.dat corresponds to the MSDs for the elements in the order provided in XDATCAR. The three columns in sacf_average.dat are the MD step, bulk SACF and shear SACF.

How-To Guides
=============

Calculate diffusion coefficient
-------------------------------
Choose option 1 will lead to the calculation of diffusion coefficient. Then, MD2D will exhaust the available time intervals to calculate the diffusion coefficient as a function of time interval. The procedure and results will be streamed on the screen, and the D-tau data will be written into the file X_D_tau.dat. 

Two integer parameters will be asked for, namely, the number of initial MD steps to drop and the number of MD steps to skip for each segment. The first is to drop out the unwanted/unequilibrated initial stage of MD. The second is to exclude the ballistic motion stage from fitting the Einstein relation. The second parameter can be determined from the log-log plot or MSD.

Plot D-tau
----------
If diffusion coefficients have been calculated by doing option 1, then choose the option 4 will plot the D-tau curve.
If step one hasn’t been done, they a warning message will prompt to do option 1 firstly.

Calculate MSD
-------------
The MSD can be calculated and plotted by choosing option 3. The plot usually sees a change of trending, and the length of the first part is the length of ballistic motion that needs to be removed from fitting the Einstein relation.

Calculate viscosities and correction to diffusion coefficient
-------------------------------------------------------------
The stress auto-correlation function, viscosities, and correction to diffusion coefficient will be calculated by choosing option 2. The results will be streamed on the screen. 
The correction to diffusion coefficient is best applied to single-component system, e.g., pure Fe liquid or molecular liquid. For complex system like silicate melt, an element- or component-wise correction according to Eq. (4) in Li and Ni, (2022) should be considered.

Plot stress auto-correlation function
-------------------------------------
The stress auto-correlation function can be plotted by choosing option 5 after finishing the calculation in option 2. Otherwise, MD2D will ask for doing option 2 firstly.

License & disclaimer
====================
MD2D Copyright (C) University of Science and Technology of China. This file is part of MD2D. 
MD2D is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation version 3 of the License. 
MD2D is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. 
You should have received a copy of the GNU General Public License along with MD2D. If not, see <http: //www.gnu.org/licenses/>. 


