计算化学公社

标题: 请教关于旋轨耦合矩阵元和系间窜越速率的计算 [打印本页]

作者
Author:
ABQTrap    时间: 2024-1-4 10:39
标题: 请教关于旋轨耦合矩阵元和系间窜越速率的计算
各位老师好,本人用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. *
复制代码







欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3