|||
考虑一个结构稳不稳定,或者通过第一性原理分子动力学获得合理初始构型,需要使用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都是不错的。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-27 06:49
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社