关于vaspkit算能带时费米能级的取法

lyf111

对于PBE，scf计算和能带计算的E-fermi有较大的的差值。我查找了我之前计算的结果 两者的差值可以达到0.4ev
有看到说scf的E-fermi准确的，取用scf中进行绘制能带图。
http://muchong.com/t-3237725-1

也看到说不能从scf中提取的。
http://bbs.vaspkit.cn/thread-7-1-1.html

另外我按第二个帖子里说的将Fermi-energy中的值修改后 运行vaspkit211功能他又自动的读取了DOSCAR中的E-fermi然后进行了覆盖 所以能带图并未改变。
也比较困惑是否会存在由于E-fermi能级的不同导致对材料的性质的误判，比如说费米能级是否穿过了价带这样。

lyf111

希望大家解答一下。
另外，我用pymatgen绘制能带图和vaspkit生成的BAND.dat绘制出来的基本一样。
然后pymatgen读取的也是能带计算的数据

vaspkit

因为能带计算只选取了布里渊区少部分的K点，显然用这部分K点计算得到物理量不能反映整个布里渊区的特性，所以PBE能带计算得到的费米能级是不准确，需要读取scf得到的费米能级。对于半导体，即使从DOSCAR读取一般不影响结论。

lyf111

vaspkit 发表于 2022-4-10 08:29
因为能带计算只选取了布里渊区少部分的K点，显然用这部分K点计算得到物理量不能反映整个布里渊区的特性，所 ...

您好，我想请问下，杂化泛函进行能带计算的时候也是选取的高对称点，为什么这个又计算的Ef是准确的呢？
另外，我想的是 我在计算前并不一定知道他是不是半导体，或者说他带隙值很小，可能会被忽略之类的，那这种情况下会有影响吗？
还有一个问题就是，我自己尝试修改Fermi-energy中的值后，在运行vaspkit211功能时，他又自动读取了Ef的值，那我该如何自己决定能带图中0对应的能量值呢？

vaspkit

因为杂化能带计算中权重不为零的部分K点是在整个布里渊区撒点得到的。后两个问题完全可以通过测试得到结论。