计算[Fe-CF3CF3COO]2+的配位前后内能的变化scf不收敛

暗黑烧酒小圆

各位老师你们好，我在用BP86/6-31g*优化[Fe-CF3CF2COO]2+配合物结构时，一直SCF不收敛，而且最后一步的能量激增。谢谢老师们的指点呀>.< />.<已尝试方法：改变自旋多重度2，4，6；scf=vshift=300/400/500; 延长0-Fe键长 均无法收敛。
%chk=F:\MAHOSHAOJIU\Gaussian data\luohe\longer C-C\Fe3+-CF3CF2COO(O-Fe longer).chk
#p opt freq 6-31g(d)  bp86


Title Card Required


2 6
C                 -6.19603421    2.07252670    0.01959729
C                 -5.79360245    2.75332280    1.34102279
C                 -4.26987636    2.64017043    1.53352040
F                 -5.56624242    2.67893666   -1.00905681
F                 -7.53176812    2.17171871   -0.14915062
F                 -5.84287317    0.77012294    0.05860681
F                 -6.42339424    2.14691283    2.36967689
F                 -6.14676349    4.05572655    1.30201326
O                 -3.56264554    3.65452751    2.25171556
O                 -3.56221866    1.52074330    0.99407301
Fe                -1.55231543    2.15472358    1.92511847

输出文件
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB-P86) =  -2025.65229963     A.U. after  129 cycles
            NFock=128  Conv=0.74D-05     -V/T= 2.0048
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 2.5000 <S**2>= 8.7567 S= 2.5011
<L.S>= 0.000000000000E+00
KE= 2.016061827326D+03 PE=-6.357393672256D+03 EE= 1.517291876606D+03
Annihilation of the first spin contaminant:
S**2 before annihilation     8.7567,   after     8.7500
Convergence failure -- run terminated.
Error termination via Lnk1e in C:\Gaussian setup\G09W\l502.exe at Fri Apr 16 17:18:48 2021.
Job cpu time:       0 days  0 hours 58 minutes 15.0 seconds.
File lengths (MBytes):  RWF=     38 Int=      0 D2E=      0 Chk=      4 Scr=      1

paramecium86

喵星大佬 发表于 2021-4-18 03:49
很明显的问题，你觉得可以存在配位这么不饱和的玩意么

合理建模

同意~ 我也觉得他这个结构即使让SCF收敛了，几何优化也挺难收敛的。

喵星大佬

很明显的问题，你觉得可以存在配位这么不饱和的玩意么

合理建模

paramecium86

http://sobereva.com/61
解决SCF不收敛问题的方法

sob大这个帖子里提到的一些方法既然已经试过 那就继续试试其它的。不改变泛函和基组的前提下，比如还有二次收敛qc没有试过可以写 scf=xqc  看起来你这个scf收敛到-5次方了，并不是那么难收敛的 加上这个可能很快就能收敛了。

至于多重度看起来应该是6重态问题不大。

njfuzjs

试了一下pbe/6-31g*,收敛是没问题的

njfuzjs

可以试试使用def2svp基组跑跑看一看