求助cp2k计算铁磁性材料磷酸铁锂LFP几何优化不收敛

流水先生

各位老师好，最近使用cp2k进行磷酸铁锂的几何优化，使用diagonalization算法，考虑了Fe元素的自旋，发现无法收敛。

通过查看帖子
CP2K计算时，材料磁性如何确定？磁矩、多重度及RELAX_MULTIPLICITY该如何设置？ - 第一性原理 (First Principle) - 计算化学公社 (keinsci.com)
CP2K中遇到SCF难收敛时的解决方法 - 第一性原理 (First Principle) - 计算化学公社 (keinsci.com)
加速SCF收敛可以尝试的一些手段 - 第一性原理 (First Principle) - 计算化学公社 (keinsci.com)
CP2K优化磷酸锰铁锂晶胞时SCF不收敛 - 第一性原理 (First Principle) - 计算化学公社 (keinsci.com)
尝试修改MAGNETIZATION  5/3.5、MULTIPLICITY、ALPHA 0.2、NBROYDEN 12添加RELAX_MULTIPLICITY 1E-4始终无法收敛，Convergence在0.01上下波动。
根据materials project给出的INCAR中 MAGMOM = 4*0.6 4*5.0 20*0.6修改Li Fe P O的MAGNETIZATION 为对应值，Convergence会小一点，在0.001上下波动，仍然不收敛。
请问该体系使用diagonalization算法有其他促使收敛的设置吗？
目前还没有尝试开启OT
请各位老师指点。

sobereva

CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

Priscilla

OT的确在某些体系中更容易收敛，但是有些体系即便用OT收敛的波函数跑对角化照样不收敛，如果设完磁矩没啥用考虑下BS，最后直接上猛药smearing开个300-600，之所以叫smearing猛药是因为确实解决scf不收敛蛮有效的，但是也有弊，具体这个弊有多大得看你体系了。另外如果只是优化早期不收敛很正常啊，减少SCF上限步数让它多跑跑可能就收敛了。