计算化学公社

标题: 请教关于nwchem中transition charge density的计算 [打印本页]

作者
Author: besphed 时间: 2018-2-15 07:22
标题: 请教关于nwchem中transition charge density的计算
各位好!
我最近在用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. 如何在输入文件中设置, 可以打印出更多的系数呢? 谢谢!

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

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

作者
Author: besphed 时间: 2018-2-15 08:43

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

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

作者
Author: sobereva 时间: 2018-2-15 20:43

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

(, 下载次数 Times of downloads: 81)

作者
Author: besphed 时间: 2018-2-16 00:07

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

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

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