文章来源:“分子动力学”公众号
链接:https://mp.weixin.qq.com/s/TGU_0H2jktEwXIkqCdG3ug
vasp计算中用到的三种赝势:模守恒赝势,超软赝势,PAW赝势(按产生顺序)。
按方法不同分为USPP(ultrasoft pesudopotential,超软赝势)和PAW(projector augmented wave,投影缀加平面波),两种方法都可以相当程度地减少过渡金属或第一行元素的每个原子所必需的平面波数量。
按交换关联函数不同分为LDA(local density approximation,局域密度近似)和GGA(generalized-gradient approximation,广义梯度近似),GGA又分为PW91(Perdew -Wang 91)和PBE(Perdew-Burke-Ernzerhof)。
求解K-S方程时电子之间的交换关联泛函取局域密度泛函LDA(Rc处电子结构当成是密度相同的均匀电子气体来算)。
广义梯度泛函GGA(考虑了密度的梯度变化),但泛函的不同也会引起贋势的不同。
paw文件夹:PAW-LDA
paw_gga文件夹:PAW-GGA-PW91
paw_pbe文件夹:PAW-GGA-PBE
pot文件夹:USPP-LDA
pot_GGA文件夹:USPP-GGA
选择某个目录进去,我们还会发现对应每种元素往往还会有多种赝势存在。这是因为根据ENMAX的大小还可以分为 Ga,Ga_s(soft),Ga_h,或者根据处理半芯态的不同还可以分为Ga,Ga_sv(s电子作为半芯态),Ga_pv的不同。
对于化合物(不同原子半径的元素混合)来说,PAW赝势比超软赝势精确度高。
US型赝势所需截至能较小,计算速度快,PAW赝势截至能通常较大,而且考虑的电子数多,计算慢,但精确度高。
在做具体的计算之前,针对需要计算的性质选取不同的贋势做下测试。对于不同贋势得到的结果,关注其最终的能量差,而不是绝对能量。