计算化学公社

标题: CASSCF如何输出旋轨耦合矩阵元 [打印本页]

作者
Author: BinWang 时间: 2026-1-30 11:05
标题: CASSCF如何输出旋轨耦合矩阵元
本帖最后由 BinWang 于 2026-1-30 11:05 编辑

各位老师，大家好！
   我用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的转化速率在皮秒量级。

作者
Author: wal 时间: 2026-1-30 11:15
本帖最后由 wal 于 2026-1-30 11:18 编辑

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

作者
Author: BinWang 时间: 2026-1-30 11:30
本帖最后由 BinWang 于 2026-1-30 11:34 编辑

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

作者
Author: wal 时间: 2026-1-30 11:45

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

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

作者
Author: BinWang 时间: 2026-1-30 11:46
本帖最后由 BinWang 于 2026-1-30 11:49 编辑

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

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

作者
Author: wal 时间: 2026-1-30 12:51

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

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

作者
Author: BinWang 时间: 2026-1-30 13:32

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

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

作者
Author: 风起~ 时间: 2026-2-1 17:04
你这轨道指定地太随意了吧
这样的结果有意义么

作者
Author: BinWang 时间: 2026-2-2 08:59

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)