 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 13 | 1/1 | 返回列表 |
| 查看: 3870 | 回復(fù): 12 | |||||
| 【獎勵】 本帖被評價8次，作者sobereva增加金幣 6.4 個 | |||||
[資源]
Multiwfn結(jié)合VMD繪制AIM拓撲分析圖 已有2人參與
|
|||||
|
Multiwfn結(jié)合VMD繪制AIM拓撲分析圖 實空間函數(shù)的拓撲分析是Multiwfn（http://multiwfn.codeplex.com）的重要功能之一，電子密度及其拉普拉斯函數(shù)的拓撲分析是AIM分析中的核心，而ELF、LOL的拓撲分析在分析成鍵方式、芳香性等問題上也有重要應(yīng)用。在Multiwfn中這些都可以通過主功能2來實現(xiàn)。關(guān)于實空間函數(shù)的拓撲分析的基本概念和在Multiwfn里的操作在此帖都有介紹：《使用Multiwfn做拓撲分析以及計算孤對電子角度》（http://hi.baidu.com/sobereva/item/c84fa5089156806dd55a11c3） 雖然Multiwfn的拓撲分析功能本身就能繪制出臨界點和拓撲路徑，也足夠滿足通常的需要，但是如果想得到更漂亮的效果，或者更自由地控制顯示方式，則可以將分析結(jié)果從Multiwfn導出并通過VMD程序來繪制。此文就介紹下過程。這里以苯酚的電子密度拓撲分析（或稱AIM拓撲分析）為例來說明。 首先我們還是先照常在Multiwfn里找出體系的各種臨界點并生成鍵徑。啟動Multiwfn，然后依次輸入 examples\PhenolDimer.wfn 2 //主功能2 2 //以原子核位置作為搜索臨界點的初猜 3 //以每兩個原子的中點作為搜索臨界點的初猜 0 //查看結(jié)果 從命令行界面中可以看到找出了61個臨界點，并且Poincare-Hopf關(guān)系已經(jīng)滿足了。同時從蹦出來的圖形窗口中也可以憑直覺看到該有的臨界點也都有了，說明臨界點都找全了。于是關(guān)閉窗口，然后選8，就把鍵徑都生成了，這即是(3,-3)和(3,-1)臨界點之間的拓撲路徑。雖然也可以再選9把(3,+1)和(3,+3)之間的拓撲路徑生成出來，但沒太大價值，這里就不做了。然后選0，可以看到臨界點和鍵徑都出現(xiàn)在圖中了，紫/橙/黃/綠小球分別代表核臨界點/鍵臨界點/環(huán)臨界點/籠臨界點。如果點擊窗口右側(cè)"Paths labels"復(fù)選框，還可以看到每條鍵徑的序號，如下所示  關(guān)閉圖形窗口，然后將臨界點和拓撲路徑都導出為.pdb文件。具體做法是選-4進入臨界點修改/導出界面，然后選6 Export CPs as pdb file in current folder，臨界點就被導出到當前目錄下的CPs.pdb里了，每個臨界點對應(yīng)此文件中的一個原子，并且如屏幕提示所示，不同元素的原子對應(yīng)不同類型臨界點，即碳/氮/氧/氟分別對應(yīng)(3,-3)/(3,-1)/(3,+1)/(3,+3)。選0退回到拓撲分析主菜單，再進入選項-5，即修改、顯示、導出以及計算拓撲路徑上的性質(zhì)的界面。然后選6 Export paths as paths.pdb file in current folder，這樣拓撲路徑就被導出到當前目錄下的path.pdb里了。實際上拓撲路徑就是由一連串緊挨著的點來表示的，整體來看就連成了曲線。在path.pdb里，這些構(gòu)成鍵徑的點對應(yīng)于一個個碳原子，每個原子所屬殘基編號就對應(yīng)于拓撲路徑的編號。 下面我們就要在VMD中繪圖了。VMD是個十分靈活、效果出色的免費的分子可視化軟件，下載地址見http://www.ks.uiuc.edu/Research/vmd/。這里用VMD1.9版。 啟動VMD，將CPs.pdb拖到主窗口VMD Main中，然后就會看到一副雖然有點藝術(shù)感，但是連線亂七八糟的圖，這是因為VMD自動判斷了原子間連接關(guān)系。我們選Graphics-Representation來改變顯示方式。這里讓(3,-3)臨界點都用紫色圓點顯示，就改成下圖中這樣，即Selected Atoms輸入name C，Drawing Method用Points，Coloring Method用ColorID并自己指定顏色，大小通過Size來調(diào)  然后把(3,-1)也顯示。點擊Create Rep新建顯示方式，輸入name N，然后設(shè)定好繪制方法、顏色和尺寸。然后再類似地顯示出(3,+1)和(3,+3)，它們用name O和name F選擇。最后繪制的效果如下。為了凸顯(3,-3)臨界點，它的size由默認的11改為了20。  上圖用的是黑色背景，臨界點位置看得很清楚，但是如果用白色背景（命令行窗口輸入color Display Background white），臨界點的顏色應(yīng)該用得深一些。 接下來顯示鍵徑。我們把paths.pdb拖進VMD Main窗口，進入Graphics-Representation后確認Selected Molecule已經(jīng)切換到paths.pdb了。將Drawing Method改為Points，Size改為1，顏色改為青色。這時圖像看起來效果已經(jīng)很好了。  我們還可以對圖像進行更進一步修改。比如，我們想讓兩個苯酚之間的鍵徑和分子內(nèi)的鍵徑區(qū)分開以便觀察。為此，我們先從本文的第一張圖找出分子間鍵徑的編號，可以看出編號是35 36 55 56 49 50。然后我們在把鍵徑對應(yīng)的selected atoms框的內(nèi)容由默認的all改為all not resid 35 36 55 56 49 50，這時由圖可見只剩下分子內(nèi)的鍵徑顯示在圖中了。然后點擊Create Rep，將新的顯示方式的所選原子設(shè)為resid 35 36 55 56 49 50，然后把顏色設(shè)為紅色。我們再做些調(diào)整，比如這里把分子內(nèi)的鍵徑加粗（即設(shè)大size），相應(yīng)地把(3,-3)和(3,-1)臨界點的size也加大一些，并且不讓(3,+3)顯示出來，最終的圖像如下所示  可見顯示效果非常不錯。依照這種方法，還可以繪制諸如ELF、LOL等實空間函數(shù)的拓撲路徑，操作過程一樣，這里就不累述了。 Multiwfn雖然內(nèi)建了很方便的可視化結(jié)果的功能，但為了靈活性考慮很多地方都留出了專門的選項用于導出數(shù)據(jù)，以便通過外部程序諸如VMD、sigmaplot等程序繪圖，以滿足一些有個性用戶的需求。使用Multiwfn時應(yīng)當多嘗試，對照著手冊相應(yīng)章節(jié)領(lǐng)會各個選項的含義。Multiwfn的靈活性與VMD的靈活性相結(jié)合使得繪制許多需要用戶高度自定義的圖像成為了可能，這是其它任何程序都無法代替的。其它的Multiwfn與VMD相結(jié)合進行繪圖的例子見《通過鍵級曲線和ELF/LOL/RDG等值面動畫研究化學反應(yīng)過程》（http://hi.baidu.com/sobereva/item/e839b2d2efd0aadf241f4015）和《使用Multiwfn結(jié)合VMD分析和繪制分子表面靜電勢分布》（http://hi.baidu.com/sobereva/item/692f9a83d8e9a5c8b071545d）。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 資源收集 | 量化軟件學習 | SOB 集錦 |
» 猜你喜歡
廣西師范大學化學與藥學學院楊秀林課題組2026年擬招收博士研究生1名
已經(jīng)有25人回復(fù)
-
26年儲能、電池方向博士申請
已經(jīng)有1人回復(fù)
-
物理化學論文潤色/翻譯怎么收費?
已經(jīng)有185人回復(fù)
-
求調(diào)劑，本科吉大，一志愿華工339
已經(jīng)有2人回復(fù)
-
意大利博洛尼亞大學招聘研究助理2名（2D材料/鋰離子與鈉離子電池方向）
已經(jīng)有1人回復(fù)
-
關(guān)于碳量子點抗菌的研究
已經(jīng)有15人回復(fù)
-
南大化學調(diào)劑
已經(jīng)有11人回復(fù)
-
0703調(diào)劑
已經(jīng)有2人回復(fù)
-
1
已經(jīng)有14人回復(fù)
-
0703 調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 一起北京看美女吧：量子化學波函數(shù)分析暨Multiwfn程序使用培訓班
已經(jīng)有18人回復(fù)
- 如何用Multiwfn軟件做出P DOS?
已經(jīng)有4人回復(fù)
- AIM鍵臨界點處電子密度拉普拉斯值符號判斷相互作用類型失敗原因的圖形分析
已經(jīng)有24人回復(fù)
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- Multiwfn的初步學習，越詳細越好。
已經(jīng)有6人回復(fù)
- MultiWFN想要生成用于NBO分析的cub時不能載入.37文件
已經(jīng)有3人回復(fù)
- 使用Multiwfn繪制NBO及相關(guān)軌道
已經(jīng)有30人回復(fù)
- 用Multiwfn、Gaussview、Molekel、VMD觀看龍蝦、盆景、骨盆、大腦
已經(jīng)有16人回復(fù)
- 引用了multiwfn的論文。現(xiàn)在multiwfn軟件已經(jīng)有19次引用率了呢
已經(jīng)有13人回復(fù)
- G09中計算AIM計算過渡金屬原子和氫原子間作用的問題
已經(jīng)有7人回復(fù)
- 使用Multiwfn做拓撲分析以及計算孤對電子角度
已經(jīng)有18人回復(fù)
- 回答一些關(guān)于Multiwfn的疑問以及未來Multiwfn的發(fā)展打算
已經(jīng)有37人回復(fù)
- 【sobereva個人文集】電子定域性的圖形分析
已經(jīng)有56人回復(fù)
2樓2013-11-05 09:51:34
4樓2013-11-05 15:43:33
6樓2013-11-05 23:39:40
7樓2013-11-06 00:54:26
9樓2013-11-07 05:44:17
10樓2013-11-07 08:13:09
11樓2015-03-12 11:13:12
12樓2015-03-12 11:13:18
13樓2015-03-12 11:13:37
簡單回復(fù)
2013-11-05 10:01
回復(fù)
VMD
ifmc12345樓
2013-11-05 18:52
回復(fù)
五星好評 
2013-11-06 08:39
回復(fù)
五星好評 
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 sobereva 的主題更新
| 13 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +3 | 煙汐憶海 2026-03-16 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|