dos与band形状及费米能级对应不一致

hf_hua

ice-wine 发表于 2022-6-7 12:00
老师您好，我遇到了和楼主同样的问题，也看了您的解释，有一点小疑问。我用pymatgen绘图，脚本如下
bs_vas ...

同学你的问题解决了吗

乘风万里

ice-wine 发表于 2022-6-7 12:00
老师您好，我遇到了和楼主同样的问题，也看了您的解释，有一点小疑问。我用pymatgen绘图，脚本如下
bs_vas ...

图1是半导体，参考卡卡老师上面的回复，非金属体系的Fermi能级没有严格的定义，常见的定义在价带顶，但是VASP数据不是简单地用pymatgen就能画出这个效果的，你得手动调。
参考你的代码，可以在读入bs_data后这么写:

if not bs_data.is_metal():

    bs_data.efermi =bs_data.get_vbm()['energy']

然后就可以画出想要的效果，DOS数据同理。

图6你的sigma跟plot range都取得太大了。

fux

ice-wine 发表于 2022-6-7 12:00
老师您好，我遇到了和楼主同样的问题，也看了您的解释，有一点小疑问。我用pymatgen绘图，脚本如下
bs_vas ...

这里感觉画的DOS是不是有点平滑太多了，细节都被吞了好多，先把能量对应到相近的范围减小平滑再分析吧。

zhangbo11011

您好，请问一下，价带顶穿过费米能级必须是金属吗？会不会有可能是P型简并半导体呢

文学

运气加持 发表于 2021-12-2 15:45
老师，我想请教一下您说的第四点，如果计算上费米能级被定义到价带顶，那么是不是无法通过计算判断一种半 ...

时隔多年，遇到了类似的问题，看到了帖子。这个问题我想谈谈想法，无缺陷的TiO2其实是本征半导体，所以算出来费米能级在价带顶。而我们常说的“TiO2是N型半导体”这句话是因为实际存在的TiO2多少会有O空位，导致其成为N型半导体。如果计算O空位TiO2的DOS，我们会发现费米能级在导带底，这是典型的N型半导体的特征。
以上为个人理解，欢迎纠错讨论。