sobereva 发表于 2020-4-23 06:47
我最近太忙，实在没工夫看，你可以问问作者有没有现成的程序，并且表示希望引用他们的文章

好的，我问问，谢谢老师！

万卷书万里路 发表于 2020-4-22 23:15
谢谢Sob老师！
终于觉得是明白了老师回复的这句话，我这几天重新计算了文献里的这个分子，按照理解，去 ...

我最近太忙，实在没工夫看，你可以问问作者有没有现成的程序，并且表示希望引用他们的文章

万卷书万里路 发表于 2020-4-3 12:25
谢谢Sob老师！我去学习下。
另外，我又仔细读了下文献，有句话说：The first term on the rhs of eq 7 i ...

Eq.7里i、a是指占据轨道和空轨道，显然和输出文件里的激发能没直接关系。Eq.7右边第一部分是需要基于以MO为基的跃迁密度矩阵连同占据轨道和空轨道的能量来计算。
每种MO跃迁都有对应的组态系数，这部分式子本质上就是MO能量差结合相应组态系数的平方（即MO跃迁对当前电子激发的贡献）做线性组合

sobereva 发表于 2020-4-3 11:57
和空穴、电子不是一码事
ρI(r)这里指的是实空间形式的跃迁密度，是下文说的
使用Multiwfn绘制跃迁密度矩 ...

谢谢Sob老师！我去学习下。
另外，我又仔细读了下文献，有句话说：The first term on the rhs of eq 7 is the single-particle contribution to the excitation energy, being an average of single-electron transition energies (ϵa – ϵi) weighted by the transition density matrix expansion elements ρaiI and ρiaI. 那是不是可以这样理解，激发能我们能从输出文件里得到，激发能减去这个单粒子对激发能的贡献（single-particle contribution to the excitation energy），就是第二部分这个积分了（也就是Eplas, I），
1. 那这里这个single-particle contribution to the excitation energy怎么理解，是激发态组成部分里占比最大的那个MO1->MO2的贡献么(比如H->L 80%)？那这个积分值就是 激发能*（1-n%）?
2.后半句解释我读不懂，an average of single-electron transition energies (ϵa – ϵi) weighted by the transition density matrix expansion elements ρaiI and ρiaI,这里是否就是要用到老师博文里写的 跃迁密度矩阵了？

和空穴、电子不是一码事
ρI(r)这里指的是实空间形式的跃迁密度，是下文说的
使用Multiwfn绘制跃迁密度矩阵和电荷转移矩阵考察电子激发特征（含视频演示）
http://sobereva.com/436（http://bbs.keinsci.com/thread-11383-1-1.html）
可以用Multiwfn导出跃迁密度函数的cube文件，然后自己写程序去算那个积分（可以参考Multiwfn里算空穴和电子之间库仑积分的代码）