计算化学公社

标题: IRC曲线正常跑完,但是有一个点掉下去了 [打印本页]

作者
Author:
radical    时间: 2021-12-20 17:25
标题: IRC曲线正常跑完,但是有一个点掉下去了
各位大神,学生找到双自由基反应的过渡态后,做irc虽然正常结束了,但图上有一点突然掉下去,输入文件和irc图如下,

(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 2)


作者
Author:
liyuanhe211    时间: 2021-12-20 19:01
负数一侧的第一个点可能与右侧的第一个点电子态不同。单重态双自由基的情况先获得正确的电子结构的chk,再读它作为初猜跑IRC更稳妥。
作者
Author:
radical    时间: 2021-12-20 22:35
liyuanhe211 发表于 2021-12-20 19:01
负数一侧的第一个点可能与右侧的第一个点电子态不同。单重态双自由基的情况先获得正确的电子结构的chk,再 ...

谢谢老师的解答,我马上来试一下
作者
Author:
sobereva    时间: 2021-12-21 18:53
电子结构特殊(如双自由基)而必须从现有过渡态chk里读初猜波函数的IRC任务,必须正、逆方向分别跑,此文明确说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400http://bbs.keinsci.com/thread-7736-1-1.html
另外,出现这种问题还有一种可能是在IRC任务中,基于自动初猜的波函数做SCF后收敛到的波函数与找过渡态任务最终得到的波函数不同,此时相当于IRC任务所在的势能面和过渡态搜索任务所在的势能面不同,这也会导致IRC异常,因为类似于违背了前述的走IRC的“任何影响势能面的设定必须严格相同”的这个前提。出现这种情况时,你会发现IRC任务第一次输出的SCF Done能量和找过渡态最后一步的SCF Done能量明显不同。为解决此问题,走IRC的时候可以用guess=read关键词,从优化过渡态的chk文件中读取最后的波函数(并且最好用forward和reverse关键词通过两个任务分别跑正向和逆向IRC),这样通常可以确保IRC任务所在的势能面和优化过渡态时相同。





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3