计算化学公社

标题: IRC不平滑求助 [打印本页]

作者
Author: 413 时间: 2019-10-23 01:43
标题: IRC不平滑求助
最近算了一个稍复杂的分子。找到了一个过渡态，只有一个虚频，从虚频对应的原子振动来看，是想找的过渡态。然后算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

作者
Author: sobereva 时间: 2019-10-23 03:23
这个IRC完全错误，基于这个讨论什么都没有意义

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

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

作者
Author: 413 时间: 2019-10-23 08:54
本帖最后由 413 于 2019-10-23 09:15 编辑

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

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

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

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

作者
Author: sobereva 时间: 2019-10-24 07:21

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

另：在使用L ...

没输出文件没法判断

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