D1，D2的给电子能力比较与A1，A2的吸电子能力比较的计算方法

yezhonghua

    关于D1，D2的给电子能力比较与A1，A2的吸电子能力比较的计算方法。
我的想法是：
     1.D1-D2（A1-A2)连在一起，进行opt+freq优化，观察HOMO-LUMO电子云分布。比较D1，D2片段对HOMO、LUMO的贡献大小。哪个片度对于HOMO贡献大，则其在D1-D2中为D型。这种方法万一出现分布差不多的情况怎么办？
    2.给电子能力考查Ip能力；吸电子能力考查EA。通过计算D1，D2 Ip大小判断其给电子能力；A1，A2的EA大小判断其吸电子能力。
这种方法是否会对于EA，Ip的计算有可能不准导致判断错误。

望大师指点一二。非常感谢！

liyuanhe211

D1, D2, A1, A2 是啥？

sobereva

虽然1、2都有一定道理，但更省事、直接的做法是，对同一个体系，接上A1、A2，比较两种情况下谁的净电荷更负。对于D1、D2的比较也是类似。

上面我说的是基团供电子、吸电子能力，但如果你说的是电子激发过程中作为电子给体、电子受体的能力，那是另外一回事，二者没直接关系。

alwens

计算电荷比看轨道更好

小范范1989

sobereva 发表于 2016-4-28 04:48
虽然1、2都有一定道理，但更省事、直接的做法是，对同一个体系，接上A1、A2，比较两种情况下谁的净电荷更负 ...

sob老师好，您说的这个，我的理解如下
1：基于基态来看，就是楼主这样说的，连接起来看电荷分布，比如说D1-D2，通过基态计算，发现D1净电荷是-0.5，那么D1就是基团吸电子。
2：如果研究的是激发过程中的donor和accpetor。那么需要在基态的基础上，做激发，看激发态和基态的电荷的差。（这个我当时在群里咨询过您）。
有很多人，对于两者的却别不清楚，所以当时我上量子化学的时候有疑问，然后在计算化学公社发帖子问过，
不知道上面的理解对不对，谢谢sob老师

yezhonghua

liyuanhe211 发表于 2016-4-28 00:32
D1, D2, A1, A2 是啥？

D1，D2就是类似咔唑，吖啶之类的给电子基团；A1，A2就是类似三嗪，吡啶等吸电子基团；这是在OLED材料体系来说的！

yezhonghua

sobereva 发表于 2016-4-28 04:48
虽然1、2都有一定道理，但更省事、直接的做法是，对同一个体系，接上A1、A2，比较两种情况下谁的净电荷更负 ...

     谢谢sober老师.
     老师的意思是D1-A，D2-A；计算两者的基态，然后比较D1-A中D1的净电荷和D2-A中的D2的净电荷大小就可以了，这样保证是一个体系下的对比。

yezhonghua

alwens 发表于 2016-4-28 08:18
计算电荷比看轨道更好

麻烦问一下，净电荷log文件中能够直接读取么？或者还是可以通过multiwfn中的方法找到。谢谢

alwens

yezhonghua 发表于 2016-4-28 14:26
麻烦问一下，净电荷log文件中能够直接读取么？或者还是可以通过multiwfn中的方法找到。谢谢

默认mulliken 是有的。

yezhonghua

sobereva 发表于 2016-4-28 04:48
虽然1、2都有一定道理，但更省事、直接的做法是，对同一个体系，接上A1、A2，比较两种情况下谁的净电荷更负 ...

sober老师，我将D1-D2链接为一个体系，然后进行基态优化，通过multiwfn里面的布居分析，得到D1的net charge 和D2的net charge 比较两个数值，为负数的片段具有吸电子的能力，也就是其供电子能力较弱。
谢谢老师！

sobereva

小范范1989 发表于 2016-4-28 08:29
sob老师好，您说的这个，我的理解如下
1：基于基态来看，就是楼主这样说的，连接起来看电荷分布，比如说 ...

对

sobereva

yezhonghua 发表于 2016-4-28 14:26
麻烦问一下，净电荷log文件中能够直接读取么？或者还是可以通过multiwfn中的方法找到。谢谢

高斯本身能算一些原子电荷，如Mulliken、CHELPG、NPA等。Multiwfn也可以算原子电荷，Mulliken、CHELPG、MK、Hirshfeld、ADCH等。默认输出Mulliken电荷。
建议用ADCH或Hirshfeld或NPA来讨论当前问题。差异详见物理化学学报中的《原子电荷计算方法的对比》一文。