计算化学公社

标题: 求助：cp2k带晶胞参数的结构优化未收敛却正常结束 [打印本页]

作者
Author: Penson 时间: 6 hour ago
标题: 求助：cp2k带晶胞参数的结构优化未收敛却正常结束
对附件中的cif结构进行带晶胞参数的结构优化，一开始用1E-6的能量精度很难能量收敛，改为1E-5后意外发现收敛速度很快，检查确认后发现结构优化正常结束，但无论是力参数还是压力参数都未达收敛要求，却认为达到收敛目标了。请大佬们帮忙判断一下问题出在哪里了

(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

作者
Author: kirSTY 时间: 4 hour ago
本帖最后由 kirSTY 于 2026-2-22 10:21 编辑

建议参考
1. CP2K对manual对L-BFGS的介绍
2. GOOGLE上对于相同问题的讨论
3. CP2K mailing list对于相同问题的讨论

其核心原因是L-BFGS判据与CP2K内部判据准则不同，程序在满足L-BFGS判据后会提前退出。可以尝试设置一个L-BFGS的额外判据，保证CP2K达到内部收敛准则之前不会因为满足L-BFGS判据而停止计算。
另外你这么多核心，但是每个电子步计算用时都很长啊，附上我自己的在7950X平台上的计算片段供你参考。

作者
Author: Penson 时间: 4 hour ago

kirSTY 发表于 2026-2-22 10:18
建议参考
1. CP2K对manual对L-BFGS的介绍
2. GOOGLE上对于相同问题的讨论

谢谢，我的体系由于扩胞后原子数相对比较多（500左右），所以每一个电子步用时不低。
再请教一下，如果我需要给L-BFGS加额外判据，是需要在.inp文件中其子关键词增加吗？或者说如果我换成CG的优化器，是不是就不存在这种问题？

作者
Author: kirSTY 时间: 4 hour ago

Penson 发表于 2026-2-22 10:42
谢谢，我的体系由于扩胞后原子数相对比较多（500左右），所以每一个电子步用时不低。
再请教一下，如果 ...

- 参考manual加入相关关键词即可。不过这个关键词判据设置应该相当苛刻，为的是避免计算提前退出。
- 更换算法应当可以从根源上避免这种情况出现

作者
Author: Penson 时间: 4 hour ago

kirSTY 发表于 2026-2-22 10:51
- 参考manual加入相关关键词即可。不过这个关键词判据设置应该相当苛刻，为的是避免计算提前退出。
- 更 ...

好的，非常感谢

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