wzkchem5 发表于 2023-5-4 15:56
最严格的方法是跑非绝热动力学；次严格的方法是用orca的VG、AHAS等不需要激发态平衡结构的方法做近似计算 ...

十分感谢老师的详细回复！在用您说的"较不严格的方法了"找一下，先减小步长做结构优化找找看，再次感谢

万卷书万里路 发表于 2023-5-3 13:31
嗯，谢谢老师！
如果以绝热计算各个三线态能级与S1之间的系间窜越可能性，那计算T3->S1的窜越速率应该怎 ...

最严格的方法是跑非绝热动力学；次严格的方法是用orca的VG、AHAS等不需要激发态平衡结构的方法做近似计算；再次严格的方法是，证明从T3势能面的Franck-Condon点可以走downhill IRC到T2/T3的锥形交叉点，以此证明T3到T2的内转换极其快，因此T3到S1的ISC量子产率可以忽略。
严格证明整个T3势能面上都没有平衡结构是极其困难的，用非绝热动力学可以在一定程度上说明，但也不是完全严格的证据。但是可以证明如果从Franck-Condon点一直走下坡会经过T2/T3锥形交叉点发生内转换到T2，这个可以间接证明T3平衡结构即使存在也很可能比较远，大部分从T3 Franck-Condon点出发的分子走不到这个平衡结构。证明这一点又分两种办法。较严格但耗时较长的方法是我刚才说的走downhill IRC，较不严格的是用较小的步长做结构优化，假如卡在了T2、T3激发能极其接近的结构附近而不能收敛（假如跟踪的是绝热态），又或者激发态成分在若干步后自动变成了T2的成分（假如跟踪的是透热态），就说明T3会直接经过T2/T3锥形交叉点发生内转换

wzkchem5 发表于 2023-5-2 19:02
并不是每个绝热态都存在一个唯一对应的透热态。所以优化T3透热态总是优化到T2是正常的，很可能代表T3透热 ...

嗯，谢谢老师！
如果以绝热计算各个三线态能级与S1之间的系间窜越可能性，那计算T3->S1的窜越速率应该怎么求得？现在找不到T3的结构，还是说用优化结果的T2透热态结构代替（这样能量求到的实际还是T2-S1吧），或者放弃绝热，用垂直激发的？
想知道遇到类似情况该怎么判断和处理。
还有，请问老师，有什么办法能证明T3没有透热态平衡结构或者找到T3的平衡结构吗，比如说要做些什么计算？

万卷书万里路 发表于 2023-5-2 11:07
谢谢老师！
仔细阅读了博文，现在明确了是T3在优化过程中滑到了T2的势能面上， 应该是文中提到的类似的 ...

并不是每个绝热态都存在一个唯一对应的透热态。所以优化T3透热态总是优化到T2是正常的，很可能代表T3透热态平衡结构本来就不存在。

万卷书万里路 发表于 2023-4-28 11:51
谢谢老师！
按照高斯的做法，我把相关态全部优化了，但是发现T3优化完（opt td=(triplets,root=3) 自旋 ...

仔细阅读http://sobereva.com/468

wzkchem5 发表于 2023-4-28 16:20
这个是溶剂模型的原因，属于老生常谈问题了。参见http://sobereva.com/314
计算跃迁速率的话，按高斯的做 ...

谢谢老师！
按照高斯的做法，我把相关态全部优化了，但是发现T3优化完（opt td=(triplets,root=3) 自旋设置1），尽管在输出文件里显示顺序是第三个激发态，但是它的能量（激发能和Total Energy这俩都）却变成了 第二高的，明显不是最高的了。
结果长这样：
Excitation energies and oscillator strengths:

Excited State   1:      Triplet-A      2.2990 eV  539.29 nm  f=0.0000  <S**2>=2.000
      80 -> 86        -0.18432
      84 -> 86        -0.10213
      85 -> 86         0.64924

Excited State   2:      Triplet-A      3.5356 eV  350.67 nm  f=0.0000  <S**2>=2.000
      78 -> 86         0.20288
      80 -> 86        -0.26589
      81 -> 86         0.33314
      81 -> 87        -0.19574
      82 -> 86         0.18394
      82 -> 87        -0.11450
      84 -> 86        -0.15179
      84 -> 87         0.15442
      85 -> 86        -0.13673
      85 -> 87        -0.18628
      85 -> 91        -0.12409

Excited State   3:      Triplet-A      2.4727 eV  501.41 nm  f=0.0000  <S**2>=2.000
      79 -> 86         0.16350
      81 -> 86        -0.27124
      81 -> 87         0.10859
      81 -> 91         0.11643
      84 -> 86        -0.49190
      84 -> 87         0.24074
      84 -> 91         0.15477
      85 -> 86        -0.10884
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -1053.46958359  

这是说 结构跑到了T2那里么，这个能量明显等于我单独优化T2的（opt td=(triplets, root=2)），还是其他原因？
怎样才能得到更准确的T3？比如我把初始结构重新打乱，再重新优化可以解决吗？  
谢谢老师！

这个是溶剂模型的原因，属于老生常谈问题了。参见http://sobereva.com/314
计算跃迁速率的话，按高斯的做法，必须有优化过的激发态结构；如果用orca，只有基态结构也可以（用VG、AHAS等方法计算）