 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 2708 | 回復(fù): 3 | ||
xiaomeichen金蟲 (小有名氣)
|
[求助]
分享gromac的基本步驟，有木有好心人告知MC模擬 用towhee基本步驟
|
|
偶然看到這個，以下是進(jìn)行分子動力學(xué)模擬的最基本步驟，來自Gromacs維基。具體的步驟可能會因為模擬目標(biāo)的不同而存在差異，這里只做基本的概述。 1. 十分清楚要什么樣的模擬系統(tǒng)性質(zhì)和現(xiàn)象，明白模擬的目標(biāo)。 2. 選擇適當(dāng)?shù)墓ぞ咭詫崿F(xiàn)模擬目標(biāo)。這需要熟知其他研究人員對相似模擬系統(tǒng)的模擬方法等資料，要好好讀別人的論文，主要包括： -模擬使用的軟件，當(dāng)然還要考慮使用的力場，有時因為力場關(guān)系，需要更改使用的軟件。 -力場描述了系統(tǒng)中原子的性質(zhì)和相互作用，選擇一個適合自己研究目標(biāo)的力場至關(guān)重要，還是需要多看論文。根據(jù)sensenbobo血淚史和別人的經(jīng)驗，模擬核酸的話，可以采用amber力場，水分子可以使用TIP3P；模擬蛋白質(zhì)的話，可以采用OPLS-AA，水分子使用TIP4P（重申一下，我本人沒有跑過測試）。 3. 獲取模擬分子的坐標(biāo)文件（PDB），或者自己生成一個。自己生成很好啊，強(qiáng)烈鼓勵，但是要十分明白有一定的危險性，要明白每一步模擬系統(tǒng)的現(xiàn)象是為什么。 4. 根據(jù)模擬軟件，生成模擬系統(tǒng)初步坐標(biāo)文件。 5. 獲取或者生成系統(tǒng)的拓?fù)湮募热鏕romacs的pdb2gmx命令，或者使用網(wǎng)絡(luò)服務(wù)器PRODRG（google一下這幾個字母），amber的leap命令和antechamber，NAMD的psfgen或者VMD。請參考相關(guān)手冊。 6. 給模擬系統(tǒng)添加一個盒子（盒子的翻譯方法我十分不喜歡，箱子顯得太大，邊界顯得太難理解，但是盒子聽起來想到蛋糕，我們是嚴(yán)肅的科學(xué)家。娃哈哈。）在盒子中添加水分子，添加離子如NA和CL之類（添油加醋，生活就是這樣進(jìn)行的）。添加任何東西到你的系統(tǒng)中，就要相應(yīng)的修改一下拓?fù)湮募�，以保持與系統(tǒng)對應(yīng)，不然會出錯哦。 7. 進(jìn)行能量最優(yōu)話。這樣做是消去系統(tǒng)中原子之間距離太近造成的高能量。必要的時候，分子可以現(xiàn)在真空中能量最優(yōu)話，然后再做溶劑環(huán)境中能量最優(yōu)話。 8. 選擇恰當(dāng)?shù)哪M參數(shù)模擬平衡模擬系統(tǒng)。需要慎重考慮力場的作用，必要的時候保持系統(tǒng)體積溫度不變（NVT系綜）進(jìn)行位置限制性平衡系統(tǒng)，然后在保持系統(tǒng)溫度壓強(qiáng)不變（NPT系綜）修正系統(tǒng)的密度。然后再選擇正確的系綜（NPT，NVT，NET，BT？？）進(jìn)行數(shù)據(jù)收集。（真煩！） 9. 跑足夠長時間的模擬，使系統(tǒng)達(dá)到理想的平衡。使用模擬軟件的工具或者其他軟件觀察系統(tǒng)的性質(zhì)，如溫度，能量瞪瞪。 10. 使用適當(dāng)?shù)膮?shù)進(jìn)行數(shù)據(jù)采集模擬，注意模擬銜接過程中各個系統(tǒng)性質(zhì)保持不變，如果原子速度之類。這里還要考慮力場的作用，同時要考慮怎么測量描述模擬目標(biāo)。 11. 模擬足夠長的時間使模擬目標(biāo)明顯出現(xiàn)（可能也出現(xiàn)跟你期望相反的結(jié)果，出現(xiàn)率應(yīng)該百分之五十，加上個人感情因素，可能會達(dá)到百分之八十。保重�。�。 12. 分析模擬結(jié)果，解釋得到結(jié)果。要好好利用視圖工具，同時要十分注意漂亮的三維圖說明不了什么問題。圖表，比較，統(tǒng)計，一個都不能少。 就這樣吧。  分割線有沒有好心人高人 高尚的厲害的好心人幫我說一說蒙特卡洛模擬 MC模擬 最好是用towhee軟件模擬的基本步驟，參照上面的那樣，我覺得他說的很清楚 謝謝謝謝~~~ |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
aken5555» 猜你喜歡
316求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿211 初試270分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
350求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有4人回復(fù)
-
接收2026碩士調(diào)劑(學(xué)碩+專碩)
已經(jīng)有6人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有5人回復(fù)
-
328求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助懂醫(yī)學(xué)的好心人
已經(jīng)有11人回復(fù)
- 大家好，有好心人幫我計算四個分子的homo和lumo軌道能級圖嗎？急需，可以支付適當(dāng)費(fèi)用
已經(jīng)有4人回復(fù)
- MATERIAL STUDIO畫高分子晶體（聚乙烯PE）
已經(jīng)有10人回復(fù)
- 如何用擬薄水鋁石制備鋁溶膠？？求具體步驟~~
已經(jīng)有4人回復(fù)
- udf 重力方向變化
已經(jīng)有17人回復(fù)
- 使用towhee進(jìn)行GEMC模擬相平衡的MC運(yùn)動定義問題..
已經(jīng)有15人回復(fù)
- 【求助】英語太難了，哪位好心人能幫我翻譯下下面關(guān)于旗袍的這段中文，翻譯成英文，謝
已經(jīng)有5人回復(fù)
- 【求助】哪位好心人幫我查查該物質(zhì)的單晶結(jié)構(gòu)有發(fā)表沒？十分感謝
已經(jīng)有7人回復(fù)
cuiyoutian88
金蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 890.1
- 散金: 27
- 帖子: 262
- 在線: 198.5小時
- 蟲號: 1193813
- 注冊: 2011-01-21
- 性別: GG
- 專業(yè): 生物化學(xué)
2樓2012-02-05 11:10:55
五絕神鷹
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 240.9
- 散金: 5
- 帖子: 54
- 在線: 22.2小時
- 蟲號: 3913553
- 注冊: 2015-06-07
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
3樓2016-12-07 17:26:55
aken5555
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3034.1
- 散金: 200
- 紅花: 1
- 帖子: 387
- 在線: 57.4小時
- 蟲號: 345700
- 注冊: 2007-04-14
- 性別: GG
- 專業(yè): 資源化工
4樓2020-02-07 16:32:20
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xiaomeichen 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 310求調(diào)劑 +4 | baibai1314 2026-03-16 | 4/200 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-16 | 3/150 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|