| 4 | 1/1 | 返回列表 |
| 查看: 3117 | 回復(fù): 3 | |||
[交流]
分子動力學(xué)(MD)模擬,模擬體系構(gòu)建經(jīng)驗總結(jié) 已有3人參與
|
|
在一個完整的分子動力學(xué)模擬中,一般包括以下幾個步驟: 1.選擇將要使用的力場,并根據(jù)模擬體系確定力場參數(shù),構(gòu)建力場文件; 2.產(chǎn)生初始構(gòu)型,搭建模擬體系; 3.模擬退火(simulated annealing); 4.平衡體系; 5.模擬數(shù)據(jù)采樣; 6.數(shù)據(jù)分析、處理。 產(chǎn)生初始構(gòu)型經(jīng)驗總結(jié) 1.產(chǎn)生初試構(gòu)型的一般步驟是: 1) 給出分子內(nèi)原子坐標(biāo); 2) 按照某種分布將分子安放到空間不同位置(可以按照某種晶格排列,對于液體、氣體,亦可隨排列); 2.分子內(nèi)原子位置排布不能偏離平衡位置太遠,否則,可能會因為能量或力發(fā)散,導(dǎo)致模擬出錯; 3.如果分子比較復(fù)雜,可以先跑一個幾皮秒的單分子模擬,再取其構(gòu)型作為平衡的分子構(gòu)型; 4.分子間位置可以偏離平衡位置,但應(yīng)注意避免原子位置重疊的情況,這樣同樣會導(dǎo)致能量發(fā)散的問題,所以,初試構(gòu)型亦可讓分子距離遠一點; 5.不要讓原子超出你給定的盒子大。 模擬退火經(jīng)驗總結(jié) 1.模擬退火過程需要根據(jù)你的要求和體系的性質(zhì)來設(shè)計,目的是加速體系達到平衡; 2.如果需要達到較好的平衡效果,則可以選取更為細致的溫度梯度,使體系溫度緩慢降低;另外,為了提高效率,高溫的模擬時間可以較短,溫度越低,所需的平衡時間越長; 3.如果進行高溫退火(沸點以上),則需要進行NVT模擬,這時第一步是確定,體系的盒子大; 4.盒子大小可以從實驗中得到;如果缺少實驗數(shù)據(jù),則需要通過NPT模擬來計算,雖然體積的平衡較快,但由于初始構(gòu)型偏離平衡態(tài)較遠,這一步定出的體積可能不準(zhǔn);可以待模擬退火結(jié)束后,在進行NVT模擬(數(shù)據(jù)采樣)前,再通過NPT模擬準(zhǔn)確的定出平衡態(tài)盒子大。 5.初始構(gòu)型給出的體系大小往往大于平衡態(tài)的體系大小,這時可以跑一個高壓(eg.10000 atm)的NPT模擬,讓體系快速壓縮; 6.在平衡過程(退火)中,建議選取berendsen熱耦,不建議選擇Nose-Hoover熱耦,后則較為細致,在遠離平衡時,可能會出現(xiàn)較大的溫度漲落; |
新蟲 (初入文壇)
新蟲 (初入文壇)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 學(xué)碩材料275調(diào)劑 +7 | 路三三 2026-03-03 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 295求調(diào)劑 +5 | 等春來, 2026-03-04 | 5/250 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考研] 085601 材料工程 320 +3 | 和樂瑤 2026-03-03 | 3/150 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 中國林科院林化所(南京)2026年招收化學(xué)/材料/環(huán)境工程等背景碩士研究生3名 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[考研] 成績276,專業(yè)代碼0856求調(diào)劑 +7 | 小陳朵 2026-03-03 | 7/350 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
[考研] 0856調(diào)劑 +10 | 劉夢微 2026-02-28 | 10/500 |
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
[基金申請]
|
Doma 2026-03-01 | 7/350 |
|
|
[考研] 調(diào)劑 +3 | 簡木ChuFront 2026-02-28 | 3/150 |
|
|
[高分子] 求環(huán)氧樹脂研發(fā)1名 +3 | 孫xc 2026-02-25 | 11/550 |
|