 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 853 | 回復(fù): 2 | |||
[求助]
能否用multiwfn計(jì)算重疊布居電荷
|
|
我又厚臉皮來(lái)了，還是關(guān)于zhou2009老師的“關(guān)于HF、H2O的電荷計(jì)算”帖子。 原帖中先是用multiwfn做出了Δρ圖，據(jù)圖分析了分子軌道的性質(zhì)，然后用GsGrid計(jì)算各個(gè)軌道的重疊布居電荷。最近multiwfn加了盆積分功能，我想試試用multiwfn代替GsGrid計(jì)算上述重疊布居電荷。(不想重新學(xué)個(gè)軟件了) 我的做法是，首先按之前sobereva專家給出的方法做Δρ圖，只是在“減掉波函數(shù)”時(shí)選擇功能5而不是4，生成格點(diǎn)文件。再選擇功能17，對(duì)內(nèi)存里的格點(diǎn)文件計(jì)算盆積分。 我用hf/sto-3g算了水分子的mo3和mo4： mo3的盆積分得到了兩個(gè)積分值為0.32的盆 mo4的盆積分得到了一個(gè)積分值為0.29，一個(gè)為-0.32的盆 zhou2009老師帖子里的計(jì)算結(jié)果為 “MO3中H上電荷為0.2649。從它的Δρ得到所謂重疊布居電荷為0.6298” “MO4的O、H1、H2三者之間的重疊布居電荷為0.3293” 這樣來(lái)看，mo3我得到的兩個(gè)盆積分加和大約就是zhou2009老師帖子里的計(jì)算結(jié)果。而mo4就不好說(shuō)了，我這里兩個(gè)盆積分值加和，正負(fù)抵消成-0.03了。不知道-0.32的盆是不是就對(duì)應(yīng)了O的電荷集中轉(zhuǎn)移到2H上的值，可不可以理解成zhou2009老師帖子里說(shuō)的mo4上的OH重疊布局全劃分給H了，那么-0.32和zhou老師給出的0.3293對(duì)應(yīng)呢？ 我自己無(wú)法準(zhǔn)確判斷，主要是因?yàn)槲矣酶咚褂?jì)算的水分子的mulliken布居分析結(jié)果和zhou2009老師帖子里的不太一致，那么就沒(méi)法和帖子里的數(shù)據(jù)直接對(duì)應(yīng)了，特來(lái)發(fā)問(wèn)，求指教。 另附我計(jì)算的水分子優(yōu)化后的坐標(biāo)和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 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 SOB 集錦 |
» 猜你喜歡
求催化作用基礎(chǔ)第三版 （甄開(kāi)吉) 課件
已經(jīng)有3人回復(fù)
-
26年材料/化學(xué)相關(guān)專業(yè)博士研究生招生 溫州大學(xué)化學(xué)與材料工程學(xué)院(26年秋季入學(xué))
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有76人回復(fù)
-
哈爾濱工業(yè)大學(xué)（深圳）綠色化學(xué)團(tuán)隊(duì)招收2026年秋季入學(xué)博士生
已經(jīng)有11人回復(fù)
-
澳大利亞昆士蘭科技大學(xué)PhD博士獎(jiǎng)學(xué)金
已經(jīng)有0人回復(fù)
-
求助電催化劑制備，怎么在碳布表面生長(zhǎng)PdCo合金
已經(jīng)有0人回復(fù)
-
求助電催化劑制備，怎么在碳布表面生長(zhǎng)PdCo合金
已經(jīng)有0人回復(fù)
-
304求調(diào)劑
已經(jīng)有3人回復(fù)
-
上海華東理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院先進(jìn)材料與制造研究組招聘公告
已經(jīng)有0人回復(fù)
-
304求調(diào)劑
已經(jīng)有1人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 分子中電荷變化種種
已經(jīng)有82人回復(fù)
- 化學(xué)鍵重疊布居數(shù)分析中的負(fù)值問(wèn)題
已經(jīng)有3人回復(fù)
- 建議multiwfn開(kāi)發(fā)GUI界面 文件打開(kāi)什么的效率太低了
已經(jīng)有10人回復(fù)
- Gaussian計(jì)算 Mulliken overlap populations 和 bond overlap population
已經(jīng)有8人回復(fù)
- 如何用Multiwfn軟件做出P DOS?
已經(jīng)有4人回復(fù)
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對(duì)論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- 使用Multiwfn繪制NBO及相關(guān)軌道
已經(jīng)有30人回復(fù)
- 【求助】NBO軌道指的是nbo電荷的分布情況嗎？
已經(jīng)有3人回復(fù)
- 用Multiwfn、Gaussview、Molekel、VMD觀看龍蝦、盆景、骨盆、大腦
已經(jīng)有16人回復(fù)
- 使用Multiwfn做拓?fù)浞治鲆约坝?jì)算孤對(duì)電子角度
已經(jīng)有18人回復(fù)
- 回答一些關(guān)于Multiwfn的疑問(wèn)以及未來(lái)Multiwfn的發(fā)展打算
已經(jīng)有37人回復(fù)
- GAUSSIAN怎么計(jì)算overlap matrix？
已經(jīng)有31人回復(fù)
- 談?wù)凪olekel做靜電勢(shì)填色等值面圖的方法及誤區(qū)
已經(jīng)有50人回復(fù)
- 關(guān)于分子間范德華作用的隨想
已經(jīng)有68人回復(fù)
- 【求助】castep計(jì)算電荷布局?jǐn)?shù)分析
已經(jīng)有9人回復(fù)
- 【zhou2009個(gè)人文集】關(guān)于電荷計(jì)算的第二組分子
已經(jīng)有4人回復(fù)
- 【zhou2009個(gè)人文集】關(guān)于HF、H2O的電荷計(jì)算
已經(jīng)有15人回復(fù)
- 【sobereva個(gè)人文集】電子定域性的圖形分析
已經(jīng)有56人回復(fù)
sobereva
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
- QC強(qiáng)帖: 17
- 應(yīng)助: 226 (大學(xué)生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時(shí)
- 蟲(chóng)號(hào): 281271
- 注冊(cè): 2006-09-23
- 專業(yè): 宇宙學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
gmy1990: 金幣+5 2013-04-01 20:13:10
klaus1987: 金幣+10, ★★★★★最佳答案 2013-04-02 09:27:15
感謝參與，應(yīng)助指數(shù) +1
gmy1990: 金幣+5 2013-04-01 20:13:10
klaus1987: 金幣+10, ★★★★★最佳答案 2013-04-02 09:27:15
|
你不需要用主功能17來(lái)做盆積分。只要你用主功能5計(jì)算完了格點(diǎn)數(shù)據(jù)，屏幕上直接就會(huì)輸出當(dāng)前格點(diǎn)數(shù)據(jù)在整個(gè)空間的積分值。 如果你按照我之前給你回復(fù)的方法，但只是將主功能4改成了主功能5，那么這個(gè)積分值就是MO2當(dāng)中O和兩個(gè)氫這兩個(gè)片段間的重疊布居值。 高斯輸出的那個(gè)矩陣是考慮了所有軌道的，而且內(nèi)容和你要分析的也不是一回事。 gsgrid早已并入Multiwfn作為主功能13了。主功能13比最后一個(gè)版本gsgrid的功能也要強(qiáng)大很多。所以沒(méi)必要學(xué)也沒(méi)必要用gsgrid。 如果你只想要得到某個(gè)MO當(dāng)中某兩個(gè)原子間的重疊布居值，那么最簡(jiǎn)單的方法就是載入fch文件，然后選7，再選5，之后選2。這時(shí)會(huì)輸出gross atomic population matrix，這和高斯輸出的Condensed to atoms矩陣是一致的，但這是考慮了所有軌道的結(jié)果。接下來(lái)選y，就會(huì)將每個(gè)軌道的貢獻(xiàn)輸出到當(dāng)前目錄的groatmdcp.txt當(dāng)中。從中可以找到例如 ************ 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當(dāng)中1和3號(hào)原子的重疊布居，那么就是0.12815040乘以2。 Multiwfn也可以獲得特定MO當(dāng)中基函數(shù)之間的重疊布居。對(duì)于Mulliken布居的分析，Multiwfn是現(xiàn)有任何程序中最強(qiáng)大的。 |
2樓2013-04-01 18:07:17
3樓2013-04-02 09:27:01
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 klaus1987 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚(yú) 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 材料專業(yè)求調(diào)劑 +3 | hanamiko 2026-03-18 | 3/150 |
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +3 | 蔡大美女 2026-03-13 | 3/150 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研]
 機(jī)械專碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|