请问gmx hbond -life输出结果中的tp(t)能说明什么结果？

UTM · 发表于 Post on 2022-6-26 21:57:54

如题，请问各位老师同学，gmx hbond -life输出结果中p(t)是氢键能维持到t时刻的概率，那tp(t)说明了什么？另外输出结果中给出的氢键平均寿命是如何计算出来的？

sobereva · 发表于 Post on 2022-6-27 05:14:28

程序显示的平均氢键寿命相当于t*p(t)的积分值

UTM · 发表于 Post on 2022-6-27 13:15:50

本帖最后由 UTM 于 2022-6-27 13:20 编辑

sobereva 发表于 2022-6-27 05:14
程序显示的平均氢键寿命相当于t*p(t)的积分值

请问是t*p(t)对什么的积分？另外单个t*p(t)的值有物理意义吗？

Frozen-Penguin · 发表于 Post on 2022-6-27 21:12:19

UTM 发表于 2022-6-27 13:15
请问是t*p(t)对什么的积分？另外单个t*p(t)的值有物理意义吗？

对时间积分，用统计学语言说这个值就是时间的期望值，也就是平均寿命，因为是期望值，单个值没有意义

UTM · 发表于 Post on 2022-6-27 21:31:00

本帖最后由 UTM 于 2022-6-27 21:32 编辑

Frozen-Penguin 发表于 2022-6-27 21:12
对时间积分，用统计学语言说这个值就是时间的期望值，也就是平均寿命，因为是期望值，单个值没有意义

t*p(t)不是已经是p(t)和t的乘积了吗？怎么还能对时间积分呢？还有p(t)的物理意义是氢键能维持到t时刻的概率吗？

Frozen-Penguin · 发表于 Post on 2022-6-27 21:39:18

UTM 发表于 2022-6-27 21:31
t*p(t)不是已经是p(t)和t的乘积了吗？怎么还能对时间积分呢？还有p(t)的物理意义是氢键能维持到t时刻的概 ...

p(t)是t的概率密度函数，就是统计不同的t出现的频率得到的结果，假设某个氢键在轨迹中形成了4次，时长分别为100ps，200ps，100ps，300ps，那么p(100)=1/4, p(200)=1/2, p(300)=1/4, 平均寿命就是100*1/4+200*1/2+300*1/4，数值的随便说的

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

[GROMACS] 请问gmx hbond -life输出结果中的tp(t)能说明什么结果？

浏览过的版块