 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 7264 | 回復: 29 | |||||||||
| 本帖產生 1 個 模擬EPI ，點擊這里進行查看 | |||||||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||||||
pigrass專家顧問 (著名寫手)
伊卡斯特的小豬
|
[交流]
CH4在MOF-5中的自擴散系數(shù)計算 已有12人參與
|
||||||||
|
【本帖重點是MS的操作】 這個帖子主要是講自擴散系數(shù)的MS操作，主要涉及forcite模塊、adsorption locator模塊。限于知識和硬件水平，帖子內計算得到的自擴散系數(shù)無參考意義。帖子最后我會單獨寫寫自己關于提高結果準確性的理解。 因為我自己對MOF材料不是很熟悉，如果關于MOF的說法有錯誤的地方希望大家看到了指出來。 CH4在MOF-5中的自擴散系數(shù)計算 1、MOF-5建模 2、MOF-5吸附CH4后的模型 3、forcite動力學計算 4、自擴散系數(shù)分析 1、MOF-5模型來自論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=1785384 選的是沒有化學鍵的一個模型，導入MS 用calculate bonds自動生成鍵【這個功能平時不常用，容易出現(xiàn)奇怪不合理成鍵，不過這次挺好用】 成鍵后觀察結構，注意苯環(huán)的畫法、Zn的4配位、對苯二甲酸中-COO-的連線方式 makeP1取消結構對稱性【否則下步的CH4會按對稱性加入而不是隨機加入】 【不知道大家對這個成鍵有沒有意見，對苯二甲酸部分，這樣畫是不是有問題，常規(guī)畫法是O-C=O吧，大家可以試試看按常規(guī)畫法能不能畫出合理結構】 新建CH4的結構文件，這個應該不用多說吧 2、激活MOF-5的結構文件，打開adsorption locator task沒得挑，quality為medium adsorbate選擇CH4結構文件，loading為10表示吸附10個CH4分子 energytab：forcefield選dreiding，旁邊more，取消勾選calculate automatically，點擊calculate，沒有錯誤表示dreiding力場可用，勾選calculate automatically。charge use current，quality medium，summation全選atom based location、properties全部默認 run 生成的std文件有10個可用的結構，選擇一個結構打開 ctrl+A，ctrl+C，新建一個.xsd，ctrl+V【或者直接用目錄下生成的.xsd文件】 單擊MOF-5上的一個原子，右鍵select fragment，display style-none隱藏MOF-5 ctrl連續(xù)選擇所有CH4分子，edit——edit sets——new——命名CH4 空白處右鍵display style-line，恢復MOF-5的顯示 3、forcite Geometry optimization more-勾選optimize cell；energy-dreiding，calculate automatically，use current，atom based 跟之前設置都一樣。優(yōu)化晶胞。 Dynamics more dynamics-NVT系綜，298K，time step0.01fs，total simulation time5.0ps，frame output every每500步輸出一個結構 thermostat-velocity scale，波動1.0K advanced-energy deviation 默認5w【time step0.01fs的設置不合理容易出現(xiàn)deviation超過w5的錯誤】 energy 同之前的設置 run 4、打開.xtd文件，包含了整個dynamic過程每500steps輸出的frame forcite——analysis——mean square displacement，sets為CH4 origin step設定最小值，length設定最大值 analyse 生成MSD.xcd和MSD.xtd，觀察xcd文件，選擇有效數(shù)據(jù)范圍，對上述的origin、length重新設置重新analyze 打開新生成的MSD.xcd選擇兩列數(shù)據(jù) tools——plot graph——scatter 2D，show best fit line——plot 生成的Scatter plot.xcd，有斜率，自擴散系數(shù)=斜率/6 【該步可以用其他數(shù)據(jù)軟件處理，只需復制所有數(shù)據(jù)到其他軟件即可】 我計算的自擴散系數(shù)只有3.5/6，遠遠小于文獻值http://www.docin.com/p-462693767.html —————————————————————————————————— 之前說過，限于知識和硬件，僅能做MS操作的演示，本帖的最后結果沒有參考意義 下面關于如何盡量做到合理數(shù)據(jù)寫寫我自己的理解 1、力場，MS自帶的力場中只是能夠自動分配MOF-5力場的就沒幾個 例如Zn的力場，dreiding分配的是zinc；UFF分配的是zinc，tetrahedral，+2 oxidation state，可以說通用力場只做通用計算，并不專門用于MOF，計算的數(shù)值可能會差點。MOF好用的力場有OPLS、Trappe等，可惜MS都不能用 http://www.gaoyang168.com/html/201103/3000701.html http://www.gaoyang168.com/html/200905/1311556.html【介紹了如何修改力場，可以試試用需要的LJ勢函數(shù)修改力場文件】 dynamic過程應該也可以用第一性原理或者量子化學計算吧，結果可能精確點，時間么，大家都懂的，有服務器會好點吧 2、dynamic的時間，依舊來源文獻http://www.docin.com/p-462693767.html。 文獻里做了6ns，前面2.5ns用于平衡結構，后面3.5ns的數(shù)據(jù)用于計算MSD和自擴散系數(shù)。我這里的5ps，打開energy.xcd看波動太大，完全沒有平衡的意思 最后附上doc版和MOF-5的未成鍵cif |
回復此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : CH4在MOF-5中的自擴散系數(shù)計算.doc
- 附件 2 : MOF.cif
2014-04-27 21:08:44, 913.5 K
2014-04-27 21:08:44, 7.38 K
» 收錄本帖的淘帖專輯推薦
 多孔材料 | MOFs之第一定理與軟件 | 動力學或熱力學計算 | 擴散系數(shù) |
