PVDF和二价钴离子结构优化SCF不收敛/无法正确优化

agone997

gaussian计算有机分子和二价钴离子相互作用时SCF无法收敛，结构偏差十分巨大。我尝试了PBEO/PBE38/B3LYP泛函、尝试了def2svp/6-31g*/lanl2dz赝势基组。均无法获得正常的结果。
十分苦恼到底是什么原因导致此类问题的，为什么？请求老师给予方向指导，感激不尽！

关键词如下：%chk=1.chk%mem=80GB
%nprocshared=28
#p opt freq b3lyp genecp em=gd3bj iop(5/37=10) scf=(maxcycle=300,fermi,novaracc,noincfock,vshift=300)



Title Card Required

2 6
C                -13.00300000  -13.17900000   -3.46500000
…………………………

F                -10.51100000   -9.49900000   -7.65000000

C H F 0
6-31G*
****
Co     0
S    1   1.00
      7.1760000              1.0000000
S    1   1.00
      2.0090000              1.0000000
S    1   1.00
      0.8055000              1.0000000
S    1   1.00
      0.1070000              1.0000000
S    1   1.00
      0.0375000              1.0000000
P    1   1.00
     21.3900000              1.0000000
P    1   1.00
      2.6500000              1.0000000
P    1   1.00
      0.8619000              1.0000000
P    1   1.00
      0.0800000              1.0000000
P    1   1.00
      0.0230000              1.0000000
D    1   1.00
     39.2500000              1.0000000
D    1   1.00
     10.7800000              1.0000000
D    1   1.00
      3.4960000              1.0000000
D    1   1.00
      1.0660000              1.0000000
D    1   1.00
      0.2606000              1.0000000
F    1   1.00
      2.7800000              1.0000000
****

CO     0
CO-ECP     2     10
d potential
  3
1    434.1076104            -10.0000000
2     78.8402758            -66.7240276
2     19.6037749            -11.7192243
s-d potential
  4
0    151.2307582              3.0000000
1    169.6790026             20.1553110
2     64.9909423            413.3762244
2     12.5076287            133.8727705
p-d potential
  4
0     92.7319924              5.0000000
1    121.1910183              6.1892353
2     48.6265909            326.3215919
2     11.9966971             73.8788970

sobereva

Co摆得那么远完全莫名其妙
没事别写IOp。本来都写了noincfock还写5/37=10完全冲突
那叫PBE0不叫PBEO

该说的下文都说了
解决SCF不收敛问题的方法
http://sobereva.com/61

并且输入文件直接上传，别就贴一部分，别人都不方便测试

agone997

sobereva 发表于 2025-3-5 06:10
Co摆得那么远完全莫名其妙
没事别写IOp。本来都写了noincfock还写5/37=10完全冲突
那叫PBE0 ...

谢sob老师指导！摆的远是我的错觉瞎尝试，IoP和noincfock冲突谢谢老师指导，PBE0PBE打字着急了，以后会严谨！提问次数太少，文件上传问题会注意！再次大谢sob老师！

Uus/pMeC6H4-/キ

确认一下，你要研究的二价钴离子是裸露着和小分子在真空气相下结合的吗？如果是溶液环境，钴离子周围应当有很多溶剂分子配位；如果是晶体表面，更是不应该取出单独的钴离子。而且单个聚合度为4的小分子与真实的PVDF长链有些区别。初始结构可以结合无机配位化学对双齿配体螯合的描述构建，自旋多重度按配体场理论确定，没有相应基础知识乱尝试只是浪费时间。

agone997

Uus/pMeC6H4-/キ 发表于 2025-3-5 10:39
确认一下，你要研究的二价钴离子是裸露着和小分子在真空气相下结合的吗？如果是溶液环境，钴离子周围应当有 ...

您的问题启发了我，我再学习学习相关知识！