 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術(shù)科研互動平臺 » 計算模擬區(qū) » 第一性原理 » MS » 【求助】請教分子動力學模擬
| 9 | 1/1 | 返回列表 |
| 查看: 4814 | 回復: 8 | |||||
liuxianlv金蟲 (正式寫手)
|
[交流]
【求助】請教分子動力學模擬 已有5人參與
|
|
本人初學MS，本實驗室我是第一個搞這個的，純粹是摸索前進，對文獻上的很多不懂特請教各位大哥大姐，謝謝各位。 分子動力學模擬 選取Fe 晶體(001)晶面為吸附表面[18], 表面體系厚度為11 層共計1584 個Fe 原子, 體系大小為3.20 nm×3.20 nm×6.43 nm. 計算過程中“凍結(jié)”表面體系中的Fe 原子, 而吸附分子保持與金屬表面自由相互作用. 模擬中使用Materials Studio 4.0 軟件包中的COMPASS 力場對體系進行優(yōu)化, 正則系綜(NVT)[19]的MD 模擬通過Discover[20] 模塊完成. 模溫度為298 K, 采用Andersen[21]恒溫器進行溫度控制. 各分子起始速度由Maxwell-Boltzmann 分布隨機產(chǎn)生. 利用周期性邊界條件和時間平均等效于系綜平均等基本假設(shè), 運用Velocity Verlet 算法[22]求解牛頓運動方程. 范德瓦爾斯和庫侖相互作用采用Charge Group 方法[23]計算. 截斷半徑為1.5 nm, 截斷距離之外的分子間相互作用能按平均密度近似方法進行校正. 模擬中時間步長為1 fs, 模擬時間為500 ps, 每1000 步記錄一次體系的軌跡信息. 所有計算均在Pentium IV 計算機上完成 分子在金屬表面的吸附能由下列公式計算得出[10]adsorption molecule surface E ＝E ＋E －Etotal (1) 其中, Eadsorption 為吸附能; Emolecule 是自由分子的能量;Esurface 是未吸附分子時金屬表面的能量; Etotal 是包含一個分子和金屬表面體系的總能量. 通過分析體系能量和溫度隨時間演化軌跡發(fā)現(xiàn): 100 ps 時體系就開始達到平衡, 后250 ps 時體系已充分平衡, 因此選擇后250 ps進行取樣平均, 求得各分子吸附能和吸附構(gòu)型結(jié)構(gòu)參數(shù)的統(tǒng)計平均值.（文獻來源《1-(2-羥乙基)-2-烷基-咪唑啉緩蝕劑緩蝕機理的理論研究》化 學 學 報2008 年第66 卷，張 軍） 我對這段模擬過程和設(shè)置細節(jié)有些不清楚，現(xiàn)摸索如下：各位幫我看看是不是這樣的步驟？ 第一步建了表面按要求。 第二步加入表面分子。 第三步  選discover 里的setup選compass直接再點calculation（是不是不用點 和 以及其他設(shè)置）。 第四步選discover 里的Dynimics，   （原文中 各分子起始速度由Maxwell-Boltzmann 分布隨機產(chǎn)生. 利用周期性邊界條件和時間平均等效于系綜平均等基本假設(shè), 運用Velocity Verlet 算法[22]求解牛頓運動方程. 范德瓦爾斯和庫侖相互作用采用Charge Group 方法[23]計算. 截斷半徑為1.5 nm, 截斷距離之外的分子間相互作用能按平均密度近似方法進行校正.是怎么設(shè)的？） 第五步點NVT下的run. Maxwell-Boltzmann 沒找到這東西，在哪設(shè)？Velocity Verlet 也沒有， 范德瓦爾斯和庫侖相互作用采用Charge Group 方法[23]計算. 截斷半徑為1.5 nm好像在Setup里有不知是不是？vdw&coulomb是指范德瓦爾斯和庫侖相互作用嗎？以及 截斷半徑為1.5 nm是不是如圖設(shè)的：  Charge Group 方法是不是指group based |
回復此樓» 收錄本帖的淘帖專輯推薦
 晶體模型計算分析 | 分子動力學 |
