计算化学公社

标题: Molcas 自洽场开壳层(UHF)收敛困难，该如何解决 [打印本页]

作者
Author: 函数与激情 时间: 2019-6-16 15:32
标题: Molcas 自洽场开壳层(UHF)收敛困难，该如何解决
我计算一分子带一个单电子 (spin=2,charge=0)，使用Molcas 8.0 计算自洽场始终过不去，而将其改成带电体系( spin=1,charge=1) 则可以较快收敛，不知道大家有何解决方法，请诸位不吝赐教，谢谢！

作者
Author: 函数与激情 时间: 2019-6-16 15:35
file:///C:/Users/Administrator/Desktop/sp190616_153449.png Manual中提到UHF的收敛情况确实不太给力

作者
Author: jiangning198511 时间: 2019-6-17 14:24
可以考虑以收敛的 spin=1,charge=1的结果为初始猜测试试

作者
Author: zjxitcc 时间: 2019-6-19 09:49
本帖最后由 zjxitcc 于 2021-2-13 00:58 编辑

我自己写了个程序（还未开源），可以转换高斯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圈就完成了。

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