24小时热门版块排行榜    

Znn3bq.jpeg
查看: 6324  |  回复: 1

xueyiran2007

铜虫 (小有名气)

[求助] MSD均方位移如何计算 已有1人参与

对MSD计算不是很清楚,这两天搜了一些相关文献,还是不太懂,很惆怅,请各位大神指点,不胜感激!

1. MSD均方位移,<[r(t)-r(0)]2>,
是指t时刻内所有相邻数据点的平方和除以粒子数N?还是指t时刻的位置与初始时刻位置距离的平方和除以粒子数N?
对于MSD-t曲线图是怎么算的不清楚。如果我要计算某一时刻t1的MSD,是要在初始时刻t0与目标时刻t1之间等分n步,分别计算各相邻步之间距离的平方和,再将n步的平方和相加,再除以系统粒子数N,就是t1时刻的MSD吗?

2. 尖括号<>表示的意思
我看书上<>指系综平均,根据遍历性原则,也就等价于时间平均,是将不同时刻系统的宏观状态相加除以系统运行时间
但MSD这里的尖括号<>指的是:粒子数N的平均
问题:分子模拟中尖括号<>代表的意思不是固定的吗?可以表示不同意思?
回复此楼

» 本帖@通知

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
月只蓝: 金币+3, 鼓励交流! 2014-05-30 18:49:40
1. 还是指t时刻的位置与初始时刻位置距离的平方和除以粒子数N?
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                              yes

2. 如果我要计算某一时刻t1的MSD,是要在初始时刻t0与目标时刻t1之间等分n步,分别计算各相邻步之间距离的平方和,再将n步的平方和相加,再除以系统粒子数N,就是t1时刻的MSD吗?
calculate for each particle, then divided by N.

You can use 'compute msd' command in lammps. For more accurate MSD calculation algorithm, please refer to
1) Frenkel and Smit, Understanding Molecular Simulations, 2nd Edition, Section 4.4, Page 91
2) Allen and Tidelsey, Computer Simulation of Liquids, Section 6.3, page 85

3.我看书上<>指系综平均,根据遍历性原则,也就等价于时间平均,是将不同时刻系统的宏观状态相加除以系统运行时间,但MSD这里的尖括号<>指的是:粒子数N的平均

MSD is a property of single particle, so you have to average it over N particles.
'<>' still means ensemble average here.
2楼2014-05-17 06:50:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xueyiran2007 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 7/350 2026-05-20 19:56 by 88817753
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +11 1234567wang 2026-05-17 13/650 2026-05-20 18:55 by tanfancy
[基金申请] 评审有感 +14 popular289 2026-05-18 23/1150 2026-05-20 17:39 by lincyb
[基金申请] 提交了我也来说说感想 +6 fummck 2026-05-20 6/300 2026-05-20 17:29 by 姜太公烤鱼
[论文投稿] Sci. Bull. 悲剧经验 +6 jyang1999 2026-05-16 7/350 2026-05-20 16:39 by feng6531
[基金申请] 今年审到国自然15份,谈谈感受 +19 国自然国社科中 2026-05-17 20/1000 2026-05-20 14:14 by 仲夏夜的星星
[基金申请] 国自然评分 +3 无名者登山 2026-05-20 4/200 2026-05-20 13:29 by vito刘
[基金申请] 河北省自然科学基金 +3 Peterchao 2026-05-18 3/150 2026-05-20 08:57 by 霸_霸
[教师之家] 上海大学实验技术岗位非升即走 +12 嘻嘻哈哈乐呵呵 2026-05-15 13/650 2026-05-20 08:34 by xli1984
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 11/550 2026-05-19 22:26 by 而立得乐
[考博] 找博士生导师 +6 小代想上岸 2026-05-15 7/350 2026-05-19 10:22 by free_fisher
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 别被青基扩招骗了!26年科研内卷才刚刚开始 +3 国自然国社科中 2026-05-14 4/200 2026-05-19 08:48 by archvillain
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
信息提示
请填处理意见