如何解决开壳层结构PSI4能量分解不收敛的问题

科研小强

各位大佬，我正在用PSI4软件进行羰基改性碳平面吸附二氧化碳或水构型的能量分解工作，但是现在出现了scf不收敛的状况，尝试多种方法后仍无法解决，求各位大佬不吝赐教，谢谢！
输入文件：
memory 140 GB

molecule dimer {
0 1
C     -5.47807794    0.29661087    2.51685673
O     -6.31698226    1.10121835    2.53109678
O     -4.63512743   -0.50535144    2.51458816
--
0 2
C     -3.21421659    0.70381178   -0.44247470
C     -2.32603417   -0.41787632   -0.36446947
C     -0.92700021   -0.20882423   -0.25926942
C     -0.39954470    1.11809081   -0.22614072
C     -1.28284886    2.23868113   -0.29959806
C     -2.68588340    2.02967518   -0.40794634
C     -0.03821398   -1.32436345   -0.18619197
C      1.00067807    1.32538594   -0.11967919
C      1.88381123    0.21064222   -0.04745725
C      1.36168805   -1.11366272   -0.08125665
C      3.29245273    0.42094348    0.05872622
C      3.81352497    1.74467167    0.09293693
C      2.92981151    2.86205451    0.02140624
C      1.52670328    2.65476209   -0.08483573
C      0.64602556    3.76864833   -0.15638370
C     -0.75894004    3.56013144   -0.26385678
C     -2.85270987   -1.75016383   -0.39503232
C     -1.96255513   -2.86410654   -0.32177557
C     -0.55935693   -2.65428352   -0.21867748
C      0.32507047   -3.76558283   -0.14638816
C      1.72998263   -3.55474780   -0.04203536
C      2.25012684   -2.23194820   -0.00900270
C      4.17448803   -0.69339168    0.13015361
C      3.65263745   -2.02044973    0.09616244
C     -2.48127917   -4.19840787   -0.34957269
C     -1.58284767   -5.28748837   -0.27752095
C     -0.20324362   -5.10051172   -0.17801839
C      0.72145457   -6.20279486   -0.10340502
C      2.62513930   -4.67277232    0.03074180
C      2.06591779   -5.99921626   -0.00441393
C      4.00239614   -4.44255804    0.13332808
H      4.68172154   -5.29726444    0.18857185
C      4.53737509   -3.14340482    0.16788911
H      0.31756201   -7.21790127   -0.12746896
H     -1.98403085   -6.30407520   -0.29907996
H      2.74959814   -6.84990916    0.05105247
C     -3.89108645   -4.38021619   -0.44663623
C     -4.74364876   -3.29798519   -0.51187603
C     -4.24595074   -1.97625649   -0.49018593
H     -4.28849742   -5.39795813   -0.46577895
H     -5.82509384   -3.42610071   -0.58062887
C     -5.19422716   -0.84589702   -0.54761300
C     -4.61263959    0.51298688   -0.54429105
C      5.58554367   -0.47738125    0.23609332
C      6.44838289   -1.62429286    0.30591639
C      5.94834794   -2.89536521    0.27334233
C      6.07803499    0.83668177    0.26798955
H      7.52551209   -1.45832800    0.38635548
H      6.62269142   -3.75356458    0.32785931
C      5.22604594    1.95025133    0.19910048
H      7.15637553    0.99751641    0.34846635
C     -3.56884196    3.15505802   -0.47932521
C     -3.02478830    4.45954802   -0.44117778
C     -1.65121174    4.68325355   -0.33639132
C     -1.08575516    6.00798087   -0.29732548
C      1.17855490    5.09972574   -0.12046378
C      0.25925680    6.20560612   -0.19481659
C      2.56267186    5.28168159   -0.01588526
H      2.96493156    6.29782768    0.01007544
C      3.45093422    4.19472802    0.05575120
H     -1.76603788    6.86137056   -0.35248417
H     -3.70316666    5.31491160   -0.49553210
H      0.66694606    7.21913685   -0.16754067
C      4.87258797    4.36960286    0.16232626
C      5.71929559    3.29942585    0.23065176
H      5.26890852    5.38765786    0.18740297
H      6.79797132    3.45502690    0.31097142
C     -5.47145497    1.63278235   -0.61391979
C     -4.96955324    2.91723952   -0.58327226
H     -6.54316509    1.44082052   -0.68172787
H     -5.64551332    3.77413852   -0.63440299
O     -6.41525200   -1.02381455   -0.57454138
}

set {
    reference rohf
    soscf_max_iter 50
    soscf true
    scf_type DF
    freeze_core True
    basis jun-cc-pVDZ
}

energy('sapt0')


报错信息：

zjxitcc

你俩的体系咋这么像http://bbs.keinsci.com/thread-35959-1-1.html
你这个羰基改性碳平面确定是二重态的自由基么。。。假设是的话，无需尝试多种方法，使用frag_guess_wfn产生PSI4 SAPT0输入文件（含坐标、基组和轨道），一次可成，例子见MOKIT线上手册https://jeanwsr.gitlab.io/mokit-doc-mdbook/chap5-3.html （SAPT在5.3.5章节）