» 猜你喜歡
生物學學碩，一志愿湖南大學，初試成績338
已經(jīng)有6人回復
-
321求調(diào)劑
已經(jīng)有7人回復
-
343求調(diào)劑
已經(jīng)有4人回復
-
305求調(diào)劑
已經(jīng)有3人回復
-
本科新能源科學與工程，一志愿華理能動285求調(diào)劑
已經(jīng)有6人回復
-
332求調(diào)劑
已經(jīng)有3人回復
-
學碩274求調(diào)劑
已經(jīng)有8人回復
-
壓國家一區(qū)線，求導師收留，有恩必謝！
已經(jīng)有7人回復
-
081200-11408-276學碩求調(diào)劑
已經(jīng)有4人回復
-
22408 359分調(diào)劑
已經(jīng)有4人回復
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教各位，分子動力學模擬NVE系統(tǒng)，不對速度進行標度可以嗎？
已經(jīng)有7人回復
- 分子動力學模擬長程庫侖勢的處理問題
已經(jīng)有8人回復
- 【求助】VASP做分子動力學模擬NVT系綜下為啥內(nèi)能竟然是升高的？
已經(jīng)有8人回復
- 請問各位熟悉分子動力學模擬的朋友，熟悉精通MS的分子動力學模擬需要多久？
已經(jīng)有17人回復
- 【討論】關(guān)于汽液兩相分子動力學模擬中粒子在汽液兩相盒子交換的編程問題
已經(jīng)有3人回復
- 【討論】大家好，加水后跑完分子動力學，在log中發(fā)現(xiàn)了問題
已經(jīng)有15人回復
- 【求助】分子動力學模擬尺度問題？
已經(jīng)有11人回復
- 【原創(chuàng)】分子對接，同源模建，分子動力學軟件介紹帖
已經(jīng)有35人回復
- 【求助】新手請教分子模擬問題
已經(jīng)有27人回復
2樓2010-06-02 21:23:40
liuxianlv
金蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 868
- 散金: 42
- 帖子: 883
- 在線: 61.4小時
- 蟲號: 839985
- 注冊: 2009-09-04
- 性別: GG
- 專業(yè): 理論和計算化學
|
可以再回答點嗎？（原文中 各分子起始速度由Maxwell-Boltzmann 分布隨機產(chǎn)生. 利用周期性邊界條件和時間平均等效于系綜平均等基本假設(shè), 運用Velocity Verlet 算法[22]求解牛頓運動方程. 范德瓦爾斯和庫侖相互作用采用Charge Group 方法[23]計算. 截斷半徑為1.5 nm, 截斷距離之外的分子間相互作用能按平均密度近似方法進行校正.是怎么設(shè)的？） 第五步點NVT下的run. Maxwell-Boltzmann 沒找到這東西，在哪設(shè)？Velocity Verlet 也沒有， 范德瓦爾斯和庫侖相互作用采用Charge Group 方法[23]計算. 截斷半徑為1.5 nm好像在Setup里有不知是不是？vdw&coulomb是指范德瓦爾斯和庫侖相互作用嗎？以及 截斷半徑為1.5 nm是不是如圖設(shè) |
3樓2010-06-02 21:43:27
043114076
木蟲 (正式寫手)
- 應助: 49 (小學生)
- 金幣: 3364.9
- 散金: 52
- 紅花: 23
- 帖子: 788
- 在線: 369.5小時
- 蟲號: 570592
- 注冊: 2008-06-07
- 性別: GG
- 專業(yè): 理論和計算化學
4樓2010-06-03 09:46:33
liuxianlv
金蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 868
- 散金: 42
- 帖子: 883
- 在線: 61.4小時
- 蟲號: 839985
- 注冊: 2009-09-04
- 性別: GG
- 專業(yè): 理論和計算化學
5樓2010-06-03 10:44:17
6樓2010-06-03 20:48:41
7樓2010-09-21 09:11:21
8樓2010-09-21 09:12:03
倚天艦
木蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 1487.1
- 散金: 4
- 紅花: 3
- 帖子: 82
- 在線: 106小時
- 蟲號: 858328
- 注冊: 2009-09-27
- 專業(yè): 生物化工與食品化工
9樓2011-09-30 18:22:43
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 liuxianlv 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +5 | 芝士雪baoo 2026-03-28 | 6/300 |
|
|
[考研] 070305高分子化學與物理 304分求調(diào)劑 +4 | c297914 2026-03-28 | 4/200 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +7 | 麻辣魷魚 2026-03-27 | 8/400 |
|
|
|
[考研] 一志愿南昌大學324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[有機交流]
 有機合成求助
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
|
[考研] 305分求調(diào)劑（食品工程） +5 | Sxy112 2026-03-21 | 7/350 |
|