计算化学公社

标题: 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: 5)

所有的输入输出文件均见上面的附件,请问各位老师:
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步数上限也是一样,因为如果遇到了能量很平缓的 ...

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




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