计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › Molcas 自洽场开壳层(UHF)收敛困难，该如何解决

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

查看数: 6795 | 评论数: 3 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

函数与激情

发布时间: 2019-6-16 15:32

正文摘要:

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

回复 Reply

zjxitcc 发表于 Post on 2019-6-19 09:49:59

本帖最后由 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圈就完成了。

jiangning198511 发表于 Post on 2019-6-17 14:24:40

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register