应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » MSD分析扩散系数的取值段
211/1返回列表
查看: 4580  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

我是砖头

银虫 (小有名气)

[交流] MSD分析扩散系数的取值段

在MS中分析熔体中某一离子的扩散系数时，得到如下曲线。扩散系数等于曲线斜率的六分之一，但是这个曲线的斜率不是唯一的，请教一下怎么取值

MSD分析扩散系数的取值段
862C.tmp.jpg
回复此楼

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

分子动力学模拟

» 本帖已获得的红花（最新10朵）

chuntao118

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-06-06 16:34:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

brucefan

专家顾问 (著名写手)
最近很多问题都和扩散系数以及热导率有关。我打算这个周末写个文章，把我计算扩散系数和热导率的程序公布出来。处理液态氩体系的程序一共就两百多行代码。

[ 发自手机版 http://muchong.com/3g ]
赞一下(2人) 回复此楼
6楼2016-06-06 18:37:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

chuntao118

新虫 (著名写手)
★
我是砖头(金币+1): 谢谢参与
送红花一朵
怎么说那？你理解的扩散系数我认为是有点点问题的。尽管看公式的话是这个意思，但是这毕竟是一个统计过程，所以，应该是多次统计后去平均值吧！你在计算是要考虑模拟系统的大小吧？也就说统计的时间不应该足够的大。因为这个均方位移必须是有界的。具体的，你可以参考《The Art of Molecular Dynamics Simulation》的第五章。有详细的讨论。至于扩散系数会达到一个稳定值。若你的数据简单拟合的话。应该5-45之间的就行，建议你好好读读这本书的第五章。
赞一下 回复此楼
2楼2016-06-06 17:00:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
★
我是砖头(金币+1): 谢谢参与
我感觉熔体还没有达到平衡状态，应该再增加弛豫时间。最后取达到平衡状态后的一段的斜率
赞一下 回复此楼
4楼2016-06-06 18:14:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)
★
我是砖头(金币+1): 谢谢参与
我还是觉得你的数据反常。一般来说，msd一开始是关联时间的二次函数。随着关联时间的增加，会过渡到一次函数。你的msd的初始行为不大对劲。所以，我怀疑你计算msd的方法有问题。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
5楼2016-06-06 18:21:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiezhou5322

银虫 (正式写手)
★
我是砖头(金币+1): 谢谢参与
引用回帖:
6楼: Originally posted by brucefan at 2016-06-06 18:37:43
最近很多问题都和扩散系数以及热导率有关。我打算这个周末写个文章，把我计算扩散系数和热导率的程序公布出来。处理液态氩体系的程序一共就两百多行代码。
...

期待楼主！
回复此楼
9楼2016-06-11 08:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

707220828

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
在计算MSD的时候出现，楼主有没有出现过不能分析 显示Mean square displacement analysis failed.
No time information is available for trajectories created by Annealed or Quenched dynamics.的情况。楼主知道怎么分析吗 谢谢
赞一下 回复此楼
13楼2016-06-13 11:47:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是砖头

银虫 (小有名气)
引用回帖:
13楼: Originally posted by 707220828 at 2016-06-13 11:47:04
在计算MSD的时候出现，楼主有没有出现过不能分析 显示Mean square displacement analysis failed.
No time information is available for trajectories created by Annealed or Quenched dynamics.的情况。楼主知道 ...

我没有出现过这种情况，不介意的话把你怎么分析的以及你分析的窗口发我看一下
赞一下 回复此楼
14楼2016-06-13 16:05:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

707220828

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
14楼: Originally posted by 我是砖头 at 2016-06-13 16:05:21
我没有出现过这种情况，不介意的话把你怎么分析的以及你分析的窗口发我看一下...

我也才接触MS，觉得好难啊，卡在这里了。我之前也咨询过 但是还是没有解决。楼主是怎么学习MS的呢。谢谢！http://muchong.com/bbs/viewthread.php?tid=10249609
赞一下 回复此楼
15楼2016-06-13 20:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是砖头

银虫 (小有名气)
引用回帖:
6楼: Originally posted by brucefan at 2016-06-06 18:37:43
最近很多问题都和扩散系数以及热导率有关。我打算这个周末写个文章，把我计算扩散系数和热导率的程序公布出来。处理液态氩体系的程序一共就两百多行代码。
...

