计算化学公社

标题: 求助cp2k分子表面吸附模型优化无法完成，结构一直没有收敛的趋势 [打印本页]

作者
Author: 满天星斗 时间: 2025-12-4 09:40
标题: 求助cp2k分子表面吸附模型优化无法完成，结构一直没有收敛的趋势
本帖最后由满天星斗于 2025-12-8 11:09 编辑

老师们，我想用cp2k优化一个分子的吸附模型。但是这个结构一直没有收敛的趋势，请大家帮忙看看是什么原因呢？输入和输出文件已经上传到附件，谢谢老师们

后续将cufoff增加到650，将步长缩小至0.02和0.01，发现几个参数仍然没有明显的收敛趋势。轨迹文件和输出文件已上传至附件

作者
Author: sobereva 时间: 2025-12-4 11:35
对于含有Na的情况用GTH赝势，CUTOFF=400不够用
加大到650还不行就尝试CG算法，或进一步减小步长上限
主动上传压缩后的轨迹文件，否则没法判断结构收敛情况

作者
Author: ljh123 时间: 2025-12-4 11:39
本帖最后由 ljh123 于 2025-12-4 11:57 编辑

确实需要调整cutoff，还有CG和DIIS可以切换下看看效果

作者
Author: 满天星斗 时间: 2025-12-4 13:33

ljh123 发表于 2025-12-4 11:39
确实需要调整cutoff，还有CG和DIIS可以切换下看看效果

谢谢老师，我增大cutoff看看。另外我这个体系SCF似乎很难收敛，之前试了很多，所以后面我SCF都用的是比较稳妥的方法，将本来Multiwfn默认的SCF收敛限提高了点。之前发现在10-6这个数量级就降得很慢了。想请教下老师建议切换DIIS的原因是什么呀？

作者
Author: 满天星斗 时间: 2025-12-4 13:33

sobereva 发表于 2025-12-4 11:35
对于含有Na的情况用GTH赝势，CUTOFF=400不够用
加大到650还不行就尝试CG算法，或进一步减小步长上限
主动 ...

谢谢老师，我增大cutoff试一试

作者
Author: ljh123 时间: 2025-12-4 14:47

满天星斗发表于 2025-12-4 13:33
谢谢老师，我增大cutoff看看。另外我这个体系SCF似乎很难收敛，之前试了很多，所以后面我SCF都用的是比较 ...

因为SCF不收敛一般试这些。我看了下你的结构，我觉得你可以试一下把材料和分子分开来分别做一下结构优化，然后再放在一起。我以前这样做SCF收敛更容易一些

作者
Author: 满天星斗 时间: 2025-12-5 10:22

ljh123 发表于 2025-12-4 14:47
因为SCF不收敛一般试这些。我看了下你的结构，我觉得你可以试一下把材料和分子分开来分别做一下结构优化 ...

好的，谢谢老师。我后面试试。分子我之前已经用Gaussian优化过，不过单纯的材料表面模型没有优化过

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