| 3 | 1/1 | 返回列表 |
| 查看: 856 | 回復: 2 | |||
[求助]
能否用multiwfn計算重疊布居電荷
|
|
我又厚臉皮來了,還是關于zhou2009老師的“關于HF、H2O的電荷計算”帖子。 原帖中先是用multiwfn做出了Δρ圖,據(jù)圖分析了分子軌道的性質,然后用GsGrid計算各個軌道的重疊布居電荷。最近multiwfn加了盆積分功能,我想試試用multiwfn代替GsGrid計算上述重疊布居電荷。(不想重新學個軟件了) 我的做法是,首先按之前sobereva專家給出的方法做Δρ圖,只是在“減掉波函數(shù)”時選擇功能5而不是4,生成格點文件。再選擇功能17,對內(nèi)存里的格點文件計算盆積分。 我用hf/sto-3g算了水分子的mo3和mo4: mo3的盆積分得到了兩個積分值為0.32的盆 mo4的盆積分得到了一個積分值為0.29,一個為-0.32的盆 zhou2009老師帖子里的計算結果為 “MO3中H上電荷為0.2649。從它的Δρ得到所謂重疊布居電荷為0.6298” “MO4的O、H1、H2三者之間的重疊布居電荷為0.3293” 這樣來看,mo3我得到的兩個盆積分加和大約就是zhou2009老師帖子里的計算結果。而mo4就不好說了,我這里兩個盆積分值加和,正負抵消成-0.03了。不知道-0.32的盆是不是就對應了O的電荷集中轉移到2H上的值,可不可以理解成zhou2009老師帖子里說的mo4上的OH重疊布局全劃分給H了,那么-0.32和zhou老師給出的0.3293對應呢? 我自己無法準確判斷,主要是因為我用高斯計算的水分子的mulliken布居分析結果和zhou2009老師帖子里的不太一致,那么就沒法和帖子里的數(shù)據(jù)直接對應了,特來發(fā)問,求指教。 另附我計算的水分子優(yōu)化后的坐標和mulliken布居分析,高斯版本是09W A02 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 8 0 0.000000 0.000000 0.127185 2 1 0 0.000000 0.758222 -0.508739 3 1 0 0.000000 -0.758222 -0.508739 --------------------------------------------------------------------- Full Mulliken population analysis: 1 2 3 4 5 1 1 O 1S 2.10788 2 2S -0.10780 2.00685 3 2PX 0.00000 0.00000 2.00000 4 2PY 0.00000 0.00000 0.00000 0.75079 5 2PZ 0.00000 0.00000 0.00000 0.00000 1.17276 6 2 H 1S -0.00112 -0.02504 0.00000 0.16111 0.11889 7 3 H 1S -0.00112 -0.02504 0.00000 0.16111 0.11889 6 7 6 2 H 1S 0.62634 7 3 H 1S -0.04537 0.62634 |
SOB 集錦 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
你不需要用主功能17來做盆積分。只要你用主功能5計算完了格點數(shù)據(jù),屏幕上直接就會輸出當前格點數(shù)據(jù)在整個空間的積分值。 如果你按照我之前給你回復的方法,但只是將主功能4改成了主功能5,那么這個積分值就是MO2當中O和兩個氫這兩個片段間的重疊布居值。 高斯輸出的那個矩陣是考慮了所有軌道的,而且內(nèi)容和你要分析的也不是一回事。 gsgrid早已并入Multiwfn作為主功能13了。主功能13比最后一個版本gsgrid的功能也要強大很多。所以沒必要學也沒必要用gsgrid。 如果你只想要得到某個MO當中某兩個原子間的重疊布居值,那么最簡單的方法就是載入fch文件,然后選7,再選5,之后選2。這時會輸出gross atomic population matrix,這和高斯輸出的Condensed to atoms矩陣是一致的,但這是考慮了所有軌道的結果。接下來選y,就會將每個軌道的貢獻輸出到當前目錄的groatmdcp.txt當中。從中可以找到例如 ************ Orbital 2 Occ: 2.00000000 ************ 1 2 3 1 1.36473981 0.12815040 0.12815040 2 0.12815040 0.05514825 0.00618105 3 0.12815040 0.00618105 0.05514825 4 1.62104062 0.18947970 0.18947970 (最后一行是相應列的加和) 如果你想要得到MO2當中1和3號原子的重疊布居,那么就是0.12815040乘以2。 Multiwfn也可以獲得特定MO當中基函數(shù)之間的重疊布居。對于Mulliken布居的分析,Multiwfn是現(xiàn)有任何程序中最強大的。 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 311求調劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研]
|
.6lL 2026-03-18 | 4/200 |
|
|
[考研] 0703化學 305求調劑 +3 | FY_yy 2026-03-14 | 3/150 |
|
|
[考研] 085601材料工程專碩求調劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 302求調劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 環(huán)境工程調劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 301求調劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 085601求調劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 070303 總分349求調劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 求老師收留調劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 復試調劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 復試調劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 26調劑/材料科學與工程/總分295/求收留 +9 | 2026調劑俠 2026-03-12 | 9/450 |
|
|
[考研] 0703化學求調劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 考研調劑 +4 | 芬達46 2026-03-12 | 4/200 |
|
|
[考研] 求調劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 085600材料與化工 309分請求調劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|