lyf111 发表于 2022-4-9 10:52:41

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

本帖最后由 lyf111 于 2022-4-9 10:56 编辑

对于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 发表于 2022-4-9 10:54:08

本帖最后由 lyf111 于 2022-4-9 11:06 编辑

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

vaspkit 发表于 2022-4-10 08:29:30

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

lyf111 发表于 2022-4-10 11:49:31

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

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

vaspkit 发表于 2022-4-10 16:10:55

因为杂化能带计算中权重不为零的部分K点是在整个布里渊区撒点得到的。后两个问题完全可以通过测试得到结论。
页: [1]
查看完整版本: 关于vaspkit算能带时费米能级的取法