计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › gmx hbond工具进行氢键分析数据解读和作图求助

gmx hbond工具进行氢键分析数据解读和作图求助

查看数: 3440 | 评论数: 3 | 收藏 Add to favorites 4

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

莫落

发布时间: 2023-10-28 22:55

正文摘要:

1、gmx hbond工具进行氢键分析分析时，得到两列数据： @ title "Hydrogen Bond Distribution" @ xaxis label "Donor - Acceptor Distance (nm)" @ ya ...

回复 Reply

sobereva 发表于 Post on 2023-11-1 08:29:10

莫落发表于 2023-10-31 21:15
感谢社长，另外我还有几点疑问麻烦社长解答：
1、hbang.xvg里的数据，前面一列是角度，后面一列应该对应 ...

1 不是加和为1，而是积分为1
hbdist类似，体现距离分布概率
没法直接乘。距离和角度不是无相关性的

2 说清楚怎么得到的
如果是我说的那种自写VMD tcl做法，脚本给出的是每一帧在(r,theta)图里的位置，转化成概率分布是之后的事（用origin之类实现，或者直接就通过VMD tcl脚本实现）

3 只能区分组间和组内的氢键。不是以分子为单位的。

莫落 发表于 Post on 2023-10-31 21:15:47

sobereva 发表于 2023-10-30 02:28
1 曲线的积分为1，是归一化的

2 自己写个VMD脚本对每一帧的结构计算对应的r和theta，散点数据弄到诸如Or ...

感谢社长，另外我还有几点疑问麻烦社长解答：
1、hbang.xvg里的数据，前面一列是角度，后面一列应该对应的是角度的统计分布概率，加和为1；那对应hbdist里的第二例是什么意义呢？如果也是概率，那我把这两列对应相乘也会得到一个矩阵，能绘制成温度图；
2、VMD得到的r和theta是分布概率吗？
3、gmx hbond可以区分分子内氢键和分子间氢键吗？VMD可以吗？
再次感谢社长

sobereva 发表于 Post on 2023-10-30 02:28:03

1 曲线的积分为1，是归一化的

2 自己写个VMD脚本对每一帧的结构计算对应的r和theta，散点数据弄到诸如Origin里，转化成概率密度分布，绘制成填色图就完了。跟gmx hbond没关系

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register