| 4 | 1/1 | 返回列表 |
| 查看: 3121 | 回復(fù): 3 | |||
可愛(ài)草莓新蟲 (初入文壇)
|
[交流]
分子動(dòng)力學(xué)(MD)模擬,模擬體系構(gòu)建經(jīng)驗(yàn)總結(jié) 已有3人參與
|
|
在一個(gè)完整的分子動(dòng)力學(xué)模擬中,一般包括以下幾個(gè)步驟: 1.選擇將要使用的力場(chǎng),并根據(jù)模擬體系確定力場(chǎng)參數(shù),構(gòu)建力場(chǎng)文件; 2.產(chǎn)生初始構(gòu)型,搭建模擬體系; 3.模擬退火(simulated annealing); 4.平衡體系; 5.模擬數(shù)據(jù)采樣; 6.數(shù)據(jù)分析、處理。 產(chǎn)生初始構(gòu)型經(jīng)驗(yàn)總結(jié) 1.產(chǎn)生初試構(gòu)型的一般步驟是: 1) 給出分子內(nèi)原子坐標(biāo); 2) 按照某種分布將分子安放到空間不同位置(可以按照某種晶格排列,對(duì)于液體、氣體,亦可隨排列); 2.分子內(nèi)原子位置排布不能偏離平衡位置太遠(yuǎn),否則,可能會(huì)因?yàn)槟芰炕蛄Πl(fā)散,導(dǎo)致模擬出錯(cuò); 3.如果分子比較復(fù)雜,可以先跑一個(gè)幾皮秒的單分子模擬,再取其構(gòu)型作為平衡的分子構(gòu)型; 4.分子間位置可以偏離平衡位置,但應(yīng)注意避免原子位置重疊的情況,這樣同樣會(huì)導(dǎo)致能量發(fā)散的問(wèn)題,所以,初試構(gòu)型亦可讓分子距離遠(yuǎn)一點(diǎn); 5.不要讓原子超出你給定的盒子大; 模擬退火經(jīng)驗(yàn)總結(jié) 1.模擬退火過(guò)程需要根據(jù)你的要求和體系的性質(zhì)來(lái)設(shè)計(jì),目的是加速體系達(dá)到平衡; 2.如果需要達(dá)到較好的平衡效果,則可以選取更為細(xì)致的溫度梯度,使體系溫度緩慢降低;另外,為了提高效率,高溫的模擬時(shí)間可以較短,溫度越低,所需的平衡時(shí)間越長(zhǎng); 3.如果進(jìn)行高溫退火(沸點(diǎn)以上),則需要進(jìn)行NVT模擬,這時(shí)第一步是確定,體系的盒子大小; 4.盒子大小可以從實(shí)驗(yàn)中得到;如果缺少實(shí)驗(yàn)數(shù)據(jù),則需要通過(guò)NPT模擬來(lái)計(jì)算,雖然體積的平衡較快,但由于初始構(gòu)型偏離平衡態(tài)較遠(yuǎn),這一步定出的體積可能不準(zhǔn);可以待模擬退火結(jié)束后,在進(jìn)行NVT模擬(數(shù)據(jù)采樣)前,再通過(guò)NPT模擬準(zhǔn)確的定出平衡態(tài)盒子大; 5.初始構(gòu)型給出的體系大小往往大于平衡態(tài)的體系大小,這時(shí)可以跑一個(gè)高壓(eg.10000 atm)的NPT模擬,讓體系快速壓縮; 6.在平衡過(guò)程(退火)中,建議選取berendsen熱耦,不建議選擇Nose-Hoover熱耦,后則較為細(xì)致,在遠(yuǎn)離平衡時(shí),可能會(huì)出現(xiàn)較大的溫度漲落; |
新蟲 (初入文壇)
新蟲 (初入文壇)
新蟲 (初入文壇)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過(guò)來(lái) +4 | 啊擺啊擺 2026-03-05 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 復(fù)試調(diào)劑 +4 | 呼呼?~+123456 2026-03-05 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 泡了個(gè)椒 2026-03-04 | 4/200 |
|
|
[考研] 江西師范大學(xué)王昆-盧章輝課題組招化學(xué)方向碩士研究生兩名 +3 | 哈哈大俠 2026-02-28 | 3/150 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-04 | 3/150 |
|
|
[考研] 0703化學(xué) 一志愿東北大學(xué) 279分求調(diào)劑 +4 | 聶樂(lè)平 2026-03-03 | 4/200 |
|
|
[考研] 本科太原理工采礦工程,求調(diào)劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考博] 26申博 +5 | north, 2026-02-28 | 5/250 |
|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
[考研] 主題 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 280分 +10 | yyqqhh 2026-03-03 | 10/500 |
|
|
[考研] 085602化學(xué)工程350,調(diào)劑,有沒(méi)有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
[考研] 26考研報(bào)考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 化工299分求調(diào)劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|