| 8 | 1/1 | 返回列表 |
| 查看: 1625 | 回復: 7 | ||
gq0502木蟲 (正式寫手)
|
[求助]
為什么計算出來的凝聚Fukui函數(shù)值有的是負的呢?
|
|
如題,計算出來的Fukui函數(shù)值有正有負,文獻中都是正的啊,大家能告訴一下我怎么計算嗎? 我用的公式就是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數(shù)值,以下是其中一個輸入文件 %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 |

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

木蟲 (正式寫手)

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

| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +3 | 超級伊昂大王 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +5 | 匹克i 2026-03-23 | 5/250 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
[考研] 材料專碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|
[考研] 070300化學求調(diào)劑 +8 | 苑豆豆 2026-03-20 | 8/400 |
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 3/150 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|
[考研] 材料學碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 化學求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學大學 ,080500材料科學與工程學碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 材料學求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 收復試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|