碳簇吸附模型的IRC结果分析求助

zihang

求助老师们一个问题，我的IRC结果如下：（研究体系：构建了碳簇吸附模型，然后在体系中引入CH3自由基进攻有机分子上的碳原子，以催化氧解离过程），目前跑CH3自由基进攻的TS结果是正常的，但是跑IRC的结果比较奇怪，曲线右端出现了两个阶段。


想找的曲线右端结构如下：


CH3自由基进攻的过渡态结果如下：（观察进攻过程没有问题）


请问为什么会出现第一张IRC结果图中右端曲线呈两段的趋势呢？且随着IRC继续跑下去，最终结果已经不是自己想找的了。
附上计算关键词：
寻找过渡态：B3LYP/6-311G** em=GD3BJ opt(ts,calcfc,noeigen) freq
跑IRC：B3LYP/6-311G** em=GD3BJ IRC(calcfc,maxpoints=150,LQA)

求助各位老师，这种情况下该怎么分析呢？也不知道为什么会出现这种情况，谢谢各位老师！

zihang

sobereva 发表于 2024-10-8 03:50
两个stage是很常见的事，出现类似肩峰说明中途有造成能量上升的因素，或者造成能量下降的动力中途已显著减 ...

谢谢sob老师！

sobereva

两个stage是很常见的事，出现类似肩峰说明中途有造成能量上升的因素，或者造成能量下降的动力中途已显著减小，但效果又没强到导致产生新的过渡态。通过研究IRC过程的电子结构改变可以进行讨论（例如结合Multiwfn做波函数分析考察原子电荷、键级变化曲线等，参考《通过键级曲线和ELF/LOL/RDG等值面动画研究化学反应过程》http://sobereva.com/200），可以在文中具体讨论反应的不同的stage的特点