这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版后一修订版两侧同时换到之后的修订记录 | ||
adf:difussion [2022/03/10 19:22] – [通过均方位移求扩散系数] liu.jun | adf:difussion [2023/10/05 18:23] – [ReaxFF-扩散系数:通过速度自相关函数或均方位移MSD,得到分析分子、原子的扩散系数] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
======ReaxFF-扩散系数:通过速度自相关函数或均方位移MSD,得到分析分子、原子的扩散系数====== | ======ReaxFF-扩散系数:通过速度自相关函数或均方位移MSD,得到分析分子、原子的扩散系数====== | ||
- | 本文虽然以ReaxFF分子动力学模拟结果为例,但实际上BAND、DFTB、MOPAC分子动力学结果,均可如此分析。 | + | 本文虽然以ReaxFF分子动力学模拟结果为例,但实际上BAND、DFTB、MOPAC分子动力学结果,均可如此分析。分子动力学计算方面有如下需要注意的地方: |
+ | - 体系已经达到平衡,否则数据就存在平均的问题 | ||
+ | - 使用自相关函数计算扩散系数,则Sample frequency要设置的非常小,例如5;如果使用MSD计算扩散系数,则采用频率可以大一些(但并非必要),例如50 | ||
计算扩散系数有两种算法: | 计算扩散系数有两种算法: | ||
行 12: | 行 14: | ||
=====通过速度自相关函数求扩散系数===== | =====通过速度自相关函数求扩散系数===== | ||
- | 常规的分子动力学计算完毕之后,在Movie中即可分析指定分子、原子的自相关函数与扩散系数: | + | 常规的分子动力学计算完毕(Step足够大,分子充分扩散,达到平衡)之后,在Movie中即可分析指定分子、原子的自相关函数与扩散系数: |
SCM - Movie - MD Properties - Autocorrelation Functions: | SCM - Movie - MD Properties - Autocorrelation Functions: | ||
行 18: | 行 20: | ||
{{ : | {{ : | ||
- | * Step,设置分析的时间段,从多少帧到多少帧 | + | * Steps,设置分析的时间段,从多少帧到多少帧 |
* Property,选择扩散系数(同时也会生成自相关函数) | * Property,选择扩散系数(同时也会生成自相关函数) | ||
* Atoms,在左侧窗口选择要分析的分子或原子,点击右侧窗口的+按钮 | * Atoms,在左侧窗口选择要分析的分子或原子,点击右侧窗口的+按钮 | ||
行 30: | 行 32: | ||
=====通过均方位移求扩散系数===== | =====通过均方位移求扩散系数===== | ||
- | AMSMovie → MD Properties → MSD,弹出窗口中,通过设置帧数的起始值控制分析MD的时间范围 | + | AMSMovie → MD Properties → MSD,弹出窗口中: |
- | 2. In Steps set 2000 - 22001 | + | * Steps,类似地通过设置帧数的起始值控制分析MD的时间范围 |
- | 3. Set Atoms to Li | + | |
- | 4. Set Max MSD Step to 5000 (corresponding 5000*1.25 fs = 6250 fs). | + | |
+ | |||
+ | 点击Generate MSD即生成MSD函数,同时也会在曲线图中显示得到的扩散系数D: | ||
+ | |||
+ | {{ : |