请问您用的是matlab编程吗？我现在在学习多元金属氧化物熔体中桥氧，非桥氧与自由氧，请问你会吗？QQ735446057，希望可以指导一下，谢谢
赞一下 回复此楼
16楼2016-06-13 21:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是砖头

银虫 (小有名气)
引用回帖:
15楼: Originally posted by 707220828 at 2016-06-13 20:07:47
我也才接触MS，觉得好难啊，卡在这里了。我之前也咨询过 但是还是没有解决。楼主是怎么学习MS的呢。谢谢！http://muchong.com/bbs/viewthread.php?tid=10249609...

我也是自学的，就是不断尝试吧，根据别人论文来学
赞一下 回复此楼
17楼2016-06-13 21:50:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by 我是砖头 at 2016-06-13 21:48:00
请问您用的是matlab编程吗？我现在在学习多元金属氧化物熔体中桥氧，非桥氧与自由氧，请问你会吗？QQ735446057，希望可以指导一下，谢谢...

Yes, I use matlab, C/C++, and CUDA (for graphics processing units)  in my research.
Sorry, I know nothing about "桥氧，非桥氧与自由氧". Molecular dynamics simulation
is a very broad subject, and I am only working on mechanical and thermal properties of
materials.
赞一下 回复此楼
18楼2016-06-13 22:21:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
5楼: Originally posted by brucefan at 2016-06-05 22:21:51
我还是觉得你的数据反常。一般来说，msd一开始是关联时间的二次函数。随着关联时间的增加，会过渡到一次函数。你的msd的初始行为不大对劲。所以，我怀疑你计算msd的方法有问题。
...

谁说MSD一开始是时间的二次函数了? MSD一开始是一个subdiffusive behavior，代表粒子走出第一个溶剂笼的过程，在subdiffusion的过程中MSD和时间的关系根本不是你说的二次函数，而是小于一次函数，就像lz给出的图那样，讨论subdiffusion的论文很多,比如这两篇，都没出现你所说的那种MSD是时间二次函数的情况
http://www.cell.com/biophysj/abstract/S0006-3495(13)05802-5
http://journals.aps.org/prl/abst ... sRevLett.115.220601
赞一下 回复此楼
19楼2016-06-14 15:32:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
19楼: Originally posted by lsloneil at 2016-06-14 15:32:45
谁说MSD一开始是时间的二次函数了? MSD一开始是一个subdiffusive behavior，代表粒子走出第一个溶剂笼的过程，在subdiffusion的过程中MSD和时间的关系根本不是你说的二次函数，而是小于一次函数，就像lz给出的图那 ...

Thanks very much for the references. I think I have to update my knowledge about diffusion. It seems that the physics discussed in these references is quite different from that in simple liquids discussed in textbooks.
赞一下 回复此楼
20楼2016-06-14 16:05:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
19楼: Originally posted by lsloneil at 2016-06-14 15:32:45
谁说MSD一开始是时间的二次函数了? MSD一开始是一个subdiffusive behavior，代表粒子走出第一个溶剂笼的过程，在subdiffusion的过程中MSD和时间的关系根本不是你说的二次函数，而是小于一次函数，就像lz给出的图那 ...

非常感谢lsloneil给出这两篇关于扩散的文献，让我了解了什么是反常扩散。我之前的回答完全是基于正常扩散理论的。但即使考虑反常扩散的可能性，楼主的数据还是有些奇怪。因为，楼主的数据显然也不能用MSD(t)～t^n (n < 1)来描述。所以，一个可能性是楼主计算的只是一个离子（从楼主的描述推断）的MSD(t)，并没有做大量的统计平均，所得数据不具备代表性。
赞一下 回复此楼
21楼2016-06-15 01:16:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
怕哦哦哦3楼
2016-06-06 17:26   回复  
我是砖头(金币+1): 谢谢参与
pingyanwe7楼
2016-06-06 19:32   回复  
我是砖头(金币+1): 谢谢参与
Jameskelly8楼
2016-06-11 01:59   回复  
我是砖头(金币+1): 谢谢参与
nono200910楼
2016-06-11 21:35   回复  
我是砖头(金币+1): 谢谢参与
star_zhang11楼
2016-06-11 21:39   回复  
我是砖头(金币+1): 谢谢参与
springburn12楼
2016-06-11 21:43   回复  
我是砖头(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 我是砖头 的主题更新
211/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 +6 晴天加油 2026-08-26 7/350 2026-08-27 07:30 by pinghe007
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +3 a089 2026-08-26 3/150 2026-08-26 23:21 by young90
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