| 查看: 612 | 回复: 0 | ||
[求助]
pmf结果中怎么处理VDW和静电
|
|
请问下面的每一步怎么做啊，实在是搞不懂，谢谢 看文章像是用adaptive biasing force (ABF)算的PMF。我查了一下，ABF方法的作者Chipot以前在NAMD mailing list 里回复了该怎样对结果进行处理来分解PMF，具体在这里 http://www.ks.uiuc.edu/Research/ ... 2009-2010/1566.html 他说 (1) loop over your trajectory and calculate the pair interactions that you are interested in, e.g. solute 1 and solute 2. (2) collect the forces at play, van der Waals and electrostatics. (3) project these forces onto your reaction coordinate, which for each configuration, should be determined independently. (4) average the projected forces. 翻译成中文大概就是 1. 读取轨迹文件，对于每一帧，计算粒子1和2之间相互作用力 2. 把1和2之间的相互作用力分解成vdw和静电两部分 3. 把计算的力投影到反应坐标上。这句话的实际意思是，你会把 反应坐标分成多个bin（比如在0到10之间，每隔0.1设一个bin），假如 你读取的这帧构型反应坐标为2.0， 就把vdw和静电力的数值放到2.0 这个bin里 4. 如此反复直到处理完这个轨迹，然后对每个bin里力的数值加和求平均，就得到 不同反应坐标值时的vdw和静电力 5. 再对这些vdw和静电力沿反应坐标积分，即得到两部分力对PMF的贡献。 我不清楚NAMD或者VMD里是否有现成的命令做这个，否则的话你得自己写个小程序来做这个。 |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有88人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有6人回复
-
散金币祈福
已经有35人回复
