计算化学公社

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

[CP2K] 加了隐式溶剂(水)的能量计算求助

[复制链接 Copy URL]

55

帖子

0

威望

691

eV
积分
746

Level 4 (黑子)

老师好,我用cp2k算加了隐式溶剂的能量,但是最后却显示这样的结果,麻烦老师帮忙看一下是哪里出了问题

202406291328541499..png (42.74 KB, 下载次数 Times of downloads: 9)

202406291328541499..png

1-3-zzzzLTA-SH-B-2Sr6-2Sr8-A-1.out

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

1-3-zzzzLTA-SH-B-2Sr6-2Sr8-A-1.inp

17.47 KB, 下载次数 Times of downloads: 8

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

2#
发表于 Post on 2024-6-30 01:20:13 | 只看该作者 Only view this author
这个体系最好用Gaussian等量化程序算,省心得多。CP2K的SCCS很容易遇上SCF难收敛

如果非要用CP2K就只能自己折腾SCF收敛相关的参数和SCCS参数了

北京科音自然科学研究中心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!

55

帖子

0

威望

691

eV
积分
746

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-6-30 20:27:18 | 只看该作者 Only view this author
sobereva 发表于 2024-6-30 01:20
这个体系最好用Gaussian等量化程序算,省心得多。CP2K的SCCS很容易遇上SCF难收敛

如果非要用CP2K就只能 ...

好的,谢谢老师,我用的是水溶液,multiwfn产生的SCCS中对应的参数应该也是水的参数,我只需要将SCF的参数修改一下

40

帖子

1

威望

732

eV
积分
792

Level 4 (黑子)

4#
发表于 Post on 2024-7-1 23:26:48 | 只看该作者 Only view this author
小可几何 发表于 2024-6-30 20:27
好的,谢谢老师,我用的是水溶液,multiwfn产生的SCCS中对应的参数应该也是水的参数,我只需要将SCF的参 ...

还有判定溶质边界的参数RHO_MAX和RHO_MIN。
但是cp2k的SCCS很奇怪,就算你调参把单点计算搞收敛了,你用相同的参数做结构优化也很容易遇到能量跳变然后不收敛。感觉cp2k开发组就没有很用心在SCCS上,andreussi的environ插件中的SCCS感觉就比较稳健,就是只能qe用。不知道是基组问题还是其他什么问题

55

帖子

0

威望

691

eV
积分
746

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2024-7-2 15:10:26 | 只看该作者 Only view this author
ultramanm87 发表于 2024-7-1 23:26
还有判定溶质边界的参数RHO_MAX和RHO_MIN。
但是cp2k的SCCS很奇怪,就算你调参把单点计算搞收敛了,你用 ...

话说可以直接用气相能量代替液相能量吗?尤其是在隐式溶剂为水的情况下

55

帖子

0

威望

691

eV
积分
746

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2024-7-2 15:14:40 | 只看该作者 Only view this author
ultramanm87 发表于 2024-7-1 23:26
还有判定溶质边界的参数RHO_MAX和RHO_MIN。
但是cp2k的SCCS很奇怪,就算你调参把单点计算搞收敛了,你用 ...

请问RHO_MAX和RHO_MIN应该如何设置,我直接用multiwfn输出的如图片这样,按道理来说应该是将MAX-MIN的间隔调的大一些,对吗?

202407021514347150..png (5.06 KB, 下载次数 Times of downloads: 11)

202407021514347150..png

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

7#
发表于 Post on 2024-7-2 19:40:37 | 只看该作者 Only view this author
小可几何 发表于 2024-7-2 15:14
请问RHO_MAX和RHO_MIN应该如何设置,我直接用multiwfn输出的如图片这样,按道理来说应该是将MAX-MIN的间 ...

参考北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)的ppt




北京科音自然科学研究中心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!

55

帖子

0

威望

691

eV
积分
746

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2024-7-2 22:07:30 | 只看该作者 Only view this author
sobereva 发表于 2024-7-2 19:40
参考北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)的ppt


...

谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 10:38 , Processed in 0.199836 second(s), 24 queries , Gzip On.

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