求助open shell singlet 与triplet（T1）旋轨耦合（SOC）常数计算

13799725819

最近，我使用ORCA程序做MECP处的旋轨耦合（SOC）常数计算，发现只能做close shell singlet（S0）与triplet（T1）的SOC常数，请问：open shell singlet 与triplet（T1）存在SOC常数吗？以及用什么程序可以做这种计算呢？

ORCA计算输入：
%pal nprocs 12 end
%maxcore 2000
! B3LYP 6-311++G(d,p) TightSCF
%TDDFT  NROOTS  5
         DOSOC   TRUE
  END
* xyz 0 1
mecp的坐标
*

计算输出：
      --------------------------------------------------------------------------------
                      CALCULATED SOCME BETWEEN TRIPLETS AND SINGLETS                  
      --------------------------------------------------------------------------------
           Root                          <T|HSO|S>  (Re, Im) cm-1                     
         T      S             Z                     X                     Y           
      --------------------------------------------------------------------------------
         1      0    (   0.00 ,    0.47 )    (   0.00 ,   32.80 )    (  -0.00 ,   15.24 )
         1      1    (   0.00 ,   -0.04 )    (   0.00 ,   -0.00 )    (  -0.00 ,    0.00 )
         1      2    (   0.00 ,   -0.02 )    (   0.00 ,    0.00 )    (  -0.00 ,    0.00 )
         1      3    (   0.00 ,    0.01 )    (   0.00 ,   -6.02 )    (  -0.00 ,   -0.16 )
         1      4    (   0.00 ,    0.02 )    (   0.00 ,   -0.00 )    (  -0.00 ,   -0.00 )
         1      5    (   0.00 ,    0.38 )    (   0.00 ,   -0.31 )    (  -0.00 ,   10.16 )
。。。

13799725819

zjxitcc 发表于 2024-2-21 15:21
虽然能看明白你的意思，不过描述是不合理的，是S0与T1的SOC常数，不是“S0的SOC常数”和“T1的SOC常数”， ...

谢谢，您表达的更清晰一些，已更正！
明白了，我去试一试。

13799725819

paramecium86 发表于 2024-2-21 15:04
ORCA本身就可以做。做SOC并不是非得局限于TD方法。比如基于QD-NEVPT2来做SOC就是可以的。

嗯嗯 谢谢 我研究一下

zjxitcc

虽然能看明白你的意思，不过描述是不合理的，是S0与T1的SOC常数，不是“S0的SOC常数”和“T1的SOC常数”，它是一个值，不是两个值，是同时指定两个电子态才对应一个值。类似地，没有“open shell singlet的SOC常数”。

你要的东西目前只能用多组态方法（如CASSCF）或多参考方法（如CASPT2,NEVPT2）算，ORCA就支持。需要先构造合理的活性空间，完成CASSCF计算，才能继续CASPT2或NEVPT2计算。

paramecium86

ORCA本身就可以做。做SOC并不是非得局限于TD方法。比如基于QD-NEVPT2来做SOC就是可以的。