ORCA IRC验证过渡态出现问题，过渡态能量比两侧低

yaoli2

先用ORCA OptTs计算出过渡态结构，这个结构是有两个未成对电子，并且算的激发态过渡态。然后跑IRC发现IRC很奇怪，过渡态能量比两侧低（下面是irc计算输出的能量），这是哪里出错了。
---------------------------------------------------------------
All gradients are in Eh/Bohr.

Step        E(Eh)      dE(kcal/mol)  max(|G|)   RMS(G)
   1     -1078.220687   -1.060432    0.001721  0.000300
   2     -1078.220527   -0.960172    0.002846  0.000468
   3     -1078.220305   -0.820978    0.003799  0.000761
   4     -1078.220166   -0.733479    0.004148  0.001123
   5     -1078.219529   -0.334079    0.005675  0.001317
   6     -1078.216440    1.604730    0.010431  0.002647
   7     -1078.211279    4.843317    0.017145  0.004226
   8     -1078.203107    9.971158    0.023670  0.005491
   9     -1078.192112   16.870779    0.031310  0.006941
  10     -1078.178488   25.419672    0.037861  0.008406
  11     -1078.162554   35.418554    0.048172  0.011266
  12     -1078.218997    0.000000    0.000664  0.000129 <= TS
  13     -1078.162173   35.657517    0.049159  0.011145
  14     -1078.177994   25.729727    0.039011  0.008379
  15     -1078.191569   17.211160    0.032448  0.006923
  16     -1078.202528   10.334253    0.024859  0.005485
  17     -1078.210658    5.232792    0.018871  0.004264
  18     -1078.215801    2.005726    0.012418  0.002661
  19     -1078.218914    0.051951    0.006472  0.001333
  20     -1078.219584   -0.368114    0.004223  0.001112
  21     -1078.219722   -0.454843    0.003588  0.000759
  22     -1078.219944   -0.594280    0.002349  0.000465
  23     -1078.220100   -0.692045    0.001541  0.000292

yaoli2

好的，谢谢老师。

wzkchem5

应该是CI-NEB优化过渡态收敛有问题
自己取第11步或第13步的结构，重新跑一个OptTS即可