Molcas 自洽场开壳层(UHF)收敛困难，该如何解决

函数与激情

我计算一分子带一个单电子 (spin=2,charge=0)，使用Molcas 8.0 计算自洽场始终过不去，而将其改成带电体系( spin=1,charge=1) 则可以较快收敛，不知道大家有何解决方法，请诸位不吝赐教，谢谢！

函数与激情

file:///C:/Users/Administrator/Desktop/sp190616_153449.png Manual中提到UHF的收敛情况确实不太给力

jiangning198511

可以考虑以收敛的 spin=1,charge=1的结果为初始猜测试试

zjxitcc

我自己写了个程序（还未开源），可以转换高斯fchk文件里的轨道给OpenMolcas作为初始，这样就可以使OpenMolcas在1-3圈左右收敛（对于RHF, UHF, CASSCF等功能）。众所周知高斯HF收敛性几乎是最为强大的。
如果你最终没解决这个问题的话，可私信找我（要求高斯使用关键词nosymm int=nobasistransform punch=gamess计算得到.fchk文件和fort.7文件）

---

更新：后面这个功能包含在我写的开源程序MOKIT中，称作fch2inporb小程序，可以从高斯.fch(k)文件 生成 OpenMolcas的输入文件.input 和 .INPORB轨道文件，传入轨道后可让OpenMolcas能量立即达到收敛值，接着跑1-2圈就完成了。