计算化学公社

标题: Mg、Bi组成的体系做CP2K的变胞优化，第一步时SCF不收敛且震荡 [打印本页]

作者
Author: Ying_Zhang 时间: 2026-2-5 19:12
标题: Mg、Bi组成的体系做CP2K的变胞优化，第一步时SCF不收敛且震荡
我的体系就是直接从Material Projects直接下载的结构，进行一个简单的扩胞，然后用CP2K进行变胞优化。结果达到优化上限，第一步也没收敛，且 convergence也没有缩小趋势，请问各位老师，如何修改输入文件，使其收敛？

作者
Author: sobereva 时间: 2026-2-5 19:39
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

作者
Author: Ying_Zhang 时间: 2026-2-5 22:53

sobereva 发表于 2026-2-5 19:39
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.ht ...

感谢sob老师回复，我尝试了一下里面提到的方法，现在正在运行。但是，想请教一下老师，我的convergce里面数值过大，跑了两步了，均在20到30左右，是否存在别的问题？且这是我拿VASP已计算结构，VAP相同结构仅需要跑1小时，结束收敛，相同配置情况下，我现在快2个小时，还没收敛。

作者
Author: yxdd98 时间: 2026-2-5 23:14

Ying_Zhang 发表于 2026-2-5 22:53
感谢sob老师回复，我尝试了一下里面提到的方法，现在正在运行。但是，想请教一下老师，我的convergce里面 ...

看了你的原始输入，你现在用的基组对Mg要求的截断能很高，大概是1000以上吧，推荐改用全电子基组试试

作者
Author: Ying_Zhang 时间: 2026-2-6 00:25

yxdd98 发表于 2026-2-5 23:14
看了你的原始输入，你现在用的基组对Mg要求的截断能很高，大概是1000以上吧，推荐改用全电子基组试试

哦哦好的我尝试一下，我之前也有这个想法，如果按照经验来说，CUTOFF得1000以上。后来我进行测试，发现700左右就行，结构现在还是不收敛，看来还是不行。感谢老师回复。

作者
Author: ultramanm87 时间: 2026-2-6 12:13
我直接拿你的输入文件跑，没有什么问题，版本cp2k2026.1 (, 下载次数 Times of downloads: 1)

作者
Author: UW_0728. 时间: 2026-2-6 13:03
本帖最后由 UW_0728. 于 2026-2-6 13:43 编辑

我这里的情况和6L的一样，CP2K v2026.1跑的，很快就收敛了。
(, 下载次数 Times of downloads: 0)

但其实2025.2和2026.1版本没有什么能够对这个造成显著影响的更改，而且我也拿2025.1验证了，结果跟上面输出的一模一样，都是很快收敛。

所以大概率是你编译或运行CP2K加载环境等方面的问题，CP2K在实际计算中没有正确利用数学库，或者你编译时使用了不兼容当前版本CP2K的数学库，或者是运行时加载的数学库跟编译时用的不严格一致，等等（但你这种情况真的实属罕见）

另外还有两个问题：
1. 单纯好奇，为什么要扩胞做变胞优化呢？
2. 考虑了k点，最好把EXTRAPOLATION USE_PREV_P的注释取消掉，让程序能够利用前一步收敛的波函数当初猜，可以大大降低优化用时。

作者
Author: Ying_Zhang 时间: 2026-2-6 16:25

UW_0728. 发表于 2026-2-6 13:03
我这里的情况和6L的一样，CP2K v2026.1跑的，很快就收敛了。

是这样的老师，扩胞做变胞优化是我想测试一下CP2K的计算效率。我这个体系在VASP已经计算过了，但是因为后续计算任务，是把相同材料的结构扩展到上千个原子去计算，VASP计算能力恐怕完成不了，所以我才又学的CP2K去完成我的目标。因为也是刚使用CAP2K，所以一是为了测试CP2K的计算效率，二是为了熟悉一下这个过程。EXTRAPOLATION USE_PREV_P这个我注释掉了，还是没计算成功

我再看看我编译环境

作者
Author: Ying_Zhang 时间: 2026-2-6 16:29

ultramanm87 发表于 2026-2-6 12:13
我直接拿你的输入文件跑，没有什么问题，版本cp2k2026.1

感谢回复，想问一下，您一共用了几个优化步完成这个结构优化的？

作者
Author: Ying_Zhang 时间: 2026-2-6 18:09

UW_0728. 发表于 2026-2-6 13:03
我这里的情况和6L的一样，CP2K v2026.1跑的，很快就收敛了。

老师，我刚才排除了一下问题，又发了一个帖子详细的说了一下，链接在这里http://bbs.keinsci.com/thread-58188-1-1.html，麻烦请您帮忙看看能否解决。

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