计算化学公社

标题: ORCA IRC验证过渡态出现问题，过渡态能量比两侧低 [打印本页]

作者
Author: yaoli2 时间: 2023-4-18 14:50
标题: ORCA IRC验证过渡态出现问题，过渡态能量比两侧低
先用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

作者
Author: wzkchem5 时间: 2023-4-18 15:35
应该是CI-NEB优化过渡态收敛有问题
自己取第11步或第13步的结构，重新跑一个OptTS即可

作者
Author: yaoli2 时间: 2023-4-18 16:22
好的，谢谢老师。

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