计算化学公社

标题: 分子吸附在晶面时结构优化难以收敛怎么解决？ [打印本页]

作者
Author: 科研好难 时间: 2023-12-2 13:24
标题: 分子吸附在晶面时结构优化难以收敛怎么解决？
各位好，我在做分子吸附在晶面（结构如下图）时，首先使用BFGS方法结合k点3，3，1、截断能600和置信半径0.05进行计算，scf可以收敛，不过结构优化一直是震荡的，后来我分别尝试了相同参数结合LBFGS和CG方法还是同样的趋势，接下来我再次将k点改为2，2，1分别结合BFGS\LBFGS和CG计算还是收敛不了，最后这次我直接不设置k点并结合了波函数外推法USE_PREV_P计算（输入文件和输出文件如下），结果还是不收敛，请教一下各位是否还有什么好的调整方案？个人确实没辙了，还望各位不吝赐教，在此先感谢各位！
注：本人所有尝试皆使用PBE泛函和DZVP-MOLOPT-SR-GTH基组，计算是在80核，256G双路服务器中计算的。

(, 下载次数 Times of downloads: 11) (, 下载次数 Times of downloads: 13)

作者
Author: sobereva 时间: 2023-12-2 14:49
改k点毫无意义
本来Z方向就是设的非周期性，开偶极校正作甚，毫无意义
分清楚SCF不收敛和几何优化不收敛，几何优化不收敛弄USE_PREV_P毫无意义

不提供轨迹文件，别人难以判断情况

作者
Author: 科研好难 时间: 2023-12-2 17:39
谢谢sob老师解答，以下是我的见解，麻烦sob老师结合我二次提交的轨迹文件再指正一下：
1、k点主要影响的是scf的收敛，结构优化不收敛时可以不随便修改k点。
2、偶极校正只在三维周期性的垂直界面方向上使用。
@sobereva

作者
Author: sobereva 时间: 2023-12-3 00:46
从最后几十帧来看被吸附的分子结构还在微调而且没有显著震荡，这种情况继续跑

k点影响计算精度，对SCF收敛影响只是间接的，即便k点不合适SCF照样可能顺利收敛。k点和几何优化能否收敛没直接关系，用个够用的就完了，几何优化不收敛时尝试改小k点完全是无用功

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