molpro的transition charge density求助

besphed · 发表于 Post on 2017-11-26 09:43:50

大家好!
我最近用molpro的casscf 和 ci做一些激发态的计算和分析.
现在, 我得到了AO基(就是Atomic orbitals basis)下的transition density matrix, 例如基态0与激发态1, 或激发态2, 3等等的 transition density matrix.
下面, 我想做一下相应的transition charge density. 请问应该怎么做?

我的想法是, 既然得到AO基下的transition density matrix, 记为矩阵形式, 其矩阵元就是 rho_{ij},其中, i,j 相应于第 i, 第j个原子轨道(AO)
而原子轨道, 可以标记为 phi_i, phi_j, 等等;
那么, transition charge density (记为sigma)就可以如下计算:
sigma(x,y,z) = \sum_{i,j} rho_{i,j} phi_i phi_j

请问, 上述想法是否正确? 还请推荐一些相关的文献.谢谢!

此外, 我还想计算transition charge current density, 又应该怎么计算? 是否有相关的文献?
我的想法如下:
j(x,y,z) =\sum_{i,j} [ rho_i \nabla rho_j - rho_j \nabla rho_i ],
其中, j(x,y,z) 即为transition charge current density; 而 \nabla 是梯度算符.

谢谢大家!

sobereva · 发表于 Post on 2017-11-26 11:21:14

sigma(x,y,z) = \sum_{i,j} rho_{i,j} phi_i phi_j式子没错

besphed · 发表于 Post on 2017-11-26 11:53:49

sobereva 发表于 2017-11-26 11:21
sigma(x,y,z) = \sum_{i,j} rho_{i,j} phi_i phi_j式子没错

谢谢sob老师!

besphed · 发表于 Post on 2017-11-26 14:16:54

sobereva 发表于 2017-11-26 11:21
sigma(x,y,z) = \sum_{i,j} rho_{i,j} phi_i phi_j式子没错

还有, 请问这个关于transition charge density 的式子, 是怎么得到的?
请问有哪些文献可以参考,学习一下?
谢谢!!

sobereva · 发表于 Post on 2017-11-26 14:42:11

这就是把希尔伯特空间下的跃迁密度矩阵变换为实空间形式（式子其实在Multiwfn手册3.21.1.1节也有）
可以看看这个文档，是Mcweeny书里的一部分，对密度矩阵讲得十分详细

Chapter5, The Electron Distribution-from Methods of molecular quantum mechanics .pdf (1.63 MB, 下载次数 Times of downloads: 9)

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Molpro] molpro的transition charge density求助