VASP使用之AIMD计算
VASP做AIMD计算时,可以通过讲体系原子的初始速度都追加到POSCAR后面,具体参考CONTCAR的格式。在静力学计算过程输出的CONTCAR文件中,原子坐标后面追加的一行行的数字0.000…其实就是每个原子所对应的速度。
AIMD计算时INCAR的关键设置参数:
IBRION=0 #AIMD计算时该参数设置为0
POTIM=0.5 #时间步长控制参数,单位是fs
NSW=1000 #控制弛豫步数
TEBEG=300 #初始温度控制开关,单位为K
TEEND=300 #弛豫最终温度控制开关,默认值与TEBEG相同
系综设置主要通过 MDALGO 和ISIF控制,VASP目前支撑的系综有NVT,NPT以及NPH系综具体设置如下表:
此外对于原子间的键长键角还可以通过输入文件ICONST来控制,具体参考官网说明: https://www. vasp.at/wiki/index.php/ ICONST
最后结果分析的主要文件:每一步的温度、能量信息在 OSZICAR 中,XDATCAR为分子动力学计算过程中的轨迹文件。
分子动力学示例:
祝大家科研愉快!!!
公众号推荐:计算运维鸟