| 5 | 1/1 | 返回列表 |
| 查看: 854 | 回復: 2 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
能否用multiwfn計算重疊布居電荷
|
|||
|
我又厚臉皮來了,還是關(guān)于zhou2009老師的“關(guān)于HF、H2O的電荷計算”帖子。 原帖中先是用multiwfn做出了Δρ圖,據(jù)圖分析了分子軌道的性質(zhì),然后用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老師帖子里的計算結(jié)果為 “MO3中H上電荷為0.2649。從它的Δρ得到所謂重疊布居電荷為0.6298” “MO4的O、H1、H2三者之間的重疊布居電荷為0.3293” 這樣來看,mo3我得到的兩個盆積分加和大約就是zhou2009老師帖子里的計算結(jié)果。而mo4就不好說了,我這里兩個盆積分值加和,正負抵消成-0.03了。不知道-0.32的盆是不是就對應(yīng)了O的電荷集中轉(zhuǎn)移到2H上的值,可不可以理解成zhou2009老師帖子里說的mo4上的OH重疊布局全劃分給H了,那么-0.32和zhou老師給出的0.3293對應(yīng)呢? 我自己無法準確判斷,主要是因為我用高斯計算的水分子的mulliken布居分析結(jié)果和zhou2009老師帖子里的不太一致,那么就沒法和帖子里的數(shù)據(jù)直接對應(yīng)了,特來發(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矩陣是一致的,但這是考慮了所有軌道的結(jié)果。接下來選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 (最后一行是相應(yīng)列的加和) 如果你想要得到MO2當中1和3號原子的重疊布居,那么就是0.12815040乘以2。 Multiwfn也可以獲得特定MO當中基函數(shù)之間的重疊布居。對于Mulliken布居的分析,Multiwfn是現(xiàn)有任何程序中最強大的。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +5 | 司空. 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +3 | 何潤采123 2026-03-18 | 3/150 |
|
|
[考研] 281求調(diào)劑(0805) +6 | 煙汐憶海 2026-03-16 | 15/750 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 一志愿985,本科211,0817化學工程與技術(shù)319求調(diào)劑 +7 | Liwangman 2026-03-15 | 7/350 |
|
|
[考研]
|
.6lL 2026-03-18 | 4/200 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 化學工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
[考研] 080500,材料學碩302分求調(diào)劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|