第一步直接看参数 以及参数简单解析
#读取文件磁性设置等等
ISTART=0 #随机生成波函数
ICHARG=2 #原子电荷密度叠加
#以上两个参数结构优化的时候 ISTART=0 ICHARG=2
#对于自洽计算可以0 2组合也可以1 1组合;对于非自洽计算1 11组合
PREC=M
ISPIN=2
LCHARG=.T.
LWAVE=.T.
#电子步
ALGO=N #这步其实是废话因为本来默认就是N 一般都用N
NELM=100
EDIFF=1E-5
ENCUT=500
#离子步
IBRION=2
NSW=100
ISIF=3
EDIFFG=-0.01
ISYM=0 #小体系可以0关掉,大体系可以2先找对称性,0优化更加充分,2速度更快
LREAL=Auto # 这个东西太迷以后再议论 但是默认值是.FALSE.
ISMEAR=-5 #金属不能用-5
SIGMA=0.05
#DFT+U
MAGMOM=32*5.0 47*2 0
LDAU=.T.
LDAUL=2 1 -1
LDAUU=5.3 0 0
LDAUJ=0 0 0
LDAUTYPE=2
LMAXMIX=4
ISMEAR
物理含义:ISMEAR这一个参数规定的费米面附近电子是否严格的遵守费米面以上没有电子,-5含义,时候表示严格遵守费米面上没有电子,0表示有半填充的能带
ISMEAR=-5 0半导体绝缘体等有禁带的东西;
ISMEAR=0能带是要跨越费米面的,金属。
总结:金属不能用-5
- 不同的值对应的是不同的展宽方法;
- 对于分子,原子的计算用0;
- 对所有体系,如果想获取更加精确能量的时候用-5,但这时候如果K点数目小于3,程序则会罢工;K 点少,半导体或者绝缘体,那么只能用 ISMEAR = 0;
- 在DOS能带计算中,使用ISMEAR= -5 用于获取精确的信息。
SIGMA
- 如果用了ISMEAR = -5; SIGMA的值可以忽略,也可以不管。(对于所有体系都适用)
- 对于金属: ISMEAR = 1 或者0,非金属: ISMEAR= 0 的时候,一般取 SIGMA = 0.10 即可,默认值是0.20。不放心的话,用0.05;
- 对于气体分子,原子体系(也就是你把分子或者原子放到一个box里面):死死记住下面组合就可以了, ISMEAR = 0; SIGMA = 0.01
这里介绍一个验证SIGMA取值正确性的方法:
grep 'entropy T' OUTCAR
用这个命令,得出的能量除以体系中原子的数目(也就是平均到每个原子上),然后再和 0.001 eV 比较。如果小于,SIGMA取值OK,如果大于,再换个小点的数值进行测试。(这里的说的entropy,是因为ISMEAR展宽这个计算方法引入的,跟我们物理化学书里面的entropy完全是两码事。后面我们会介绍如何通过频率分析计算物理化学里面的Entropy。千万别把这个数值当成频率所对应的entropy。)