求助Molpro做cas不收敛解决方法

tongxin8045

初用Molpro，计算含有重金属原子的cas单点计算，提示不收敛，请大侠们帮忙看看，万分感谢！输入命令和最后的结果提示见截屏。
PS：第一次来论坛，有不规范之处，请指正。

wzkchem5

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

beefly

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

zjxitcc

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

1. %mem=64GB
   
2. %nprocshared=32
   
3. #p CASSCF/def2TZVP
   
4. 
   
5. mokit{CASSCF_prog=Molpro}
   
6. 
   
7. 0 1
   
8. O      -0.23497692     0.90193619    -0.068688
   
9. H       1.26502308     0.90193619    -0.068688
   
10. H      -0.73568721    2.31589843    -0.068688

复制代码

tongxin8045

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

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

tongxin8045

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

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

tongxin8045

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

...

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

zjxitcc

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

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