| 9 | 1/1 | 返回列表 |
| 查看: 3002 | 回復(fù): 8 | ||||
[求助]
急求。PMF收斂問題求助 已有1人參與
|
| 各位蟲友好,我是用MAMD模擬生物跨膜體系,利用ABF方法算出的PMF。但是關(guān)于PMF的部分,審稿人給出這樣的意見“The PMF results do not have any errors associated and no probe of the convergence of these simulations is shown in any part of the manuscript. This analysis should be included in a new version of the paper.” 請問其中的“ the convergence of these simulations”是指什么呢,應(yīng)該怎么回答呢? 多謝多謝 |
分子動力學(xué) |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
|
再多跑一段時間估計我的時間就太緊張了,我這個體系一個PMF下來少說要一個多星期。我還想請教一個問題,就是我在NAMD官網(wǎng)的很多自由能教程里面看到說:通過各個分窗口的.count文件或者.grad文件也可以判斷PMF是否收斂。如果PMF收斂了,則各個窗口的counts should increase uniformly;蛘呤莗lot the forces stored in the *.grad files for all windows, checking for continuity across the window boundaries。但是我還是不太明白這是什么意思呢?“增加一致”和“窗口連續(xù)性”....... |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
|
我沒用過NAMD,也沒用過Adaptive Biasing Force (ABF)計算過PMF,以下僅為我個人的解釋,可能不太嚴謹。 根據(jù)我的理解,ABF跟metadynamics有些類似,也是一種增強抽樣(enhanced sampling)的方法。其實現(xiàn)方法是在模擬過程中施加biasing force,逐漸把粒子從勢能面上的坑中拉出來。當勢能面上的坑被填平后,粒子相當于在一個沿反應(yīng)坐標自由能處處相等的勢能面上運動,這時采樣可以認為收斂了,就可以根據(jù)biasing force求出PMF。 我想你關(guān)于收斂的判據(jù)可能是來自這個文檔: http://www.ks.uiuc.edu/Training/ ... nel/channel-tut.pdf 里面提到三種判斷收斂的方法 1. One way to assess convergence is to plot sequential PMFs and observe any changes as the simulation progresses 這就是我之前回帖提到的判據(jù),多跑一段時間,觀察PMF的變化。 2. You can also plot the counts, which measure the number of times each bin has been sampled. 這條也不難理解,當采樣收斂后,粒子在一個自由能處處相等的勢能面上運動,那么在反應(yīng)坐標為任何值處出現(xiàn)的概率是均等的,所以假如你把體系處于各個反應(yīng)坐標(each bin along the reaction coordinate)的頻率畫出來的話,得到的是一個均勻的分布。 3. As another way to judge convergence, you can plot the forces stored in the *.grad files for all windows, checking for continuity across the window boundaries. 這也不是太難理解。為了方便采樣我們把整個反應(yīng)坐標劃分成了幾個窗口。因為PMF是連續(xù)函數(shù),所對應(yīng)的力(PMF的導(dǎo)數(shù))應(yīng)該也是連續(xù)函數(shù),那么當采樣收斂后,這個力在每相鄰兩個窗口的交界處應(yīng)該是連續(xù)的。 感覺你對這種方法的原理還是不太了解。建議你有時間的話,閱讀一下相關(guān)的文獻,深入了解一下ABF的原理。 |


| 9 | 1/1 | 返回列表 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|