IRC不平滑求助

413

最近算了一个稍复杂的分子。找到了一个过渡态，只有一个虚频，从虚频对应的原子振动来看，是想找的过渡态。然后算IRC。                                                                     LQA, RCFC，forward方向能够回到反应物。
但是reverse方向只跑出了两个点，这两个点前后坐标的变化看，确实往生成物方向走了，只是相应键长变化的幅度较小，只有0.0几埃。（拿出这个构型做opt，奇怪的是优化回了反应物）。于是尝试调整stepsize，从默认的10，调整到了30，得到了不光滑的IRC曲线，如下图。针对现有结果有几点疑问请教大家。
1. 从能量看，是上升的？？？
2. 为什么irc曲线如此波折？
3. 为什么拿最开始stepsize=10时跑出的第二个点直接opt，分子会回到反应物？可是这个irc跑的时候，第二个点确实是在往产物走啊，只是幅度有点小
4. 仅从irc跑出两点，第二点在往产物走（再次说明，幅度有点小），是否就可以说明这个过渡态找的对？还需不需要跑出完美的irc曲线？应该怎么调整才能得到完美的irc曲线？



下图说明，LQA，stepsize=30

sobereva

这个IRC完全错误，基于这个讨论什么都没有意义

跑不成功都是此文里提及的问题没有注意到
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

如果说不清楚就把优化过渡态和IRC任务的输出文件都上传

413

sobereva 发表于 2019-10-23 03:23
这个IRC完全错误，基于这个讨论什么都没有意义

跑不成功都是此文里提及的问题没有注意到

请教下，我调整了stepsize=5得到了上边类似波折的曲线。请问我这个跑的IRC错在了哪里呢？

另：在使用LQA，stepsize=10时得到了较为合理的IRC曲线，但是只跑出了3个点，如下图。第三个点往产物走了，原子位置变化的幅度较小。
问题是我拿第三个点的直接opt又回到了反应物。这可咋整

sobereva

413 发表于 2019-10-23 08:54
请教下，我调整了stepsize=5得到了上边类似波折的曲线。请问我这个跑的IRC错在了哪里呢？

另：在使用L ...

没输出文件没法判断