VASP的输出文件中如何看是否收敛

流年2022

卡开发发 发表于 2021-5-27 10:47
看下OSZICAR好了。

是因为力收敛标准尚未达成？

卡开发发

流年2022 发表于 2021-5-27 12:00
已经上传

你这是构型优化没收敛。构型优化看的不是EDIFF而是EDIIFG，而且EDIFFG设置负值是参考力收敛。这个力参考的是OUTCAR当中3N个原子的力中最大的那个要小于EDIFFG才能达到收敛。如果安装了Henckelman的NEB，也会有FORCES的字段直接给出最大力。EDIFFG>0才是以力收敛作为参考，此时应该参考自洽迭代后的d E，不过一般不太建议用能量作为收敛标准。

另外你的第一步构型优化中的自洽场计算一开始就没收敛，导致力和能量都有偏差，最终构型优化难收敛也就不足为奇了。

流年2022

卡开发发 发表于 2021-5-27 12:11
你这是构型优化没收敛。构型优化看的不是EDIFF而是EDIIFG，而且EDIFFG设置负值是参考力收敛。这个力参考 ...

您的意思是EDIFFG设置的有问题？

卡开发发

流年2022 发表于 2021-5-27 13:56
这样的话，请问应该如何去解决？

1、确定模型合理、磁矩初始猜测合理。
2、确定参数合理，该做测试的做测试。
3、SCF难收敛的话适当通过AMIX和AMIX_mag之类参数来减小振荡。

卡开发发

流年2022 发表于 2021-5-27 14:00
您的意思是EDIFFG设置的有问题？

不是，和这个问题没关系。

流年2022

卡开发发 发表于 2021-5-27 14:01
1、确定模型合理、磁矩初始猜测合理。
2、确定参数合理，该做测试的做测试。
3、SCF难收敛的话适当通过 ...

还有两个问题1 您是如何确定第一步就没有收敛？通过dE和d eps嘛？
2 您说的AMIX和AMIX_mag是要添加到INCAR中？

卡开发发

流年2022 发表于 2021-5-27 14:11
还有两个问题1 您是如何确定第一步就没有收敛？通过dE和d eps嘛？
2 您说的AMIX和AMIX_mag是要添加到INC ...

到了NELM最大步数60不你的dE和d eps都没有达到设置的EDIFF。你需要明确一些观念，就是什么是自洽场（电子步）什么是构型优化（离子步），他们分别在做些什么，两者的关系是什么。

流年2022

卡开发发 发表于 2021-5-27 14:16
到了NELM最大步数60不你的dE和d eps都没有达到设置的EDIFF。你需要明确一些观念，就是什么是自洽场（电子 ...

好的，谢谢您的回复

乐平

卡开发发 发表于 2021-5-26 09:28
什么收敛，自洽场（电子步）还是构型优化（离子步）？
如果是电子步的收敛，需要看dE和d eps两个。

请问 d eps 里的 eps 具体是什么意思？

卡开发发

乐平 发表于 2022-9-29 16:21
请问 d eps 里的 eps 具体是什么意思？

应该是占据态本征值的总变化吧。

乐平

卡开发发 发表于 2022-9-29 17:20
应该是占据态本征值的总变化吧。

谢谢回复

我在 VAPS wiki 网页的  OSZICAR 里看到这样的说法链接如下
https://www.vasp.at/wiki/index.php/OSZICAR





d eps the change in the bandstructure energy

d eps 是能带能量的变化。

您说的“占据态本征值的总变化”是不是可以认为是价带能量的变化呢？

卡开发发

乐平 发表于 2022-9-29 18:37
谢谢回复

我在 VAPS wiki 网页的  OSZICAR 里看到这样的说法链接如下

应该不包含空轨道的，去翻了下源码，求和部分应该是用到了占据信息。说价带也不大合适，因为可能还包含一些不在价带的占据轨道（满带）。

乐平

卡开发发 发表于 2022-9-29 18:55
应该不包含空轨道的，去翻了下源码，求和部分应该是用到了占据信息。说价带也不大合适，因为可能还包含一 ...

感谢讨论！

这样说来，应该是所有的电子占据轨道（占据的能带）的能量的变化量。
eps 也许是 electron populated states 的缩写？

卡开发发

乐平 发表于 2022-9-29 20:00
感谢讨论！

这样说来，应该是所有的电子占据轨道（占据的能带）的能量的变化量。

我觉得也可能是epsilon的缩写。