计算化学公社

标题: CP2K计算molden文件不收敛问题 [打印本页]

作者
Author: zhifanfan 时间: 2025-1-3 14:13
标题: CP2K计算molden文件不收敛问题
本帖最后由 zhifanfan 于 2025-1-3 16:25 编辑

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

作者
Author: zhifanfan 时间: 2025-1-3 16:27
在这个过程中，使用过6-31G*和dzvp基组，都不能收敛，对于构建的SiO2模型

作者
Author: sobereva 时间: 2025-1-3 21:35
看清楚板块再发贴，CP2K不是量子化学程序，甭发到量子化学版块。这次给你移动了，以后注意

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

作者
Author: sobereva 时间: 2025-1-3 21:36
我不知道你想算什么，如果是当体相晶体算，盒子尺寸不对，都不是三维周期性。如果当slab算，真空区太小，并且底端的O应当加氢饱和
不要把DZVP-MOLOPT-SR-GTH简写为dzvp，有严重歧义

作者
Author: zhifanfan 时间: 2025-1-4 09:54

sobereva 发表于 2025-1-3 21:35
看清楚板块再发贴，CP2K不是量子化学程序，甭发到量子化学版块。这次给你移动了，以后注意

有别人回复之 ...

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

作者
Author: zhifanfan 时间: 2025-1-4 09:59
本帖最后由 zhifanfan 于 2025-1-4 10:00 编辑

sobereva 发表于 2025-1-3 21:36
我不知道你想算什么，如果是当体相晶体算，盒子尺寸不对，都不是三维周期性。如果当slab算，真空区太小，并 ...

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

作者
Author: sobereva 时间: 2025-1-5 01:37

zhifanfan 发表于 2025-1-4 09:59
卢老师，这个是我用GView调整好的晶胞，感觉就是不对，放到cp2k里计算也报错，而我用这个方式构建的晶胞 ...

重新判断成键
O加上氢进行饱和
问报错时始终主动说清楚什么报错

作者
Author: zhifanfan 时间: 2025-1-5 09:44

sobereva 发表于 2025-1-5 01:37
重新判断成键
O加上氢进行饱和
问报错时始终主动说清楚什么报错

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

作者
Author: sobereva 时间: 2025-1-5 21:27

zhifanfan 发表于 2025-1-5 09:44
谢谢卢老师，您说的氧上加氢进行饱和我有点不是很明白，同时，若是在氧上加氢饱和了，在最表层的Si不也是 ...

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

怎么处理模型取决于你要研究什么，把这最重要的前提先交代清楚

作者
Author: zhifanfan 时间: 2025-1-5 21:47

sobereva 发表于 2025-1-5 21:27
氧不饱和会有奇怪的悬键，SCF难收敛或收敛得巨恶心

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)