| 5 | 1/1 | 返回列表 |
| 查看: 3258 | 回復: 5 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
jianjunjiang鐵桿木蟲 (正式寫手)
|
[求助]
如何計算均方位移(MSD)中原子在不同時刻的位移差
|
|||
|
各位大俠: 在計算均方位移(MSD)時,要計算出不同時刻某個粒子的位移(比如說ri(t)-ri(0)),ri(t)和ri(0)是粒子在不同時刻的位置, 請問這種情況下,是不是也要利用最小鏡像原理計算出ri(t)-ri(0)相應模的最小值。(請注意,不同時刻模擬盒子的大小不一樣) |
鐵桿木蟲 (正式寫手)
專家顧問 (著名寫手)
![]() |
專家經驗: +218 |
|
"有一點要特別注意的是,在模擬流體(氣體和液體)時,要定義兩套坐標,對其中一套坐標施加一個周期邊界條件(即每當某個粒子跑出模擬盒子,就把它拉回來;對于固體,這是沒有必要的),對另一套坐標不用施加周期邊界條件(即不用將跑出盒子的粒子拉回盒子,而是任其自由跑動,即允許體系“擴散”)。求力的時候,一定要用那套施加了周期邊界條件的坐標(不然求的力是不對的);而計算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所保存的坐標是“自由”坐標,不需要施加周期邊界條件的。 |
鐵桿木蟲 (正式寫手)
專家顧問 (著名寫手)
![]() |
專家經驗: +218 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 319求調劑 +14 | 太容易1018 2026-04-01 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 265求調劑 +6 | 梁梁校校 2026-04-01 | 6/300 |
|
|
[考研] 化學308分調劑 +14 | 你好明天你好 2026-03-30 | 15/750 |
|
|
[考研] 085600 295分求調劑 +16 | W55j 2026-03-30 | 19/950 |
|
|
[考研] 材料化工340求調劑 +5 | jhx777 2026-03-30 | 5/250 |
|
|
[考研] 一志愿廈門大學材料工程專碩354找調劑。! +7 | 貝唄鋇鋇 2026-03-30 | 7/350 |
|
|
[考研] 070305高分子化學與物理 304分求調劑 +14 | c297914 2026-03-28 | 14/700 |
|
|
[考研] 一志愿211,335分,0856,求調劑院校和導師 +14 | 傾____蕭 2026-03-27 | 15/750 |
|
|
[考研] 一志愿北京科技,085601總分305求調劑 +9 | 半生瓜! 2026-04-01 | 11/550 |
|
|
[考研] 一志愿安徽大學計算機科學與技術學碩,331分求調劑 +5 | 蔣昌鵬qtj 2026-04-01 | 5/250 |
|
|
[考研] 307分求調劑 +14 | (o~o) 2026-03-31 | 15/750 |
|
|
[考研] 生物與醫(yī)藥考研調劑 +5 | 鐵憨憨123425 2026-03-31 | 5/250 |
|
|
[考研] 0710生物學考研調劑 +3 | 李多米lee. 2026-03-27 | 4/200 |
|
|
[考研] 0855機械初試280求調劑 +3 | kazenotori 2026-03-31 | 3/150 |
|
|
[考研] 【調劑】一志愿廈大生物與醫(yī)藥調劑 +3 | Echo蝦米 2026-03-31 | 3/150 |
|
|
[考研] 352分-085602-一志愿985 +6 | 海納百川Ly 2026-03-29 | 6/300 |
|
|
[考研] 346求調劑 一志愿070303有機化學 +11 | 蘿卜燉青菜 2026-03-28 | 12/600 |
|
|
[考研] 318求調劑 +10 | 陳晨79 2026-03-30 | 10/500 |
|
|
[考研] 085602 化學工程專碩 340分求調劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 292求調劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|