 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 7227 | 回復(fù): 29 | |||||||||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ，點(diǎn)擊這里進(jìn)行查看 | |||||||||
pigrass專家顧問(wèn) (著名寫手)
伊卡斯特的小豬
|
[交流]
CH4在MOF-5中的自擴(kuò)散系數(shù)計(jì)算 已有12人參與
|
||||||||
|
【本帖重點(diǎn)是MS的操作】 這個(gè)帖子主要是講自擴(kuò)散系數(shù)的MS操作，主要涉及forcite模塊、adsorption locator模塊。限于知識(shí)和硬件水平，帖子內(nèi)計(jì)算得到的自擴(kuò)散系數(shù)無(wú)參考意義。帖子最后我會(huì)單獨(dú)寫寫自己關(guān)于提高結(jié)果準(zhǔn)確性的理解。 因?yàn)槲易约簩?duì)MOF材料不是很熟悉，如果關(guān)于MOF的說(shuō)法有錯(cuò)誤的地方希望大家看到了指出來(lái)。 CH4在MOF-5中的自擴(kuò)散系數(shù)計(jì)算 1、MOF-5建模 2、MOF-5吸附CH4后的模型 3、forcite動(dòng)力學(xué)計(jì)算 4、自擴(kuò)散系數(shù)分析 1、MOF-5模型來(lái)自論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=1785384 選的是沒有化學(xué)鍵的一個(gè)模型，導(dǎo)入MS 用calculate bonds自動(dòng)生成鍵【這個(gè)功能平時(shí)不常用，容易出現(xiàn)奇怪不合理成鍵，不過(guò)這次挺好用】 成鍵后觀察結(jié)構(gòu)，注意苯環(huán)的畫法、Zn的4配位、對(duì)苯二甲酸中-COO-的連線方式 makeP1取消結(jié)構(gòu)對(duì)稱性【否則下步的CH4會(huì)按對(duì)稱性加入而不是隨機(jī)加入】 【不知道大家對(duì)這個(gè)成鍵有沒有意見，對(duì)苯二甲酸部分，這樣畫是不是有問(wèn)題，常規(guī)畫法是O-C=O吧，大家可以試試看按常規(guī)畫法能不能畫出合理結(jié)構(gòu)】 新建CH4的結(jié)構(gòu)文件，這個(gè)應(yīng)該不用多說(shuō)吧 2、激活MOF-5的結(jié)構(gòu)文件，打開adsorption locator task沒得挑，quality為medium adsorbate選擇CH4結(jié)構(gòu)文件，loading為10表示吸附10個(gè)CH4分子 energytab：forcefield選dreiding，旁邊more，取消勾選calculate automatically，點(diǎn)擊calculate，沒有錯(cuò)誤表示dreiding力場(chǎng)可用，勾選calculate automatically。charge use current，quality medium，summation全選atom based location、properties全部默認(rèn) run 生成的std文件有10個(gè)可用的結(jié)構(gòu)，選擇一個(gè)結(jié)構(gòu)打開 ctrl+A，ctrl+C，新建一個(gè).xsd，ctrl+V【或者直接用目錄下生成的.xsd文件】 單擊MOF-5上的一個(gè)原子，右鍵select fragment，display style-none隱藏MOF-5 ctrl連續(xù)選擇所有CH4分子，edit——edit sets——new——命名CH4 空白處右鍵display style-line，恢復(fù)MOF-5的顯示 3、forcite Geometry optimization more-勾選optimize cell；energy-dreiding，calculate automatically，use current，atom based 跟之前設(shè)置都一樣。優(yōu)化晶胞。 Dynamics more dynamics-NVT系綜，298K，time step0.01fs，total simulation time5.0ps，frame output every每500步輸出一個(gè)結(jié)構(gòu) thermostat-velocity scale，波動(dòng)1.0K advanced-energy deviation 默認(rèn)5w【time step0.01fs的設(shè)置不合理容易出現(xiàn)deviation超過(guò)w5的錯(cuò)誤】 energy 同之前的設(shè)置 run 4、打開.xtd文件，包含了整個(gè)dynamic過(guò)程每500steps輸出的frame forcite——analysis——mean square displacement，sets為CH4 origin step設(shè)定最小值，length設(shè)定最大值 analyse 生成MSD.xcd和MSD.xtd，觀察xcd文件，選擇有效數(shù)據(jù)范圍，對(duì)上述的origin、length重新設(shè)置重新analyze 打開新生成的MSD.xcd選擇兩列數(shù)據(jù) tools——plot graph——scatter 2D，show best fit line——plot 生成的Scatter plot.xcd，有斜率，自擴(kuò)散系數(shù)=斜率/6 【該步可以用其他數(shù)據(jù)軟件處理，只需復(fù)制所有數(shù)據(jù)到其他軟件即可】 我計(jì)算的自擴(kuò)散系數(shù)只有3.5/6，遠(yuǎn)遠(yuǎn)小于文獻(xiàn)值http://www.docin.com/p-462693767.html —————————————————————————————————— 之前說(shuō)過(guò)，限于知識(shí)和硬件，僅能做MS操作的演示，本帖的最后結(jié)果沒有參考意義 下面關(guān)于如何盡量做到合理數(shù)據(jù)寫寫我自己的理解 1、力場(chǎng)，MS自帶的力場(chǎng)中只是能夠自動(dòng)分配MOF-5力場(chǎng)的就沒幾個(gè) 例如Zn的力場(chǎng)，dreiding分配的是zinc；UFF分配的是zinc，tetrahedral，+2 oxidation state，可以說(shuō)通用力場(chǎng)只做通用計(jì)算，并不專門用于MOF，計(jì)算的數(shù)值可能會(huì)差點(diǎn)。MOF好用的力場(chǎng)有OPLS、Trappe等，可惜MS都不能用 http://www.gaoyang168.com/html/201103/3000701.html http://www.gaoyang168.com/html/200905/1311556.html【介紹了如何修改力場(chǎng)，可以試試用需要的LJ勢(shì)函數(shù)修改力場(chǎng)文件】 dynamic過(guò)程應(yīng)該也可以用第一性原理或者量子化學(xué)計(jì)算吧，結(jié)果可能精確點(diǎn)，時(shí)間么，大家都懂的，有服務(wù)器會(huì)好點(diǎn)吧 2、dynamic的時(shí)間，依舊來(lái)源文獻(xiàn)http://www.docin.com/p-462693767.html。 文獻(xiàn)里做了6ns，前面2.5ns用于平衡結(jié)構(gòu)，后面3.5ns的數(shù)據(jù)用于計(jì)算MSD和自擴(kuò)散系數(shù)。我這里的5ps，打開energy.xcd看波動(dòng)太大，完全沒有平衡的意思 最后附上doc版和MOF-5的未成鍵cif |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : CH4在MOF-5中的自擴(kuò)散系數(shù)計(jì)算.doc
- 附件 2 : MOF.cif
2014-04-27 21:08:44, 913.5 K
2014-04-27 21:08:44, 7.38 K
» 收錄本帖的淘帖專輯推薦
 多孔材料 | MOFs之第一定理與軟件 | 動(dòng)力學(xué)或熱力學(xué)計(jì)算 | 擴(kuò)散系數(shù) |
