请教关于旋轨耦合矩阵元和系间窜越速率的计算

ABQTrap

各位老师好，本人用CASSCF进行了基态单重态和基态三重态的最低能量交叉点（MECP）的结构优化，现在想看看单重态-三重态之间的系间窜越速率（ISC rate）。目前看到论坛里讨论较多的是用Gaussian或ORCA的TDDFT计算旋轨耦合矩阵元，Sob老师也写了详细的教程，但用CASSCF计算SOC的较少。作为新手，想请问几个问题：
（1）本人算的是基态单重态和基态三重态的MECP，旋轨耦合矩阵元如果用TDDFT算合理吗？因为不涉及激发态。还是最好在CASSCF级别下计算（ORCA可以算）？
（2）Sob老师在《使用Gaussian+PySOC在TDDFT下计算旋轨耦合矩阵元》提到“计算旋轨耦合矩阵元还被用于考察系间窜越速率（正比于相应两个态之间旋轨耦合矩阵元的模方）”，请问旋轨耦合矩阵元的模方就是实部和虚部的平方和吗？另外如何用旋轨耦合矩阵元计算系间窜越速率呢？多大的系间窜越速率可以认为系间窜越容易发生呢？
（3）用ORCA的CASSCF计算SOC（输入文件见下），得到的旋轨耦合矩阵元格式与TDDFT不同（见图片），新手不太理解表格中字符的含义，如Block、S、MS的含义？谢谢！

1. ! TZVP miniprint tightSCF pal3 nopop
   
2. %casscf
   
3. nel 6
   
4. norb 4
   
5. mult 1,3
   
6. nroots 3,3
   
7. rel
   
8. dosoc true
   
9. printlevel 3
   
10. end
    
11. end
    
12. * xyz 0 1
    
13. C                  0.00000000    0.00000000   -0.52513500
    
14. H                  0.00000000    0.93987900   -1.11261300
    
15. H                  0.00000000   -0.93987900   -1.11261300
    
16. O                  0.00000000    0.00000000    0.67200400
    
17. *

复制代码