计算化学公社

标题: 福井函数相关问题：q（N）、q（N+1）、q（N-1）代表的是什么 [打印本页]

作者
Author: chiara 时间: 2019-9-5 11:24
标题: 福井函数相关问题：q（N）、q（N+1）、q（N-1）代表的是什么
本帖最后由 chiara 于 2019-9-5 11:25 编辑

图1、2两张图是我计算的福井函数相关的数据，q（N）、q（N+1）、q（N-1）代表的是什么，是不是就是电子的分布情况啊，是NPA吗，还是muliken电荷呀，是否就是图3中的charge（0）、charge（+）、charge（-），求老师指点

作者
Author: sobereva 时间: 2019-9-5 11:35
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484（http://bbs.keinsci.com/thread-13328-1-1.html）

文中已经说得非常明确了

上述原子指数都可以写为基于原子电荷计算的形式。根据一些文章的测试，Hirshfeld电荷非常适合这种目的（不建议用NPA等其它电荷），所以Multiwfn自动计算它们都是用的Hirshfeld电荷。

显然既不是NPA也不是Mulliken电荷

仔细把Multiwfn手册4.7.3节的例子看了自然就明白诸如q是什么

更多信息看
概念密度泛函综述和重要文献合集 http://bbs.keinsci.com/thread-384-1-1.html

作者
Author: chiara 时间: 2019-9-5 11:49

sobereva 发表于 2019-9-5 11:35
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484（http://bbs.kein ...

谢谢

作者
Author: Changwenxu 时间: 2019-10-22 14:42

sobereva 发表于 2019-9-5 11:35
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484（http://bbs.kein ...

sob老师，我也有问题请教一下您

如下是我按Multiwfn "Calculate various quantitative indices"功能计算的结果
Atom    q(N) q(N+1) q(N-1)    f-       f+       f0       CDD
   1(C ) 0.0490  -0.0496 0.1485 0.0995 0.0986 0.0991  -0.0009
   2(C )  -0.1195  -0.2772 0.0725 0.1920 0.1577 0.1749  -0.0342
   3(H ) 0.0387  -0.0265 0.1025 0.0638 0.0652 0.0645 0.0015
   4(H ) 0.0363  -0.0154 0.0889 0.0526 0.0517 0.0522  -0.0009
   5(C )  -0.0843  -0.1119  -0.0558 0.0285 0.0275 0.0280  -0.0010
   6(H ) 0.0376 0.0024 0.0764 0.0388 0.0352 0.0370  -0.0036
   7(H ) 0.0375 0.0090 0.0651 0.0276 0.0285 0.0281 0.0009
   8(H ) 0.0396 0.0032 0.0786 0.0390 0.0364 0.0377  -0.0026
   9(N )  -0.1579  -0.2312  -0.0473 0.1105 0.0734 0.0920  -0.0372
10(C ) 0.1624 0.0277 0.2347 0.0723 0.1347 0.1035 0.0624
11(N )  -0.1538  -0.2603  -0.0174 0.1363 0.1066 0.1215  -0.0298
12(C )  -0.0146  -0.0674 0.0227 0.0374 0.0528 0.0451 0.0155
13(H ) 0.0451 0.0045 0.0800 0.0349 0.0406 0.0378 0.0057
14(H ) 0.0441 0.0052 0.0745 0.0304 0.0389 0.0346 0.0085
15(H ) 0.0399  -0.0121 0.0761 0.0362 0.0521 0.0441 0.0159
Condensed local electrophilicity/nucleophilicity index (e*eV)
   Atom             Electrophilicity       Nucleophilicity
   1(C )                0.05062                0.30843
   2(C )                0.08097                0.59492
   3(H )                0.03348                0.19763
   4(H )                0.02653                0.16310
   5(C )                0.01414                0.08838
   6(H )                0.01808                0.12027
   7(H )                0.01464                0.08566
   8(H )                0.01868                0.12080
   9(N )                0.03766                0.34253
