sobereva 发表于 2024-10-10 16:59 好的,我去研究一下 |
把此文认真看了 使用Gaussian做镧系金属配合物的量子化学计算 http://sobereva.com/581(http://bbs.keinsci.com/thread-21274-1-1.html) |
wjc404 发表于 2024-10-10 15:45 好的,谢谢您,我了解一下 |
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收敛问题” 这条路子的前提是体系结构/电荷/自旋多重度合理。如果前提不成立,则应先思考前提条件。 |
我之前算钐的SmI2(THF)n的时候遇到过收敛性的问题,后来用ORCA软件的TRAH可以解决。(不过当时算到单电子还原卤代烃的过渡态时TRAH也失败了,没有再继续了) |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-27 09:10 , Processed in 0.718642 second(s), 26 queries , Gzip On.