分裂成两个自由基反应的过渡态跑IRC时发现两边都是原料，怎么得不到自由基产物？

TQcomputation

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



#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

sobereva

对称破缺单重态的IRC计算应当forward和reverse分别跑，两次都读取之前找过渡态得到的chk里的波函数当初猜

TQcomputation

现在能跑出游离基产物了！非常感谢！！！