计算化学公社

标题: cp2k做GEO_OPT过程中SCF不收敛 [打印本页]

作者
Author: cuihongqiang 时间: 2021-3-10 09:12
标题: cp2k做GEO_OPT过程中SCF不收敛
我从官网上找到的一个例子用我自己的体系做，结果收敛不起来，想请教下是什么原因，谢谢

作者
Author: sobereva 时间: 2021-3-10 10:12
我不知道你算的什么体系，如果就是算个水分子，直接按下文的做法产生个优化的输入文件绝对能顺利跑出结果
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

作者
Author: wth1219 时间: 2021-3-10 12:55
20.0*20.0*20.0的体系，LDA/DZVP水平（水分子为什么要用LDA计算？），CUTOFF也很低（对于水分子而言过低了。），怎么会这么慢的？

作者
Author: cuihongqiang 时间: 2021-3-10 14:22

wth1219 发表于 2021-3-10 12:55
20.0*20.0*20.0的体系，LDA/DZVP水平（水分子为什么要用LDA计算？），CUTOFF也很低（对于水分子而言过低了 ...

基组的话是我在cp2k官网上找的例子，cutoff这个的选取标准是什么，望指点

作者
Author: sobereva 时间: 2021-3-11 02:35

cuihongqiang 发表于 2021-3-10 14:22
基组的话是我在cp2k官网上找的例子，cutoff这个的选取标准是什么，望指点

别盲目效仿官网上的例子，有的例子明显有不合理、不理想的地方，你的输入文件里有一大把不当的地方。用Multiwfn产生CP2K输入文件能让你少陷好多坑

CUTOFF问题看https://www.cp2k.org/howto:converging_cutoff
懒得做测试就拿相应基组里最大指数乘以40

作者
Author: cuihongqiang 时间: 2021-3-12 16:55

sobereva 发表于 2021-3-11 02:35
别盲目效仿官网上的例子，有的例子明显有不合理、不理想的地方，你的输入文件里有一大把不当的地方。用Mu ...

谢谢提醒

作者
Author: wth1219 时间: 2021-3-12 21:09

cuihongqiang 发表于 2021-3-10 14:22
基组的话是我在cp2k官网上找的例子，cutoff这个的选取标准是什么，望指点

赞同卢老师的回复。

CUTOFF精确测试：https://www.cp2k.org/howto:converging_cutoff
CUTOFF经验方法：https://www.cp2k.org/faq:cutoff

对于纯水体系，早期常用280 Ry，现在一般400 Ry。

作者
Author: sobereva 时间: 2021-3-14 03:02

jitou11 发表于 2021-3-13 14:16
这个SCF有一个问题，我也遇到了，我是用Multiwfn产生的输入文件。
问题是它SCF的stepsize都是一样的，他 ...

SCF过程并没有什么固定的stepsize，下一步波函数是什么样是通过诸如DIIS、OT等算法确定的。
倘若有个固定的步长，必须要求这一步和上一步波函数变化了多少，那SCF永远收敛不了

作者
Author: ma455173220 时间: 2021-12-7 06:29
因为 CP2K 默认的 MXING 方法（DIRECT_P_MIXING）不好用，改成BROYDEN_MIXING一般没啥问题。

作者
Author: chen0201 时间: 2023-3-13 09:44
请问对于一百多个原子并且含水分子的体系，aximum number of steps of inner SCF 和Maximum number of steps of outer SCF设置多少合适？

作者
Author: ykhuang 时间: 2023-3-13 10:50

chen0201 发表于 2023-3-13 09:44
请问对于一百多个原子并且含水分子的体系，aximum number of steps of inner SCF 和Maximum number of step ...

OT情况，30，5，一般就够

作者
Author: chen0201 时间: 2023-3-13 16:00

ykhuang 发表于 2023-3-13 10:50
OT情况，30，5，一般就够

请问这个设置依据是什么呢

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