本文虽然以ReaxFF分子动力学模拟结果为例,但实际上基于BAND(DFT)、DFTB、MOPAC、经典力场、机器学习势的分子动力学结果,均可如此分析。分子动力学计算方面有如下需要注意的地方:
计算扩散系数有两种算法:
其中关于自相关函数的定义参考手册:
常规的分子动力学计算完毕(Step足够大,分子充分扩散,达到平衡)之后,在Movie中即可分析指定分子、原子的自相关函数与扩散系数:
SCM - Movie - MD Properties - Autocorrelation Functions:
生成了ACF与扩散系数,这里扩散系数为4.085 × 10-8 m2/s
注意,扩散系数是一个严重受统计样本影响的量,所以建议多个样本求均值。