(一)对于不知道系统是绝缘体、半导体还是金属,使用ISMEAR=0 和SIGMA=0.03-0.05;(二)半导体或...
PREC、EDIFF、ENCUT根据自己的计算需求来定,ISTART= 0。ISMEAR的设置需要注意一下,对于半导体和绝缘体设置为0,对于金属设置为1,同时SIGMA设置一个合适的值,笔者一般取SIGMA = 0.01进行计算。如果不知道如何取SIGMA值,可以查找侯柱峰老师的测试脚本进行测试。如果需要获得芯能级、真空能级,则需要再INCAR中添加一...
熵(entropy)越小越好,选择 entropy T*S EENTRO 值中最小的那个所对应的 SIGMA 。( SIGMA should be as large as possible keeping the difference between the free energy and the total energy (i.e. the term 'entropy T*S') in the OUTCAR file negligible (1meV/atom). ) (l六)注意: 1、当 k...
一般情况N取1和2就好,而且大多数情况,N=1和2给出的结果很接近。 注意:进行任何的静态计算或态密度计算,且k点数目大于4时,取ISMEAR=-5;当由于原胞较大而k点数目较少(小于4个)时,取ISMEAR=0,并设置一个合适的SIGMA值。另外对半导体或绝缘体的计算(不论是静态还是结构优化),取ISMEAR=-5:当体系呈现金属性...
vasp说明里提到,SIGMA的取值一般使总能量与Energy without entropy的差值小于1meV就可以了。我算ZnO体系,SIGMA=0.01,上述两个能量的差值仍然大于1meV,而且SIGMA趋于0时两个能量的差值依然大于1meV。这是什么原因? 返回小木虫查看更多分享至: 更多 今日热帖求助...
SIGMA = 0.05 这两个参数决定在做布里渊区积分时,如何计算分布函数。ISMEAR常用的选择为0、1、2和-5。其中ISMEAR=0表示使用Gaussian展宽,一般用于半导体或绝缘体,同时设置展宽大小SIGMA为一个较小的值,如 0.05 eV。ISMEAR=1或者2表示Methfessel-Paxton方法,一般用于金属体系,同时可设置一个较大的SIGMA值,如0.2eV...
SIGMA = 0.2 ALGO = FastKPOINTS mesh auto 0 G 1 1 1 0 0 0POSCAR(由MS或者MedeA画出结构图,然后导出cif文件到 VESTA,再由VESTA导出POSCAR文件到vasp)POTCAR (赝势文件库) 2、记录计算结果,保存文件 能量:建议采用本子记下来,注明体系及重要参数。
ISMEAR用来确定如何或用何种方法来设置每个波函数的部分占有数fnk。在采用有限温度方法设置fnk时,smearing方法中的smearing宽度σ。它们的默认值分别为ISMEAR =1,SIGMA =0.2。可能的取值为-5,-4,-3,-2,-1,0,(N表示正整数),一般很少用-2和-3。 ISMEAR=-5,表示采用Blöchl修正的四面体方法。
一般情况下ISMEAR取0,同时取一个合适的SIGMA就不会出问题,笔者一般取ISMEAR= 0 ; SIGMA = 0.01进行计算。如果不知道如何取SIGMA值,可以查找侯柱峰老师的测试脚本进行测试。ALGO可以设置为Normal,Fast和Very_Fast,设置为Normal表示用blocked Davidson算法进行电子优化;设置为Very_Fast则表示用RMM-DIIS算法进行...
5、kmesh, SIGMA收敛问题还跟kmesh及SIGMA(当使用ISMEAR不等于-5 和-4时)的设置有关。要达到同样的精度,较小的SIGMA则需要较大的kmesh;而且,当SIGMA较小时,若KPOINTS不够多,也会出现难收敛的情况。以上便是VASP计算过程中的一些注意事项的总结,在这里给大家一一列出,当然实际的问题还要大家根据实际情况进行具体...