| 查看: 1451 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
gq0502木虫 (正式写手)
|
[求助]
为什么计算出来的凝聚Fukui函数值有的是负的呢?
|
|
|
如题,计算出来的Fukui函数值有正有负,文献中都是正的啊,大家能告诉一下我怎么计算吗? 我用的公式就是f+=q(N+1)-q(N) f-=q(N)-q(N-1) q用的是Mulliken atomic charges值，求问这样对吗？为什么计算的值有正有负呢？求大神解答，不然下一步走不下去了……  1 1 C 0.463589 2 C 0.611716 3 C -0.908433 4 C -0.894842 5 C -0.430052 6 C -0.271609 7 N -0.336673 8 N -0.159266 9 C -0.238646 10 H 0.152731 11 H 0.168085 12 H 0.153209 13 H 0.154597 14 H 0.374051 15 H 0.161543 Sum of Mulliken charges= -1.00000 N电子体系 Mulliken atomic charges: 1 1 C 0.475274 2 C 0.582663 3 C -0.729210 4 C -0.739263 5 C -0.389324 6 C -0.296695 7 N -0.378213 8 N -0.162042 9 C 0.119814 10 H 0.213386 11 H 0.229152 12 H 0.217081 13 H 0.216587 14 H 0.425287 15 H 0.215502 Sum of Mulliken charges= 0.00000 N+1电子体系 Mulliken atomic charges: 1 1 C 0.572753 2 C 0.669263 3 C -0.697888 4 C -0.738609 5 C -0.366257 6 C -0.185471 7 N -0.304112 8 N -0.081522 9 C 0.162207 10 H 0.289698 11 H 0.303989 12 H 0.292944 13 H 0.304758 14 H 0.490514 15 H 0.287734 Sum of Mulliken charges= 1.00000 N-1电子体系 Mulliken atomic charges: 1 1 C 0.463589 2 C 0.611716 3 C -0.908433 4 C -0.894842 5 C -0.430052 6 C -0.271609 7 N -0.336673 8 N -0.159266 9 C -0.238646 10 H 0.152731 11 H 0.168085 12 H 0.153209 13 H 0.154597 14 H 0.374051 15 H 0.161543 Sum of Mulliken charges= -1.00000 |
回复此楼» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
物理化学论文润色/翻译怎么收费?
已经有132人回复
-
已知面上挂了，小学校实验方向无平台，感觉太难了
已经有26人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
Fukui函数的计算
已经有5人回复- FUKUI 函数和BADER电荷的联系？
已经有3人回复
- fukui函数是动力学性质还是热力学性质
已经有6人回复
- 【求助】如何求Fukui函数？
已经有11人回复
- 【求助】Fukui函数的计算
已经有4人回复
gq0502
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 7933.6
- 散金: 265
- 红花: 6
- 帖子: 893
- 在线: 525.2小时
- 虫号: 1422748
- 注册: 2011-09-29
- 性别: GG
- 专业: 化学工程及工业化学
5楼2013-11-24 21:48:11
gq0502
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 7933.6
- 散金: 265
- 红花: 6
- 帖子: 893
- 在线: 525.2小时
- 虫号: 1422748
- 注册: 2011-09-29
- 性别: GG
- 专业: 化学工程及工业化学
|
附：是计算单点能的log文件读的Mulliken charges数值，以下是其中一个输入文件 %nproc=4 %mem=600MB %chk=bs.chk #P b3lyp/6-311+g(d) pop=full Title Card Required 0 1 C -0.25786200 0.69021200 -0.00000100 C -0.26676300 -0.72396000 0.00000000 C 0.92329500 1.43213700 0.00000000 C 0.94330000 -1.42618300 -0.00000100 C 2.11349600 0.71399100 0.00000300 C 2.12329100 -0.69464800 0.00000200 N -1.59713200 1.04409600 -0.00001500 N -1.57251600 -1.19805500 -0.00000400 C -2.31812700 -0.12823000 0.00001200 H 0.92244000 2.51749400 -0.00000500 H 0.94520300 -2.51068200 -0.00000300 H 3.05666900 1.25053900 0.00000400 H 3.07587500 -1.21425500 0.00000100 H -1.97674000 1.97701500 0.00003700 H -3.39969400 -0.12232000 0.00001400 |
2楼2013-11-24 19:54:28
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2013-11-25 06:14:07
gq0502: 金币+25, ★★★★★最佳答案 2013-11-25 16:17:22
感谢参与，应助指数 +1
gmy1990: 金币+3 2013-11-25 06:14:07
gq0502: 金币+25, ★★★★★最佳答案 2013-11-25 16:17:22
|
你算的是简缩福井函数(condensed Fukui function)而非福井函数。 福井函数是连续的函数，由于轨道的弛豫，尽管积分值为1，但是空间分布上却可以有正有负，甚至有的化学现象还必须通过负值部分才能说明，这在一定程度上也体现出福井函数相较前线轨道理论的优势（前线轨道理论完全忽略了轨道弛豫），更多的相关讨论见JPCA, 111, 10017-10019、Chemical Reactivity Theory-A Density Functional View一书的第18章。相对应的简缩福井函数也因此可以有负值出现。出现负值既可能是计算原子电荷的方法不合适（你所用的Mulliken电荷几乎是最不适合计算简缩福井函数的方法，用Hirshfeld明显更好，负值相对来说不容易出现），也可能表现的正是实际情况，即本来就该是负的，是有物理意义的。 |
3楼2013-11-24 21:40:47
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
4楼2013-11-24 21:43:33