| 24小時(shí)熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1343 | 回復(fù): 6 | ||
| 本帖產(chǎn)生 2 個(gè) 模擬EPI ,點(diǎn)擊這里進(jìn)行查看 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
yuelingxiao新蟲 (小有名氣)
|
[求助]
DPD方法初始條件的實(shí)現(xiàn)
|
|
|
小弟DPD初學(xué)者,請(qǐng)教其初始條件是如何實(shí)現(xiàn)的。 在分子動(dòng)力學(xué)中,常用的初始條件可以選為:(1)令初始位置在差分劃分網(wǎng)格的格子上,初始速度則從玻爾茲曼分布隨機(jī)抽樣得到。(2)令初始位置隨機(jī)地偏離差分劃分網(wǎng)格的格子,初始速度為零。(3)令初始位置隨機(jī)地偏離差分劃分網(wǎng)格的格子,初始速度也是從玻爾茲曼分布隨機(jī)抽樣得到。 關(guān)于第一條:是不是還要保證抽樣得到的速度滿足能量均分定理,也即總動(dòng)能應(yīng)該為(3/2)NkT.還是只要是從玻爾茲曼分布抽樣得出的速度分布,這個(gè)速度分布就會(huì)滿足能量均分定理 關(guān)于第二條:倘使這么給定的話,其初始速度為都為0,是不是意味著系統(tǒng)溫度也為0呢,這樣的話不是不符合常理了嗎? 關(guān)于第三條,同理 DPD的初始條件與分子動(dòng)力學(xué)的初始條件是一樣的嗎? |
新蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
|
根據(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的偽隨機(jī)數(shù),(一般我們還保證體系的總動(dòng)量為0,為此我們需求得體系質(zhì)心的動(dòng)量,然后在每個(gè)粒子上除去質(zhì)心的動(dòng)量),然后rescale總動(dòng)能,使其等于$3/2 N k_B T$。這樣速度的初始化就完成了。 |
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 282求調(diào)劑 +17 | ycy1201 2026-04-01 | 19/950 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑推薦 材料 304 +21 | 荷包蛋hyj 2026-03-26 | 21/1050 |
|
|
[考研]
|
廈大化工 2026-04-01 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 5/250 |
|
|
[考研] 求調(diào)劑 +5 | 零八# 2026-03-27 | 5/250 |
|
|
[考研] 一志愿中農(nóng)0710生物學(xué),微生物方向總分338求調(diào)劑 +3 | 柒xxxx. 2026-03-26 | 3/150 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +6 | 篤行致遠(yuǎn). 2026-03-28 | 6/300 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +7 | 12.27 2026-03-27 | 7/350 |
|
|
[考研] 復(fù)試調(diào)劑 +7 | 雙馬尾痞老板2 2026-03-31 | 7/350 |
|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 8/400 |
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊(duì)招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|