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

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

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

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

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

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都取得太大了。

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

同学你的问题解决了吗

老师您好，我遇到了和楼主同样的问题，也看了您的解释，有一点小疑问。我用pymatgen绘图，脚本如下
bs_vasprun = Vasprun("./band/vasprun.xml",parse_projected_eigen=True)
bs_data = bs_vasprun.get_band_structure(line_mode=True, force_hybrid_mode=True)

dos_vasprun=Vasprun("./dos/vasprun.xml")
dos_data = bs_vasprun.complete_dos

plt_1 = BSDOSPlotter(bs_projection=None, dos_projection=None)
plt_1.get_plot(bs=bs_data,dos=dos_data)

plt_6 = DosPlotter(stack=False,sigma=0.5)
plt_6.add_dos('total dos',dos=dos_data)

理论上绘图读取的dos数据都一样，为何在图1中费米能级在DOS图中间，而在图6中费米能级进入了价带。感谢

珊珊来迟 发表于 2021-2-27 18:36
谢谢您的耐心解答！对于展宽问题，我在计算dos时都用的ISMEAR=-5,计算能带时都采用的ISMEAR=0。另外还有 ...

萌新求问，一般band用ISMEAR=0计算时，此时涉及展宽，而如果dos用布洛赫修正四面体法，就不涉及展宽的困扰，此时如果band的K点取值没有问题，那么费米能级处的情况，应该是dos结果更准一点吧。

卡开发发 发表于 2021-12-3 01:04
可能不大好用这种方式判断。

嗯嗯，谢谢老师

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

可能不大好用这种方式判断。

卡开发发 发表于 2021-2-27 08:30
1、可以看下不同的k对应的不同能带在这个窄峰位置能级的对应，造成这样问题的原因我暂时不是很明确。
2、 ...

老师，我想请教一下您说的第四点，如果计算上费米能级被定义到价带顶，那么是不是无法通过计算判断一种半导体是N型还是P型呢？我计算TiO2得到的费米能级就在价带顶附近，但是资料显示TiO2是N型半导体。

珊珊来迟 发表于 2021-2-28 11:04
老师我看了一下，dos对应能级在不同k点下占据数是一样的，而band出现了不一样的情况，主要集中在高对称点 ...

1、在Fermi能级附近从占据数2（非自旋极化的情况）经过分数过渡到0这种也都是展宽引起的，无需过于担心。这也是对你前面给的图占据数分数情形的解释。DOS没有是因为你的DOS计算用了四面体的方法做的。
2、就是说，我们先确定窄峰在DOS上的能级范围，然后对照OUTCAR不同的k点中同样的能级范围看看这个窄峰是由哪些能级组成，然后在能带当中也根据这个能量范围确定能带的标号画出来，看看是否这些能带能够确实与DOS一致。
3、虽然不是必须，只是免得审稿人问起来没办法解释，不如提前把顾虑解决掉，也顺带检查下计算中是否因为某些设置因粗心而犯错。

卡开发发 发表于 2021-2-27 22:57
208-211的分数占据是因为你使用了展宽引起的，应该不是这个问题引起的那个局域峰和能带对不上。我觉得可 ...

老师我看了一下，dos对应能级在不同k点下占据数是一样的，而band出现了不一样的情况，主要集中在高对称点Y以及R右边的几个点，他们的占据数为大小不等的分数，而其他均为0. 但还是不太懂这个结果怎么分析。。这怎么看出窄峰是否为窄带的贡献呢？希望老师再帮学生指点迷津
另外如果可确定窄峰不是窄带的贡献，就说明dos和band确实是这种不对应的情况，但具体为什么不对应还得再找原因对吗

珊珊来迟 发表于 2021-2-27 22:30
感谢老师！学到很多！您第一点说到了能级占据数的情况，我查看了dos和band的能带占据数，发现有几条能带 ...

208-211的分数占据是因为你使用了展宽引起的，应该不是这个问题引起的那个局域峰和能带对不上。我觉得可能需要排查的是那个窄峰出现的位置所对应的能级在能带计算和DOS计算不同k的情况下，两者有所不同。如果确实从DOS上看窄峰也不是窄带的贡献，那应该就没啥问题了。