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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » MSD分析扩散系数的取值段
北京石油化工学院2026年研究生招生接收调剂公告
211/1返回列表
查看: 4325  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿生物与医药，296分，求调剂 +11 66鹿 2026-04-03 13/650 2026-04-06 21:45 by 52305043001
[考研] 专硕304找调剂，一线城市最好 +3 李lsl李 2026-04-05 3/150 2026-04-06 12:16 by ffffjjjj
[考研] 求助071001调剂！！！ +4 黄守松 2026-04-05 5/250 2026-04-06 10:55 by 1028907439
[考研] 372分，材料与化工，一志愿湖南大学，求调剂 +3 蓝笺片 2026-04-01 3/150 2026-04-06 09:04 by 无际的草原
[考研] 求调剂 +11 熊二想上岸 2026-04-04 11/550 2026-04-05 22:21 by 醉翁wl
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-02 11/550 2026-04-04 23:37 by 永字号
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 材料调剂 +12 一样YWY 2026-04-02 13/650 2026-04-04 20:49 by 蓝云思雨
[考研] 309求调剂 +6 刘刘刘1231 2026-04-02 7/350 2026-04-04 13:41 by liucky
[考研] 387求调剂 +4 爱吃片豆土 2026-04-03 5/250 2026-04-04 08:10 by 岸上的一条鱼
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 学硕化学工程与技术，一志愿中国海洋大学320+求调剂 +8 披星河 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[考研] 324求调剂 +5 想上学求调 2026-04-01 6/300 2026-04-02 10:16 by sanrepian
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 【求调剂】新能源材料本科，一志愿211，初试321 +6 求调剂学校， 2026-04-02 6/300 2026-04-02 09:41 by 晴空210210
[考研] 省双一流重点一本大学招收调剂 +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
信息提示
关闭
请填处理意见
关闭