计算化学公社

标题: PVDF和二价钴离子结构优化SCF不收敛/无法正确优化 [打印本页]

作者
Author: agone997 时间: 2025-3-5 00:11
标题: PVDF和二价钴离子结构优化SCF不收敛/无法正确优化
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

作者
Author: sobereva 时间: 2025-3-5 06:10
Co摆得那么远完全莫名其妙
没事别写IOp。本来都写了noincfock还写5/37=10完全冲突
那叫PBE0不叫PBEO

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

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

作者
Author: agone997 时间: 2025-3-5 09:08

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

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

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

作者
Author: agone997 时间: 2025-3-5 23:11

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

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

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