 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 1336 | 回復(fù): 6 | ||
| 本帖產(chǎn)生 2 個(gè) 模擬EPI ，點(diǎn)擊這里進(jì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é)的初始條件是一樣的嗎？ |
回復(fù)此樓» 本帖@通知
» 猜你喜歡
291求調(diào)劑
已經(jīng)有15人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有17人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有5人回復(fù)
-
070300一志愿211，312分求調(diào)劑院校
已經(jīng)有6人回復(fù)
-
物理學(xué)調(diào)劑
已經(jīng)有3人回復(fù)
-
08工科，295，接受跨專業(yè)調(diào)劑
已經(jīng)有5人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有12人回復(fù)
-
本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿西交288化工專碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
308求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 國(guó)外做多因素條件的實(shí)驗(yàn)中一般選擇什么方法？
已經(jīng)有8人回復(fù)
- 急求 一篇 專利文獻(xiàn) ！�。�！一種非平行文本條件下的語音轉(zhuǎn)換方法
已經(jīng)有2人回復(fù)
- 【500金幣本版活動(dòng)】“經(jīng)驗(yàn)力場(chǎng)分子模擬+第一性原理計(jì)算問題討論互助及資源共享”
已經(jīng)有82人回復(fù)
- 【求助】求能夠氧化二級(jí)醇為酮的催化劑制備方法及反應(yīng)條件
已經(jīng)有12人回復(fù)
- 【求助】高真空條件下甲醇回收方法
已經(jīng)有19人回復(fù)
- 【求助】由間硝基苯甲酸合成間氨基苯甲酸的方法和條件
已經(jīng)有13人回復(fù)
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時(shí)
- 蟲號(hào): 74894
- 注冊(cè): 2005-06-15
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
2樓2011-10-22 13:19:52
yuelingxiao
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 82.6
- 散金: 7
- 紅花: 2
- 帖子: 77
- 在線: 40小時(shí)
- 蟲號(hào): 1241850
- 注冊(cè): 2011-03-22
- 專業(yè): 傳熱傳質(zhì)學(xué)
3樓2011-10-22 22:06:22
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時(shí)
- 蟲號(hào): 74894
- 注冊(cè): 2005-06-15
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★
ghcacj(金幣+4, 模擬EPI+1): 謝謝 2011-10-23 14:42:52
ghcacj(金幣+4, 模擬EPI+1): 謝謝 2011-10-23 14:42:52
|
根據(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$。這樣速度的初始化就完成了。 |
4樓2011-10-22 22:42:44
yuelingxiao
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 82.6
- 散金: 7
- 紅花: 2
- 帖子: 77
- 在線: 40小時(shí)
- 蟲號(hào): 1241850
- 注冊(cè): 2011-03-22
- 專業(yè): 傳熱傳質(zhì)學(xué)
5樓2011-10-23 10:16:56
yuelingxiao
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 82.6
- 散金: 7
- 紅花: 2
- 帖子: 77
- 在線: 40小時(shí)
- 蟲號(hào): 1241850
- 注冊(cè): 2011-03-22
- 專業(yè): 傳熱傳質(zhì)學(xué)
6樓2011-10-23 13:22:05
yahoohoo
鐵桿木蟲 (著名寫手)
- 模擬EPI: 10
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.55
- 金幣: 7632.4
- 散金: 251
- 紅花: 4
- 帖子: 1176
- 在線: 167小時(shí)
- 蟲號(hào): 74894
- 注冊(cè): 2005-06-15
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
7樓2011-10-23 15:17:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yuelingxiao 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +6 | 吃的不少 2026-03-30 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +4 | 睿08 2026-03-27 | 4/200 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +7 | 蘿卜燉青菜 2026-03-28 | 8/400 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +8 | 26要上岸 2026-03-28 | 8/400 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|