请教关于nwchem中transition charge density的计算

besphed

各位好!
我最近在用nwchem中的tddft模块做激发态的计算. 输入文件中设置如下:
TDDFT
CIS        
NROOTS  100
notriplet
transden
civecs
END
然后得到了如下的输出文件: (以下为第一个激发态的信息)
----------------------------------------------------------------------------
Root   1 singlet a              0.111500118 a.u.                3.0341 eV
----------------------------------------------------------------------------
    Transition Moments    X  0.00000   Y -0.00000   Z  0.00000
    Transition Moments   XX  0.00000  XY -0.00000  XZ  0.00000
    Transition Moments   YY -0.00000  YZ -0.00000  ZZ  0.00000
    Dipole Oscillator Strength                    0.0000000000
    Electric Quadrupole                           0.0000000000
    Magnetic Dipole                               0.0000000000
    Total Oscillator Strength                     0.0000000000

    Occ.   26  a   ---  Virt.   38  a       -0.05015
    Occ.   27  a   ---  Virt.   37  a        0.05015
    Occ.   28  a   ---  Virt.   35  a        0.11617
    Occ.   29  a   ---  Virt.   34  a       -0.11617
    Occ.   30  a   ---  Virt.   33  a       -0.68902
    Occ.   31  a   ---  Virt.   32  a        0.68902
当然, 还有其他激发的信息, 输出格式是类似的, 这里没有贴出来.

我想请教一下, 如何得到任意两个态(可以是基态和激发态, 也可以是任意两个激发态)之间的transition charge density?
我对以上的输出文件的分析, 例如从Occ. 和 Virt. 可以知道相应于该激发态的configuration state function, 从而构建相应于该激发态的波函数Phi_i(r).
那么,任意两个激发态之间的transition charge density, 就是Phi_i(r) Phi_j(r), 不知道我的理解对不对? 请指教!
当然, 激发态与基态之间的transition charge density, 就是Phi_i(r) Phi_0(r), 其中, Phi_0(r)是基态波函数.
不过以上都需要自己编程实现. 不知道nwchem有没有更直接,简便的方法,给出与transition charge density有关的信息? 请指教!

另外, 我还想请教, 如何让输出文件打印出更多的configuration state function 的系数, 上面的例子中, 只打印出了较大的一些系数, -0.05015, 0.05015, ..., 0.68902. 如何在输入文件中设置, 可以打印出更多的系数呢? 谢谢!

sobereva

如果不是非得用nwchem不可，建议用Gaussian或ORCA，基于其输出，Multiwfn直接就可以很方便地给出跃迁密度（矩阵形式以及实空间函数形式，还可以得到跃迁电荷）。用NWChem的话，本身其轨道波函数信息导出来就不方便，自己写程序搞很麻烦。

基于TDDFT计算跃迁密度的公式在Multiwfn手册3.21.1.1节里有。

besphed

sobereva 发表于 2018-2-15 07:52
如果不是非得用nwchem不可，建议用Gaussian或ORCA，基于其输出，Multiwfn直接就可以很方便地给出跃迁密度（ ...

谢谢sob老师!
在Multiwfn手册3.21.1.1节里, Theory 3: Transition density matrix and transition density给出了基态与激发态之间的跃迁密度矩阵. 对任意两个激发态之间的跃迁密度矩阵呢? 不能直接用如下的式子了吧?
谢谢!

sobereva

besphed 发表于 2018-2-15 08:43
谢谢sob老师!
在Multiwfn手册3.21.1.1节里, Theory 3: Transition density matrix and transition densi ...

besphed

sobereva 发表于 2018-2-15 20:43

感谢sob老师! 祝您新年快乐,万事如意!