计算化学公社

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

[CP2K] CP2K计算molden文件不收敛问题

[复制链接 Copy URL]

21

帖子

0

威望

97

eV
积分
118

Level 2 能力者

本帖最后由 zhifanfan 于 2025-1-3 16:25 编辑

卢老师您好,各位老师好,我在做cp2k产生molden文件,参考了您的651博文,同时想要做类似于716博文的分析,这里建了一个模型,发现无论怎么样使用对角化方法均不能收敛,OT方法很快能够收敛,体系共200左右原子,而您的716博文文件在我的服务器上几秒钟即可完成收敛,查看设置和您的inp文件基本相似,还请卢老师指点,需要如何修改。现附上输入文件以及模型。

SiO2.zip

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

输入源文件

O2Si.inp

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

更新后的输入文件

21

帖子

0

威望

97

eV
积分
118

Level 2 能力者

2#
 楼主 Author| 发表于 Post on 2025-1-3 16:27:06 | 只看该作者 Only view this author
在这个过程中,使用过6-31G*和dzvp基组,都不能收敛,对于构建的SiO2模型

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

3#
发表于 Post on 2025-1-3 21:35:08 | 只看该作者 Only view this author
看清楚板块再发贴,CP2K不是量子化学程序,甭发到量子化学版块。这次给你移动了,以后注意

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了。
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

4#
发表于 Post on 2025-1-3 21:36:09 | 只看该作者 Only view this author
我不知道你想算什么,如果是当体相晶体算,盒子尺寸不对,都不是三维周期性。如果当slab算,真空区太小,并且底端的O应当加氢饱和
不要把DZVP-MOLOPT-SR-GTH简写为dzvp,有严重歧义
北京科音自然科学研究中心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

21

帖子

0

威望

97

eV
积分
118

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2025-1-4 09:54:39 | 只看该作者 Only view this author
sobereva 发表于 2025-1-3 21:35
看清楚板块再发贴,CP2K不是量子化学程序,甭发到量子化学版块。这次给你移动了,以后注意

有别人回复之 ...

收到,卢老师,我下次注意,就是我想实现您的716博文http://sobereva.com/716这种使用Multiwfn结合CP2K对周期性体系做电荷分解分析(CDA),用您给的文件,NaCl_CO文件,发现您构建的体系是能够实现收敛的,而且您的体系能够实现周期性的同时还能实现留白,也就是很大的空白地方没有原子,而我构建的SiO2体系实现了类似您那种晶胞边界的扩展后,就不再是三维周期性体系了,我也不知道该怎么操作,我也是想计算某个分子的加入从而对SiO2表面电荷的转移的影响

21

帖子

0

威望

97

eV
积分
118

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2025-1-4 09:59:03 | 只看该作者 Only view this author
本帖最后由 zhifanfan 于 2025-1-4 10:00 编辑
sobereva 发表于 2025-1-3 21:36
我不知道你想算什么,如果是当体相晶体算,盒子尺寸不对,都不是三维周期性。如果当slab算,真空区太小,并 ...

卢老师,这个是我用GView调整好的晶胞,感觉就是不对,放到cp2k里计算也报错,E:\SiO2.jpg而我用这个方式构建的晶胞E:\2.jpg,这个就可以跑通

SiO2.jpg (55.4 KB, 下载次数 Times of downloads: 2)

调整晶胞边界

调整晶胞边界

2.jpg (129.45 KB, 下载次数 Times of downloads: 1)

未调整

未调整

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

7#
发表于 Post on 2025-1-5 01:37:25 | 只看该作者 Only view this author
zhifanfan 发表于 2025-1-4 09:59
卢老师,这个是我用GView调整好的晶胞,感觉就是不对,放到cp2k里计算也报错,而我用这个方式构建的晶胞 ...

重新判断成键
O加上氢进行饱和
问报错时始终主动说清楚什么报错
北京科音自然科学研究中心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

21

帖子

0

威望

97

eV
积分
118

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2025-1-5 09:44:51 | 只看该作者 Only view this author
sobereva 发表于 2025-1-5 01:37
重新判断成键
O加上氢进行饱和
问报错时始终主动说清楚什么报错

谢谢卢老师,您说的氧上加氢进行饱和我有点不是很明白,同时,若是在氧上加氢饱和了,在最表层的Si不也是存在未满配的情况吗,这个操作确实不会做了。

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

9#
发表于 Post on 2025-1-5 21:27:13 | 只看该作者 Only view this author
zhifanfan 发表于 2025-1-5 09:44
谢谢卢老师,您说的氧上加氢进行饱和我有点不是很明白,同时,若是在氧上加氢饱和了,在最表层的Si不也是 ...

氧不饱和会有奇怪的悬键,SCF难收敛或收敛得巨恶心

怎么处理模型取决于你要研究什么,把这最重要的前提先交代清楚
北京科音自然科学研究中心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

21

帖子

0

威望

97

eV
积分
118

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2025-1-5 21:47:22 | 只看该作者 Only view this author
sobereva 发表于 2025-1-5 21:27
氧不饱和会有奇怪的悬键,SCF难收敛或收敛得巨恶心

怎么处理模型取决于你要研究什么,把这最重要的前 ...

好的,辛苦老师了,我打算把最近的问题重新总结一下,发一个求助稿,第一次求助,有很多地方有疏漏,还请老师原谅

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 13:04 , Processed in 0.192786 second(s), 24 queries , Gzip On.

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