Calculating diffusivity: Difference between revisions
Jump to navigation
Jump to search
(Created page with "<P ALIGN="CENTER"> <FONT SIZE="+3" color="darkred"><STRONG> Diffusivity calculation in MD++</STRONG></font></P> <DIV> <P ALIGN="CENTER"><STRONG>Yanming Wang</STRONG></P> </DIV...") |
|||
| (9 intermediate revisions by the same user not shown) | |||
| Line 4: | Line 4: | ||
<DIV> |
<DIV> |
||
<P ALIGN="CENTER"><STRONG>Yanming Wang</STRONG></P> |
<P ALIGN="CENTER"><STRONG>Yanming Wang</STRONG></P> |
||
</DIV> |
|||
<P ALIGN="CENTER"><STRONG>Sep 2015</STRONG></P> |
<P ALIGN="CENTER"><STRONG>Sep 2015</STRONG></P> |
||
</DIV> |
|||
This tutorial describes how to compute diffusivity in MD++ based on the random walk model. Here as an example, we will provide MD++ Tcl and Matlab scripts for calculating Si and Au diffusivities in a Au-Si binary liquid alloy system. |
|||
== Theoretical background == |
|||
We can assume the motion of each individual atom follows the random walk model. Thus we can find the following relation between the mean square displacement (MSD) <math> \sigma </math> and the diffusivity <math>D</math>. |
|||
{|border="0" align="center" |
|||
|<math> D = \frac{\sigma}{6t}</math> |
|||
|} |
|||
This is the definition of MSD <math> \sigma = <r(t) - r_0> </math>, which can be easily obtained from the trajectories of the atoms in MD. |
|||
Latest revision as of 08:21, 17 September 2015
Diffusivity calculation in MD++
Yanming Wang
Sep 2015
This tutorial describes how to compute diffusivity in MD++ based on the random walk model. Here as an example, we will provide MD++ Tcl and Matlab scripts for calculating Si and Au diffusivities in a Au-Si binary liquid alloy system.
Theoretical background
We can assume the motion of each individual atom follows the random walk model. Thus we can find the following relation between the mean square displacement (MSD) and the diffusivity .
This is the definition of MSD , which can be easily obtained from the trajectories of the atoms in MD.
