计算化学公社

 找回密码 Forget password
 注册 Register
Views: 10883|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[Multiwfn使用咨询] 福井函数相关问题:q(N)、q(N+1)、q(N-1)代表的是什么

[复制链接 Copy URL]

14

帖子

0

威望

81

eV
积分
95

Level 2 能力者

本帖最后由 chiara 于 2019-9-5 11:25 编辑

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

3.png (59.36 KB, 下载次数 Times of downloads: 63)

3

3

2.png (49.44 KB, 下载次数 Times of downloads: 83)

2

2

1.png (36.16 KB, 下载次数 Times of downloads: 58)

1

1

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

2#
发表于 Post on 2019-9-5 11:35:10 | 只看该作者 Only view this author
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484http://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



北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

81

eV
积分
95

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2019-9-5 11:49:17 | 只看该作者 Only view this author
sobereva 发表于 2019-9-5 11:35
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484(http://bbs.kein ...

谢谢

41

帖子

0

威望

109

eV
积分
150

Level 3 能力者

4#
发表于 Post on 2019-10-22 14:42:38 | 只看该作者 Only view this author
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。又和上面都矛盾。
请您指点。

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

5#
发表于 Post on 2019-10-23 04:47:19 | 只看该作者 Only view this author
Changwenxu 发表于 2019-10-22 14:42
sob老师,我也有问题请教一下您
如下是我按Multiwfn "Calculate various quantitative indices"功能 ...

当前情况不容易下结论
考虑到CDD和Hirshfeld电荷结论相同,都体现C10比C2更容易发生亲核反应,应当认为C10发生的概率更高一些
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

41

帖子

0

威望

109

eV
积分
150

Level 3 能力者

6#
发表于 Post on 2019-10-23 09:42:26 | 只看该作者 Only view this author
本帖最后由 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.简缩局部亲电亲核指数何时适用?

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

7#
发表于 Post on 2019-10-24 07:15:03 | 只看该作者 Only view this author
Changwenxu 发表于 2019-10-23 09:42
好的,谢谢老师。但还是几点我没理解到位,打扰您了
1.Hirshfeld 电荷,是用来算FF和DD的数据,直接来看 ...

1 看的是q(N),即你这个体系的中性状态下Hirshfeld电荷
2 被亲核试剂进攻
3 讨论当前体系各个原子亲电性、亲核性时候用
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

41

帖子

0

威望

109

eV
积分
150

Level 3 能力者

8#
发表于 Post on 2019-10-24 09:13:22 | 只看该作者 Only view this author
sobereva 发表于 2019-10-24 07:15
1 看的是q(N),即你这个体系的中性状态下Hirshfeld电荷
2 被亲核试剂进攻
3 讨论当前体系各个原子亲电 ...

好的,谢谢您

1

帖子

0

威望

17

eV
积分
18

Level 1 能力者

9#
发表于 Post on 2022-6-8 11:09:31 | 只看该作者 Only view this author
q(N)、q(N+1)、q(N-1)代表的是什么?q(n)是零价电荷,q(N+1)是正一价的吗?q(N-1)是负一价的?

1万

帖子

0

威望

9870

eV
积分
22110

Level 6 (一方通行)

10#
发表于 Post on 2022-6-8 14:30:29 | 只看该作者 Only view this author
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代表电子数,就可以知道为什么是这样的。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
lonemen + 5 谢谢

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

11#
发表于 Post on 2022-6-17 19:47:17 | 只看该作者 Only view this author
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/484http://bbs.keinsci.com/thread-13328-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

245

帖子

0

威望

2723

eV
积分
2968

Level 5 (御坂)

12#
发表于 Post on 2022-6-17 20:52:13 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-8 14:30
q(N)是你感兴趣的价态的电荷,这个价态不一定是0价的。q(N+1)是比你感兴趣的价态低一价的电荷(注意不是 ...

多谢多谢!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 17:36 , Processed in 0.228899 second(s), 30 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list