| 5 | 1/1 | 返回列表 |
| 查看: 1196 | 回復(fù): 4 | ||||
懶羊羊Q金蟲 (小有名氣)
|
[求助]
discover模擬擴(kuò)散行為 已有1人參與
|
|
在用discover模擬物質(zhì)的擴(kuò)散行為時(shí),有沒(méi)有人知道擴(kuò)散位移是怎樣得到的? 我說(shuō)的不是均方位移(MSD),而是指隨著時(shí)間的變化,分子相對(duì)于原點(diǎn)的擴(kuò)散位移。希望懂的人指點(diǎn)一二。 擴(kuò)散位移圖.gif |
金蟲 (小有名氣)
專家顧問(wèn) (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
對(duì)于各向同性系統(tǒng),并且處于平衡態(tài),粒子相對(duì)于t=0起始點(diǎn)的平均位移為0.因?yàn)檫@時(shí)粒子的運(yùn)動(dòng)就是個(gè)隨機(jī)游走,隨機(jī)游走的平均位移為0,但MSD隨時(shí)間成正比。 如果你的系統(tǒng)不處于平衡態(tài),或者各向不同性,則情況很復(fù)雜。 說(shuō)白了,這個(gè)問(wèn)題就是求解系統(tǒng)的Fokker-Planck Equation,得到粒子在時(shí)間t處于x的概率P(x,t),然后算平均值就可以。對(duì)于各向同性下的自由擴(kuò)散,P(x,t)是個(gè)中心在原點(diǎn)的高斯分布,所以平均位移為0,但MSD(方差)隨時(shí)間成正比。 |
金蟲 (小有名氣)
專家顧問(wèn) (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
我前面講的可能有點(diǎn)問(wèn)題,我的大致意思是在各向同性的平衡態(tài)中的自由擴(kuò)散,如果粒子在t=0時(shí)坐標(biāo)為x0,t=t1時(shí)坐標(biāo)為x1,那么x1的平均值是等于x0的,但是(x1-x0)^2(也就是MSD)的平均值是與時(shí)間成正比的。 如果你的體系在某個(gè)方向上是非isotropic的,那么x1的平均值可能不等于x0的,但是(x1-x0)^2還是會(huì)隨時(shí)間增加的,直到粒子碰到了confinement(比如盒子的墻壁)為止。 所以你把MSD開(kāi)方對(duì)時(shí)間作圖發(fā)現(xiàn)隨時(shí)間增加是正常的。而你貼的圖片里的鋸齒形圖形,大致的趨勢(shì)也是隨時(shí)間逐步增加的。至于你得到的MSD開(kāi)方后的曲線不是鋸齒形的,是因?yàn)镸SD是多條軌跡的平均,而你圖片中只給出了一條軌跡,對(duì)多條軌跡求平均后鋸齒形的噪音就會(huì)消除不少。 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 8/400 |
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 085600材料與化工 +4 | 安全上岸! 2026-03-16 | 4/200 |
|
|
[考研] 283求調(diào)劑 +3 | 聽(tīng)風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[基金申請(qǐng)]
有必要更換申報(bào)口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 314求調(diào)劑 +7 | 無(wú)懈可擊的巨人 2026-03-12 | 7/350 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|