VASP干货 | 如何使用VASP做第一性原理分子动力学计算？

考虑一个结构稳不稳定，或者通过第一性原理分子动力学获得合理初始构型，需要使用VASP进行第一性原理分子动力学计算。

首先，结构的晶胞应该足够大，包含几十个原子，这样温度才会波动比较小。

其次，KPOINTS采用比较少的撒点，因为晶胞比较大。采用1 1 1的采样就可以。

0  0  0

第三，在INCAR里设置参数。降低EDIFF的精度，一般10-4, 10-5 eV/atom就很高了。加入分子动力学的参数：

#md

IBRION=0

NSW=100

POTIM=3

第四，注意INCAR里电子结构计算的一些参数设置，有些系统计算出来OSZICAR里温度T一直为0，经过摸索才发现与ISMEAR, SIGMA等有关：

ISYM=0

第五，从OSZICAR里提取温度和能量。

grep “T= ” OSZICAR | awk ‘{print $1 ” ” $3 ” ” $5}’ > T_E.txt

第六，一些摸索的经验。POTIM不宜取得太大，否则会导致非常大的温度T，使得程序停掉。尝试过POTIM=5，温度infinite。POTIM取1.5，2都是不错的。

本文转自胡涛科学网博客，转载目的在于知识分享,如有侵权问题，请联系我们删除。链链接地址：https://blog.sciencenet.cn/blog-907836-1237593.html

👇👇扫描二维码，立即报名👇👇