请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 270|回复: 4

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

[复制链接]

24

帖子

0

威望

120

eV
积分
144

Level 2 能力者

发表于 2018-2-15 07:22:32 | 显示全部楼层 |阅读模式
各位好!
我最近在用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. 如何在输入文件中设置, 可以打印出更多的系数呢? 谢谢!







1万

帖子

25

威望

1万

eV
积分
30943

管理员

公社社长

发表于 2018-2-15 07:52:57 | 显示全部楼层
如果不是非得用nwchem不可,建议用Gaussian或ORCA,基于其输出,Multiwfn直接就可以很方便地给出跃迁密度(矩阵形式以及实空间函数形式,还可以得到跃迁电荷)。用NWChem的话,本身其轨道波函数信息导出来就不方便,自己写程序搞很麻烦。

基于TDDFT计算跃迁密度的公式在Multiwfn手册3.21.1.1节里有。
北京科音自然科学研究中心:http://www.keinsci.com  不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训
思想家公社的门口Blog:http://sobereva.com  发布大量原创计算化学相关博文
Multiwfn量子化学波函数分析程序主页:http://sobereva.com/multiwfn
计算化学公社论坛:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

24

帖子

0

威望

120

eV
积分
144

Level 2 能力者

 楼主| 发表于 2018-2-15 08:43:32 | 显示全部楼层
sobereva 发表于 2018-2-15 07:52
如果不是非得用nwchem不可,建议用Gaussian或ORCA,基于其输出,Multiwfn直接就可以很方便地给出跃迁密度( ...

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

tdm

tdm

1万

帖子

25

威望

1万

eV
积分
30943

管理员

公社社长

发表于 2018-2-15 20:43:36 | 显示全部楼层
besphed 发表于 2018-2-15 08:43
谢谢sob老师!
在Multiwfn手册3.21.1.1节里, Theory 3: Transition density matrix and transition densi ...

无标题.png
北京科音自然科学研究中心:http://www.keinsci.com  不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训
思想家公社的门口Blog:http://sobereva.com  发布大量原创计算化学相关博文
Multiwfn量子化学波函数分析程序主页:http://sobereva.com/multiwfn
计算化学公社论坛:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

24

帖子

0

威望

120

eV
积分
144

Level 2 能力者

 楼主| 发表于 2018-2-16 00:07:53 | 显示全部楼层

感谢sob老师! 祝您新年快乐,万事如意!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-5-24 20:02 , Processed in 0.175188 second(s), 26 queries .

快速回复 返回顶部 返回列表