||| |
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的参数。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-7-30 08:50
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社