24小时热门版块排行榜    

Znn3bq.jpeg
查看: 557  |  回复: 0

晨曦1981225

新虫 (初入文坛)

[求助] 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里是否有现成的命令做这个,否则的话你得自己写个小程序来做这个。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 晨曦1981225 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 化工调剂求导师收留!一志愿失利,踏实肯干,有植物提取科研经历 +20 yzyzx 2026-04-09 21/1050 2026-04-12 00:12 by 小小小小啦啦啦
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 288求调剂,一志愿华南理工大学071005 +18 ioodiiij 2026-04-08 18/900 2026-04-11 20:25 by liyun12321
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 化学308分求调剂 +22 你好明天你好 2026-04-07 24/1200 2026-04-11 11:14 by ChemPharm
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 计算机类求调剂,22408-274分 +7 上岸de小虫 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
[考研] 085404,285分求调剂 +12 薇薇考研 2026-04-07 14/700 2026-04-09 23:10 by parmtree
[考研] 367求调剂 +10 hffQAQ 2026-04-09 10/500 2026-04-09 18:06 by lijunpoly
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +23 熊哥xtk 2026-04-07 26/1300 2026-04-09 17:17 by 18446523
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
信息提示
请填处理意见