求助：Guassian模拟含金属原子的有机分子吸附金属离子的SCF不收敛

qooqooli

问题如下： 做的是一个含有Cu原子的有机小分子吸附Fe3+的吸附能计算，但一直不收敛，将自旋多重度也改为1,3,5,9都试过了，但是还是没有收敛。（试过不含Cu的有机分子吸附Fe3+，能够收敛，自旋多重度为6），跪求大神们指点一下。。。谢谢！ 输入文件如下：

%chk=C:\Users\LQ\Desktop\Cu-N-OH-Fe-2.chk
#p opt freq ub3lyp/genecp em=gd3bj

Title Card Required

3 7
C                  0.89170961   -2.81718925    0.09704860
C                  1.59774193   -1.55500635    0.06645061
C                  0.86564761   -0.37040321    0.00084535
C                 -0.56444479   -0.44954479   -0.02038004
C                 -0.46924386   -2.85181729    0.06140128
C                  1.40037761    0.96328825   -0.00150652
C                 -1.37973342    0.69969243   -0.01748503
C                 -0.81794604    2.04005040    0.00467730
C                  0.56486797    2.13673919   -0.00161621
C                 -1.71929716    3.13897358    0.03148208
H                 -1.33172743    4.15431894    0.05168446
C                 -3.08475804    2.90781007    0.04691639
C                 -3.56428183    1.58528021    0.03616957
H                  1.45014608   -3.74722534    0.10640544
H                  2.67861876   -1.56419054    0.10887680
H                 -1.03961132   -3.77341206    0.03146459
H                  1.04145448    3.11553042   -0.01042348
H                 -3.80080186    3.72258003    0.05667731
H                 -4.63309082    1.38601041    0.05346759
N                 -1.19537962   -1.66170711   -0.01273507
N                 -2.74390687    0.52370221   -0.02777699
O                  2.74524519    1.13045774   -0.00084452
H                  2.95862465    2.08598709    0.00706687
Cu                -3.10887766   -1.40896163   -0.04316111
Fe                 5.23550176    0.07059312   -0.02643180

Fe Cu 0
lanl2dz
****
C H O N 0
6-31G*
****

Fe Cu 0
lanl2dz

构建分子模型如下：

snljty

SCF不收敛先别opt freq，先把单点算出来再说。
LanL2DZ太糙，换LanL2TZ。
具体的解决自洽场收敛的方法看卢老师的博文
《解决SCF不收敛问题的方法》
http://sobereva.com/61

qooqooli

snljty 发表于 2020-7-23 13:51
SCF不收敛先别opt freq，先把单点算出来再说。
LanL2DZ太糙，换LanL2TZ。
具体的解决自洽场收敛的方法看 ...

谢谢您，好的我再试试，请问下这个体系自旋多重度应该设定多少才好？7是正确的吗？

冰释之川

qooqooli 发表于 2020-7-23 14:29
谢谢您，好的我再试试，请问下这个体系自旋多重度应该设定多少才好？7是正确的吗？

哪个自旋多重度下最稳定，需要一个个试试

qooqooli

冰释之川 发表于 2020-7-23 15:44
哪个自旋多重度下最稳定，需要一个个试试

关键我的结构优化还没收敛，那我优化的时候应该先选择那个自旋多重度呢？还是只要先解决SCF不收敛，这个自旋多重度的设定不影响收敛？刚刚试了http://sobereva.com/61里面的能级移动法SCF=vshift=500，但还是没有收敛，现在在使用SCF=XQC,如果这个方法也不行，请问还有什么方法解决？

thanhtam

你这个铁离子是准备吸附到哪个位置啊，现在放的位置有点远

qooqooli

snljty 发表于 2020-7-23 13:51
SCF不收敛先别opt freq，先把单点算出来再说。
LanL2DZ太糙，换LanL2TZ。
具体的解决自洽场收敛的方法看 ...

上一个帖子发的是有机分子吸附Fe3+的结构优化的时候，用过(1)SCF=vshift, (2)SCF=conver放宽收敛，(3)SCF=XQC都试过了，还是没有收敛。

qooqooli

thanhtam 发表于 2020-7-23 17:50
你这个铁离子是准备吸附到哪个位置啊，现在放的位置有点远

吸附到O原子的附近，现在距离是2.5A左右。PS:不含Cu原子的有机小分子吸附Fe3+离子能够收敛，Fe3+——O距离是2.3A左右。

thanhtam

qooqooli 发表于 2020-7-23 17:50
上一个帖子发的是有机分子吸附Fe3+的结构优化的时候，用过(1)SCF=vshift, (2)SCF=conver放宽收敛，(3)SCF ...

并不是只能用一个，你可以多管齐下

qooqooli

thanhtam 发表于 2020-7-23 17:59
并不是只能用一个，你可以多管齐下

那SCF收敛请问下这个跟自旋多重度关系大吗？因为现在这个体系的自旋多重度用的是7，Fe3+是5个单电子，Cu原子是一个单电子，总体系自旋多重度是7吗？

sobereva

qooqooli 发表于 2020-7-23 17:42
关键我的结构优化还没收敛，那我优化的时候应该先选择那个自旋多重度呢？还是只要先解决SCF不收敛，这个 ...

我总是纳闷，明明博文里写了那么多方法，最顶端还用红字强调，怎么总有人就用其中几个发现没收敛就无视其它的？

sobereva

qooqooli 发表于 2020-7-23 18:04
那SCF收敛请问下这个跟自旋多重度关系大吗？因为现在这个体系的自旋多重度用的是7，Fe3+是5个单电子，Cu ...

我博文里已经明确强调可能有很大关系

qooqooli

sobereva 发表于 2020-7-24 07:16
我总是纳闷，明明博文里写了那么多方法，最顶端还用红字强调，怎么总有人就用其中几个发现没收敛就无视其 ...

不好意思，sob老师，我自己再仔细阅读。谢谢您！