» 本帖已獲得的紅花（最新10朵）
» 猜你喜歡
一志愿北化085600材料專碩275|有文章專利｜求調劑
已經有6人回復
-
一志愿華理，數(shù)一英一285求A區(qū)調劑
已經有12人回復
-
289求調劑
已經有11人回復
-
食品工程專碩一志愿中海洋309求調劑
已經有7人回復
-
生物學學碩，一志愿湖南大學，初試成績338
已經有6人回復
-
321求調劑
已經有7人回復
-
343求調劑
已經有4人回復
-
305求調劑
已經有3人回復
-
本科新能源科學與工程，一志愿華理能動285求調劑
已經有6人回復
-
332求調劑
已經有3人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- MS中sorption模塊使用自定義力場出現(xiàn)的問題
已經有26人回復
- 計算MOF的摩爾磁化率中分子的量是多少？
已經有8人回復
- 求助“CH4 四步反應簡化機理”計算甲烷燃燒反應
已經有3人回復
- CH4在Ir（111）表面的解離吸附過渡態(tài)計算：Henkelman的PRL文章
已經有10人回復
- 【討論】MS中sorption模塊的MC計算
已經有24人回復
doz0803
金蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 843.5
- 帖子: 13
- 在線: 35.4小時
- 蟲號: 2882160
- 注冊: 2013-12-19
- 性別: GG
- 專業(yè): 高分子材料結構與性能
9樓2014-08-07 20:49:36
lusostty
金蟲 (小有名氣)
- 應助: 6 (幼兒園)
- 金幣: 1430.2
- 紅花: 3
- 帖子: 141
- 在線: 84.5小時
- 蟲號: 1318378
- 注冊: 2011-06-08
- 性別: GG
- 專業(yè): 環(huán)境化工
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
送紅花一朵
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
送紅花一朵
雪中送炭吶  我是做分子篩吸附氣體的動力學模擬，剛接觸，看了樓主的分享，覺得受益匪淺！不過我在模擬過程中遇到幾個問題： 1 您的第三步中，勾選optimize cell，會不會改變晶胞參數(shù)？ 2 還有MSD的分析中，origin step、length的設置問題。如果我設置成最小值和最大值，出現(xiàn)的是一條直線；但是如果我取后面一段時間（origin step不設置成最小值1）分析，則會出現(xiàn)一條波動很大的線。這是為什么？如何解決呢？ 3 跑完動力學后，我看了下動畫，發(fā)現(xiàn)小分子會跑出骨架外（如圖3），這是否合理？看到有些帖子說改變display style，但我覺得這并未改變本質。 謝謝 3.PNG |
2樓2014-04-28 08:50:47
pigrass
專家顧問 (著名寫手)
伊卡斯特的小豬
-
專家經驗: +503 - 模擬EPI: 4
- 應助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時
- 蟲號: 1463681
- 注冊: 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結構與性能
- 管轄: 分子模擬
|
1 您的第三步中，勾選optimize cell，會不會改變晶胞參數(shù)？ 會 2 還有MSD的分析中，origin step、length的設置問題。如果我設置成最小值和最大值，出現(xiàn)的是一條直線；但是如果我取后面一段時間（origin step不設置成最小值1）分析，則會出現(xiàn)一條波動很大的線。這是為什么？如何解決呢？ origin取最小值，你以為看到的是最小值，你拿鼠標滾輪放大看看圖，或者在其他數(shù)據(jù)處理軟件看。 這就是為什么dynamic時間要盡量長，MSD圖的跨度也要盡量長的原因吧，可以將某段大波動的影響降到最低 我見過一張MSD圖是這樣的，你覺得斜率是取哪一段？ 3 跑完動力學后，我看了下動畫，發(fā)現(xiàn)小分子會跑出骨架外（如圖3），這是否合理？看到有些帖子說改變display style，但我覺得這并未改變本質。 晶胞沒有跑出晶胞外的說法，晶胞是無窮大空間的簡化，你看到MS里一個晶胞，它代表整個背景、無窮景深處都是這種結構。小分子“跑”出骨架外？你以為旁邊就沒骨架了？ |
» 本帖已獲得的紅花（最新10朵）
3樓2014-04-28 09:04:32
4樓2014-04-28 09:18:49
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學0703 調劑 306分 一志愿211 +4 | 26要上岸 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085600 286分 材料求調劑 +7 | 麻辣魷魚 2026-03-27 | 8/400 |
|
|
|
[考研] 083000學碩274求調劑 +7 | Li李魚 2026-03-26 | 7/350 |
|
|
|
[考研] 086000調劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 265求調劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 0703化學338求調劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 316求調劑 +5 | 江辭666 2026-03-26 | 5/250 |
|
|
|
[考研] 求調劑 一志愿 本科 北科大 化學 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研] 材料考研求調劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 086000生物與醫(yī)藥292求調劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
|
[考研] 機械學碩總分317求調劑！�。。� +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 0854AI CV方向招收調劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
|
[考研] 293求調劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 300分，材料，求調劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 材料專碩找調劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 求調劑一志愿武漢理工大學材料工程（085601） +5 | WW.' 2026-03-23 | 7/350 |
|
|
|
[考研] 284求調劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[考研] 275求調劑 +6 | shansx 2026-03-22 | 8/400 |
|