VASP_AIMD+VASPKIT计算含温力学性质

材料的力学性质一直是DFT计算的重要方向，笔者在以往已有针对于静态结构的力学性质诸如弹性常数的相关计算，同时可通过VASPKIT借助相关方程导出力学性能。

bash+vasp+vaspkit能量应变计算弹性常数

vaspkit计算弹性常数的对称性指定

vaspkit计算弹性常数脚本补充

本文依据VASPKIT的EXAMPLE，简单介绍使用AIMD方法和VASPKIT软件计算含温力学性能。

材料选择为面心立方结构的Al为模型，空间群为 Fm-3m。

注意计算力学相关性质需要使用晶体学晶胞，而非原胞。

Generated by VASPKIT code 1.000000    4.0389299999999997    0.0000000000000000    0.0000000000000000    0.0000000000000000    4.0389299999999997    0.0000000000000000    0.0000000000000000    0.0000000000000000    4.0389299999999997   Al     4Direct    0.0000000000000000    0.0000000000000000    0.0000000000000000    Al001    0.0000000000000000    0.5000000000000000    0.5000000000000000    Al002    0.5000000000000000    0.0000000000000000    0.5000000000000000    Al003    0.5000000000000000    0.5000000000000000    0.0000000000000000    Al004

    在计算AIMD时首先进行扩胞，这里使用2×2×2的超胞结构进行下一步计算。

    手动生成VASPKIT计算所需要的设置文件INPUT.in，内容如下，参数详细说明附在后面，可根据需要自行更改。

1                                                              ! 1 for prep-rocessing, 3 for finite-temperature MD3D                                                             ! 2D for slab, 3D for bulk4