在VASP(Vienna Ab initio Simulation Package)中,判断计算是否收敛主要依据以下两个参数:
EDIFF:
控制电子步收敛,即前后能量差值。默认值为1E-4,建议设置为1E-5。
EDIFFG:
控制离子步收敛。如果以力作为收敛标准,则EDIFFG为负值,一般取-0.01到-0.05之间;如果以能量作为收敛标准,则EDIFFG为正值,一般取0.0001到0.001之间。
除了上述参数外,还可以通过以下方式判断收敛:
检查OUTCAR文件:
OUTCAR文件包含了每个离子步的详细计算结果,包括原子所受的力和能量变化。通过观察这些数据,可以判断计算是否收敛。特别是,可以检查能量差值是否小于预设的EDIFF值。
检查OSZICAR文件:
OSZICAR文件包含了每个单点计算的能量信息。通过观察OSZICAR文件中的能量变化,可以判断电子自洽迭代是否收敛。通常,能量差值小于EDIFF即可认为收敛。
检查最大迭代步数:
在VASP中,可以设置最大迭代步数(NELM参数)。当达到最大迭代步数时,计算会自动停止。因此,可以通过设置合理的NELM值来控制计算过程,确保计算在合理步数内完成。
示例代码
```bash
选择合适的超胞大小
ncore = 3
nspin = 1
nsym = 1
选择合适的k点网格
nk1 = 12
nk2 = 12
nk3 = 12
运行计算
vasp -incar 'NSW=0,IBRION=2' -kpoints 'automatic' -ncore $ncore -nspin $nspin -nsym $nsym -nk1 $nk1 -nk2 $nk2 -nk3 $nk3
分析结果
检查OUTCAR文件中的能量变化
grep F= OUTCAR | awk '{print $1,$5}' > conv.dat
将conv.dat文件复制到Origin中直观查看能量收敛情况
```
通过上述步骤和工具,可以有效地判断VASP计算是否收敛,并确保计算结果的准确性和可靠性。