10(C )                0.06914                0.22415
11(N )                0.05470                0.42251
12(C )                0.02710                0.11575
13(H )                0.02085                0.10826
14(H )                0.01996                0.09419
15(H )                0.02673                0.11218
分析福井f+,C2比C10更易被亲核试剂进攻。而CDD(我理解的是Δf)分析得C10最易被亲核进攻，而C2最易被亲电进攻。这两者是否矛盾？
然后后面的简缩局部亲电亲核指数，结论又给出C2无论是亲电还是亲核能力都优于C10。又和上面都矛盾。
请您指点。

作者
Author: sobereva 时间: 2019-10-23 04:47

Changwenxu 发表于 2019-10-22 14:42
sob老师，我也有问题请教一下您
如下是我按Multiwfn "Calculate various quantitative indices"功能 ...

当前情况不容易下结论
考虑到CDD和Hirshfeld电荷结论相同，都体现C10比C2更容易发生亲核反应，应当认为C10发生的概率更高一些

作者
Author: Changwenxu 时间: 2019-10-23 09:42
本帖最后由 Changwenxu 于 2019-10-23 17:20 编辑

sobereva 发表于 2019-10-23 04:47
当前情况不容易下结论
考虑到CDD和Hirshfeld电荷结论相同，都体现C10比C2更容易发生亲核反应，应当认为C ...

好的，谢谢老师。但还是几点我没理解到位，打扰您了
1.Hirshfeld 电荷，是用来算FF和DD的数据，直接来看是看q(N+1)值更小者更易被亲核进攻？
2.如您所说“更易发生亲核反应”是指更易被亲核进攻还是更易去亲核进攻？
3.简缩局部亲电亲核指数何时适用？

作者
Author: sobereva 时间: 2019-10-24 07:15

Changwenxu 发表于 2019-10-23 09:42
好的，谢谢老师。但还是几点我没理解到位，打扰您了
1.Hirshfeld 电荷，是用来算FF和DD的数据，直接来看 ...

1 看的是q(N)，即你这个体系的中性状态下Hirshfeld电荷
2 被亲核试剂进攻
3 讨论当前体系各个原子亲电性、亲核性时候用

作者
Author: Changwenxu 时间: 2019-10-24 09:13

sobereva 发表于 2019-10-24 07:15
1 看的是q(N)，即你这个体系的中性状态下Hirshfeld电荷
2 被亲核试剂进攻
3 讨论当前体系各个原子亲电 ...

好的，谢谢您

作者
Author: Cassie 时间: 2022-6-8 11:09
q（N）、q（N+1）、q（N-1）代表的是什么？q（n）是零价电荷，q（N+1）是正一价的吗？q（N-1）是负一价的？

作者
Author: wzkchem5 时间: 2022-6-8 14:30

Cassie 发表于 2022-6-8 04:09
q（N）、q（N+1）、q（N-1）代表的是什么？q（n）是零价电荷，q（N+1）是正一价的吗？q（N-1）是负一价的？

q(N)是你感兴趣的价态的电荷，这个价态不一定是0价的。q(N+1)是比你感兴趣的价态低一价的电荷（注意不是高一价），q(N-1)是比你感兴趣的价态高一价的电荷。
只要知道这里的N代表电子数，就可以知道为什么是这样的。

作者
Author: sobereva 时间: 2022-6-17 19:47

Cassie 发表于 2022-6-8 11:09
q（N）、q（N+1）、q（N-1）代表的是什么？q（n）是零价电荷，q（N+1）是正一价的吗？q（N-1）是负一价的？

把此文最后一节看了，专门说得非常详细
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484（http://bbs.keinsci.com/thread-13328-1-1.html）

作者
Author: lonemen 时间: 2022-6-17 20:52

wzkchem5 发表于 2022-6-8 14:30
q(N)是你感兴趣的价态的电荷，这个价态不一定是0价的。q(N+1)是比你感兴趣的价态低一价的电荷（注意不是 ...

多谢多谢！

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