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

珊珊来迟

如图一所示，计算态密度时导带有一个尖峰，然而能带的导带是较为分散的一条能带，并没有与态密度导带对应的局域态出现（计算态密度时kpoints已做过收敛性测试，dos的形状应该是没有什么问题）。猜想可能的原因是band只计算了特殊路径的能带，但这里已经考虑了布里渊区中给出的所有路径（图二）。1.想从band中找到与dos导带处所出现的这个局域的峰所对应的能带，但不知该如何找到？或者说并没有必要去找，出现这种band与dos不符的情况是正常的？

2、另外还有一个问题困惑已久：在画band和dos时（用vaspkit处理的数据），会将费米能级平移至零点位置。由于band和dos的费米能级不一致，会出现对不齐的情况，一般是会将band或dos做一个平移。但请问应该按照band的费米能级平移还是按照dos的呢？（个人感觉是按照dos平移，因为band计算的是一些特殊的点，费米能级不准）
3、比如我计算的体系（小分子单晶，三维体系，含C,S,N原子），band费米能级位于带隙中间，而dos的费米能级进入了价带。此时应该将band按照dos的费米能级平移吗？另外，对于本征半导体，vasp计算出来的费米能级不是在禁带中间位置吗？为何dos显示费米能级进入了价带呢？如果是dos的这种费米能级穿过价带的情况，岂不是会定性得出体系具有半金属性的结论？
4.还有的程序会直接将价带顶设置为费米能级，一直不太理解，如果说费米能级是可以人为随意移动的，那我们根据费米能级与能带的相对位置来定性判断体系是否有金属性/半金属性还有意义吗？再者，如果都将费米能级设置到价带顶，对于掺杂的体系，如何判断掺杂之后造成的费米能级向导带/价带移动的情况呢？

汪杰

LY3353311 发表于 2022-1-26 12:47
**** 作者被禁止或删除 内容自动屏蔽 ****

费米能级的计算应该是自洽最准确，自洽的时候，K点设的密集一点，就可以算的比较准确的。能带计算肯定不如自洽准的

汪杰

卡开发发 发表于 2025-9-9 19:39
对DFT+U，使用LMAXMIX官网有描述，多看看官网可能可以学习到很多。

明白了老师，以后还是要多看看官网

卡开发发

汪杰 发表于 2025-9-9 17:25
问题解决了，非常感谢卡老师的帮助。这下又学到一个坑

对DFT+U，使用LMAXMIX官网有描述，多看看官网可能可以学习到很多。

汪杰

卡开发发 发表于 2025-9-9 10:27
我就是拿你的输入加了LMAXMIX=4做了scf，然后在把CHGCAR和WAVECAR拷贝到你算band的目录，然后把INCAR也加 ...

问题解决了，非常感谢卡老师的帮助。这下又学到一个坑

卡开发发

汪杰 发表于 2025-9-9 09:24
卡老师，请问您是用什么软件计算的呢？方便的话，可以将输出给我看一眼吗。我实在不知道问题出在哪里。。 ...

我就是拿你的输入加了LMAXMIX=4做了scf，然后在把CHGCAR和WAVECAR拷贝到你算band的目录，然后把INCAR也加上了LMAXMIX，如果还不行就私下讨论吧。

汪杰

卡开发发 发表于 2025-9-9 00:41
非自洽0.78左右，自洽和你差不多。

卡老师，请问您是用什么软件计算的呢？方便的话，可以将输出给我看一眼吗。我实在不知道问题出在哪里。。。

卡开发发

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

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

汪杰

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

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

卡开发发

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

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

汪杰

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

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

卡开发发

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

你试试把LMAXMIX设置到4。

汪杰

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

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

卡开发发

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

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

汪杰

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

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