计算化学公社

标题: molpro的transition charge density求助 [打印本页]

作者
Author: besphed 时间: 2017-11-26 09:43
标题: molpro的transition charge density求助
大家好!
我最近用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 是梯度算符.

谢谢大家!

作者
Author: sobereva 时间: 2017-11-26 11:21
sigma(x,y,z) = \sum_{i,j} rho_{i,j} phi_i phi_j式子没错

作者
Author: besphed 时间: 2017-11-26 11:53

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

谢谢sob老师!

作者
Author: besphed 时间: 2017-11-26 14:16

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

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

作者
Author: sobereva 时间: 2017-11-26 14:42
这就是把希尔伯特空间下的跃迁密度矩阵变换为实空间形式（式子其实在Multiwfn手册3.21.1.1节也有）
可以看看这个文档，是Mcweeny书里的一部分，对密度矩阵讲得十分详细
(, 下载次数 Times of downloads: 9)

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