Gaussian计算Tb镧系化合物SCF不收敛

DUAN112233 · 发表于 Post on 2024-10-10 15:11:42

本帖最后由 DUAN112233 于 2024-10-10 15:09 编辑

按照社长在解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)中写的做了测试，体系现在还是无法收敛，上传的文件是最近的一次计算参数。麻烦问一下还有什么设置能帮助Tb镧系化合物收敛吗？

wjc404 · 发表于 Post on 2024-10-10 15:45:28

我之前算钐的SmI2(THF)n的时候遇到过收敛性的问题，后来用ORCA软件的TRAH可以解决。（不过当时算到单电子还原卤代烃的过渡态时TRAH也失败了，没有再继续了）

zjxitcc · 发表于 Post on 2024-10-10 16:01:23

本帖最后由 zjxitcc 于 2024-10-10 16:03 编辑

（1）假设该体系基态为单重态，或者你不管基态/激发态你只想研究体系是最低单重态时的情况，那这不仅仅是收敛与否的问题，你这个体系看上去不是闭壳层的；如果先用闭壳层计算，然后依靠stable=opt发现不稳定、自动优化至开壳层单重态，会花费巨量的时间，因为 RKS解与稳定的UKS解可能异常遥远，还很容易收敛到错误的高能SCF解。合理的做法是使用片段组合波函数构建初始猜测，一开始就想好、设置好不同Tb之间的单电子耦合情况；即使不用片段组合波函数，也应该至少考虑构建对称破缺初猜。建议先仔细阅读这几篇再做这类计算
《用Gaussian做UHF计算》
《谈谈Gaussian软件中的guess=mix》
《片段组合波函数实例1.双原子分子》
谈谈片段组合波函数与自旋极化单重态http://sobereva.com/82

（2）假设该体系基态不是单重态，或者你是随便写的自旋多重度、还没想明白基态自旋多重度，那你很大概率做了错误的计算，自然也不需要关心收敛与否。

另外，“这个SCF不收敛” -> “我就去解决SCF收敛问题” 这条路子的前提是体系结构/电荷/自旋多重度合理。如果前提不成立，则应先思考前提条件。

DUAN112233 · 发表于 Post on 2024-10-10 16:33:23

zjxitcc 发表于 2024-10-10 16:01
（1）假设该体系基态为单重态，或者你不管基态/激发态你只想研究体系是最低单重态时的情况，那这不仅仅是 ...

好的，谢谢您的指导，我去学习一下

DUAN112233 · 发表于 Post on 2024-10-10 16:34:17

wjc404 发表于 2024-10-10 15:45
我之前算钐的SmI2(THF)n的时候遇到过收敛性的问题，后来用ORCA软件的TRAH可以解决。（不过当时算到单电子还 ...

好的，谢谢您，我了解一下

sobereva · 发表于 Post on 2024-10-10 16:59:22

把此文认真看了
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

DUAN112233 · 发表于 Post on 2024-10-10 17:16:49

sobereva 发表于 2024-10-10 16:59
把此文认真看了
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.kein ...

好的，我去研究一下

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

[新手求助] Gaussian计算Tb镧系化合物SCF不收敛