| 查看: 3400 | 回复: 5 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
jianjunjiang铁杆木虫 (正式写手)
|
[求助]
如何计算均方位移（MSD）中原子在不同时刻的位移差
|
|||
|
各位大侠： 在计算均方位移（MSD）时，要计算出不同时刻某个粒子的位移（比如说ri（t）-ri(0)），ri（t）和ri(0)是粒子在不同时刻的位置， 请问这种情况下，是不是也要利用最小镜像原理计算出ri（t）-ri(0)相应模的最小值。（请注意，不同时刻模拟盒子的大小不一样） |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有78人回复
-
filecode=后面第一个是大写字母
已经有8人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有13人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
jianjunjiang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8343.2
- 散金: 60
- 红花: 1
- 帖子: 932
- 在线: 252.1小时
- 虫号: 509147
- 注册: 2008-02-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2016-08-01 08:14:13
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 (幼儿园)
- 金币: 8343.2
- 散金: 60
- 红花: 1
- 帖子: 932
- 在线: 252.1小时
- 虫号: 509147
- 注册: 2008-02-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2016-07-31 23:23:53
jianjunjiang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8343.2
- 散金: 60
- 红花: 1
- 帖子: 932
- 在线: 252.1小时
- 虫号: 509147
- 注册: 2008-02-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2016-07-31 23:37:12
25