计算化学公社

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

[CP2K] 求助:大体系PBE0计算单点能不收敛

[复制链接 Copy URL]

46

帖子

0

威望

378

eV
积分
424

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
大家好,我们想要计算吸附在金表面的纳米带的电子结构,使用杂化泛函PBE0,基组为DZVP-MOLOPT-SR-GTH,利用PBE得到restart文件作为初猜。输入文件借助Multiwfn生成
但是在计算的时候收敛值反而越变越大?是我的输入文件有问题吗?



6-40-gfn1-neworder-3au-less-box-move-PBE0-ADMM.inp

54.09 KB, 下载次数 Times of downloads: 11

out.txt

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

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2025-2-15 07:25:52 | 只看该作者 Only view this author
不要用PRINT_LEVEL HIGH,输出信息看着费劲死了,还会浪费耗时
能考虑的做法这里该说的都说了
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665http://bbs.keinsci.com/thread-37196-1-1.html
输入文件设置没明显问题,但Y方向只有20埃,对于不考虑k点的情况来说计算导体不合适
北京科音自然科学研究中心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

46

帖子

0

威望

378

eV
积分
424

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-2-15 09:31:33 | 只看该作者 Only view this author
sobereva 发表于 2025-2-15 07:25
不要用PRINT_LEVEL HIGH,输出信息看着费劲死了,还会浪费耗时
能考虑的做法这里该说的都说了
CP2K中遇到 ...

老师,对于这种体系,cp2k的一些细节我把握不住,要是我用VASP-PBE0去算,会好点吗

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

4#
发表于 Post on 2025-2-15 10:55:24 | 只看该作者 Only view this author
zhouchen 发表于 2025-2-15 09:31
老师,对于这种体系,cp2k的一些细节我把握不住,要是我用VASP-PBE0去算,会好点吗

VASP哪算得动这么大体系
北京科音自然科学研究中心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

894

帖子

3

威望

1684

eV
积分
2638

Level 5 (御坂)

傻傻的木瓜

5#
发表于 Post on 2025-2-15 20:14:29 | 只看该作者 Only view this author
zhouchen 发表于 2025-2-15 09:31
老师,对于这种体系,cp2k的一些细节我把握不住,要是我用VASP-PBE0去算,会好点吗

平时PRINT_LEVEL LOW(最多MEDIUM)的输出有什么细节是你能把握、比较关心但又遗漏的么?

另外我也好奇一下,为啥这个模型里纳米带只有一半吸附在金表面,另一半直立在真空中,顶上还有一个若干金原子的团簇……
√546=23.36664289109

46

帖子

0

威望

378

eV
积分
424

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-2-16 18:35:37 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-2-15 20:14
平时PRINT_LEVEL LOW(最多MEDIUM)的输出有什么细节是你能把握、比较关心但又遗漏的么?

另外我也好 ...

我看cp2k用了PRINT_LEVEL HIGH后会产生ELF.cube文件,我为了得到这个文件,所以用的HIGH
关于我们的体系是这样的,顶上的金原子是金锥电极,金锥电极将完全吸附在平板电极上的纳米带给提拉起来,我给出的文件结构只是其中的一个过程结构

894

帖子

3

威望

1684

eV
积分
2638

Level 5 (御坂)

傻傻的木瓜

7#
发表于 Post on 2025-2-16 22:38:23 | 只看该作者 Only view this author
zhouchen 发表于 2025-2-16 18:35
我看cp2k用了PRINT_LEVEL HIGH后会产生ELF.cube文件,我为了得到这个文件,所以用的HIGH
关于我们的体系 ...

完全没必要只为这个目的用PRINT_LEVEL HIGH。在Multiwfn创建CP2K输入文件时直接用选项-3 Set exporting cube file的子选项2就能产生相应文件,或者产生的molden载入Multiwfn计算格点数据。这与平时计算细节的把握关系不大。

感觉应该考虑一下两个金电极附近的外电场,这种外势也能通过函数式或者.cube文件在CP2K里设置,不过写起来可能比较复杂,而且盒子z方向尺寸也得大到足以消除相邻镜像间电极多余的相互作用。
√546=23.36664289109

46

帖子

0

威望

378

eV
积分
424

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2025-3-6 09:29:30 | 只看该作者 Only view this author
sobereva 发表于 2025-2-15 07:25
不要用PRINT_LEVEL HIGH,输出信息看着费劲死了,还会浪费耗时
能考虑的做法这里该说的都说了
CP2K中遇到 ...

老师,我尝试了很久,对角化还是没法收敛,但是OT可以收敛,我也尝试了OT波函数作为初猜,但是还是收敛振荡。我这个体系OT算的结果能用吗

328

帖子

3

威望

1759

eV
积分
2147

Level 5 (御坂)

9#
发表于 Post on 2025-3-9 21:18:28 | 只看该作者 Only view this author
zhouchen 发表于 2025-3-6 09:29
老师,我尝试了很久,对角化还是没法收敛,但是OT可以收敛,我也尝试了OT波函数作为初猜,但是还是收敛振 ...

如果目的是得到ELF,那对角化和OT都行。可以把这俩理解成不同的scf收敛方法,只要都到了足够严的收敛限,就都可以认为是准的。

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 16:42 , Processed in 0.152421 second(s), 23 queries , Gzip On.

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