质子迁移的过渡态IRC路径两端的结构异常

lyt0226

各位老师，大家好，

在计算一个CH2CCH2分子次序解离的路径中，其第一步解离成H+C3H3-2+后，C3H3-2+要进一步解离成C+C2H3，所以要先发生质子迁移，在寻找对应的过渡态时虚频的振动方向合理，IRC正常结束，但是IRC跑向两端的结构却类似，如图“product”，IRC中质子迁移到一半后又返回原来连接的碳原子上。请问可能是什么问题。

具体文件见附件，谢谢各位老师！

TS3的输入文件：
#p opt=(cartesian,calcfc,ts,noeigen) m062x/def2TZVPP freq

TS3

2 2
C                  1.36696300    0.00005800    0.00018300
H                  0.81183586   -0.16021930   -1.20563336
C                  0.09097800   -0.00015500    0.00046300
C                 -1.26126900    0.00005300   -0.00016300
H                 -1.82556600   -0.95831300   -0.00048700
H                 -1.82542300    0.95849600   -0.00048600



TS3-IRC的输入文件：
#p irc(LQA,calcfc,maxpoint=50,stepsize=15) m062x/def2tzvpp pop=npa

TS3-IRC

2 2
C                  1.51056200   -0.16804500   -0.05243400
H                  1.13464500    1.05881200    0.37938100
C                  0.16565200    0.03442700   -0.00802800
C                 -1.26272400   -0.02381700   -0.00098500
H                 -1.86111200    0.87537200   -0.24860400
H                 -1.75447300   -0.98957700    0.23790200

sobereva

如果用的g09，把优化和IRC都加上int=ultrafine。IRC里calcfc改成calcall，stepsize设成20

lyt0226

sobereva 发表于 2021-4-29 15:52
如果用的g09，把优化和IRC都加上int=ultrafine。IRC里calcfc改成calcall，stepsize设成20

谢谢sob老师。

但是我更改输入文件以后还是一样的结果，还有另外一个碎裂成两个碎片的过渡态，也有同样的问题，IRC的一端是碎片分开，但是另一端是靠近然后又分开。