24小时热门版块排行榜    

查看: 3071  |  回复: 11

jolly1987

金虫 (小有名气)

[求助] 在DPD中如何计算获得均方末端距?

请问在DPD中如何计算获得均方末端距?在Analysis里有看到End-to-end distribution,帮助软件说:“End-point distribution: returns a histogram of the end-to-end distance of a molecule for each of the frames selected. The histogram is accumulated over all molecules and all steps up until the previous frame. Also returned is the corresponding radial distribution, which is the histogram divided by the total bin count and the volume of each bin."
没太明白什么意思,能用来计算均方末端距吗?望高人指点~谢谢!
回复此楼

» 收录本帖的淘帖专辑推荐

xuexijisuan

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yahoohoo

铁杆木虫 (著名写手)

★ ★ ★
ghcacj(金币+3): 谢谢 2011-09-05 18:00:51
引用回帖:
5楼: Originally posted by jolly1987 at 2011-09-05 16:59:55:
谢谢大侠指!
Distribution是概率密度,是吗?我一开始迷惑的是,这两幅图应该是对同一结果的不同表达方式,那么他们的图不应该是一样的吗?但有时输出的结果分布曲线不同。
另外,左图的 P(R_e)太小,积分得到 ...

在使用很多bins的条件下,图1应该是图2归一化的结果。
完整的计算均方末端矩的公式是 \langle R_e^2 \rangle = {\int P(R_e) R_e^2 d R_e} \over {\int P(R_e) R_e^2 d R_e}。不明白为什么你用两个图计算的结果不同。
6楼2011-09-05 17:31:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yahoohoo

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★
jolly1987(金币+1): 2011-09-05 16:22:07
ghcacj(金币+3): 谢谢 2011-09-05 18:00:30
引用回帖:
1楼: Originally posted by jolly1987 at 2011-09-05 14:40:55:
请问在DPD中如何计算获得均方末端距?在Analysis里有看到End-to-end distribution,帮助软件说:“End-point distribution: returns a histogram of the end-to-end distance of a molecule for each of the frame ...

这段话的意思是:先建立末端矩的分布,也就是P(R_e)。那么均方末端矩$\sqrt \langle R_e^2 \rangle = \int P(R_e) R_e^2 d R_e$。
2楼2011-09-05 16:06:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jolly1987

金虫 (小有名气)

引用回帖:
2楼: Originally posted by yahoohoo at 2011-09-05 16:06:11:
这段话的意思是:先建立末端矩的分布,也就是P(R_e)。那么均方末端矩$\sqrt \langle R_e^2 \rangle = \int P(R_e) R_e^2 d R_e$。

不好意思,请问“$\sqrt \langle R_e^2 \rangle = \int P(R_e) R_e^2 d R_e$“是什么?
请问图中Distribution的单位是什么?Bin Count是什么?单位是什么?谢谢!

3楼2011-09-05 16:21:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
jolly1987(金币+1): 2011-09-05 16:52:29
ghcacj(金币+4): 谢谢 2011-09-05 18:00:41
引用回帖:
3楼: Originally posted by jolly1987 at 2011-09-05 16:21:39:
不好意思,请问“$\sqrt \langle R_e^2 \rangle = \int P(R_e) R_e^2 d R_e$“是什么?
请问图中Distribution的单位是什么?Bin Count是什么?单位是什么?谢谢!

“$\sqrt \langle R_e^2 \rangle = \int P(R_e) R_e^2 d R_e$“是求均方末端矩的公式啊

自己思考思考distribution表示什么,那你就该知道单位是什么。在分子模拟中我们一般用离散的histogram来近似连续的distribution,那么我们需要用到bins,例如Re的范围在1-10,我们可以用100个bins,第一个bin对应的Re范围就是[1,1.1]。顾名思义,bin count就是落在特定bin中的样本数目。
4楼2011-09-05 16:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jolly1987

金虫 (小有名气)

引用回帖:
4楼: Originally posted by yahoohoo at 2011-09-05 16:27:05:
“$\sqrt \langle R_e^2 \rangle = \int P(R_e) R_e^2 d R_e$“是求均方末端矩的公式啊

自己思考思考distribution表示什么,那你就该知道单位是什么。在分子模拟中我们一般用离散的histogram来近似连续的dist ...

谢谢大侠指!
Distribution是概率密度,是吗?我一开始迷惑的是,这两幅图应该是对同一结果的不同表达方式,那么他们的图不应该是一样的吗?但有时输出的结果分布曲线不同。
另外,左图的 P(R_e)太小,积分得到的均方末端矩很小,不太符合实际;将右图转换成概率密度分布图后积分的结果还比较靠谱~
5楼2011-09-05 16:59:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jolly1987

金虫 (小有名气)

引用回帖:
6楼: Originally posted by yahoohoo at 2011-09-05 17:31:51:
在使用很多bins的条件下,图1应该是图2归一化的结果。
完整的计算均方末端矩的公式是 \langle R_e^2 \rangle = {\int P(R_e) R_e^2 d R_e} \over {\int P(R_e) R_e^2 d R_e}。不明白为什么你用两个图计算的结果 ...

明白了,谢谢!
7楼2011-09-05 18:54:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jolly1987

金虫 (小有名气)

引用回帖:
6楼: Originally posted by yahoohoo at 2011-09-05 17:31:51:
在使用很多bins的条件下,图1应该是图2归一化的结果。
完整的计算均方末端矩的公式是 \langle R_e^2 \rangle = {\int P(R_e) R_e^2 d R_e} \over {\int P(R_e) R_e^2 d R_e}。不明白为什么你用两个图计算的结果 ...

不好意思,大侠,刚学习MS,想再问一句,\langle R_e^2 \rangle = {\int P(R_e) R_e^2 d R_e} \over {\int P(R_e) R_e^2 d R_e}这个公式是可以直接写到MS里,然后给出结果的吗?谢谢!
8楼2011-09-05 22:31:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

引用回帖:
8楼: Originally posted by jolly1987 at 2011-09-05 22:31:07:
不好意思,大侠,刚学习MS,想再问一句,\langle R_e^2 \rangle = {\int P(R_e) R_e^2 d R_e} \over {\int P(R_e) R_e^2 d R_e}这个公式是可以直接写到MS里,然后给出结果的吗?谢谢!

我没用ms。公式是latex格式的,
9楼2011-09-06 00:18:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jolly1987

金虫 (小有名气)

引用回帖:
9楼: Originally posted by yahoohoo at 2011-09-06 00:18:11:
我没用ms。公式是latex格式的,

哦~好吧,非常感谢!
10楼2011-09-06 08:24:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jolly1987 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 6/300 2026-08-25 00:11 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见