计算化学公社

标题: 求助：能带图显示为导体，DOS图显示为半导体 [打印本页]

作者
Author: 每天进步一点点 时间: 2024-5-20 11:03
标题: 求助：能带图显示为导体，DOS图显示为半导体
本帖最后由每天进步一点点于 2024-5-20 13:29 编辑

各位老师好，我第一次用VASP计算，态密度和能带一直不对应，能带图显示为导体，DOS显示为半导体，这是怎么回事呢。以下分别是DOS和能带的的INCAR及结果。希望各位老师帮我看一下，谢谢！

作者
Author: 每天进步一点点 时间: 2024-5-20 11:07
而且计算过很多次结果都差不多，不知道哪里出了问题

作者
Author: Eudaimonia 时间: 2024-5-20 12:33
重新复习能带和DOS的关系，这俩图是对应的，都是典型的半导体
确定Fermi能级是否正确被程序读取

作者
Author: logzzz 时间: 2024-5-20 12:54
你这个DOS，不也是金属相么。哪里不对应了。

作者
Author: 每天进步一点点 时间: 2024-5-20 13:07

Eudaimonia 发表于 2024-5-20 12:33
重新复习能带和DOS的关系，这俩图是对应的，都是典型的半导体
确定Fermi能级是否正确被程序读取

谢谢老师，您的意思是因为费米能级在能带图中的位置不对才导致无法对应吗？我在静态自洽中得到的的费米能级只有0.07左右，计算band时vaspkit读取的费米能级是0.25ev，我是在vaspkit得到的数据基础上加上这两个的差值绘制图像的。
另外想请问老师，这个能带图是怎么看出来和dos对应的呢，我只能看出能带图0.8-1.8eV之间的带隙对应dos图0附近的直线，但是能带图上还有一个宽度在1.5eV左右的能带和3.5eV左右的带隙，这些在dos图上没有显示呀，dos图从0往左几乎全部有值，只在-2.5eV左侧有一个非常小的带隙。
学生比较愚钝，问题有点多，请老师指教！

作者
Author: 每天进步一点点 时间: 2024-5-20 13:10

logzzz 发表于 2024-5-20 12:54
你这个DOS，不也是金属相么。哪里不对应了。

老师我这个dos不是在0右侧有一段直线吗，虽然等于0的时候dos不为0，但是我觉得那个应该是展宽造成的，所以我判断dos图表示的是半导体。主要是我的图大体上看能带和带隙不对应，我不理解为什么dos图在0左侧几乎全部有值，但是我的能带图里没有宽度那么大的能带呀

作者
Author: 每天进步一点点 时间: 2024-5-20 13:28

Eudaimonia 发表于 2024-5-20 12:33
重新复习能带和DOS的关系，这俩图是对应的，都是典型的半导体
确定Fermi能级是否正确被程序读取

老师我上一条回复可能没表达清楚，更新了一张图片在主楼，我的意思是图里红色部分对应的是半导体的带隙吗，如果两个图片对应为什么能带图的绿色部分比较窄，dos图（没截全）的绿色部分对应的能量这么宽呢

作者
Author: 白皛白 时间: 2024-5-20 15:12
本帖最后由白皛白于 2024-5-20 15:18 编辑

请问K点是怎样设置的呢？如果K点数量不多的话，建议提高K点数量再计算band与dos不知结果会不会有所不同。另外在提高K的同时，使用ISMEAR=-5尝试一下？与此同时再提高些EDIFF的精度，比如1E-6或者1E-7试试。
在这个帖子里问题与楼主的类似，卡卡老师进行了细致的回答，楼主也可以参考下dos与band形状及费米能级对应不一致 - 第一性原理 (First Principle) - 计算化学公社 (keinsci.com)

作者
Author: logzzz 时间: 2024-5-20 15:22

每天进步一点点发表于 2024-5-20 13:10
老师我这个dos不是在0右侧有一段直线吗，虽然等于0的时候dos不为0，但是我觉得那个应该是展宽造成的，所 ...

我大概明白你是啥问题了。你是想问，为什么能带为0的地方投影出了态密度吧。你INCAR只有非自洽的部分，我觉得可以把NEDOS和NBANDS的数设置的再大些。还有你对费米能级的一顿处理，我没看太懂。但费米能级的处理可能也不太对。

作者
Author: 每天进步一点点 时间: 2024-5-20 15:37

白皛白发表于 2024-5-20 15:12
请问K点是怎样设置的呢？如果K点数量不多的话，建议提高K点数量再计算band与dos不知结果会不会有所不同。另 ...

老师我的体系有二百多个原子，态密度k点设置的是5*5*8，也用ISMEAR=-5算过，不过形状和主楼的图差不多，就是峰值更尖锐了一些。谢谢您，我再看一下这个帖子！

作者
Author: 每天进步一点点 时间: 2024-5-20 15:44

logzzz 发表于 2024-5-20 15:22
我大概明白你是啥问题了。你是想问，为什么能带为0的地方投影出了态密度吧。你INCAR只有非自洽的部分，我 ...

好的谢谢老师！关于费米能级，是因为我查到费米能级要取静态自洽得到的值（0.07eV），然后能带计算我用vaspkit生成数据时，vaspkit把Ef设置到0eV了，但是vaspkit读取的Ef是能带计算所得到的新的费米能级（在我这个例子里是0.25eV），相当于vaspkit多减了0.18eV，所以作图的时候我给所有能量值都加了个0.18，我也不知道这个操作对不对

欢迎光临计算化学公社 (http://bbs.keinsci.com/)