VASP达到收敛标准，但计算仍在继续

种花家的蓝兔

各位老师前辈好，我在OUTCAR文件中看到TOTAL-FORCE (eV/Angst)下边的数值已经小于|EDIFFG|，以及dE<EDIFF，这应该是已经达到收敛标准了吧，为啥计算仍在计算没有停止呢？望各位老师不吝赐教

jing@lll

grep RMS OUTCAR 看一下

种花家的蓝兔

jing@lll 发表于 2022-9-28 20:17
grep RMS OUTCAR 看一下

啥也没有

含光君

是不是这原子没收敛来着

种花家的蓝兔

含光君 发表于 2022-9-28 21:54
是不是这原子没收敛来着

看力收敛，是看上边这些数值<|EDIFFG|吗，如果是的话，我的EDIFF设置的是-0.02，您画的那一行满足小于0.02吧

含光君

种花家的蓝兔 发表于 2022-9-28 22:04
看力收敛，是看上边这些数值

VASPwiki（https://www.vasp.at/wiki/index.php/EDIFFG）上没有说明是每个方向的力还是合力，我个人理解是每个离子所受合力达到收敛标准才行。以及，在VASPwiki上一个实例（https://www.vasp.at/wiki/index.php/H2O）上出现如下表述：

从表达式看是对每一个离子判定收敛条件，所以我想大概是要达到合力收敛才行？以上内容仅为个人理解，如有不正确之处欢迎大佬指正。

卡开发发

种花家的蓝兔 发表于 2022-9-28 20:32
啥也没有

不要看RMS的Force，应该看Max的Force，但只有装了Henckelman他们的VTST插件这两个值才会打印到OUTCAR。

卡开发发

含光君 发表于 2022-9-28 22:43
VASPwiki（https://www.vasp.at/wiki/index.php/EDIFFG）上没有说明是每个方向的力还是合力，我个人理解 ...

应该是看每一个分量的，也就是N个forces（F=sqrt(Fx^2+Fy^2+Fz^2)）当中的最大的那一个分量的绝对值要≤-EDIFFG。LZ给出的力信息不完善其实没办法判断，另外有固定的原子也应该被排除。

种花家的蓝兔

卡开发发 发表于 2022-9-29 08:02
应该是看每一个分量的，也就是3N个forces（N个离子*3坐标）当中的最大的那一个分量的绝对值要≤-EDIFFG。 ...

那老师可以帮忙看看吗，我这个应该是达到了收敛标准，为啥计算仍在继续呢

卡开发发

种花家的蓝兔 发表于 2022-9-29 08:24
那老师可以帮忙看看吗，我这个应该是达到了收敛标准，为啥计算仍在继续呢

你不如发完整OUTCAR上来？

卡开发发

种花家的蓝兔 发表于 2022-9-29 08:24
那老师可以帮忙看看吗，我这个应该是达到了收敛标准，为啥计算仍在继续呢

我去翻了下源码，应该6楼说法是对的，应该是看N个原子的|F|=sqrt(Fx^2+Fy^2+Fz^2)最大的小于EDIFFG才行。

种花家的蓝兔

卡开发发 发表于 2022-9-29 08:55
我去翻了下源码，应该6楼说法是对的，应该是看N个原子的|F|=sqrt(Fx^2+Fy^2+Fz^2)最大的小于EDIFFG才行。

学到了学到了，谢谢老师的回复

种花家的蓝兔

含光君 发表于 2022-9-28 22:43
VASPwiki（https://www.vasp.at/wiki/index.php/EDIFFG）上没有说明是每个方向的力还是合力，我个人理解 ...

谢谢您的回复

含光君

卡开发发 发表于 2022-9-29 08:55
我去翻了下源码，应该6楼说法是对的，应该是看N个原子的|F|=sqrt(Fx^2+Fy^2+Fz^2)最大的小于EDIFFG才行。

谢谢发发老师~

wsz

vasp的收敛标准就是看所有非固定原子中受力最大的原子的受力，三个方向的模长
之前写过一个脚本，对于没有添加VTST的可以用这个查看受力，但是写的比较粗糙，需要手动修改一下体系原子数（第24行），没有去除固定的原子（需要的话可以自己修改一下）python2和3都可以，只用了内置的math库