计算化学公社

标题: IRC的生成物经过优化后生成完全不同的结构 [打印本页]

作者
Author: Vicky_Sugar 时间: 2025-3-12 16:08
标题: IRC的生成物经过优化后生成完全不同的结构
本帖最后由 Vicky_Sugar 于 2025-3-12 16:08 编辑

附件中的TS为搜索到的过渡态，有且仅有一个虚频，且振动方向符合预期。
对该过渡态使用相同的泛函和基组进行IRC（IRC关键词为irc=(calcall,LQA,maxpoints=200,stepsize=10)）。
此后取产物方向的末端结构进行优化，也使用的是相同的泛函和基组，且采用了小步长（优化的关键词为opt=(calcall,maxstep=5) freq和opt=(calcall,maxstep=1) freq），但是发现优化过程中结构发生了很大的变化。
(, 下载次数 Times of downloads: 8)

所有的输入输出文件均见上面的附件，请问各位老师：
1）这个基元反应的生成物应该是哪一种？
2）为什么IRC的末端不能跑出结构优化后的产物？
感谢各位老师的时间！

作者
Author: sobereva 时间: 2025-3-13 07:49
IRC任务结束的位置未必和极小点足够接近，哪怕你用很大的IRC步数上限也是一样，因为如果遇到了能量很平缓的地方往往会自动结束任务。因此如果在此之后实际的IRC曲线还有很长一段距离，做几何优化自然还会导致结构变化很多。实际极小点是什么就是什么（最好对比IRC最后一步和几何优化第一步的电子能量，确保一致，否则说明收敛到了不同的波函数，此时没可比性）。

作者
Author: Vicky_Sugar 时间: 2025-3-13 10:40

sobereva 发表于 2025-3-13 07:49
IRC任务结束的位置未必和极小点足够接近，哪怕你用很大的IRC步数上限也是一样，因为如果遇到了能量很平缓的 ...

好的我明白啦，感谢卢老师详细耐心的解答！

作者
Author: leduo 时间: 2025-12-29 14:46

sobereva 发表于 2025-3-13 07:49
IRC任务结束的位置未必和极小点足够接近，哪怕你用很大的IRC步数上限也是一样，因为如果遇到了能量很平缓的 ...

老师您好，我遇到一个情况是我找到的过渡态虚频为-30多，振动方向是正确的，但是就是IRC计算尝试了您IRC常见问题帖子中的各种方法但是就是无法走到中间体结构，拿IRC最后一步去优化也无法优化成中间体。是否也是因为这个过渡态附近能量平缓所以任务自动结束了呢？我对这个结构也进行了SCAN，如下所示 (, 下载次数 Times of downloads: 3)

作者
Author: sobereva 时间: 2025-12-30 14:53

leduo 发表于 2025-12-29 14:46
老师您好，我遇到一个情况是我找到的过渡态虚频为-30多，振动方向是正确的，但是就是IRC计算尝试了您IRC ...

IRC末端结构总能优化出一个结构来，如果这个结构不是你期望的结构（且优化用了较小的步长上限），说明和你期望的结构对应的极小点之间还隔着一个或多个过渡态

作者
Author: leduo 时间: 2025-12-31 14:19

sobereva 发表于 2025-12-30 14:53
IRC末端结构总能优化出一个结构来，如果这个结构不是你期望的结构（且优化用了较小的步长上限），说明和 ...

好的谢谢老师，我将stepsize调整到40之后可以拿到加成的中间体结构

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