» 本帖已獲得的紅花（最新10朵）
» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有9人回復(fù)
-
293求調(diào)劑
已經(jīng)有14人回復(fù)
-
0703化學(xué)336分求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有9人回復(fù)
-
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請(qǐng)
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- MS中sorption模塊使用自定義力場(chǎng)出現(xiàn)的問(wèn)題
已經(jīng)有26人回復(fù)
- 計(jì)算MOF的摩爾磁化率中分子的量是多少？
已經(jīng)有8人回復(fù)
- 求助“CH4 四步反應(yīng)簡(jiǎn)化機(jī)理”計(jì)算甲烷燃燒反應(yīng)
已經(jīng)有3人回復(fù)
- CH4在Ir（111）表面的解離吸附過(guò)渡態(tài)計(jì)算：Henkelman的PRL文章
已經(jīng)有10人回復(fù)
- 【討論】MS中sorption模塊的MC計(jì)算
已經(jīng)有24人回復(fù)
pigrass
專家顧問(wèn) (著名寫手)
伊卡斯特的小豬
-
專家經(jīng)驗(yàn): +503 - 模擬EPI: 4
- 應(yīng)助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時(shí)
- 蟲號(hào): 1463681
- 注冊(cè): 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
- 管轄: 分子模擬
21樓2014-11-06 19:36:22
lusostty
金蟲 (小有名氣)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1430.2
- 紅花: 3
- 帖子: 141
- 在線: 84.5小時(shí)
- 蟲號(hào): 1318378
- 注冊(cè): 2011-06-08
- 性別: GG
- 專業(yè): 環(huán)境化工
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
送紅花一朵
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
送紅花一朵
雪中送炭吶  我是做分子篩吸附氣體的動(dòng)力學(xué)模擬，剛接觸，看了樓主的分享，覺得受益匪淺！不過(guò)我在模擬過(guò)程中遇到幾個(gè)問(wèn)題： 1 您的第三步中，勾選optimize cell，會(huì)不會(huì)改變晶胞參數(shù)？ 2 還有MSD的分析中，origin step、length的設(shè)置問(wèn)題。如果我設(shè)置成最小值和最大值，出現(xiàn)的是一條直線；但是如果我取后面一段時(shí)間（origin step不設(shè)置成最小值1）分析，則會(huì)出現(xiàn)一條波動(dòng)很大的線。這是為什么？如何解決呢？ 3 跑完動(dòng)力學(xué)后，我看了下動(dòng)畫，發(fā)現(xiàn)小分子會(huì)跑出骨架外（如圖3），這是否合理？看到有些帖子說(shuō)改變display style，但我覺得這并未改變本質(zhì)。 謝謝 3.PNG |
2樓2014-04-28 08:50:47
pigrass
專家顧問(wèn) (著名寫手)
伊卡斯特的小豬
-
專家經(jīng)驗(yàn): +503 - 模擬EPI: 4
- 應(yīng)助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時(shí)
- 蟲號(hào): 1463681
- 注冊(cè): 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
- 管轄: 分子模擬
|
1 您的第三步中，勾選optimize cell，會(huì)不會(huì)改變晶胞參數(shù)？ 會(huì) 2 還有MSD的分析中，origin step、length的設(shè)置問(wèn)題。如果我設(shè)置成最小值和最大值，出現(xiàn)的是一條直線；但是如果我取后面一段時(shí)間（origin step不設(shè)置成最小值1）分析，則會(huì)出現(xiàn)一條波動(dòng)很大的線。這是為什么？如何解決呢？ origin取最小值，你以為看到的是最小值，你拿鼠標(biāo)滾輪放大看看圖，或者在其他數(shù)據(jù)處理軟件看。 這就是為什么dynamic時(shí)間要盡量長(zhǎng)，MSD圖的跨度也要盡量長(zhǎng)的原因吧，可以將某段大波動(dòng)的影響降到最低 我見過(guò)一張MSD圖是這樣的，你覺得斜率是取哪一段？ 3 跑完動(dòng)力學(xué)后，我看了下動(dòng)畫，發(fā)現(xiàn)小分子會(huì)跑出骨架外（如圖3），這是否合理？看到有些帖子說(shuō)改變display style，但我覺得這并未改變本質(zhì)。 晶胞沒有跑出晶胞外的說(shuō)法，晶胞是無(wú)窮大空間的簡(jiǎn)化，你看到MS里一個(gè)晶胞，它代表整個(gè)背景、無(wú)窮景深處都是這種結(jié)構(gòu)。小分子“跑”出骨架外？你以為旁邊就沒骨架了？ |
» 本帖已獲得的紅花（最新10朵）
3樓2014-04-28 09:04:32
故事是這樣的
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 125
- 帖子: 90
- 在線: 74.1小時(shí)
- 蟲號(hào): 2533243
- 注冊(cè): 2013-07-05
4樓2014-04-28 09:18:49
pigrass
專家顧問(wèn) (著名寫手)
伊卡斯特的小豬
-
專家經(jīng)驗(yàn): +503 - 模擬EPI: 4
- 應(yīng)助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時(shí)
- 蟲號(hào): 1463681
- 注冊(cè): 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
- 管轄: 分子模擬
5樓2014-04-28 09:26:10
lusostty
金蟲 (小有名氣)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1430.2
- 紅花: 3
- 帖子: 141
- 在線: 84.5小時(shí)
- 蟲號(hào): 1318378
- 注冊(cè): 2011-06-08
- 性別: GG
- 專業(yè): 環(huán)境化工
6樓2014-04-28 09:28:43
pigrass
專家顧問(wèn) (著名寫手)
伊卡斯特的小豬
-
專家經(jīng)驗(yàn): +503 - 模擬EPI: 4
- 應(yīng)助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時(shí)
- 蟲號(hào): 1463681
- 注冊(cè): 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
- 管轄: 分子模擬
7樓2014-04-28 10:38:19
8樓2014-04-29 00:17:24
doz0803
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 843.5
- 帖子: 13
- 在線: 35.4小時(shí)
- 蟲號(hào): 2882160
- 注冊(cè): 2013-12-19
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
9樓2014-08-07 20:49:36
pigrass
專家顧問(wèn) (著名寫手)
伊卡斯特的小豬
-
專家經(jīng)驗(yàn): +503 - 模擬EPI: 4
- 應(yīng)助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時(shí)
- 蟲號(hào): 1463681
- 注冊(cè): 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
- 管轄: 分子模擬
10樓2014-08-08 11:25:12
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 pigrass 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +6 | yy要上岸呀 2026-03-17 | 6/300 |
|
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考博] 2026年博士申請(qǐng) +3 | QwQwQW10 2026-03-11 | 3/150 |
|