有关镧系团簇的结构优化中，SCF无法收敛

比尔摩斯 · 发表于 Post on 2024-1-26 21:12:38

想求助各位老师有关镧系团簇的结构优化中，SCF无法收敛怎么解决。
该结构有214个原子，结构式为：[(nBuSn)8Er(μ2-O)(μ2-OH)5(μ3-O)2(μ4-O)2(NE)2(IDA)2(EG)2]，其中μn表示桥配体与周围n个金属原子配位，Bu为丁烷，NE为5-将冰片烯-2，3-二羧酸，IDA为亚氨基二乙酸，EG为乙二醇。输入文件我已上传。

Sn8Er1-opt4.gjf (12.78 KB, 下载次数 Times of downloads: 9)
写了如下关键词：#p opt freq em=gd3bj int=fine pbe1pbe/genecp scf=(novaracc,noincfock,maxcyc=200)
所用基组分别为：C、O、N、H：6-31g*，Sn：def2svp，Er：MWB58
也曾试过放宽收敛限conver=6，更换为b3lyp基组，使用大核赝势MWB57，小核赝势MWB28（SDD），自定义def2赝势等方法，在使用大核赝势时只能优化3-6步，结构向内收缩甚至有些原子会发生重合，最后也报错SCF不收敛；使用小核赝势时，第一步在200次cycle内SCF都无法收敛，产生报错。希望各位老师可以提供意见，感激不尽！！！

sobereva · 发表于 Post on 2024-1-26 22:40:52

b3lyp叫泛函

参考下文，该说的都说了
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）
解决SCF不收敛问题的方法
http://sobereva.com/61

体系中的明显硬伤必须解决，诸如这里明显非常奇怪

写int=fine完全多余

比尔摩斯 · 发表于 Post on 2024-1-27 13:08:18

sobereva 发表于 2024-1-26 22:40
b3lyp叫泛函

参考下文，该说的都说了

好的，多谢sob老师指出错误！以后会注意的！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 有关镧系团簇的结构优化中，SCF无法收敛