| 查看: 2866 | 回复: 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”是指什么呢，应该怎么回答呢？ 多谢多谢 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学 |
» 猜你喜欢
心脉受损
已经有8人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
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