计算化学公社

标题: Gaussian计算Tb镧系化合物SCF不收敛 [打印本页]

作者
Author: DUAN112233 时间: 2024-10-10 15:11
标题: Gaussian计算Tb镧系化合物SCF不收敛
本帖最后由 DUAN112233 于 2024-10-10 15:09 编辑

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

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

作者
Author: zjxitcc 时间: 2024-10-10 16:01
本帖最后由 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收敛问题” 这条路子的前提是体系结构/电荷/自旋多重度合理。如果前提不成立，则应先思考前提条件。

作者
Author: DUAN112233 时间: 2024-10-10 16:33

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

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

作者
Author: DUAN112233 时间: 2024-10-10 16:34

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

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

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

作者
Author: DUAN112233 时间: 2024-10-10 17:16

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

好的，我去研究一下

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