标签:
VASP手册对于GGA在INCAR中的设置给了较为明确的规定:
GGA = 91 —PE —RP — PS— AM Default – , XC type is chosen according to POTCAR
也就是缺省的情况下,交换关联势的选取就是根据POTCAR来确定的,这在OUTCAR中显示为:
Exchange correlation treatment: GGA = -- GGA type LEXCH = 8 internal setting for exchange type VOSKOWN = 0 Vosko Wilk Nusair interpolation LHFCALC = F Hartree Fock is set to LHFONE = F Hartree Fock one center treatment AEXX = 0.0000 exact exchange contribution
而在相应的POTCAR中,较为关键的参数就是:LEXCH = PE,如果我们在INCAR中设置了GGA=PE,那么在相应的OUTCAR中就会出现:
Exchange correlation treatment: GGA = PE GGA type LEXCH = 8 internal setting for exchange type VOSKOWN= 0 Vosko Wilk Nusair interpolation LHFCALC = F Hartree Fock is set to LHFONE = F Hartree Fock one center treatment AEXX = 0.0000 exact exchange contribution
这是很显然的事实。在VASP手册解释为:This tag was added to perform GGA calculation with pseudopotentials generated with conventional LDA reference configurations.The tag is named GGA.并且给出了具体解释:
91 Perdew -Wang 91 PE Perdew-Burke-Ernzerhof RP revised Perdew-Burke-Ernzerhof AM AM05 (Ref. [49, 50], VASP tests see Ref. [51]) PS Perdew-Burke-Ernzerhof revised for solids (PBEsol, see Ref. [52])
论坛上对于GGA-tag的设置的讨论出现了很多说法,例如:
“先看POTCAR中LEXCH后面的值是什么,如果是CA,那么INCAR中就不用设置GGA了;如果是91,那么就在 INCAR 中输入GGA=91;如果是PE,那么就设置GGA=PE。”
既然VASP完全可以采用默认值,而且即使使用GGA,也必须和POTCAR相一致,那么要这个GGA岂不是多余吗?
非也,GGA这个标签设置在INCAR里面实际上就是为了确保我们计算所采用的势函数与计算文件夹中的POTCAR一致,如果不一致,那么就会报错。如果不设置GGA这个标签的话,我们原本想采用PBE这个交换关联势,结果copy到计算文件中的POTCAR是PBEsol,那么计算过程中也不会有任何提示,分析结果的时候,我们就会把PBEsol这个交换关联泛函误认为是PBE交换泛函。
参考链接:
[1] VASP the GUIDE(Vienna, May 30, 2014)
[2] VASP新手入门十个简单问题: http://blog.chinaunix.net/uid-10577171-id-2947532.html
[3] 如何使用LDA和GGA近似:http://emuch.net/bbs/viewthread.php?tid=2598315
[4] VASP5.2 LDA,GGA对应哪个赝势啊?http://emuch.net/bbs/viewthread.php?tid=5332301
标签:
原文地址:http://www.cnblogs.com/panscience/p/4694588.html