计算化学公社

标题: 求助Molpro做cas不收敛解决方法 [打印本页]

作者
Author: tongxin8045 时间: 2021-8-26 17:00
标题: 求助Molpro做cas不收敛解决方法
初用Molpro，计算含有重金属原子的cas单点计算，提示不收敛，请大侠们帮忙看看，万分感谢！输入命令和最后的结果提示见截屏。
PS：第一次来论坛，有不规范之处，请指正。

作者
Author: wzkchem5 时间: 2021-8-26 18:38
先把CAS迭代最后一步的活性轨道作图出来看，看是否符合你的预期。如果不符合预期，调整轨道使其符合预期以后再接着算；如果符合预期，那么从最后一步的轨道开始，加一点level shift再试；如果不知道轨道理应是什么样子的，先查文献看类似体系的活性轨道是什么样子的，然后按照文献设置

作者
Author: beefly 时间: 2021-8-26 18:48
molpro的casscf对初始轨道的要求很高，uhf不是好的选择。建议用闭壳层的阴、阳离子做rhf，或者小活性空间的casscf。
为了避免不同初始轨道的芯轨道不同导致能量差异，需要把正式casscf计算中的frozen轨道设为0

作者
Author: zjxitcc 时间: 2021-8-26 19:31
本帖最后由 zjxitcc 于 2021-8-26 19:34 编辑

不妨试试自动确定活性空间《自动做多参考态计算的程序MOKIT》。可以自动调用Molpro做CAS计算，举个例子

%mem=64GB
%nprocshared=32
#p CASSCF/def2TZVP
mokit{CASSCF_prog=Molpro}
0 1
O -0.23497692 0.90193619 -0.068688
H 1.26502308 0.90193619 -0.068688
H -0.73568721 2.31589843 -0.068688

复制代码

作者
Author: tongxin8045 时间: 2021-8-27 09:35

wzkchem5 发表于 2021-8-26 18:38
先把CAS迭代最后一步的活性轨道作图出来看，看是否符合你的预期。如果不符合预期，调整轨道使其符合预期以 ...

谢谢，发现现在不收敛的还是uhf那一步，目前使用了rhf,还是不收敛。

作者
Author: tongxin8045 时间: 2021-8-27 09:37

beefly 发表于 2021-8-26 18:48
molpro的casscf对初始轨道的要求很高，uhf不是好的选择。建议用闭壳层的阴、阳离子做rhf，或者小活性空间的 ...

谢谢，咨询下，体系中有重原子，现在没有使用赝势，是否需要考虑赝势呢？如果考虑了赝势，总的电子数是不是就要扣除赝势的核数？

作者
Author: tongxin8045 时间: 2021-8-27 09:38

zjxitcc 发表于 2021-8-26 19:31
不妨试试自动确定活性空间《自动做多参考态计算的程序MOKIT》。可以自动调用Molpro做CAS计算，举个例子

...

谢谢，我先了解下，初用molpro，好多设置还没有整太明白

作者
Author: zjxitcc 时间: 2021-8-27 12:54

tongxin8045 发表于 2021-8-27 09:38
谢谢，我先了解下，初用molpro，好多设置还没有整太明白

MOKIT自动调用molpro做计算的好处就是：你不需要懂molpro；你也不需要手动看、选活性轨道。

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