计算化学公社

 找回密码 Forget password
 注册 Register
Views: 14324|回复 Reply: 12
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] cp2k做GEO_OPT过程中SCF不收敛

[复制链接 Copy URL]

67

帖子

0

威望

1283

eV
积分
1350

Level 4 (黑子)

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

捕获.PNG (62.28 KB, 下载次数 Times of downloads: 44)

捕获.PNG

test.inp

1.44 KB, 下载次数 Times of downloads: 61

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

2#
发表于 Post on 2021-3-10 10:12:44 | 只看该作者 Only view this author
我不知道你算的什么体系,如果就是算个水分子,直接按下文的做法产生个优化的输入文件绝对能顺利跑出结果
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587http://bbs.keinsci.com/thread-21668-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

106

帖子

0

威望

3567

eV
积分
3673

Level 5 (御坂)

被催的人

3#
发表于 Post on 2021-3-10 12:55:27 | 只看该作者 Only view this author
20.0*20.0*20.0的体系,LDA/DZVP水平(水分子为什么要用LDA计算?),CUTOFF也很低(对于水分子而言过低了。),怎么会这么慢的?

67

帖子

0

威望

1283

eV
积分
1350

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2021-3-10 14:22:24 | 只看该作者 Only view this author
wth1219 发表于 2021-3-10 12:55
20.0*20.0*20.0的体系,LDA/DZVP水平(水分子为什么要用LDA计算?),CUTOFF也很低(对于水分子而言过低了 ...

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

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

5#
发表于 Post on 2021-3-11 02:35:48 | 只看该作者 Only view this author
cuihongqiang 发表于 2021-3-10 14:22
基组的话是我在cp2k官网上找的例子,cutoff这个的选取标准是什么,望指点

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

CUTOFF问题看https://www.cp2k.org/howto:converging_cutoff
懒得做测试就拿相应基组里最大指数乘以40
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

67

帖子

0

威望

1283

eV
积分
1350

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2021-3-12 16:55:37 | 只看该作者 Only view this author
sobereva 发表于 2021-3-11 02:35
别盲目效仿官网上的例子,有的例子明显有不合理、不理想的地方,你的输入文件里有一大把不当的地方。用Mu ...

谢谢提醒

106

帖子

0

威望

3567

eV
积分
3673

Level 5 (御坂)

被催的人

7#
发表于 Post on 2021-3-12 21:09:33 | 只看该作者 Only view this author
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。

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

8#
发表于 Post on 2021-3-14 03:02:33 | 只看该作者 Only view this author
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

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

39

帖子

3

威望

551

eV
积分
650

Level 4 (黑子)

9#
发表于 Post on 2021-12-7 06:29:42 | 只看该作者 Only view this author
因为 CP2K 默认的 MXING 方法(DIRECT_P_MIXING)不好用,改成BROYDEN_MIXING一般没啥问题。

181

帖子

0

威望

731

eV
积分
912

Level 4 (黑子)

10#
发表于 Post on 2023-3-13 09:44:29 | 只看该作者 Only view this author
请问对于一百多个原子并且含水分子的体系,aximum number of steps of inner SCF 和Maximum number of steps of outer SCF设置多少合适?

14

帖子

1

威望

144

eV
积分
178

Level 3 能力者

11#
发表于 Post on 2023-3-13 10:50:09 来自手机 | 只看该作者 Only view this author
chen0201 发表于 2023-3-13 09:44
请问对于一百多个原子并且含水分子的体系,aximum number of steps of inner SCF 和Maximum number of step ...

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

181

帖子

0

威望

731

eV
积分
912

Level 4 (黑子)

12#
发表于 Post on 2023-3-13 16:00:12 | 只看该作者 Only view this author
ykhuang 发表于 2023-3-13 10:50
OT情况,30,5,一般就够

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

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 23:57 , Processed in 0.279637 second(s), 24 queries , Gzip On.

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