| 查看: 3253 | 回复: 8 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
急求。PMF收敛问题求助 已有1人参与
|
||||
| 各位虫友好，我是用MAMD模拟生物跨膜体系，利用ABF方法算出的PMF。但是关于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”是指什么呢，应该怎么回答呢？ 多谢多谢 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学 |
» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有21人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有8人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
filecode
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
|
我没用过NAMD，也没用过Adaptive Biasing Force (ABF)计算过PMF，以下仅为我个人的解释,可能不太严谨。 根据我的理解，ABF跟metadynamics有些类似，也是一种增强抽样(enhanced sampling）的方法。其实现方法是在模拟过程中施加biasing force，逐渐把粒子从势能面上的坑中拉出来。当势能面上的坑被填平后，粒子相当于在一个沿反应坐标自由能处处相等的势能面上运动，这时采样可以认为收敛了，就可以根据biasing force求出PMF。 我想你关于收敛的判据可能是来自这个文档： 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 这就是我之前回帖提到的判据，多跑一段时间，观察PMF的变化。 2. You can also plot the counts, which measure the number of times each bin has been sampled. 这条也不难理解，当采样收敛后，粒子在一个自由能处处相等的势能面上运动，那么在反应坐标为任何值处出现的概率是均等的，所以假如你把体系处于各个反应坐标(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. 这也不是太难理解。为了方便采样我们把整个反应坐标划分成了几个窗口。因为PMF是连续函数，所对应的力（PMF的导数)应该也是连续函数，那么当采样收敛后，这个力在每相邻两个窗口的交界处应该是连续的。 感觉你对这种方法的原理还是不太了解。建议你有时间的话，阅读一下相关的文献，深入了解一下ABF的原理。 |
yylcx6楼2015-03-24 09:12:39
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2015-03-19 11:14:13
3楼2015-03-19 14:37:55
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2015-03-20 07:49:25
