计算化学公社

 找回密码 Forget password
 注册 Register

cp2k做GEO_OPT过程中SCF不收敛

查看数: 14331 | 评论数: 12 | 收藏 Add to favorites 7
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-3-10 09:12

正文摘要:

我从官网上找到的一个例子用我自己的体系做,结果收敛不起来,想请教下是什么原因,谢谢

回复 Reply

chen0201 发表于 Post on 2023-3-13 16:00:12
ykhuang 发表于 2023-3-13 10:50
OT情况,30,5,一般就够

请问这个设置依据是什么呢
ykhuang 发表于 Post on 2023-3-13 10:50:09
chen0201 发表于 2023-3-13 09:44
请问对于一百多个原子并且含水分子的体系,aximum number of steps of inner SCF 和Maximum number of step ...

OT情况,30,5,一般就够
chen0201 发表于 Post on 2023-3-13 09:44:29
请问对于一百多个原子并且含水分子的体系,aximum number of steps of inner SCF 和Maximum number of steps of outer SCF设置多少合适?
ma455173220 发表于 Post on 2021-12-7 06:29:42
因为 CP2K 默认的 MXING 方法(DIRECT_P_MIXING)不好用,改成BROYDEN_MIXING一般没啥问题。
sobereva 发表于 Post on 2021-3-14 03:02:33
jitou11 发表于 2021-3-13 14:16
这个SCF有一个问题,我也遇到了,我是用Multiwfn产生的输入文件。
问题是它SCF的stepsize都是一样的,他 ...

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

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
jitou11 + 5 谢谢

查看全部评分 View all ratings

wth1219 发表于 Post on 2021-3-12 21:09:33
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。
cuihongqiang 发表于 Post on 2021-3-12 16:55:37
sobereva 发表于 2021-3-11 02:35
别盲目效仿官网上的例子,有的例子明显有不合理、不理想的地方,你的输入文件里有一大把不当的地方。用Mu ...

谢谢提醒
sobereva 发表于 Post on 2021-3-11 02:35:48
cuihongqiang 发表于 2021-3-10 14:22
基组的话是我在cp2k官网上找的例子,cutoff这个的选取标准是什么,望指点

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

CUTOFF问题看https://www.cp2k.org/howto:converging_cutoff
懒得做测试就拿相应基组里最大指数乘以40
cuihongqiang 发表于 Post on 2021-3-10 14:22:24
wth1219 发表于 2021-3-10 12:55
20.0*20.0*20.0的体系,LDA/DZVP水平(水分子为什么要用LDA计算?),CUTOFF也很低(对于水分子而言过低了 ...

基组的话是我在cp2k官网上找的例子,cutoff这个的选取标准是什么,望指点
wth1219 发表于 Post on 2021-3-10 12:55:27
20.0*20.0*20.0的体系,LDA/DZVP水平(水分子为什么要用LDA计算?),CUTOFF也很低(对于水分子而言过低了。),怎么会这么慢的?
sobereva 发表于 Post on 2021-3-10 10:12:44
我不知道你算的什么体系,如果就是算个水分子,直接按下文的做法产生个优化的输入文件绝对能顺利跑出结果
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587http://bbs.keinsci.com/thread-21668-1-1.html

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 06:08 , Processed in 0.191780 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list