计算化学公社

标题: Cu(111)吸附铀酰离子的k点与cutoff收敛性测试值的选择问题 [打印本页]

作者
Author: pikachuupup 时间: 2025-9-5 16:11
标题: Cu(111)吸附铀酰离子的k点与cutoff收敛性测试值的选择问题
本帖最后由 pikachuupup 于 2025-9-6 14:22 编辑

请问各位老师，我想利用cp2k做一下关于Cu的111晶面吸附铀酰离子UO22+吸附能的计算，首先优化了Cu的原胞，然后切面扩胞并固定底层原子后优化超胞，现在来到了表面上面负载UO22+这一步，有一些疑问想请教老师们：

1、我建好模之后用脚本进行k点收敛性测试，发现当达到7,7,1时才满足1 meV/atom的收敛标准，目前超胞已经挺多原子了（111个），请问是我的输入文件有什么硬伤吗，还是正常现象呢？

2、另外我又进行了cutoff的收敛性测试，刚开始在rel_cutoff=60的情况下计算，结果为CUTOFF-1.txt，发现从550到900都不满足1E-6 Hatree/atom的收敛标准，后续又在rel_cutoff=70的情况下进行了第二次和第三次测试，结果分别为CUTOFF-2.txt和CUTOFF-3.txt，但是也没有达到收敛标准，且呈现波动趋势。当前这个k点和cutoff已经很耗时了，请问接下来还要继续加大cutoff进行测试吗？

作者
Author: sobereva 时间: 2025-9-6 15:56
注意检查各个输出文件里的能量收敛情况，如果能量收敛精度还不够充分时就被判断为了SCF收敛，则收敛性测试没有意义
没必要教条式地做收敛性测试。你的情况4,4,1 k点结合CUTOFF 500完全够了。而且实际上更应该测试的是吸附能的收敛性，而不是整个体系能量的收敛性。

作者
Author: pikachuupup 时间: 2025-9-6 17:27

sobereva 发表于 2025-9-6 15:56
注意检查各个输出文件里的能量收敛情况，如果能量收敛精度还不够充分时就被判断为了SCF收敛，则收敛性测试 ...

明白了，谢谢sob老师~

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