| 查看: 3026 | 回复: 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”是指什么呢，应该怎么回答呢？ 多谢多谢 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学 |
» 猜你喜欢
291 求调剂
已经有8人回复
-
化学工程085602 305分求调剂
已经有11人回复
-
材料学硕297已过四六级求调剂推荐
已经有12人回复
-
一志愿北京化工大学材料与化工 264分各科过A区国家线
已经有6人回复
-
085600 材料与化工 329分求调剂
已经有9人回复
-
一志愿厦门大学化学学硕307求调剂
已经有3人回复
-
315分求调剂
已经有4人回复
-
化学调剂一志愿上海交通大学336分-本科上海211
已经有4人回复
-
总分293求调剂
已经有7人回复
-
环境专硕324分求调剂推荐
已经有3人回复
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
