计算化学公社

标题: 使用Multiwfn手撕电荷转移积分（基于FCD方法） [打印本页]

作者
Author: Puying 时间: 2020-1-7 20:03
标题: 使用Multiwfn手撕电荷转移积分（基于FCD方法）
本帖最后由 Puying 于 2020-1-7 22:54 编辑

首先看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凑出来算的，不保证正确。。。错了的话原谅我，【捂脸】

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

作者
Author: Puying 时间: 2020-1-7 22:44

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

社长威武！！！

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

作者
Author: Puying 时间: 2021-11-6 20:17

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

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

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

作者
Author: 虎王 时间: 2021-11-7 16:20

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

手撕是什么专业词汇吗？

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

作者
Author: Puying 时间: 2021-11-19 23:44
本帖最后由 Puying 于 2021-11-19 23:48 编辑

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

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

作者
Author: RAL 时间: 2022-7-22 16:11
我最近正好也需要算这个积分，对着帖子和原文手撕了一下，觉得△q12好像应该是sqrt(97.10%*2.81%)+sqrt(2.90%*97.19%)=0.33307

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