IRC 过程中在 TS 附近出现能量尖峰与梯度异常，这种 IRC 是否可信？

HuiYr

各位老师好，目前我利用Gaussian计算一个金属催化的反应，

图一是最近搜索的一个过渡态的结果，有且仅有一个虚频，且是我想要的；图二是根据该结果跑IRC，利用#p IRC(calcfc,lqa,maxpoints=60) b3lyp/genecp nosymm empiricaldispersion=gd3bj关键词得到的结果。左右两边分别是TS前后的两个中间体，这样的IRC结果是我没有见过的，想要请教各位老师，① 这样的结果可信吗？② 得到这样的结果会是什么原因呢？③ 会是由于我没有进行guess=read的原因吗？

附：计算的是Fe(Ⅲ)的二重态，发现S**2=0.8787，偏离0.75>10%。④ 请问会是自旋污染导致这样的结果吗？

HuiYr

pikachuupup 发表于 2025-12-24 18:32
检验一下TS的波函数稳定性，然后用guess=read读取稳定波函数，并用forward和reverse两个方向分开跑IRC

#p stable=opt b3lyp/genecp nosymm guess=check empiricaldispersion=gd3bj geom=check 利用该关键词检查了波函数，发现波函数是稳定的。 （The wavefunction is stable under the perturbations considered.
The wavefunction is already stable. ）那我重新读check，然后分开两边跑一下irc。谢谢你~

pikachuupup

检验一下TS的波函数稳定性，然后用guess=read读取稳定波函数，并用forward和reverse两个方向分开跑IRC