1) ISMEAR:

A) 对于半导体和绝缘体体系,ISMEAR的值取绝对不能大于0, 一般用0,

B)对所有体系,更加精确的时候用-5,

C)对于所有K点数目小于3的计算,不能用ISMEAR = -5。会报错。

D) K 点少,还是半导体或者绝缘体,那么只能用 ISMEAR = 0。

E) 对于金属来说,ISMEAR的取值一般为1。

F) 保守地说,ISMEAR = 0(Gaussian Smearing) 可以满足大部分的体系(金属,导体,半导体,分子),

G) 在DOS能带计算中,使用ISMEAR= -5 用于获取精确的信息。

2) SIGMA:

SIGMA的取值和ISMEAR息息相关:

A) 如果用了ISMEAR = -5 (对于所有体系),SIGMA的值可以忽略,也可以不管(VASP会自动略过)

B)对于金属: ISMEAR =1 或者0,非金属: ISMEAR= 0 的时候,一般取0.10 或者采用默认值 0.20 即可。

不过很多人都不放心,为了保险想测试一番。标准是: SIGMA 的取值要保证OUTCAR 中的 entropy T*S 这一项,平均到每个原子上,要小于 1-2 meV。

检查的时候可以用这个命令: grep 'entropy T'  OUTCAR   (这里有两个单引号,不要输入的时候漏掉) 得出的能量除以体系中原子的数目,然后再和 0.001 eV 比较,如果小于,SIGMA取值OK,如果大于,再换个测试。

C) 对于分子,原子体系(也就是你把分子或者原子放到一个box里面),必须要用很小的值,这里我们用的0.01。

3) IBRION:

一般来说，优化结构的时候有3个选择：

IBRION=3：你的初始结构很差的时候；

IBRION=2：共轭梯度算法，很可靠的一个选择，一般来说都用它。

IBRION=1：用于小范围内稳定结构的搜索。

如果你的体系遇到结构不收敛的时候，可以尝试下换下IBRION的参数。