计算化学公社

标题: 求助：计算OER吸附OOH时结构总是不收敛，不知是结构问题还是INCAR参数设置不合理 [打印本页]

作者
Author: 晨晨晨嘛 时间: 2023-3-17 10:13
标题: 求助：计算OER吸附OOH时结构总是不收敛，不知是结构问题还是INCAR参数设置不合理
在吸附OOH后，OOH结构发生很大变化，一直不收敛。
kpoints设置的是2*2*1
d E和能量都达到了所设置的要求，但是计算一直没有停下来，不知道是哪里出现了问题，请各位老师同学指教

作者
Author: 晨晨晨嘛 时间: 2023-3-17 10:29
或者各位老师，OOH的各个键长键角的合理参数范围是多少，请指教

作者
Author: 含光君 时间: 2023-3-17 11:35
EDIFF太大了可能，可以设到1E-8，你这个INCAR相当于一个很粗糙的SCF精度配合了一个比较高的opt精度

作者
Author: 晨晨晨嘛 时间: 2023-3-18 09:23

含光君发表于 2023-3-17 11:35
EDIFF太大了可能，可以设到1E-8，你这个INCAR相当于一个很粗糙的SCF精度配合了一个比较高的opt精度

老师您看我的CONTCAR，ooh已经跑的没有样子了，是哪里出现了问题呢

作者
Author: 876449830 时间: 2023-3-18 10:43

晨晨晨嘛发表于 2023-3-18 09:23
老师您看我的CONTCAR，ooh已经跑的没有样子了，是哪里出现了问题呢

不明白为什么要减小EDIFF，你设置的EDIFFG=-0.01就是以力的收敛为标准的，可以查看一下RMS是多少，如果在0.01附近震荡，可以选择改IBRION=1，或者把POTIM改小一点。

作者
Author: 含光君 时间: 2023-3-18 12:14

876449830 发表于 2023-3-18 10:43
不明白为什么要减小EDIFF，你设置的EDIFFG=-0.01就是以力的收敛为标准的，可以查看一下RMS是多少，如果在 ...

力是能量对坐标的导数，你觉得力和能量没有关系么

作者
Author: 含光君 时间: 2023-3-18 12:20

晨晨晨嘛发表于 2023-3-18 09:23
老师您看我的CONTCAR，ooh已经跑的没有样子了，是哪里出现了问题呢

5L说的方法可以一试，此外我之前也遇到过类似情况，当时重新摆了摆起始结构就好了。

作者
Author: 晨晨晨嘛 时间: 2023-3-20 14:43

876449830 发表于 2023-3-18 10:43
不明白为什么要减小EDIFF，你设置的EDIFFG=-0.01就是以力的收敛为标准的，可以查看一下RMS是多少，如果在 ...

好的谢谢您！我试一下！

作者
Author: 晨晨晨嘛 时间: 2023-3-20 14:44

含光君发表于 2023-3-18 12:20
5L说的方法可以一试，此外我之前也遇到过类似情况，当时重新摆了摆起始结构就好了。

好的，我也试一下其他位置能否成功，谢谢您！

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