计算化学公社

 找回密码 Forget password
 注册 Register

gmx hbond计算氢键时,针对氢键数目解读的疑惑

查看数: 348 | 评论数: 1 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-11-17 12:33

正文摘要:

感谢各位费心点击帖子 我在使用gmx hbond工具进行氢键分析分析时(指令:gmx hbond -s md.tpr -f mdnoPBC.xtc -n index.ndx -num hbnum_A_B.xvg;选择index.ndx文件中指定的A与B链),得到三列数据:第1列是时 ...

回复 Reply

KazusaT 发表于 Post on 2025-11-17 21:49:25
本帖最后由 KazusaT 于 2025-11-17 21:50 编辑

1. 如果我的理解没错,应该第一列是时间,第二列是满足角度+满足距离的氢键数目,第三列是满足距离+不满足角度的氢键数目,也就是说二、三列相加是只考虑距离的氢键数目,这个值可以通过-hbr指定
2. 不对,你没有看完这个hbond.ndx,你看的这个是"donors_hydrogens_A_chain",是A链的供体-氢,例如第一行1 2 1 3 1 4指的是A链上,1原子可以通过1-2、1-3、1-4作为供体形成氢键(显然这应该是一个-NH3+),你想找的那个A链和B链之间的Donor–Hydrogen–Acceptor在这个ndx最后,应该是[ hbonds_A_chain-B_chain ]部分,你完整读一下这个ndx(你这个可能非常长,你可以随便画两个短肽链稍微模拟一下得到一个比较短的文件看一看)还可以找到记录A链氢键受体的[ acceptors_A_chain ],记录B链的[donors_hydrogens_B_chain]、[ acceptors_B_chain ]

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 16:42 , Processed in 0.174613 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list