| 7 | 1/1 | 返回列表 |
| 查看: 1335 | 回復: 6 | ||
| 本帖產(chǎn)生 2 個 模擬EPI ,點擊這里進行查看 | ||
yuelingxiao新蟲 (小有名氣)
|
[求助]
DPD方法初始條件的實現(xiàn)
|
|
|
小弟DPD初學者,請教其初始條件是如何實現(xiàn)的。 在分子動力學中,常用的初始條件可以選為:(1)令初始位置在差分劃分網(wǎng)格的格子上,初始速度則從玻爾茲曼分布隨機抽樣得到。(2)令初始位置隨機地偏離差分劃分網(wǎng)格的格子,初始速度為零。(3)令初始位置隨機地偏離差分劃分網(wǎng)格的格子,初始速度也是從玻爾茲曼分布隨機抽樣得到。 關于第一條:是不是還要保證抽樣得到的速度滿足能量均分定理,也即總動能應該為(3/2)NkT.還是只要是從玻爾茲曼分布抽樣得出的速度分布,這個速度分布就會滿足能量均分定理 關于第二條:倘使這么給定的話,其初始速度為都為0,是不是意味著系統(tǒng)溫度也為0呢,這樣的話不是不符合常理了嗎? 關于第三條,同理 DPD的初始條件與分子動力學的初始條件是一樣的嗎? |
鐵桿木蟲 (著名寫手)
新蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
|
根據(jù)equipartition theorem,平均速率 $\bar{v_i^2} = \frac{k_B T}{m}$。Maxwell-Boltzmann分布可以寫成 $p(v_i) = \frac{1}{\sqrt{2\pi \sigma^2}}e^{\frac{(v_i-0)^2}{2\sigma^2}}$,其中$\simga^2 = \bar{v_i^2} = \frac{k_B T}{m}$。這表示$v_i$滿足Gaussian distribution(平均值為0,variance為\frac{k_B T}{m})。我們需要產(chǎn)生滿足Guassian distribution的平均值為0、variance為1的偽隨機數(shù),(一般我們還保證體系的總動量為0,為此我們需求得體系質(zhì)心的動量,然后在每個粒子上除去質(zhì)心的動量),然后rescale總動能,使其等于$3/2 N k_B T$。這樣速度的初始化就完成了。 |
新蟲 (小有名氣)
新蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿食品科學與工程083200求調(diào)劑 +4 | XQTJZ 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機化學 +7 | 蘿卜燉青菜 2026-03-28 | 8/400 |
|
|
[考研] 359求調(diào)劑 +5 | 王了個楠 2026-03-25 | 5/250 |
|
|
[考研] 一志愿:西北大學,英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
[考研] 289求調(diào)劑 +16 | 新時代材料 2026-03-27 | 16/800 |
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[基金申請] 面上5B能上會嗎? +7 | redcom 2026-03-29 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +8 | 楒桉 2026-03-28 | 8/400 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 085600,材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學ing 2026-03-28 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個好覺_16 2026-03-24 | 4/200 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|