求助关于Q-chem计算激发态间的非绝热耦合

shyshy

想请教一下~我用Q-chem做T1态和T2态间的非绝热耦合计算，输出了一下内容：
             CIS Derivative Couplings
---------------------------------------------------
              between states 1 and 2
get_w_dc time:  CPU 1.33 s  wall 0.06 s
AOints 129 time:  CPU 3.26 s  wall 0.14 s
DFT time  CPU 4.68 s  wall 0.19 s
Ei = -78.429144
Ej = -78.288937
Ej-Ei = 0.140207
---------------------------------------------------
CIS derivative coupling without ETF
  Atom         X              Y              Z
---------------------------------------------------
   1       0.000000      -0.000000      -0.000000
   2       0.000000       0.000000      -0.260430
   3       0.000000      -0.000000       0.260430
   4      -0.000000       0.000000      -0.000000
   5      -0.000000      -0.000000      -0.260430
   6      -0.000000       0.000000       0.260430
---------------------------------------------------
---------------------------------------------------
CIS Force Matrix Element
  Atom         X              Y              Z
---------------------------------------------------
   1       0.000000      -0.000000      -0.000000
   2       0.000000       0.000000      -0.032662
   3       0.000000      -0.000000       0.032662
   4      -0.000000       0.000000      -0.000000
   5      -0.000000      -0.000000      -0.032662
   6      -0.000000       0.000000       0.032662
---------------------------------------------------
---------------------------------------------------
CIS derivative coupling with ETF
  Atom         X              Y              Z
---------------------------------------------------
   1       0.000000      -0.000000      -0.000000
   2       0.000000       0.000000      -0.232954
   3       0.000000      -0.000000       0.232954
   4      -0.000000       0.000000      -0.000000
   5      -0.000000      -0.000000      -0.232954
   6      -0.000000       0.000000       0.232954
---------------------------------------------------



非绝热耦合元是波函数相对核坐标的一阶导数
请问应该如何计算才是耦合元的值呢？ 是Σ(xi2+yi2+zi2)^2 对所有原子坐标平方求和开根号，还是(Σxi)^2+(Σyi)^2+(Σzi)^2先对所有x、y、z的值进行加和再平方后开根号呢？

Mikasa

nonadiabatic coupling matrix本身就不是scalar

Siennasan

你好，我也有相同的问题，不知道是否已经解决