计算化学公社

标题: 请问gmx hbond -life输出结果中的tp(t)能说明什么结果？ [打印本页]

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

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

作者
Author: UTM 时间: 2022-6-27 13:15
本帖最后由 UTM 于 2022-6-27 13:20 编辑

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

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

作者
Author: Frozen-Penguin 时间: 2022-6-27 21:12

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

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

作者
Author: UTM 时间: 2022-6-27 21:31
本帖最后由 UTM 于 2022-6-27 21:32 编辑

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

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

作者
Author: Frozen-Penguin 时间: 2022-6-27 21:39

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，数值的随便说的

欢迎光临计算化学公社 (http://bbs.keinsci.com/)