计算化学公社

标题: Fe配合物计算SCF不收敛 [打印本页]

作者
Author: woshiblake 时间: 2023-6-19 15:12
标题: Fe配合物计算SCF不收敛
各位高手，
   我在做Fe配合物opt任务时碰到一个怪问题：用含碳数少的烃基的配体时优化能成功，配体加一个碳就死活SCF不收敛，按sob老师的方法基本都试过，不成功。输入文件都为：

%mem=80GB
%cpu=0-52
%gpucpu=0=41
#p opt ub3lyp/genecp

+2jia

1 2

坐标

C P O H 0
6-311g
****
Fe 0
LanL2DZ
****

Fe 0
LanL2DZ

不收敛的输出文件的描述：
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -1495.20438922    A.U. after  129 cycles
         NFock=128  Conv=0.69D-06    -V/T= 2.0397
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.9412 S= 0.5914
<L.S>=  0.00000000000
KE= 1.438150043554D+03 PE=-6.738123840125D+03 EE= 2.186539215631D+03
Annihilation of the first spin contaminant:
S**2 before annihilation    0.9412, after    0.7578
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/omnisky/g16/l502.exe at Sun Jun 18 21:00:11 2023.
Job cpu time:    0 days  1 hours 14 minutes  7.8 seconds.
Elapsed time:    0 days  0 hours  1 minutes 24.2 seconds.
File lengths (MBytes):  RWF=    39 Int=    0 D2E=    0 Chk=    2 Scr=    1

作者
Author: sobereva 时间: 2023-6-19 20:34
不带极化函数，胡算瞎算
跟碳数没直接关系
该说的都说了
解决SCF不收敛问题的方法
http://sobereva.com/61

作者
Author: PoorChaos 时间: 2023-6-21 11:55
6-311g先改为6-31g*再说

作者
Author: wzkchem5 时间: 2023-6-21 13:43
此外自旋多重度设错了，建议仔细学一下配体场理论，学会怎么正确判断这种配合物的自旋多重度

作者
Author: woshiblake 时间: 2023-6-28 11:18
选这个基组不是我的最终目的。我的配体是含多个碳原子的，最初用带极化加弥散的稍大基组优化时不收敛，所以我从少碳原子开始逐步往上加，等加到我的目的碳数了，再用带极化加弥散的稍大基组优化。然而，4个碳时能收敛，再加2个碳就出现我帖子里描述的问题。高手们的建议我会采纳，好好学习。

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