计算化学公社

标题: 求助，使用对角化方法，对使用OT进行结构优化后的结构求molden文件scf不收敛 [打印本页]

作者
Author: HYX1221 时间: 2026-4-26 13:29
标题: 求助，使用对角化方法，对使用OT进行结构优化后的结构求molden文件scf不收敛
本帖最后由 HYX1221 于 2026-4-27 08:21 编辑

如题，本人研究的体系为镁铝尖晶石表面对铁原子和硫化锰的吸附，使用OT进行结构优化，后因要绘制DOS图，并且必须考虑空轨道，所以必须使用对角化来求molden文件绘制DOS图。于是使用OT收敛后的结构，并使用收敛后的波函数作为初猜进行计算。
但是再过程中遇到了scf难收敛的问题。
自己尝试过一些解决scf难收敛的办法：

改&MIXING中的ALPHA为0.1、NBROYDEN为16
提升电子温度（默认的300K提升到500K）
先放宽收敛精度，待收敛后逐步缩小收敛精度进行续算（这个过程中，最多收敛到1e-03，然后就降不下去了）

并且还尝试了用OT收敛后的结构、用收敛后的波函数作为初猜，用对角化再次进行结构优化，scf没有收敛趋势。
想请教一下还能使用什么方法来尝试。
附件中展示的为镁铝尖晶石表面稳定吸附Fe原子后在其顶位吸附Mn原子的case，我所有涉及Mn、S原子吸附的case都涉及以上问题。

作者
Author: reid 时间: 2026-5-3 23:19
这是cp2k使用对角化的痛点，经常性的碰上不收敛

作者
Author: 喝酸奶的巧乐兹 时间: 2026-5-7 00:01
开始时先忽略scf不收敛而继续优化几步看看，看初始优化的几步会不会让convergence变小

作者
Author: HYX1221 时间: 2026-5-7 19:33

喝酸奶的巧乐兹发表于 2026-5-7 00:01
开始时先忽略scf不收敛而继续优化几步看看，看初始优化的几步会不会让convergence变小

无论是直接求molden还是用对角化进行结构优化，covergence均最多只能缩小到1e-03。

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