计算化学公社

标题: 分裂成两个自由基反应的过渡态跑IRC时发现两边都是原料，怎么得不到自由基产物？ [打印本页]

作者
Author: TQcomputation 时间: 2025-1-10 17:44
标题: 分裂成两个自由基反应的过渡态跑IRC时发现两边都是原料，怎么得不到自由基产物？
本帖最后由 TQcomputation 于 2025-1-10 17:48 编辑

如下反应 (, 下载次数 Times of downloads: 5)
通过自旋极化得到过渡态，该过渡态符合反应预期的运动方式，且只有一个虚频-289.17。但是用这个过渡态跑IRC，不论我多重度怎么设置得到的都是正态分布曲线，原料与产物相同，得不到分裂成自由基的产物。请问大佬，是不是我的过渡态找错了？或者是我IRC设置错误？

(, 下载次数 Times of downloads: 5) (, 下载次数 Times of downloads: 4)

#p irc=(calcfc,maxpoint=100,lqa) ub3lyp/6-31g(d) guess=read em=gd3bj

tsIRC

0 1 0 1 0 2 0 -2
N(Fragment=1)    -0.33073400 -0.46528800 -0.02238000
N(Fragment=1)    0.33080900 0.46540900 0.02257400
C(Fragment=2)    2.44865700 -0.04237100 0.00334900
C(Fragment=2)    2.61190100 -0.15800200 -1.48298300
H(Fragment=2)    2.34683600 0.77727700 -1.98919200
H(Fragment=2)    3.65666200 -0.39131400 -1.75006400
H(Fragment=2)    1.98566500 -0.96059600 -1.88985300
C(Fragment=2)    2.44903300 -1.32075200 0.78047700
H(Fragment=2)    3.39606500 -1.86983100 0.65097700
H(Fragment=2)    2.31095800 -1.14465200 1.85271800
H(Fragment=2)    1.64237900 -1.98231300 0.43573100
C(Fragment=2)    2.97558300 1.19672600 0.65819100
H(Fragment=2)    4.07888300 1.20969100 0.66193400
H(Fragment=2)    2.64100600 2.09535600 0.12730200
H(Fragment=2)    2.64478800 1.27068100 1.70000300
C(Fragment=3)    -2.44860800 0.04237600 -0.00336500
C(Fragment=3)    -2.97528600 -1.19622200 -0.65936000
H(Fragment=3)    -2.64094400 -2.09525800 -0.12901400
H(Fragment=3)    -4.07858400 -1.20915700 -0.66356800
H(Fragment=3)    -2.64406600 -1.26939500 -1.70109300
C(Fragment=3)    -2.44892400 1.32137900 -0.77946300
H(Fragment=3)    -3.39597900 1.87032800 -0.64961900
H(Fragment=3)    -1.64231800 1.98268200 -0.43410700
H(Fragment=3)    -2.31071700 1.14615100 -1.85182900
C(Fragment=3)    -2.61237900 0.15678300 1.48300500
H(Fragment=3)    -1.98645700 0.95918200 1.89074200
H(Fragment=3)    -3.65728300 0.38965600 1.74991400
H(Fragment=3)    -2.34728500 -0.77884600 1.98854900

作者
Author: sobereva 时间: 2025-1-10 23:39
对称破缺单重态的IRC计算应当forward和reverse分别跑，两次都读取之前找过渡态得到的chk里的波函数当初猜

作者
Author: TQcomputation 时间: 2025-1-11 14:04
现在能跑出游离基产物了！非常感谢！！！

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