电子结构特殊（如双自由基）而必须从现有过渡态chk里读初猜波函数的IRC任务，必须正、逆方向分别跑，此文明确说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://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任务所在的势能面和优化过渡态时相同。

liyuanhe211 发表于 2021-12-20 19:01
负数一侧的第一个点可能与右侧的第一个点电子态不同。单重态双自由基的情况先获得正确的电子结构的chk，再 ...

谢谢老师的解答，我马上来试一下

负数一侧的第一个点可能与右侧的第一个点电子态不同。单重态双自由基的情况先获得正确的电子结构的chk，再读它作为初猜跑IRC更稳妥。