计算化学公社

标题: 质子迁移的过渡态IRC路径两端的结构异常 [打印本页]

作者
Author: lyt0226 时间: 2021-4-29 15:32
标题: 质子迁移的过渡态IRC路径两端的结构异常
本帖最后由 lyt0226 于 2021-4-29 15:31 编辑

各位老师，大家好，

在计算一个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

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

作者
Author: lyt0226 时间: 2021-4-29 20:01

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

谢谢sob老师。

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)