计算化学公社

标题: VASP计算能带出现负占据 [打印本页]

作者
Author: 小书童 时间: 2017-6-29 08:44
标题: VASP计算能带出现负占据
VASP用ISMEAR=-5进行静态计算，结果在OUTCAR中发现能带出现了负占据，如图所示，不明白为什么会出现负占据，请各位老师指点一下迷津，先谢过！

作者
Author: 卡开发发 时间: 2017-6-29 08:50
http://bbs.keinsci.com/forum.php ... &extra=page%3D1

作者
Author: 小书童 时间: 2017-6-29 09:00

卡开发发发表于 2017-6-29 08:50
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=6050&extra=page%3D1

谢谢老师，我算的是石墨烯二维材料的掺杂一个金属Fe，这是说明我的体系不适合用ISMEAR=-5做静态计算吗？

作者
Author: 卡开发发 时间: 2017-6-29 10:25

小书童发表于 2017-6-29 09:00
谢谢老师，我算的是石墨烯二维材料的掺杂一个金属Fe，这是说明我的体系不适合用ISMEAR=-5做静态计算吗？

可能这个体系磁矩不太局域，所以实际可能需要很大的k才有可能得到好的结果，为了避免这个问题可以考虑使用ismear>=0。

作者
Author: 小书童 时间: 2017-6-29 15:28

卡开发发发表于 2017-6-29 10:25
可能这个体系磁矩不太局域，所以实际可能需要很大的k才有可能得到好的结果，为了避免这个问题可以考虑使 ...

嗯嗯，多谢老师指点

作者
Author: 小书童 时间: 2017-7-1 16:44

卡开发发发表于 2017-6-29 10:25
可能这个体系磁矩不太局域，所以实际可能需要很大的k才有可能得到好的结果，为了避免这个问题可以考虑使 ...

老师，我对比了有能带负占据（ISMEAR=-5）和没有负占据（ISMEAR=0）的能量值，发现几乎一样，说明这个能带的负占据对能量计算是没有什么影响的，那么对能带和DOS的绘制会有影响吗？（我对能带和DOS态密度绘制还不太懂，所以暂时没法通过自己绘制对比检验），还请老师指点一下。这个能带的负占据根据您的经验，主要是对那些性质影响比较大呢？

作者
Author: 卡开发发 时间: 2017-7-1 17:28

小书童发表于 2017-7-1 16:44
老师，我对比了有能带负占据（ISMEAR=-5）和没有负占据（ISMEAR=0）的能量值，发现几乎一样，说明这个能 ...

占据数是负的（有可能还同时出现占据数大于1的能态），原则上在DOS上表现也就是对峰的贡献也是负的，能带会是什么样的图像我不是很清楚。至于其他性质涉及到电子结构的性质我认为应该也会受到影响，包括force等也会有影响，本来也不建议金属性的体系用ismear=-5来做，可能需要非常大的k点才能得到定性正确的结果。

作者
Author: 小书童 时间: 2017-7-1 18:59
本帖最后由小书童于 2017-7-1 19:02 编辑

卡开发发发表于 2017-7-1 17:28
占据数是负的（有可能还同时出现占据数大于1的能态），原则上在DOS上表现也就是对峰的贡献也是负的，能带 ...

嗯嗯，谢谢老师的讨论，看来-5还真的不是和金属体系，我的体系里面就一个金属，都会出错。我刚刚查阅了官网对ISMEAR的说明，看到了这句话，如图。意思是-5不能够用于含金属的体系求DOS和能量。

作者
Author: 卡开发发 时间: 2017-7-1 20:03

小书童发表于 2017-7-1 18:59
嗯嗯，谢谢老师的讨论，看来-5还真的不是和金属体系，我的体系里面就一个金属，都会出错。我刚刚查阅了官 ...

也不是这个意思。

ismear=-5不适合用于金属体系的构型优化。另外如上面反复强调的问题，即便用于电子结构计算，如果体系具有金属性（和是否包含金属原子其实没有直接关系），电子并不是很局域的情况，可能需要很大的k点才能够正确描述。

作者
Author: 小书童 时间: 2017-7-1 21:44
本帖最后由小书童于 2017-7-1 21:46 编辑

卡开发发发表于 2017-7-1 20:03
也不是这个意思。

ismear=-5不适合用于金属体系的构型优化。另外如上面反复强调的问题，即便用于电子 ...

嗯嗯，我是静态计算的时候发现问题的，我的K点已经给了15*15*1，还有这个问题。那我就结构优化和静态计算都改成ISMEAR=0计算，这样应该是没问题的吧？谢谢老师指点

作者
Author: 卡开发发 时间: 2017-7-2 03:16
本帖最后由卡开发发于 2017-7-2 03:30 编辑

小书童发表于 2017-7-1 21:44
嗯嗯，我是静态计算的时候发现问题的，我的K点已经给了15*15*1，还有这个问题。那我就结构优化和静态计算 ...

没问题，静态计算的时候展宽小一点就可以了，也可以得到比较可靠的结果。

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