计算化学公社

标题: DOS计算，报错提示为换ISMEAR=1or2 [打印本页]

作者
Author: lll000 时间: 2024-6-22 21:10
标题: DOS计算，报错提示为换ISMEAR=1or2
我是做高熵合金磁性体系的，但是这个DOS计算中，一直提示ISMEAR要更换，有无大佬来解释一下

作者
Author: lll000 时间: 2024-6-22 21:31
本帖最后由 lll000 于 2024-6-23 11:17 编辑

目前找到一个错误ISTART

作者
Author: 乐平 时间: 2024-6-22 21:57
本帖最后由乐平于 2024-6-22 22:01 编辑

最关键的问题是，你计算合金之类的金属体系不能用 ISMEAR = -5。应该用 ISMEAR = 1 或者 ISMEAR = 2，配合 SIGMA = 0.2，最好测试一下。

无论是计算态密度（DOS），能带（Band Structure），甚至最开始几何优化的时候，金属体系应该用 ISMEAR > 0 的参数。小于等于 0 适用于半导体和绝缘体。详细说明看 VASP wiki 手册关于 ISMEAR 的说明。
https://www.vasp.at/wiki/index.php/ISMEAR

另外，DOS 是非自洽计算，ISIF 是做几何优化，根本没必要写 ISIF

作者
Author: 卡开发发 时间: 2024-6-22 22:40
对DOS计算来说那个Warning不用管[参考]，

For the calculations of the DOS and very accurate total-energy calculations (no relaxation in metals), use the tetrahedron method (ISMEAR=-5).

这纯粹是因为自洽或非自洽结束后VASP仍然进行了force和stress的计算，只不过force和stress的计算用ISMEAR=-5亦即四面体积分的方法对金属来说不太可靠。但你只需要非自洽后电子结构信息，力不准这事情完全可以忽略。但需要注意的是，金属使用ISMEAR=-5的情况不管是电子结构还是涉及到力或应力计算，对k网格尺寸的依赖都比较大，建议慎重对k网格尺寸进行测试。

作者
Author: lll000 时间: 2024-6-23 11:02

乐平发表于 2024-6-22 21:57
最关键的问题是，你计算合金之类的金属体系不能用 ISMEAR = -5。应该用 ISMEAR = 1 或者 ISMEAR = 2，配合 ...

好的，谢谢大佬的回答，我回去测试一下更改ISMEAR

作者
Author: lll000 时间: 2024-6-23 11:07

卡开发发发表于 2024-6-22 22:40
对DOS计算来说那个Warning不用管[参考]，

这纯粹是因为自洽或非自洽结束后VASP仍然进行了force和stress ...

感谢卡卡老师的回复，您的建议也是更改一下ISMEAR吗？原谅我知识浅薄，听不太懂您的意思

作者
Author: lll000 时间: 2024-6-23 11:09
我怀疑是我从一开始的结构优化，静态自洽的计算有问题

，还劳驾各位大佬帮忙指正一下错误

作者
Author: lll000 时间: 2024-6-23 11:25
纯计算小白，还请各位大佬指点一二，做的是高熵合金磁性体系 (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

作者
Author: 卡开发发 时间: 2024-6-23 13:45

lll000 发表于 2024-6-23 11:07
感谢卡卡老师的回复，您的建议也是更改一下ISMEAR吗？原谅我知识浅薄，听不太懂您的意思

就是如果你只是算态密度可以不改，但如果不改的话建议把k网格尺寸好好测试一下。

作者
Author: lll000 时间: 2024-6-23 21:40

卡开发发发表于 2024-6-23 13:45
就是如果你只是算态密度可以不改，但如果不改的话建议把k网格尺寸好好测试一下。

卡老师，请问一下DOS计算ISMEAR一定要等于-5吗，如果我不更改K点基础上（因为是vaspkit,自动生成的，我不太会测试和改），能否改ISMEAR为1or2，这样算出来的DOS结果是否合理

作者
Author: 卡开发发 时间: 2024-6-23 22:05

lll000 发表于 2024-6-23 21:40
卡老师，请问一下DOS计算ISMEAR一定要等于-5吗，如果我不更改K点基础上（因为是vaspkit,自动生成的，我不 ...

不一定非要ISMEAR=-5，ISMEAR=-5的好处在于对SIGMA的设置没什么依赖，但k网格需要的比较密。实际上对金属来说ISMEAR=-1~2也都可以用，只要SIGMA设置合适即可。

作者
Author: lll000 时间: 2024-6-27 19:02

卡开发发发表于 2024-6-23 22:05
不一定非要ISMEAR=-5，ISMEAR=-5的好处在于对SIGMA的设置没什么依赖，但k网格需要的比较密。实际上对金属 ...

卡老师，我看很多教程上，都是设置SIGMA小了以后，每个原子的能量小于1mev，就满足条件，我测试过了，ISMEAR=1,SIGMA = 0.01 算出来每个原子上能量是小于1毫伏，比起ISMEAR=-5,测试k点网络，跟这个ISMEAR=1or2，测试SIGMA,对于我这个体系计算DOS来说，您觉得哪种方法更计算结果更准确

作者
Author: 卡开发发 时间: 2024-6-27 19:12

lll000 发表于 2024-6-27 19:02
卡老师，我看很多教程上，都是设置SIGMA小了以后，每个原子的能量小于1mev，就满足条件，我测试过了，ISM ...

道理上肯定是ISMEAR=-5取足够的k点精确，尤其是对计算能量来说。但只是说起来，到底什么时候一定要要求能量或DOS计算这么精确呢，这种场景其实反而遇到的不多，但这个说法仅供参考。

作者
Author: Stars 时间: 2024-6-30 10:49
官方的建议是，对于金属，涉及结构改变都不要用ISMEAR=-5。此外楼上已经说的很清楚了，对于金属设置ISMEAR=1时可以取大的σ和小的k点，不信可以看电子熵，ISMEAR=1时也很小，至少比0时候小得多

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