计算化学公社

标题: 关于Cr@C60的不同自旋多重度的构型优化SCF不收敛求助 [打印本页]

作者
Author: JTDJZS 时间: 2025-7-3 09:57
标题: 关于Cr@C60的不同自旋多重度的构型优化SCF不收敛求助
本帖最后由 JTDJZS 于 2025-7-4 16:57 编辑

老师们好，我在Gaussian中对Cr@C60不同自旋多重度的构型进行几何优化，使用的是TPSS泛函，和赝势基组（6-31g*~SDD）,分别考虑了单、三、五、七重态，目前是五重态能量最低，因此在对七重态进行优化，但是尝试了解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元帖子中的SCF(novaracc,noincfock) int=acc2e=12，SCF=vshift=x，SCF=conver=6，SCF=Fermi，并且在b3lyp泛函下可以得到收敛的几何优化结果，以这个时候的chk文件作为初猜，依旧不收敛，又尝试SCF=xqc，stable=opt，得到了收敛的波函数，但是用该波函数进行优化，也依旧不收敛。我还尝试调整了Cr在C60中的几何位置，都会向碳笼中心移动，但是也都会在靠近中心时不收敛，观察Optimization，发现前面的优化步骤能量都还在下降，但是一到最后一步就回突然上升并不收敛，又拿倒数第二步的结构继续进行优化，同样不收敛，以下为stable=opt，和以此为初猜的输入输出文件，想请问各位老师还有什么建议，谢谢。

作者
Author: wzkchem5 时间: 2025-7-3 10:33
尝试用ORCA的TRAH方法收敛波函数，再把波函数用Multiwfn或MOKIT传给高斯

作者
Author: sobereva 时间: 2025-7-5 08:18
本来就不推荐用TPSS。纯泛函gap偏低，比杂化泛函往往更难收敛。而且对比自旋态能量差的目的，明显更建议用TPSSh。

如果TPSSh也不收敛，若发现PBE0对七重态优化能收敛，读取它的优化完的结构和波函数当做TPSSh优化的初猜。

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