针对周期性体系的分子动力学模拟(包含计算所需输入文件)
通过vasp进行分子动力学的计算需要进行以下具体步骤:
1 选择一个足够大的晶胞制作成POSCAR,原子数越多越好(100个以上),一般k点较小值即可,可以是一个k点,也就是使用Gamma点。
2 可以将计算得到的CONTCAR复制到POSCAR,CONTCAR可能提供初始速度。如果没有初始速度,动力学计算会提供一个随机初始速度,这样得计算并没有问题,但是应该注意,由于初始速度的随机性,不同计算获取的运动轨迹或者相关参数会有较小出入。
3 INCAR文件的关键设置:
IBRION=0 #启动分子动力学计算
POTIM=1 #分子动力学步长,一般1fs
NSW=10000 #设置动力学步数