关于“跃迁势”

肖富贵

最近着手做跃迁电荷（transition charge），Multiwfn可以计算Mulliken transition charge，但此电荷并不完美，因此想自己拟合。
拟合分两步，第一步生成静电势；第二步则是用电荷拟合第一步得到的静电势。
第一步用到的公式如下：

等式左边：𝜙为跃迁势（我自己不知道”跃迁势“这种叫法是否合适），B代表B分子，b代表基态，b'代表激发态。
等式右边：积分里边的𝜌是one-particle density，其定义如下：

等式右边第二项一般等于0，因为我们一般对激发态-基态之间的耦合感兴趣。

我想问的是，这种“跃迁势” 可否由Multiwfn求出来呢？

跪谢跪谢跪谢

肖富贵

不好意思，帖子中有两处乱码。第一处是phi，第二处是rho。从word中直接粘贴过来就乱码了，抱歉。

sobereva

Multiwfn当前的代码直接给不出来跃迁势，但是改改代码能得到。
function.f90里的function eleesp是计算电子密度对静电势的贡献的。计算的时候会循环每一对儿高斯函数(GTF)，循环中会看到注释：
!Now calc "term"=<psi(iprim)|1/r_Z|psi(jprim)>
也就是说下面的term变量就是<GTF_i|1/r_Z|GTF_j>这个积分。Multiwfn利用这个积分，乘上轨道对GTF的展开系数，乘上轨道占据数再累加，得到对静电势的贡献。对于获得跃迁势，在展开系数这块动动手脚改成对应跃迁密度矩阵的情况。

肖富贵

sobereva 发表于 2016-1-7 00:03
Multiwfn当前的代码直接给不出来跃迁势，但是改改代码能得到。
function.f90里的function eleesp是计算电 ...

跪谢跪谢跪谢.

肖富贵

sobereva 发表于 2016-1-7 00:03
Multiwfn当前的代码直接给不出来跃迁势，但是改改代码能得到。
function.f90里的function eleesp是计算电 ...

Sob前辈，您好. 我看了您的帖子后，试着重新对此公式做了调整，不知道对与不对。请问您可不可以帮忙看一下.
因为不能直接从word中粘贴，所以就截图了。请见谅。

sobereva

肖富贵 发表于 2016-1-18 01:56
Sob前辈，您好. 我看了您的帖子后，试着重新对此公式做了调整，不知道对与不对。请问您可不可以帮忙看一 ...

正确

肖富贵

sobereva 发表于 2016-1-18 06:55
正确

跪谢跪谢跪谢Sob前辈.

肖富贵

sobereva 发表于 2016-1-18 06:55
正确

Sob前辈，您好.我试着修改Multiwfn中的某一行代码，来计算跃迁势. 修改的内容如下：


请问Sob前辈，这样的修改合理吗？

sobereva

肖富贵 发表于 2016-1-19 19:21
Sob前辈，您好.我试着修改Multiwfn中的某一行代码，来计算跃迁势. 修改的内容如下：

合理，但你得把公式里的N乘上去

肖富贵

sobereva 发表于 2016-1-20 03:47
合理，但你得把公式里的N乘上去

跪谢跪谢跪谢Sob前辈.

肖富贵

sobereva 发表于 2016-1-20 03:47
合理，但你得把公式里的N乘上去

Sob前辈，非常不好意思，再请问您一个问题：
因为wfn文件没有虚轨道的信息，所以需要用fch文件。请问我可不可以直接用fch文件来计算跃迁势，而不需要把fch文件转换为wfn文件？

sobereva

肖富贵 发表于 2016-1-20 18:56
Sob前辈，非常不好意思，再请问您一个问题：
因为wfn文件没有虚轨道的信息，所以需要用fch文件。请问我 ...

可以

肖富贵

sobereva 发表于 2016-1-20 20:06
可以

多谢多谢多谢Sob前辈，非常不好意思总是麻烦您。