24小时热门版块排行榜    

查看: 3168  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ustbmars

铜虫 (小有名气)

[交流] 【求助】如何计算平均滞留时间(residence time) 已有3人参与

分析MD模拟轨迹, 想比较氧离子在不同元素离子周围的平均滞留时间, 即residence time, 该如何统计? 请有经验的虫虫们不吝赐教啊. 如果有code可供参考那就更好了. 先谢谢了!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+2):谢谢 2010-12-23 10:29:21
引用回帖:
Originally posted by ustbmars at 2010-12-21 23:45:53:
谢谢, 我更倾向于第一种方法, 意义比较明确.



第一种方法所谓“意义明确“其实只是操作起来更为直观罢了。第二种方法也是可行的,所以建议使用两种方法,然后对比下结果,如果一致那就很happy了不是
4楼2010-12-22 00:30:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

yahoohoo

铁杆木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+3):thank you 2010-12-21 21:22:49
引用回帖:
Originally posted by ustbmars at 2010-12-21 05:47:38:
分析MD模拟轨迹, 想比较氧离子在不同元素离子周围的平均滞留时间, 即residence time, 该如何统计? 请有经验的虫虫们不吝赐教啊. 如果有code可供参考那就更好了. 先谢谢了!

有两种可能的方法,一种是记录氧原子与某种元素离子的距离随时间的信息(trajectory),然后根据人为设定的标准(例如两者距离在$r_L$与$r_U$之间定义为相互作用的状态)从trajectory中提取滞留时间,然后就能得到该时间的分布,最后便可以确定平均滞留时间;另一种是计算距离的时间相关函数,从而确定特征的相关时间。这个时间可以看作是平均滞留时间。后一种思路你也可以参考Gromacs中氢键的指认算法。
2楼2010-12-21 20:21:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbmars

铜虫 (小有名气)

谢谢, 我更倾向于第一种方法, 意义比较明确.
引用回帖:
Originally posted by yahoohoo at 2010-12-21 20:21:16:


有两种可能的方法,一种是记录氧原子与某种元素离子的距离随时间的信息(trajectory),然后根据人为设定的标准(例如两者距离在$r_L$与$r_U$之间定义为相互作用的状态)从trajectory中提取滞留时间,然后就能 ...

3楼2010-12-21 23:45:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿吉你行

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
1197886楼: Originally posted by yahoohoo at 2010-12-21 20:21:16
有两种可能的方法,一种是记录氧原子与某种元素离子的距离随时间的信息(trajectory),然后根据人为设定的标准(例如两者距离在$r_L$与$r_U$之间定义为相互作用的状态)从trajectory中提取滞留时间,然后就能得 ...

“分析MD模拟轨迹, 想比较氧离子在不同元素离子周围的平均滞留时间, 即residence time, 该如何统计?”
你的回答“有两种可能的方法,一种是记录氧原子与某种元素离子的距离随时间的信息(trajectory),然后根据人为设定的标准(例如两者距离在$r_L$与$r_U$之间定义为相互作用的状态)从trajectory中提取滞留时间,然后就能得到该时间的分布,最后便可以确定平均滞留时间;另一种是计算距离的时间相关函数,从而确定特征的相关时间。这个时间可以看作是平均滞留时间。后一种思路你也可以参考Gromacs中氢键的指认算法。”
想请问下第一种方法具体怎么弄?我是用DL_POLY算的MD,渴望您的回复
5楼2013-03-06 19:46:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 04:21 by 4wMiSEwB6436
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +11 Tide man 2026-08-10 12/600 2026-08-15 02:12 by home3163
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 7/350 2026-08-15 00:40 by 4wMiSEwB6436
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见