Mg、Bi组成的体系做CP2K的变胞优化，第一步时SCF不收敛且震荡

Ying_Zhang

我的体系就是直接从Material Projects直接下载的结构，进行一个简单的扩胞，然后用CP2K进行变胞优化。结果达到优化上限，第一步也没收敛，且 convergence也没有缩小趋势，请问各位老师，如何修改输入文件，使其收敛？

sobereva

CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

Ying_Zhang

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个小时，还没收敛。

yxdd98

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

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

Ying_Zhang

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

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

ultramanm87

我直接拿你的输入文件跑，没有什么问题，版本cp2k2026.1

UW_0728.

我这里的情况和6L的一样，CP2K v2026.1跑的，很快就收敛了。


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

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

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

Ying_Zhang

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

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

Ying_Zhang

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

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

Ying_Zhang

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

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