| 8 | 1/1 | 返回列表 |
| 查看: 1626 | 回復: 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 |

木蟲 (正式寫手)
|
附:是計算單點能的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 |

至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
你算的是簡縮福井函數(condensed Fukui function)而非福井函數。 福井函數是連續(xù)的函數,由于軌道的弛豫,盡管積分值為1,但是空間分布上卻可以有正有負,甚至有的化學現象還必須通過負值部分才能說明,這在一定程度上也體現出福井函數相較前線軌道理論的優(yōu)勢(前線軌道理論完全忽略了軌道弛豫),更多的相關討論見JPCA, 111, 10017-10019、Chemical Reactivity Theory-A Density Functional View一書的第18章。相對應的簡縮福井函數也因此可以有負值出現。出現負值既可能是計算原子電荷的方法不合適(你所用的Mulliken電荷幾乎是最不適合計算簡縮福井函數的方法,用Hirshfeld明顯更好,負值相對來說不容易出現),也可能表現的正是實際情況,即本來就該是負的,是有物理意義的。 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
木蟲 (正式寫手)

木蟲 (正式寫手)

至尊木蟲 (著名寫手)
本人已永久離開小木蟲
木蟲 (正式寫手)

| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工328分調劑 +4 | 。,。,。,。i 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研]
|
pk3725069 2026-03-19 | 16/800 |
|
|
[考研] 0703化學求調劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
[考研] 考研化學308分求調劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
[考研] 一志愿中國石油大學(華東) 本科齊魯工業(yè)大學 +4 | 石能偉 2026-03-17 | 4/200 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 一志愿華中農業(yè)071010,總分320求調劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研] 305分求調劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 一志愿深大,0703化學,總分302,求調劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 336求調劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 0805材料320求調劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 268求調劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
[考研] 265求調劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 22 350 本科985求調劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 304求調劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 288求調劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|