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

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

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

KPOINTS(Gamma-only)
0
Gamma
1  1  1
0  0  0

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

#md
IBRION=0
NSW=100
NWRITE=0
TEBEG=300
TEEND=1300
SMASS=3
NBLOCK=50
POTIM=3

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

ISMEAR = 0 ; SIGMA = 0.05
ALGO=Fast
#POTIM = 0.10
#Wavefunction and charge
LWAVE = FALSE
LCHARG = FALSE
#Target Pressure
#PSTRESS = 0.01
LREAL=Auto
NCORE=4
LVDW=.TRUE.
ISYM=0

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

grep "T= " OSZICAR | awk '{print $1 " " $3 " " $5}' > T_E.txt

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

转载本文请联系原作者获取授权，同时请注明本文来自胡涛科学网博客。

链接地址：https://wap.sciencenet.cn/blog-907836-1237593.html?mobile=1

收藏

分享