| 查看: 3411 | 回复: 5 | ||||
jianjunjiang铁杆木虫 (正式写手)
|
[求助]
如何计算均方位移（MSD）中原子在不同时刻的位移差
|
|
各位大侠： 在计算均方位移（MSD）时，要计算出不同时刻某个粒子的位移（比如说ri（t）-ri(0)），ri（t）和ri(0)是粒子在不同时刻的位置， 请问这种情况下，是不是也要利用最小镜像原理计算出ri（t）-ri(0)相应模的最小值。（请注意，不同时刻模拟盒子的大小不一样） |
回复此楼» 猜你喜欢
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
-
2026年叶企孙基金
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有13人回复
-
基金系统什么内容也没有
已经有9人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有11人回复
-
怎么查啊
已经有6人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
jianjunjiang: 金币+5, ★★★★★最佳答案 2016-07-31 23:23:51
感谢参与，应助指数 +1
jianjunjiang: 金币+5, ★★★★★最佳答案 2016-07-31 23:23:51
|
"有一点要特别注意的是，在模拟流体（气体和液体）时，要定义两套坐标，对其中一套坐标施加一个周期边界条件（即每当某个粒子跑出模拟盒子，就把它拉回来；对于固体，这是没有必要的），对另一套坐标不用施加周期边界条件（即不用将跑出盒子的粒子拉回盒子，而是任其自由跑动，即允许体系“扩散”）。求力的时候，一定要用那套施加了周期边界条件的坐标（不然求的力是不对的）；而计算MSD时，则用那套没有施加周期边界条件的坐标（不然计算的MSD是错的）。" 上面一段话取自于我的博文（详见我的名为《用分子动力学模拟计算扩散系数》的博文：http://blog.sciencenet.cn/home.php?mod=space&uid=3102863&do=blog&id=991891）。我的博文中提供了完整的C语言源代码、matlab后处理脚本、以及算例。 再回到你的问题“ 在计算均方位移（MSD）时，要计算出不同时刻某个粒子的位移（比如说ri（t）-ri(0)），ri（t）和ri(0)是粒子在不同时刻的位置， 请问这种情况下，是不是也要利用最小镜像原理计算出ri（t）-ri(0)相应模的最小值。（请注意，不同时刻模拟盒子的大小不一样）” 答案是：不需要，因为为计算MSD所保存的坐标是“自由”坐标，不需要施加周期边界条件的。 |
2楼2016-07-31 23:11:09
jianjunjiang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8344.2
- 散金: 60
- 红花: 1
- 帖子: 932
- 在线: 252.1小时
- 虫号: 509147
- 注册: 2008-02-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2016-07-31 23:23:53
jianjunjiang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8344.2
- 散金: 60
- 红花: 1
- 帖子: 932
- 在线: 252.1小时
- 虫号: 509147
- 注册: 2008-02-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2016-07-31 23:37:12
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
5楼2016-08-01 03:29:41
jianjunjiang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8344.2
- 散金: 60
- 红花: 1
- 帖子: 932
- 在线: 252.1小时
- 虫号: 509147
- 注册: 2008-02-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2016-08-01 08:14:13
10