计算化学公社

标题: 求助：Guassian模拟含金属原子的有机分子吸附金属离子的SCF不收敛 [打印本页]

作者
Author: qooqooli 时间: 2020-7-23 11:45
标题: 求助：Guassian模拟含金属原子的有机分子吸附金属离子的SCF不收敛

问题如下：
做的是一个含有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

构建分子模型如下：
(, 下载次数 Times of downloads: 20)

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

作者
Author: qooqooli 时间: 2020-7-23 14:29

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

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

作者
Author: 冰释之川 时间: 2020-7-23 15:44

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

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

作者
Author: qooqooli 时间: 2020-7-23 17:42

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

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

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

作者
Author: qooqooli 时间: 2020-7-23 17:50

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

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

作者
Author: qooqooli 时间: 2020-7-23 17:55

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

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

作者
Author: thanhtam 时间: 2020-7-23 17:59

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

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

作者
Author: qooqooli 时间: 2020-7-23 18:04

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

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

作者
Author: sobereva 时间: 2020-7-24 07:16

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

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

作者
Author: sobereva 时间: 2020-7-24 07:17

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

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

作者
Author: qooqooli 时间: 2020-7-25 15:55

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

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

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