求助CP2K计算RuO₂ (110) 表面吸附OER中间体时的SCF收敛问题

蹇子铭

卢老师您好，我最近在用 CP2K 算 RuO₂ (110) 表面吸附 OER 中间体的时候，遇到一些 SCF 收敛的问题，想请教您看看是不是我设置有问题。我这个体系是 RuO₂ (110) 表面，截了三层，计算的时候固定了最下面一层，放开上面两层。用的是对角化方法，输入文件是通过 Multiwfn 生成的。对于光面催化剂（*）和吸附 O 的中间体（*O），SCF 都能正常收敛。但是一旦涉及到含 H 的中间体，比如 *OH 和 *OOH，SCF 就特别慢，而且经常震荡收敛不了。我试过把 MAX_SCF 调到 500，Alpha 设成 0.1，优化器也换成了 CG，效果还是不太理想。因为体系是导体，我就没用 OT 算法。现在这一步算下来已经花了一天多时间，相比之前用 VASP 算同样的模型，基本一天之内就能收完，感觉 CP2K 并未在效率上体现出优势。所以想请教您看看，我现在的设置有没有不合理的地方？或者有没有什么提高收敛效率的建议？