计算化学公社

标题: 有关镧系团簇的结构优化中，SCF无法收敛 [打印本页]

作者
Author: 比尔摩斯 时间: 2024-1-26 21:12
标题: 有关镧系团簇的结构优化中，SCF无法收敛
想求助各位老师有关镧系团簇的结构优化中，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为乙二醇。输入文件我已上传。 (, 下载次数 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都无法收敛，产生报错。希望各位老师可以提供意见，感激不尽！！！

作者
Author: sobereva 时间: 2024-1-26 22:40
b3lyp叫泛函

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

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

(, 下载次数 Times of downloads: 9)

写int=fine完全多余

作者
Author: 比尔摩斯 时间: 2024-1-27 13:08

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

参考下文，该说的都说了

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

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