 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 16621 | 回復(fù): 44 | |||||||||||||||
| 本帖產(chǎn)生 9 個 模擬EPI ，點擊這里進行查看 | |||||||||||||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||||||||||||
[交流]
分子模擬時間越長越好嗎？
|
|||||||||||||||
| 如題，在對蛋白質(zhì)做分子動力學模擬的時候，我用Gromacs先進行了3ns的MD運算，然后又進行了2ns，對接結(jié)果發(fā)現(xiàn)后者不如前者，但通過蛋白質(zhì)評價發(fā)現(xiàn)后者的各項指標要好一些，請各位高手討論一下MD多長時間比較好呢？ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 分子模擬 | 模型 - 模擬 - 統(tǒng)計 - 熱力學 | MDs-Gromacs | 計算 |
| gromacs | 蛋白模擬與改造 | 分子模擬 | 精華帖轉(zhuǎn)載 |
| 分子模擬 | PBS作業(yè)排隊系統(tǒng) | 計算化學 | 分子動力學模擬 |
» 本帖已獲得的紅花（最新10朵）
ChemiAndy» 猜你喜歡
285求調(diào)劑
已經(jīng)有5人回復(fù)
-
309求調(diào)劑
已經(jīng)有19人回復(fù)
-
08工科275分求調(diào)劑
已經(jīng)有12人回復(fù)
-
esi高被引論文是不是能對中標有所加分和幫助呢
已經(jīng)有3人回復(fù)
-
求調(diào)劑推薦 材料 304
已經(jīng)有19人回復(fù)
-
265求調(diào)劑
已經(jīng)有11人回復(fù)
-
292分，材料與化工，申請調(diào)劑
已經(jīng)有9人回復(fù)
-
調(diào)劑
已經(jīng)有5人回復(fù)
-
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑
已經(jīng)有12人回復(fù)
-
295材料工程專碩求調(diào)劑
已經(jīng)有13人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 對接的分子用于分子動力學模擬的問題
已經(jīng)有12人回復(fù)
- 分子模擬中MSD的問題討論
已經(jīng)有25人回復(fù)
- 擴散系數(shù)的計算及體系平衡的判斷
已經(jīng)有4人回復(fù)
- 分子動力學的時間
已經(jīng)有9人回復(fù)
- gromacs結(jié)果分析怎樣分析氫鍵能和結(jié)合能
已經(jīng)有12人回復(fù)
- LAMMPS中ReaxFF力場是怎么回事？
已經(jīng)有17人回復(fù)
- 蛋白質(zhì)構(gòu)象模擬
已經(jīng)有12人回復(fù)
- 關(guān)于材料的熱力學穩(wěn)定性 結(jié)構(gòu)穩(wěn)定性
已經(jīng)有16人回復(fù)
- 金屬功函數(shù)的測量方法有幾種
已經(jīng)有5人回復(fù)
- 什么動力學模擬軟件可以模擬化學反應(yīng)過程
已經(jīng)有13人回復(fù)
- 如何評價誤差和變化？
已經(jīng)有10人回復(fù)
- 做離子液體的分子模擬，用哪個分子模擬的軟件好呢？
已經(jīng)有3人回復(fù)
- gromacs 中g(shù)_msd問題
已經(jīng)有8人回復(fù)
- 有沒有人用Amber做過膜蛋白動力學模擬
已經(jīng)有14人回復(fù)
- 生測數(shù)據(jù)的標準差大于平均值，怎么辦？
已經(jīng)有7人回復(fù)
- AMBER中qm/mm的問題
已經(jīng)有5人回復(fù)
- 【討論】分子模擬還是同源建模
已經(jīng)有16人回復(fù)
- 【求助】msd算法
已經(jīng)有10人回復(fù)
- 【求助】NAMD和Gromacs哪個簡單？
已經(jīng)有5人回復(fù)
- 【求助】關(guān)于用分子動力學來做光譜的問題
已經(jīng)有14人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
原子層沉積裝備18817872921
+1/99
-
沈陽理工大學化學工程學碩和專碩調(diào)劑招生
+1/84
-
海南醫(yī)科大學（海南省醫(yī)學科學院）李東偉課題組招聘研究助理
+2/70
-
大連工業(yè)大學紡織與材料工程學院超臨界流體技術(shù)課題組 研究生調(diào)劑
+2/42
-
2026年武漢紡織大學材料學院鄧波教授團隊接收調(diào)劑（碩士生4名）
+1/36
-
新疆碳納米材料制備與應(yīng)用重點實驗室2026年招聘啟事
+1/34
-
南通大學國家級人才團隊招收衛(wèi)生檢驗、醫(yī)學檢驗相關(guān)方向調(diào)劑研究生
+2/24
-
廣西師范大學分子探針小組招有機化學、生物與醫(yī)藥專業(yè)碩士生2-3名
+1/19
-
天津理工大學功能晶體研究院（晶體材料全國重點實驗室）杰青團隊招收2026年博士研究生
+1/18
-
江西師大化學工程學院接收2026調(diào)劑生
+1/17
-
浙江寧波 生物與醫(yī)藥 接受生物學、園藝學、作物學、農(nóng)業(yè)調(diào)劑
+1/13
-
中國農(nóng)業(yè)大學資源與環(huán)境學院 - 招聘杰出人才（教授-專技四級）、青年研究員、優(yōu)秀人才
+1/8
-
海南大學海洋科學學院王慧團隊碩士、博士、博士后及科研骨干招聘啟事
+1/8
-
福建師范大學環(huán)境微生物技術(shù)課題組接受2026級碩士調(diào)劑。計劃招收2-3名。
+1/8
-
南京醫(yī)科大學-譚磊課題組-誠招博士后
+1/7
-
北京石油化工學院機器人應(yīng)用技術(shù)研究所2026年招收調(diào)劑研究生！
+1/5
-
金屬材料先進連接技術(shù)
+1/2
-
燕山大學人工智能學院（軟件學院）“知識發(fā)現(xiàn)與推薦”團隊碩士招生/調(diào)劑
+1/2
-
杭州電子科技大學 磁電功能材料創(chuàng)新團隊 招收碩士調(diào)劑生
+1/1
-
俄克拉荷馬大學2026秋季博士招生 (材料, 能源, 機械)
+1/1
41樓2012-05-05 02:22:40
2樓2011-06-30 17:08:06
3樓2011-06-30 17:34:47
★
奮斗1s(金幣+1):謝謝參與
奮斗1s(金幣+10): 很詳細，謝謝這位兄弟的指點呀！ 2011-07-01 08:21:38
奮斗1s(金幣+1):謝謝參與
奮斗1s(金幣+10): 很詳細，謝謝這位兄弟的指點呀！ 2011-07-01 08:21:38
|
平衡態(tài)模擬的時間長度取決于所感興趣的性質(zhì)的弛豫周期。獲得感興趣性質(zhì)的的統(tǒng)計平均值，需要在不同的弛豫周期采樣3-10次。比如，鍵的振動，振動一次10fs算一個周期，你需要采樣至少50fs。分子擴散的數(shù)量級在ns，中等肽鏈loop-closing需要10ns數(shù)量級，形成alpha螺旋的弛豫周期在200ns量級，beta結(jié)構(gòu)在1-10微秒，等等。 弛豫周期的長度取決于某種運動的持續(xù)多久以后就與先前的狀態(tài)無關(guān)了，或者忘記了。所有的運動都有多個狀態(tài)，這個運動所要穿越的自由能壘，決定了該運動的弛豫周期。這就是為什么我們可以模擬一個過程（如蛋白質(zhì)折疊）的平均特征時間，根據(jù)阿倫尼烏斯公式，反推自由能壘。 然而，平均特征時間為200ns的過程，并不意味著你跑200ns就一定能模擬到該過程的發(fā)生。它能否向你要的最終態(tài)跑去，也是個概率問題�？赡芎芸炫苓^去，也可能很慢，甚至完全跑不過去，即使能壘不高。因為你只有少數(shù)幾個分子阿。而實驗觀測的是10^23數(shù)量級的分子的統(tǒng)計結(jié)果。不過一般來說，采樣時間取特征時間的10倍以上，都可能獲得感興趣的性質(zhì)的良好的平均值。 那么，如何知道一個過程或某個性質(zhì)的平均馳豫時間呢？以體積弛豫為例，如果你plot整個體積數(shù)據(jù)隨時間變化的曲線，你會看到體積在平均值上下浮動。它受到盒子內(nèi)分子相互作用，振動，擴散的內(nèi)在驅(qū)動，有多種頻率的波動。瞬時的波動是受分子振動驅(qū)動的的，而緩慢的，長周期的上下浮動是受到分子聚合，擴散驅(qū)動的，但是不容易看出來。怎么估計這種緩慢浮動的周期呢？有幾種方法估計： 一是計算相關(guān)函數(shù)，看看相關(guān)函數(shù)中，從任一時刻開始，體積數(shù)據(jù)之間需要多久時間就不再有相互的關(guān)聯(lián)了(相關(guān)函數(shù)值趨近于0); 這段時間長度，就是特征弛豫時間；該方法適合與考察擴散運動的時間特征，計算自擴散系數(shù)，黏度等。 第二種方法是block avergage方法，看看多大的block計算得到的均方偏差RMSD與你取的block size無關(guān)。做法是把你的trajectory分成很多等長的小段，每一段稱為block， 每段長度稱為block size， 有的稱window size（e.g., 5ps），對于每一種block size的取法（比如，5ps, 10ps, 15ps, 20ps....），你計算每個block的平均值average，根據(jù)所有的block average你可以得到它們的均方偏差RMSD，Plot你的RMSD隨block size變化的曲線，你會發(fā)現(xiàn)，曲線逐漸升高，但在某個size(比如20ps)之后成為平臺區(qū)，與window size無關(guān)了（再之后會劇烈升降，進入混沌狀態(tài)），那么該性質(zhì)的特征弛豫時間為20ps。而20ps block得到的RMSD就是你對該過程或性質(zhì)的相對誤差。（很多人用全部數(shù)據(jù)的MSD做為誤差，是不嚴密的。因為所謂獨立測量，必須是測量結(jié)果之間沒有相關(guān)性。我們MD模擬和實驗測量一樣，也要從不相關(guān)的各個block中間取值阿兄弟們） 第三種方法，計算波動譜。我們知道振動有振動譜，那波動自然也有波動譜。比如體積波動，將體積數(shù)據(jù)進行自相關(guān)計算，并做傅立葉變化，可以得到一個譜圖。上面有一系列峰，在不同的頻率上。不同頻率即反映了波動的時間長短，而峰面積顯示了與該頻率的運動的相關(guān)性。面積最大的峰自然就是其特征峰了。根據(jù)其頻率可以計算其運動的特征時間。（紅外光譜就是這么計算的，它保留整個盒子的dipole，然后分解為很多不同頻率的運動，不同頻率的運動對應(yīng)不同模式的振動，它們對盒子的總dipole的變化有貢獻，傅立葉一分解就出來了）所以，只要是某種運動對某種性質(zhì)有貢獻，都可以通過對該性質(zhì)的自相關(guān)函數(shù)-傅立葉變化分解出其特征頻率，繼而得到其特征時間。 PS，上面我用到馳豫時間relaxation time這個概念，它其實指的是體系受到特定能量的激發(fā)以后，逐漸衰退decay到正常狀態(tài)所需的時間。因為這個時間是與受到激發(fā)的那種運動的運動周期直接相對應(yīng)的，所以該體系沒有受激發(fā)的時候，那種運動的周期也可以稱relaxation time弛豫時間。 再次回到樓主提出的問題上，模擬時間是否越長越好�？梢钥隙ǖ恼f是的，但是沒有必要。且不說成本，機時，能量，地球，綠色啥的問題，單說分子動力學模擬本身，是揭示相關(guān)運動的時間特征，從而定量估算內(nèi)在驅(qū)動的能量特征。 搞到蛋疼才夠酷阿兄弟  。[ Last edited by ChemiAndy on 2011-6-30 at 14:19 ] 需要說明的是，上面的討論均是指采樣階段，隱含著模擬已達到平衡的要求。而審稿人要求模擬時間加長的，一般是因為擔心體系并沒有達到平衡，或者說擔心體系停留在一個局部最小的亞平衡態(tài)。亞平衡態(tài)體系可能在幾百個納秒之后突然跳到另外一個亞平衡態(tài)。體系在多個亞平衡態(tài)之間跳躍是正常的。這樣的多態(tài)體系很難定義其平衡態(tài)，因為很難預(yù)知究竟有多少個亞態(tài)，也不知道亞態(tài)之間的能壘有多高，所以唯一的辦法就是盡可能的延長模擬時間。關(guān)于多態(tài)體系平衡的問題，囿于計算能力和研究周期限制，文獻中較少涉及。 |
8樓2011-07-01 02:24:58
簡單回復(fù)
zhang73110711樓
2011-07-01 07:56
回復(fù)
奮斗1s(金幣+1):謝謝參與
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學308分調(diào)劑 +12 | 你好明天你好 2026-03-30 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 342求調(diào)劑 +10 | Mary Keen 2026-03-28 | 10/500 |
|
|
|
[考研] 找調(diào)劑 +4 | 楚喬喬 2026-04-01 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +3 | 汪�。�！ 2026-03-31 | 6/300 |
|
|
|
[考研] 一志愿華南師范361分，化學求調(diào)劑 +4 | Nicole88888 2026-04-01 | 4/200 |
|
|
|
[考研] 085900土木水利336分求調(diào)劑 +3 | Zhangjiangj 2026-03-31 | 5/250 |
|
|
|
[考研] 311（085601）求調(diào)劑 +12 | liziyeyeye 2026-03-28 | 13/650 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-29 | 7/350 |
|
|
|
[考研] 一志愿西交大080500材料學碩349 +6 | jqx1258 2026-03-31 | 7/350 |
|
|
|
[考研] 070300化學354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
|
[考研] 08開頭看過來�。。� +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +10 | 陳晨79 2026-03-30 | 10/500 |
|
|
|
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學，材料工程（085601） +6 | 小小機靈蟲 2026-03-29 | 12/600 |
|
|
|
[考研]
 物理學調(diào)劑
+4
|
小羊36 2026-03-30 | 4/200 |
|
|
|
[考研] 085601一志愿西北工業(yè)大學初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研]
085602化學工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|