| 查看: 513 | 回复: 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里是否有现成的命令做这个，否则的话你得自己写个小程序来做这个。 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有226人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复