计算化学公社

标题: 高斯寻找过渡态几何不收敛（已经尝试过社长帖子的大部分方法）该怎么办？ [打印本页]

作者
Author: Vincentears 时间: 2020-12-5 11:01
标题: 高斯寻找过渡态几何不收敛（已经尝试过社长帖子的大部分方法）该怎么办？
学生在寻找五价铁和小分子污染的过渡态中，遇到了几何不收敛的问题，基于社长的帖子http://sobereva.com/164在之前的过渡态寻找过程中得到了极大的帮助，但是在这计算里面，尝试了许多方法（包括gdiis，maxstep=1都用过了，笛卡尔坐标等....）都无法收敛到我想要的的过渡态（铁氧中的氧被小分子污染物拔走），因此来发第一次贴求助，恳求老师的帮忙。先谢过大家~

作者
Author: 冰释之川 时间: 2020-12-5 11:48
先用opt=loose 试试看吧

作者
Author: liyuanhe211 时间: 2020-12-5 12:36
本帖最后由 liyuanhe211 于 2020-12-5 12:40 编辑

含铁的注意检查一下每一步的电子结构对不对，我遇到过优化过程中它震荡的两个结构其实是在两个不同的电子态之间不断切换。

既然都能接受calcall的耗时了，可以看看每一步的虚频是不是对的，你的虚频数值比较小。

可以做做柔性扫描，不只是键长，还包括从什么方向接近之类的。

作者
Author: wzkchem5 时间: 2020-12-5 12:37
没准这个反应就是没有过渡态的，活化能为0。
可以做个柔性扫描，看看是不是真的没有势垒

作者
Author: 冰释之川 时间: 2020-12-5 13:05

liyuanhe211 发表于 2020-12-5 12:36
含铁的注意检查一下每一步的电子结构对不对，我遇到过优化过程中它震荡的两个结构其实是在两个不同的电子态 ...

我在想是不是应该对波函数稳定性检测一下……

作者
Author: Vincentears 时间: 2020-12-5 13:24

liyuanhe211 发表于 2020-12-5 12:36
含铁的注意检查一下每一步的电子结构对不对，我遇到过优化过程中它震荡的两个结构其实是在两个不同的电子态 ...

学生在气相中能找出来，虚频也合适，一加隐性水就跑不出来了（虽然水的影响肯定会很大，但总是觉得气相中能找出来在水中应该也能找出来）。谢谢老师，我用柔性扫描试试

作者
Author: Vincentears 时间: 2020-12-5 13:24

wzkchem5 发表于 2020-12-5 12:37
没准这个反应就是没有过渡态的，活化能为0。
可以做个柔性扫描，看看是不是真的没有势垒

谢谢您的回复，学生待会试试

作者
Author: Vincentears 时间: 2020-12-5 13:30

冰释之川发表于 2020-12-5 11:48
先用opt=loose 试试看吧

直接上loose吗，学生先试一下，谢谢您~

作者
Author: sobereva 时间: 2020-12-6 14:56

Vincentears 发表于 2020-12-5 13:24
学生在气相中能找出来，虚频也合适，一加隐性水就跑不出来了（虽然水的影响肯定会很大，但总是觉得气相中 ...

对于反应物或产物带显著电荷的情况，极性溶剂环境可能极为显著地影响势能面，甚至决定过渡态能否出现。例如
(, 下载次数 Times of downloads: 57)

作者
Author: wzkchem5 时间: 2020-12-6 17:10

冰释之川发表于 2020-12-5 13:05
我在想是不是应该对波函数稳定性检测一下……

同意，不过波函数稳定性检测也不够，含铁体系经常有两个以上的稳定波函数。。。

作者
Author: Vincentears 时间: 2020-12-10 10:52

liyuanhe211 发表于 2020-12-5 12:36
含铁的注意检查一下每一步的电子结构对不对，我遇到过优化过程中它震荡的两个结构其实是在两个不同的电子态 ...

老师，再请教一下您，如何确定优化有没有收敛到正确的电子态？

作者
Author: wzkchem5 时间: 2020-12-10 12:09

Vincentears 发表于 2020-12-10 10:52
老师，再请教一下您，如何确定优化有没有收敛到正确的电子态？

你做计算之前必然是本来就知道正确的电子态是哪个的，比如铁是高自旋还是低自旋，其他的一些带自旋的分子和铁是铁磁耦合还是反铁磁耦合。所以一般看一下自旋布居是不是和你的预期相符就行了，有时还需要作出前线轨道的图。
如果不知道哪个电子态是正确的，那就穷举出所有可能的电子态，分别找到这些电子态，取能量最低的那个。

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