发布时间 : 星期三 文章VASP中电子态密度计算的流程更新完毕开始阅读846d8ed0240c844769eaee09
#ISIF=2 (OR 3) NPAR=10
POTCAR 文件直接在势库中拷贝
POSCAR文件 Al 4.05
1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4
Selective Dynamic Direct
1.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T
KPOINT 文件
Automatic generation 0
Mohkorst Pack 21 21 21
0.0 0.0 0.0
VASP中文手册:在进行能带和DOS计算时,ISMEAR 不能使用N阶MP方法。因为MP方法在空轨道上有负的占据,所以求得的能带和DOS是不正确的。但是从其它地方看到---“提示:在计算能带结构时,采ISMEAR = 0或1对结果的影响非常小,可以认为是一样的。但是不能采用ISMEAR = -5 或-4。”ISMEAR到底多少?
计算能带:ICHARG = 11 导体的话,用ISMEAR=1; 半导体或绝缘体,用ISMEAR=0 。 计算 DOS: ICHARG = 11 ISMEAR = -5
计算的时候,金属可用0、1,非金属不要大过0,体材料可用-4、-5(面的话就用-1、0吧?)
参考http://210.45.66.120/doc/vasp5.0/node224.html
设置完成后进行计算,计算完后,得到包含了态密度值的DOSCAR文件,采用split_dos对态密度文件DOSCAR进行分割,得到总态密度DOS0,各个原子的分波态密度DOS1,DOS2……。另外在运行split_dos程序对DOSCAR文件分割时,要保证当前目录下有对应的OUTCAR和POSCAR文件。
分割后的DOS0,DOS1…等文件的能量值是以费米能级作为能量参考零点。DOS0的第一列数据是能量值,单位为eV;第二列数据是总态密度的值,单位 State/eV.unit cell;第三列数据是总态密度的积分值,也就是电子数,单位为electrons。DOS1是第一个原子的分波态密度值,其中的第一列数据是能量值,单位为eV;第二、三、四列数据分别对应于s、p、d态的分波态密度值,单位为State/eV.atom。其他的DOS文件与DOS1类似。