VASP做纯碳氢分子几何优化，力的收敛如何判断？

ZYXA

各位老师好，
我最近刚刚接触VASP计算，向各位老师求教。
请问：
（1）力的收敛标准是三个方向同时满足小于EDIFFG吗？如果不是，如何求解呢？
（2）优化成功是必须同时满足能量收敛和力收敛标准吗？
（3）如果出现能量收敛，但力未收敛是否需要进一步调整优化？
（3）几何优化成功与否和log文件中“reached required accuracy - stopping structural energy minimisation”无关，对吗？

附上INCAR、OSZICAR、OUTCAR相关截图。
谢谢~

ZYXA

截图如下

sobereva

ZYXA 发表于 2023-12-21 20:02
截图如下

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

ZYXA

好的，老师。以后一定注意。

白皛白

1,是指最大原子受力的模。
2,是否结束结构优化只有最大原子受力这个收敛标准。关于“其中有一些离子步没有收敛但是最终最大原子受力达到了收敛标准能否算作成功的结构优化”这个问题我个人也不是很明白，我自己做的时候是会让所有离子步都收敛的，同门们一般也是如此认为，望有老师可以答疑解惑。
3,只有力达到收敛限才会结束几何结构的弛豫。
4,是最大原子受力达到收敛限的标志，是否每个SCF步骤都收敛了需要检查一下。
其实我觉得能量收敛和力收敛这种说法很容易让新手迷惑，我当时也迷惑了很久。这个能量收敛实际上应该就是SCF收敛，由EDIFF这个参数决定。力收敛有一个脚本可以直接把最大原子受力的模输出出来。一个是针对离子步的，一个是针对全局的。
注：该脚本来自知乎用户“叠加态”

卡开发发

白皛白 发表于 2023-12-22 10:56
1,是指最大原子受力的模。
2,是否结束结构优化只有最大原子受力这个收敛标准。关于“其中有一些离子步没有 ...

2、其实也不是必须，但通常建议把SCF都做收敛会比较稳妥，尤其是自旋极化情形下，如果上一步自洽场并未收敛，可能后续步骤按照这个电子密度进行初猜可能磁矩都是错的，另外那一步的力可能也是错的。这个问题以往有类似讨论，不收敛继续做下去类似于Gaussian的iop(5/13=1)。
4、Henkelman他们组的那个VTST tools装过之后其实可以直接grep FORCES OUTCAR来看，甚至Stress的最大分量也能获取到，在NEB xxx的一个block下面，这样应该会更简单。LZ指的是能量收敛可能是SCF的自洽场收敛，这种说法很容易造成误解。VASP的优化算法当中的收敛标准要么是只能力收敛要么只能是能量收敛，即EDIFFG按照负或者正来取，其他程序不一定，甚至还会有原子位移的收敛标准。

ZYXA

白皛白 发表于 2023-12-22 10:56
1,是指最大原子受力的模。
2,是否结束结构优化只有最大原子受力这个收敛标准。关于“其中有一些离子步没有 ...

谢谢指导~，祝科研顺利~

ZYXA

卡开发发 发表于 2023-12-22 11:39
2、其实也不是必须，但通常建议把SCF都做收敛会比较稳妥，尤其是自旋极化情形下，如果上一步自洽场并未收 ...

谢谢指导~
也就是说EDIFFG小于零时: 关注力收敛情况(EDIFFG; log. 文件中reached required accuracy)和电子步SCF收敛(dE)情况(关注EDIFF)，对吗？

卡开发发

ZYXA 发表于 2023-12-22 14:12
谢谢指导~
也就是说EDIFFG小于零时: 关注力收敛情况(EDIFFG; log. 文件中reached required accuracy)和 ...

对。实际上EDIFF收敛包括了dE和d eps。

ZYXA

卡开发发 发表于 2023-12-22 17:11
对。实际上EDIFF收敛包括了dE和d eps。

谢谢~，祝一切顺利~