计算化学公社

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

[CP2K] 求助cp2k优化钙钛矿结构SCF不收敛稳定在0.00000147

[复制链接 Copy URL]

4

帖子

0

威望

59

eV
积分
63

Level 2 能力者

本帖最后由 dyj_ 于 2024-11-7 15:31 编辑

按照文献,使用guassview建模,导出gjf文件,使用multiwfn生成cp2k的输入文件(输入文件按照课程结构优化示例7 Pd(100)晶面对苯分子的吸附设置的,没有冻结原子,使用OT,CUTOFF 400 REL_CUTOFF 55),参考CP2K中遇到SCF难收敛时的解决方法博文进行修改,Convergence一直维持在0.00000147没有变化,刚开始学cp2k,所以不清楚是建模有问题还是参数设置有问题。

文章图.png (36.17 KB, 下载次数 Times of downloads: 0)

文章图.png

pvk-opt.inp

64.48 KB, 下载次数 Times of downloads: 7

输入文件

pvk-opt.out

11.48 KB, 下载次数 Times of downloads: 2

输出文件

pvk.gjf

62.84 KB, 下载次数 Times of downloads: 1

模型

53

帖子

0

威望

197

eV
积分
250

Level 3 能力者

2#
发表于 Post on 2024-11-7 14:54:22 | 只看该作者 Only view this author
可以尝试修改参数中的PRECONDITIONER和MINIMIZER
改为以下参数后计算是逐渐收敛的,如图
PRECONDITIONER FULL_SINGLE_INVERSE
MINIMIZER CG

202411071452169610..png (83.84 KB, 下载次数 Times of downloads: 0)

202411071452169610..png

4

帖子

0

威望

59

eV
积分
63

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-11-7 15:22:10 | 只看该作者 Only view this author
1372133091 发表于 2024-11-7 14:54
可以尝试修改参数中的PRECONDITIONER和MINIMIZER
改为以下参数后计算是逐渐收敛的,如图
PRECONDITIONER  ...

谢谢,我试过这个MINIMIZER CG;我试试这个PRECONDITIONER FULL_SINGLE_INVERSE。

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

4#
发表于 Post on 2024-11-9 07:15:44 | 只看该作者 Only view this author
0.00000147就算不错了,可以SCF收敛限适当放宽点,诸如3E-6。等优化到极小点结构附近SCF会更容易收敛,可以到时候再读取之前的波函数当初猜结合1E-6继续优化。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

4

帖子

0

威望

59

eV
积分
63

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2024-11-9 08:36:39 | 只看该作者 Only view this author
sobereva 发表于 2024-11-9 07:15
0.00000147就算不错了,可以SCF收敛限适当放宽点,诸如3E-6。等优化到极小点结构附近SCF会更容易收敛,可以 ...

感谢社长的指导。修改结构后SCF收敛,但是结构优化跑了38步没有收敛,因为刚开始学习,不太清楚结构优化多少步收敛是比较正常的。然后28步的时候,三个yes,37步一个yes,这是正常的么?

28步.png (161.33 KB, 下载次数 Times of downloads: 0)

28步.png

37步.png (113.86 KB, 下载次数 Times of downloads: 0)

37步.png

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

6#
发表于 Post on 2024-11-9 08:52:43 | 只看该作者 Only view this author
dyj_ 发表于 2024-11-9 08:36
感谢社长的指导。修改结构后SCF收敛,但是结构优化跑了38步没有收敛,因为刚开始学习,不太清楚结构优化 ...

大几十步才收敛是家常便饭
中途偶尔震荡导致YES暂时失去也是家常便饭
要有耐心
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

4

帖子

0

威望

59

eV
积分
63

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2024-11-9 09:13:13 | 只看该作者 Only view this author
sobereva 发表于 2024-11-9 08:52
大几十步才收敛是家常便饭
中途偶尔震荡导致YES暂时失去也是家常便饭
要有耐心

好的,谢谢社长指导。

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 06:52 , Processed in 0.173242 second(s), 24 queries , Gzip On.

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