计算化学公社

标题: 求助：大体系PBE0计算单点能不收敛 [打印本页]

作者
Author: zhouchen 时间: 2025-2-14 19:25
标题: 求助：大体系PBE0计算单点能不收敛
大家好，我们想要计算吸附在金表面的纳米带的电子结构，使用杂化泛函PBE0，基组为DZVP-MOLOPT-SR-GTH，利用PBE得到restart文件作为初猜。输入文件借助Multiwfn生成
但是在计算的时候收敛值反而越变越大？是我的输入文件有问题吗？

作者
Author: sobereva 时间: 2025-2-15 07:25
不要用PRINT_LEVEL HIGH，输出信息看着费劲死了，还会浪费耗时
能考虑的做法这里该说的都说了
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）
输入文件设置没明显问题，但Y方向只有20埃，对于不考虑k点的情况来说计算导体不合适

作者
Author: zhouchen 时间: 2025-2-15 09:31

sobereva 发表于 2025-2-15 07:25
不要用PRINT_LEVEL HIGH，输出信息看着费劲死了，还会浪费耗时
能考虑的做法这里该说的都说了
CP2K中遇到 ...

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

作者
Author: sobereva 时间: 2025-2-15 10:55

zhouchen 发表于 2025-2-15 09:31
老师，对于这种体系，cp2k的一些细节我把握不住，要是我用VASP-PBE0去算，会好点吗

VASP哪算得动这么大体系

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-2-15 20:14

zhouchen 发表于 2025-2-15 09:31
老师，对于这种体系，cp2k的一些细节我把握不住，要是我用VASP-PBE0去算，会好点吗

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

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

作者
Author: zhouchen 时间: 2025-2-16 18:35

Uus/pMeC6H4-/キ发表于 2025-2-15 20:14
平时PRINT_LEVEL LOW（最多MEDIUM）的输出有什么细节是你能把握、比较关心但又遗漏的么？

另外我也好 ...

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

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-2-16 22:38

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方向尺寸也得大到足以消除相邻镜像间电极多余的相互作用。

作者
Author: zhouchen 时间: 2025-3-6 09:29

sobereva 发表于 2025-2-15 07:25
不要用PRINT_LEVEL HIGH，输出信息看着费劲死了，还会浪费耗时
能考虑的做法这里该说的都说了
CP2K中遇到 ...

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

作者
Author: Stardust0831 时间: 2025-3-9 21:18

zhouchen 发表于 2025-3-6 09:29
老师，我尝试了很久，对角化还是没法收敛，但是OT可以收敛，我也尝试了OT波函数作为初猜，但是还是收敛振 ...

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

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