使用Multiwfn手撕电荷转移积分（基于FCD方法）

Puying

首先看FCD的公式  原文链接： https://doi.org/10.1063/1.1502255
H=（E2-E1)*|△q12)|  /  sqrt[   (△q1-△q2)^2  + 4*△q12^2 ]

比如计算HOMO和HOMO-1轨道之间的电荷转移积分。
E1和E2分别是轨道能量，E2-E1为轨道能量差值。可以通过Multiwfn查看轨道得到轨道能量
△q12为transition charge difference。
△q1-△q2为adiabatic charge difference。
以上两个值可以通过Multiwfn的功能8的轨道分解分析得到

举例（C13H28烷烃二聚体）：
                     donor            acceptor
HOMO           97.10%          2.90%
HOMO-1        2.81%            97.19%

以上分别为HOMO和HOMO-1轨道在两个分子上的分布占比，则：
△q1-△q2=(97.10%-2.90%)-(2.81%-97.19%)=1.8858
△q12=sqrt(97.10%*2.90%)+sqrt(2.81%*97.19%)=0.33297
由于E2-E1=0.000149 hartree=0.0040545 eV
则 H=0.0040545*0.33297/sqrt(1.8858^2+4*0.33297^2)=0.675meV     用site energy correction 方法计算的是0.716meV。我也试了其他两三个例子，结果也都差不多

以上是我读了FCD原文加上Multiwfn凑出来算的，不保证正确。。。错了的话原谅我，【捂脸】

sobereva

Multiwfn里其实原本直接就有FCD和GMH计算功能，只不过代码完成99%后玩去了，于是一直没显示在主界面里也没写到手册里... 后来把这部分代码挪到一个新坑里了，之后有时间的时候会发布

Puying

sobereva 发表于 2020-1-7 22:14
Multiwfn里其实原本直接就有FCD和GMH计算功能，只不过代码完成99%后玩去了，于是一直没显示在主界面里也没 ...

社长威武！！！

pizaipizai@5989

老哥，我读了你贴的原文献，看你那意思是说把N聚体的结构扒成一个结构文件，然后一次计算一下单点，最后再mwf提取你说的这些参量，最后扔到公式里手撕一下。是这个意思吗？
我还有个想法，对应空穴的传输就用中性N聚体，要改为研究电子的传输我们给N聚体加上一个电荷，然后自旋多重度改为1，再去算单点，再用mwf分析最高能量那个轨道的参量，再手撕一下,

Puying

pizaipizai@5989 发表于 2021-10-25 19:31
老哥，我读了你贴的原文献，看你那意思是说把N聚体的结构扒成一个结构文件，然后一次计算一下单点，最后再m ...

N聚体不晓得，2聚体是可以的。如果是2聚体的话，算一个2聚体的能量就好了，不需要分开算两个分子的单点。如果两个分子一模一样，那么大概率homo和homo-1 是分布在两个分子上面的。

空穴算homo和homo-1，电子算lumo和lumo+1。这个方法不需要加电子

虎王

pizaipizai@5989 发表于 2021-10-25 19:31
老哥，我读了你贴的原文献，看你那意思是说把N聚体的结构扒成一个结构文件，然后一次计算一下单点，最后再m ...

手撕是什么专业词汇吗？

MAX-0309

C:\Users\Muyakui\Desktop
你好，我读了一下这个公式的部分，这里的δq1 和 δq2应该是对应轨道上的donor和acceptor的电荷差值（charges difference），楼主这里直接用分子的轨道占比相减，轨道差分和电荷差分应该是不对等的吧，盼告知@sobereva sob老师如果有时间也请帮忙看一下，谢谢

Puying

虎王 发表于 2021-11-7 16:20
手撕是什么专业词汇吗？

没有没有，我自己瞎叫的。手撕比较形象，就是自己拿手算。

RAL

我最近正好也需要算这个积分，对着帖子和原文手撕了一下，觉得△q12好像应该是sqrt(97.10%*2.81%)+sqrt(2.90%*97.19%)=0.33307