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型半导体的特征。
以上为个人理解，欢迎纠错讨论。

汪杰

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

卡老师，您提到“band计算所包含的信息并不是只有高对称点和路径上的，而是会读取自洽得到的密度和本征值去计算其余高对称点的本征值”
vasp中计算能带时，是根据KPOINTS的高对称点路径来计算的，那么其所得到的结果不就是被KPOINTS中的高对称点路径所约束吗？我最近就遇到自洽计算与非自洽计算得到的费米能级出现了显著差异，差距达到了2eV...

卡开发发

汪杰 发表于 2025-9-6 16:42
卡老师，您提到“band计算所包含的信息并不是只有高对称点和路径上的，而是会读取自洽得到的密度和本征值 ...

你尝试过增加自洽阶段k点数量吗？

汪杰

卡开发发 发表于 2025-9-6 17:30
你尝试过增加自洽阶段k点数量吗？

自洽所使用的k点密度应该是相当可靠了，是10*10*15，而我晶胞参数才4.8*4.8*3.3，且我体系还是绝缘体性质。所以这个k点数目怎么说都是足够的了吧。。

卡开发发

汪杰 发表于 2025-9-6 21:03
自洽所使用的k点密度应该是相当可靠了，是10*10*15，而我晶胞参数才4.8*4.8*3.3，且我体系还是绝缘体性质 ...

你可以给一些完整的计算信息看看。

汪杰

卡开发发 发表于 2025-9-6 23:55
你可以给一些完整的计算信息看看。

卡老师，相关的计算文件和输出已经上传至附件

卡开发发

汪杰 发表于 2025-9-7 12:52
卡老师，相关的计算文件和输出已经上传至附件

你试试把LMAXMIX设置到4。

汪杰

卡开发发 发表于 2025-9-8 01:15
你试试把LMAXMIX设置到4。

加上LMAXMIX=4的计算结果已经在附件，费米能级计算出来仍然是3 eV多。scf计算的就没有贴上去了，结果与原来一样

卡开发发

汪杰 发表于 2025-9-8 11:21
加上LMAXMIX=4的计算结果已经在附件，费米能级计算出来仍然是3 eV多。scf计算的就没有贴上去了，结果与原 ...

你传上来的是自洽计算而不是能带的结果。我用你数据测试过，差距应该没这么大。

汪杰

卡开发发 发表于 2025-9-8 16:31
你传上来的是自洽计算而不是能带的结果。我用你数据测试过，差距应该没这么大。

哦对的，我上传到的是自洽的结果，这应该是我弄混了。
但是卡老师，您计算结果，差距是多少呢？
如果完全按照我的输入文件，自洽计算大概是0.6，非自洽计算大概是3。我用的是vasp.5.4.4

卡开发发

汪杰 发表于 2025-9-8 21:57
哦对的，我上传到的是自洽的结果，这应该是我弄混了。
但是卡老师，您计算结果，差距是多少呢？
如果完 ...

非自洽0.78左右，自洽和你差不多。