求助CO在类似石墨结构上脱附IRC过程的产物优化到了反应物

wch990122

各位老师好，本人在分析一个脱附过程时遇到了问题，具体是CO在类似石墨结构上脱附，如图1。
首先根据感觉进行过渡态的判断，但始终得不到正确的过渡态。于是对键长做了柔性扫描，发现存在最大值，如图2，于是将最大值结构设置成初猜的过渡态，得到的过渡态有且至于一个虚频，并且振动方向也是正确的，进一步进行IRC分析，因为默认步长始终是走了几步就停止，按照相关帖子将步长放大到30，得到IRC曲线如图3，反应物和产物都是我想要的，但是对产物进行进一步的优化时，却得到了反应物，不知为何所以特来求救。
扫描，过渡态，IRC的输入文件如下：

%nprocshared=4
%mem=8GB
%chk=scan2.chk
#p opt=modredundant ub3lyp/6-31g(d) nosymm

.......
12 15 s 20 0.1

%nprocshared=4
%mem=8GB
%chk=TS3.chk
#p opt=(TS,calcfc,noeigen) freq ub3lyp/6-31g(d) em=gd3


%nprocshared=4
%mem=8GB
%chk=IRC3.chk
#p irc=(calcfc,maxpoint=150,stepsize=30,lqa) ub3lyp/6-31g(d)


图片：
  

sobereva

IRC加上recalc=3，然后再对末端优化。而且优化时opt加上maxstep=5,notrust

还有这问题的话，改用M06-2X重新算一遍

wch990122