计算化学公社

标题: VASP的输出文件中如何看是否收敛 [打印本页]

作者
Author: 流年2022 时间: 2021-5-26 09:22
标题: VASP的输出文件中如何看是否收敛
dE值不断变小，是否可以证明在收敛？

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

作者
Author: 流年2022 时间: 2021-5-26 09:43

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

能量

作者
Author: 卡开发发 时间: 2021-5-26 10:03

流年2022 发表于 2021-5-26 09:43
能量

如果你指的是单点能，那就是我前面说的自洽场的收敛标准参考。当dE和d eps都达到设定阈值，能量的自洽计算就会跳出循环。

作者
Author: 流年2022 时间: 2021-5-26 10:31

卡开发发发表于 2021-5-26 10:03
如果你指的是单点能，那就是我前面说的自洽场的收敛标准参考。当dE和d eps都达到设定阈值，能量的自洽计 ...

比如在高斯计算中，可以通过具体数值或者能量的变化知道计算是否收敛或者结构在往我理想的构型跑。不知道在vasp中是否有这样的输出可以清楚地把每一步的能量一起呈现出来？

作者
Author: 卡开发发 时间: 2021-5-26 10:37

流年2022 发表于 2021-5-26 10:31
比如在高斯计算中，可以通过具体数值或者能量的变化知道计算是否收敛或者结构在往我理想的构型跑。不知道 ...

$ grep "energy without entropy" OUTCAR | awk '{print $7}'

复制代码

作者
Author: 流年2022 时间: 2021-5-26 11:09

卡开发发发表于 2021-5-26 10:37

非常感谢！

作者
Author: 流年2022 时间: 2021-5-27 08:01

卡开发发发表于 2021-5-26 10:37

请问已经算了130步，deps已经为0，dE达到E-9，还没有结束，是不收敛吗

作者
Author: 卡开发发 时间: 2021-5-27 09:25

流年2022 发表于 2021-5-27 08:01
请问已经算了130步，deps已经为0，dE达到E-9，还没有结束，是不收敛吗

你可能需要上传一下你的输入输出。

作者
Author: 流年2022 时间: 2021-5-27 09:46

卡开发发发表于 2021-5-27 09:25
你可能需要上传一下你的输入输出。

INCAR 和 OUTCAR 文件？

作者
Author: 卡开发发 时间: 2021-5-27 09:56

流年2022 发表于 2021-5-27 09:46
INCAR 和 OUTCAR 文件？

对。

作者
Author: 流年2022 时间: 2021-5-27 10:11
本帖最后由流年2022 于 2021-6-9 20:21 编辑

这是INCAR文件 OUTCAR文件太大

作者
Author: 流年2022 时间: 2021-5-27 10:12

卡开发发发表于 2021-5-27 09:56
对。

您好已经上传

作者
Author: 卡开发发 时间: 2021-5-27 10:47

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

看下OSZICAR好了。

作者
Author: 流年2022 时间: 2021-5-27 12:00
本帖最后由流年2022 于 2021-6-9 20:21 编辑

已经上传

作者
Author: 流年2022 时间: 2021-5-27 12:00

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

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

作者
Author: 卡开发发 时间: 2021-5-27 12:11

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

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

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

作者
Author: 流年2022 时间: 2021-5-27 14:00

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

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

作者
Author: 卡开发发 时间: 2021-5-27 14:01

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

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

作者
Author: 卡开发发 时间: 2021-5-27 14:03

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

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

作者
Author: 流年2022 时间: 2021-5-27 14:11

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

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

作者
Author: 卡开发发 时间: 2021-5-27 14:16

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

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

作者
Author: 流年2022 时间: 2021-5-27 14:31

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

好的，谢谢您的回复

作者
Author: 乐平 时间: 2022-9-29 16:21

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

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

作者
Author: 卡开发发 时间: 2022-9-29 17:20

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

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

作者
Author: 乐平 时间: 2022-9-29 18:37
本帖最后由乐平于 2022-9-29 18:39 编辑

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

谢谢回复

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

(, 下载次数 Times of downloads: 34)

d eps the change in the bandstructure energy

d eps 是能带能量的变化。

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

作者
Author: 卡开发发 时间: 2022-9-29 18:55
本帖最后由卡开发发于 2022-9-29 19:00 编辑

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

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

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

作者
Author: 乐平 时间: 2022-9-29 20:00

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

感谢讨论！

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

作者
Author: 卡开发发 时间: 2022-9-29 22:13

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

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

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

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