CASSCF如何输出旋轨耦合矩阵元

BinWang

各位老师，大家好！
       我用CASSCF计算一个双自由基的旋轨耦合矩阵元，发现输出里面没有给出相应的旋轨耦合矩阵元，这个该怎么解决。我只从输出里面读出如下内容：
*************************************
Doing QDPT with ONLY SOC!
*************************************

Lowest eigenvalue of the SOC matrix:  -20674.88409994 Eh
Energy stabilization:    -0.00109 cm-1
Eigenvalues:     cm-1         eV      Boltzmann populations at T =  300.000 K
    0:            0.00        0.0000       3.01e-01
    1:            0.00        0.0000       3.01e-01
    2:            0.00        0.0000       3.01e-01
    3:          238.82        0.0296       9.59e-02
    4:        21118.90        2.6184       3.10e-45
    5:        21139.48        2.6210       2.81e-45
    6:        21139.48        2.6210       2.81e-45
    7:        21139.48        2.6210       2.81e-45
    8:        24349.28        3.0189       5.80e-52
    9:        24349.28        3.0189       5.80e-52
   10:        24349.28        3.0189       5.80e-52
   11:        24434.71        3.0295       3.85e-52
   12:        34204.31        4.2408       1.73e-72
   13:        34204.31        4.2408       1.73e-72
   14:        34204.31        4.2408       1.73e-72
   15:        38934.56        4.8273       2.42e-82

The threshold for printing is 0.0100
Eigenvectors:
                          Weight      Real          Image    : Block Root    Spin   Ms
  STATE   0:       0.0000
                          0.118033    -0.094585    -0.330282 :     0    0    1    1
                          0.764518    -0.000653    -0.874367 :     0    0    1    0
                          0.117449    -0.095855     0.329030 :     0    0    1   -1
  STATE   1:       0.0001
                          0.470895     0.521926    -0.445519 :     0    0    1    1
                          0.049917     0.002136     0.223411 :     0    0    1    0
                          0.479188     0.526176     0.449808 :     0    0    1   -1
  STATE   2:       0.0001
                          0.411073    -0.557405    -0.316816 :     0    0    1    1
                          0.185565     0.104046     0.418018 :     0    0    1    0
                          0.403362    -0.347635     0.531519 :     0    0    1   -1
  STATE   3:     238.8186
                          0.999999    -1.000000    -0.000000 :     1    0    0    0
  STATE   4:   21118.8989
                          0.999998    -0.999999    -0.000016 :     1    1    0    0
  STATE   5:   21139.4772
                          0.166596     0.042725    -0.405920 :     0    1    1    1
                          0.665972    -0.000090    -0.816071 :     0    1    1    0
                          0.167432     0.045925     0.406599 :     0    1    1   -1
  STATE   6:   21139.4773
                          0.495365    -0.679558    -0.183210 :     0    1    1    1
                          0.011753    -0.000404     0.108412 :     0    1    1    0
                          0.492882    -0.677878     0.182658 :     0    1    1   -1
  STATE   7:   21139.4773
                          0.338038     0.191637    -0.548920 :     0    1    1    1
                          0.322275    -0.000258     0.567692 :     0    1    1    0
                          0.339686     0.195028     0.549226 :     0    1    1   -1
  STATE   8:   24349.2824
                          0.170956     0.002224     0.413462 :     0    2    1    1
                          0.658406    -0.000147     0.811422 :     0    2    1    0
                          0.170623     0.002962    -0.413055 :     0    2    1   -1
  STATE   9:   24349.2836
                          0.369732    -0.348823    -0.498051 :     0    2    1    1
                          0.272774     0.029596     0.521439 :     0    2    1    0
                          0.357494    -0.291900     0.521812 :     0    2    1   -1




我的输入文件是：

! X2C X2C-SVPall NoIter TightSCF  pal32   # 第1行：计算指令
%maxcore 10000                         # 第5行：单行指令（无end！）

!MORead
%MOInp "SOC_1.gbw"                 # 第7行

%casscf
   nel 8
   norb 8
   mult 3,1
   nroots 4,4
  # weights 0.5, 0.5
#  weights 1,1,1,1, 1,1,1,1
   maxiter 800
   rel
   dosoc true
   end
   PrintLevel 4
end

# %cpcm smd true
#       SMDsolvent "water"
# end
%basis                                 # 第12行：基组设置
    NewGTO Ir "X2C-TZVPall" end          # 第13行
end                                    # 第14行
*xyzfile 1 1 s.xyz                   # 第15行：分子坐标文件


另外，通过输出文件给出的信息来推断T0-S1（S1和T0分别为单重态与三重态基态，T0比S1仅仅稳定1.2kcal/mol，两者的最稳定构象极其相似，总的差距不过0.38 anstrong）的旋轨耦合常数几乎为0，这个怎么解决？我们通过实验推测S0道T0的转化速率在皮秒量级。

wal

rel里要加printlevel
哎不对啊 这不是打印了嘛

BinWang

wal 发表于 2026-1-30 11:15
rel里要加printlevel
哎不对啊 这不是打印了嘛

没打印吧？没有类似这种形式的输出：Spin-Orbit Matrix Elements (cm^-1)
Between states of different multiplicity:
S0 (0.000 cm-1)  <->  T1 (xxxx cm-1)
   Cartesian components:
     X:   xx.xx cm-1
     Y:   xx.xx cm-1
     Z:   xx.xx cm-1
   Total SOC:   xx.xx cm-1

wal

BinWang 发表于 2026-1-30 11:30
没打印吧？没有类似这种形式的输出：Spin-Orbit Matrix Elements (cm^-1)
Between states of different  ...

orca能打印这种格式的输出吗？我好像没见到，你可以在rel里加大printlevel试试，我之前看手册只给到2了，继续加大可能给出

BinWang

wal 发表于 2026-1-30 11:45
orca能打印这种格式的输出吗？我好像没见到，你可以在rel里加大printlevel试试，我之前看手册只给到2了， ...

我已经设置到PrintLevel 4了，没有直接给出旋轨耦合常数，只给出了旋轨耦合常数相关的一些光谱数据，DFT方法里面有这种格式的旋轨耦合常数。

wal

BinWang 发表于 2026-1-30 11:46
我已经设置到PrintLevel 4了，没有直接给出旋轨耦合常数，只给出了旋轨耦合常数相关的一些光谱数据，DFT ...

我看你贴出来的输入文件写的printlevel并不在rel里。

BinWang

wal 发表于 2026-1-30 12:51
我看你贴出来的输入文件写的printlevel并不在rel里。

，非常感谢，我重算一遍。

风起~

你这轨道指定地太随意了吧
这样的结果有意义么

BinWang

风起~ 发表于 2026-2-1 17:04
你这轨道指定地太随意了吧
这样的结果有意义么

非常感谢老师的回复，我这个读入的是已经收敛好的基态CASSCF波函数，轨道没什么问题。