 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3256 | 回復(fù): 5 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
jianjunjiang鐵桿木蟲 (正式寫手)
|
[求助]
如何計(jì)算均方位移（MSD）中原子在不同時(shí)刻的位移差
|
||||
|
各位大俠： 在計(jì)算均方位移（MSD）時(shí)，要計(jì)算出不同時(shí)刻某個(gè)粒子的位移（比如說ri（t）-ri(0)），ri（t）和ri(0)是粒子在不同時(shí)刻的位置， 請問這種情況下，是不是也要利用最小鏡像原理計(jì)算出ri（t）-ri(0)相應(yīng)模的最小值。（請注意，不同時(shí)刻模擬盒子的大小不一樣） |
回復(fù)此樓» 猜你喜歡
805600專碩材料與化工348分求調(diào)劑
已經(jīng)有6人回復(fù)
-
重慶大學(xué)材料與化工085600，初試370+，求求調(diào)劑建議
已經(jīng)有5人回復(fù)
-
309求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北交大材料工程，總分358
已經(jīng)有4人回復(fù)
-
302求調(diào)劑一志愿北航070300，本科鄭大化學(xué)
已經(jīng)有11人回復(fù)
-
化學(xué)工程專碩324分，一志愿中國礦業(yè)大學(xué)求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料求調(diào)劑
已經(jīng)有11人回復(fù)
-
【求調(diào)劑】新能源材料本科，一志愿211，初試321
已經(jīng)有4人回復(fù)
-
336求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料求調(diào)劑
已經(jīng)有8人回復(fù)
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時(shí)
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
★ ★ ★
jianjunjiang: 金幣+3, ★★★很有幫助 2016-08-01 08:14:10
jianjunjiang: 金幣+3, ★★★很有幫助 2016-08-01 08:14:10
|
你還是沒有明白我講的要點(diǎn)。在計(jì)算這種距離時(shí)，要用一套新的坐標(biāo)。這套新的坐標(biāo)是不對其作周期邊界條件束縛的。如果用這樣的坐標(biāo)，就不會出現(xiàn)水分子中兩個(gè)氫原子坐標(biāo)一個(gè)在0.1一個(gè)在0.9的情形。你要想象，有兩套坐標(biāo)，一套是用來計(jì)算力的，一套是單獨(dú)用來計(jì)算均方位移的。你去把我給的代碼看明白就會清楚些。 [ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ] |
5樓2016-08-01 03:29:41
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時(shí)
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
jianjunjiang: 金幣+5, ★★★★★最佳答案 2016-07-31 23:23:51
感謝參與，應(yīng)助指數(shù) +1
jianjunjiang: 金幣+5, ★★★★★最佳答案 2016-07-31 23:23:51
|
"有一點(diǎn)要特別注意的是，在模擬流體（氣體和液體）時(shí)，要定義兩套坐標(biāo)，對其中一套坐標(biāo)施加一個(gè)周期邊界條件（即每當(dāng)某個(gè)粒子跑出模擬盒子，就把它拉回來；對于固體，這是沒有必要的），對另一套坐標(biāo)不用施加周期邊界條件（即不用將跑出盒子的粒子拉回盒子，而是任其自由跑動，即允許體系“擴(kuò)散”）。求力的時(shí)候，一定要用那套施加了周期邊界條件的坐標(biāo)（不然求的力是不對的）；而計(jì)算MSD時(shí)，則用那套沒有施加周期邊界條件的坐標(biāo)（不然計(jì)算的MSD是錯的）。" 上面一段話取自于我的博文（詳見我的名為《用分子動力學(xué)模擬計(jì)算擴(kuò)散系數(shù)》的博文：http://blog.sciencenet.cn/home.php?mod=space&uid=3102863&do=blog&id=991891）。我的博文中提供了完整的C語言源代碼、matlab后處理腳本、以及算例。 再回到你的問題“ 在計(jì)算均方位移（MSD）時(shí)，要計(jì)算出不同時(shí)刻某個(gè)粒子的位移（比如說ri（t）-ri(0)），ri（t）和ri(0)是粒子在不同時(shí)刻的位置， 請問這種情況下，是不是也要利用最小鏡像原理計(jì)算出ri（t）-ri(0)相應(yīng)模的最小值。（請注意，不同時(shí)刻模擬盒子的大小不一樣）” 答案是：不需要，因?yàn)闉橛?jì)算MSD所保存的坐標(biāo)是“自由”坐標(biāo)，不需要施加周期邊界條件的。 |
2樓2016-07-31 23:11:09
jianjunjiang
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 8680.7
- 散金: 60
- 紅花: 1
- 帖子: 926
- 在線: 250.9小時(shí)
- 蟲號: 509147
- 注冊: 2008-02-22
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2016-07-31 23:23:53
jianjunjiang
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 8680.7
- 散金: 60
- 紅花: 1
- 帖子: 926
- 在線: 250.9小時(shí)
- 蟲號: 509147
- 注冊: 2008-02-22
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
不過我還有個(gè)問題，假設(shè)我們要計(jì)算兩個(gè)氫原子中心的MSD（為了簡單起見，只考慮x方向），假設(shè)t=0時(shí)，兩個(gè)氫原子坐標(biāo)是0.1L和0.9L（L為t=0時(shí)的盒子在x方向長度），那么t=0他們的中心坐標(biāo)為（0.1L+L+0.9L）/2=L或者（0.9L-L+0.1L）/2=0,經(jīng)過dt時(shí)間后，假設(shè)兩個(gè)氫原子坐標(biāo)是0.2L和0.3L（為了簡單起見，假設(shè)盒子長度不變），那么此時(shí)兩個(gè)氫原子中心坐標(biāo)為0.25L�，F(xiàn)在問題就出現(xiàn)了，這兩個(gè)時(shí)刻（t=0和t=dt)氫原子中心坐標(biāo)的位移是0.25L-L還是0.25L-0呢？ 期待您的答復(fù)，謝謝！ |
4樓2016-07-31 23:37:12
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +7 | 雙馬尾痞老板2 2026-04-01 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 070300一志愿211，312分求調(diào)劑院校 +14 | 小黃鴨寶 2026-03-30 | 14/700 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +13 | 小黃鴨寶 2026-03-30 | 13/650 |
|
|
[考研] 362求調(diào)劑 +13 | 西南交材料專碩3 2026-03-31 | 13/650 |
|
|
|
[考研] 085600，材料與化工321分，求調(diào)劑 +11 | 大饞小子 2026-03-27 | 11/550 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +8 | YYYYYNNNNN 2026-03-26 | 10/500 |
|
|
|
[考研] 一志愿中國科學(xué)院大學(xué)265求調(diào)劑 +8 | 恬淡ye 2026-03-31 | 9/450 |
|
|
|
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +15 | Doleres 2026-03-31 | 16/800 |
|
|
|
[考研] 282求調(diào)劑 +6 | 呼吸都是減肥 2026-04-01 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-29 | 7/350 |
|
|
|
[考研] 080200學(xué)碩，機(jī)械工程專業(yè)277分，求帶走！ +4 | 瓶子PZ 2026-03-31 | 4/200 |
|
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-31 | 4/200 |
|
|
|
[考研] 282求調(diào)劑 +3 | 呼吸都是減肥 2026-03-31 | 3/150 |
|
|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
|
[考研] 求佛 +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|