本帖產(chǎn)生 1 個(gè) 模擬EPI ，點(diǎn)擊這里進(jìn)行查看

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

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é)

[交流] 【在線答疑】經(jīng)典粒子體系的Monte Carlo 模擬之基礎(chǔ)篇已有29人參與

應(yīng)版主lei0736之盛情相邀，特開此貼，雖名為分享，實(shí)為討論，希望與初學(xué)者以及從事MC研究的同行一起交流。第一篇關(guān)于粒子體系的MC模擬基礎(chǔ)，欲將涉及以下幾個(gè)概念與問題：

概念:
正則系綜、相空間、概率密度、細(xì)致平衡(Detailed Balance)、Metropolis重要性抽樣

問題:
(1) 為什么用MC？
(2) MC能做什么？不能做什么？
(3) 怎么做MC？

討論中打算在適當(dāng)之處添加一些偽代碼(c++)以供理解。

歡迎大家多提寶貴意見。

[ Last edited by lei0736 on 2009-12-21 at 18:04 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

仿真建模與計(jì)算	第一性原理	多孔材料	分子模擬入門學(xué)習(xí)資料心得
動(dòng)力學(xué)或熱力學(xué)計(jì)算	網(wǎng)頁

» 本帖已獲得的紅花（最新10朵）

zhang_jaj

» 猜你喜歡

生醫(yī)工0831調(diào)劑求推薦已經(jīng)有4人回復(fù)
調(diào)劑已經(jīng)有4人回復(fù)
一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! 已經(jīng)有9人回復(fù)
求調(diào)劑已經(jīng)有8人回復(fù)
086000生物與醫(yī)藥初試274求調(diào)劑已經(jīng)有3人回復(fù)
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行已經(jīng)有18人回復(fù)
318求調(diào)劑已經(jīng)有9人回復(fù)
08工科，295，接受跨專業(yè)調(diào)劑已經(jīng)有4人回復(fù)
307分求調(diào)劑已經(jīng)有4人回復(fù)
085601 329分調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

Monte carlo模擬已經(jīng)有5人回復(fù)
Monte Carlo 模擬時(shí)圖是用什么軟件做的已經(jīng)有10人回復(fù)
【求助】是不是monte carlo 模擬只適用于高分子方向的？已經(jīng)有6人回復(fù)
【求助】求如何學(xué)習(xí)monte carlo仿真已經(jīng)有7人回復(fù)
【求助】用monte carlo法做分子模擬，用什么編程語言實(shí)現(xiàn)比較好？已經(jīng)有25人回復(fù)
【資源】Introduction to Monte Carlo Method 已經(jīng)有7人回復(fù)
【討論】大家對(duì)未來GPU在monte carlo中的應(yīng)用持有什么樣的態(tài)度(CUDA)？已經(jīng)有15人回復(fù)
【轉(zhuǎn)帖】一道面試題的monte carlo實(shí)現(xiàn) 已經(jīng)有8人回復(fù)
【求助】想用Monte Carlo方法預(yù)測(cè)晶體磁性轉(zhuǎn)變溫度應(yīng)該用什么程序已經(jīng)有5人回復(fù)
【求助】有用Monte Carlo 編過有關(guān)Fluid中粒子運(yùn)動(dòng)的程序已經(jīng)有6人回復(fù)
【求助】求Monte Carlo 關(guān)于擴(kuò)散系數(shù)計(jì)算的源代碼已經(jīng)有5人回復(fù)

1樓 2009-12-21 16:53:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

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é)

★ ★ ★
lei0736(金幣+3,VIP+0):謝謝 12-22 21:44

引用回帖:

Originally posted by lei0736 at 2009-12-22 10:01:
在MD模擬里有平衡和非平衡的MD
那么有沒有平衡MC和非平衡MC之分？
他們之間有哪些區(qū)別與聯(lián)系？是否非平衡MC不需要滿足細(xì)致平衡條件？
如果有區(qū)別他們各自應(yīng)用的領(lǐng)域特點(diǎn) 又是什么？

沒聽說過非平衡MC，所以我不確定有沒有。不過有所謂的kinetic MC，有人用來模擬硬球體系的結(jié)晶成核速率。至于是否需要滿足細(xì)致平衡條件，我想應(yīng)該還是要滿足的。

贊一下(11人)

回復(fù)此樓

7樓2009-12-22 20:52:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 74 個(gè)回答

lei0736

榮譽(yù)版主 (職業(yè)作家)

模擬EPI: 5
應(yīng)助: 17 (小學(xué)生)
貴賓: 10.32
金幣: 19299
散金: 5527
紅花: 30
沙發(fā): 1
帖子: 4768
在線: 1085.3小時(shí)
蟲號(hào): 62582
注冊(cè): 2005-03-17
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 量子化學(xué)

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

祝賀開貼
順便問下細(xì)致平衡條件的重要性如何
以及如何在模擬中注意哪些事項(xiàng)以保證細(xì)致平衡條件

贊一下(1人)

回復(fù)此樓

2樓2009-12-21 18:55:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

qzhaosdu

金蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1401.6
散金: 2052
紅花: 7
帖子: 1969
在線: 263.5小時(shí)
蟲號(hào): 493351
注冊(cè): 2008-01-10
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 量子化學(xué)

★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
lei0736(金幣+1):歡迎參與討論 2010-01-30 17:00

請(qǐng)問用MC做表面活性劑自組裝和高分子稀溶液自組裝的可靠性？
看過很多文章，貌似用MC做同樣的體系和用SCFT,DPD做出來的結(jié)論大不相同。謝謝！

贊一下(2人)

回復(fù)此樓

Anewday,anewhour,anewminute,anewpeople.

3樓2009-12-21 21:13:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

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é)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金幣+10,VIP+0):呵呵值得推薦 12-22 09:57
zeoliters(金幣+5,VIP+0):謝謝！ 12-25 16:25

引用回帖:

Originally posted by lei0736 at 2009-12-21 18:55:
祝賀開貼
順便問下細(xì)致平衡條件的重要性如何
以及如何在模擬中注意哪些事項(xiàng)以保證細(xì)致平衡條件

Detailed balance 用公式表示為
$P(\mathrm{o}) \alpha(\mathrm{o}\to\mathrm{n}) P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})= P(\mathrm{n}) \alpha(\mathrm{n}\to\mathrm{o}) P_\mathrm{acc}(\mathrm{n}\to\mathrm{o})$
其中， $P(\mathrm{o})$ 表示微觀態(tài)或微觀構(gòu)象 $\mathrm{o}$ 的概率(或幾率)密度分布，在NVT系綜中 $P(\mathrm{o})=\frac{e^{-\beta U(\mathrm{o})}}{\int\mathrmibjahkw\Gamma e^{-\beta U(\Gamma)}}$ ，即Boltzmann分布。 $\alpha(\mathrm{o}\to\mathrm{n})$ 表示從構(gòu)象 $\mathrm{o}$ 到構(gòu)象 $\mathrm{n}$ 所設(shè)計(jì)進(jìn)行的嘗試運(yùn)動(dòng)的幾率，而 $P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})$ 則是接受從構(gòu)象 $\mathrm{o}$ 到構(gòu)象 $\mathrm{n}$ 的嘗試運(yùn)動(dòng)的幾率，因此 $\alpha(\mathrm{o}\to\mathrm{n}) P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})$ 表示從構(gòu)象 $\mathrm{o}$ 到構(gòu)象 $\mathrm{n}$ 的轉(zhuǎn)移幾率。用簡(jiǎn)單的語言描述，Detailed balance 是指從構(gòu)象 $\mathrm{o}$ 到構(gòu)象 $\mathrm{n}$ 的流量(flow)等于從構(gòu)象 $\mathrm{n}$ 到構(gòu)象 $\mathrm{o}$ 的流量，這也就是所謂的微觀可逆性。

Detailed balance 的重要性在于 (1) 保證在有限的MC步數(shù)內(nèi)達(dá)到所需的平衡概率分布[c.f. M.E.J. Newman & G.T. Barkema, Monte Carlo Methods in Statistical Physics, Oxford]，這一點(diǎn)有嚴(yán)格的數(shù)學(xué)證明，我們只需要知道這個(gè)結(jié)論就可以了； (2) 達(dá)到平衡后平衡就不再被破壞，這意味著根據(jù)MC算法產(chǎn)生的構(gòu)象將滿足平衡概率分布。

滿足Detailed balance 公式的 $P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})$ 解有很多，使用最為廣泛的當(dāng)屬M(fèi)etropolis 算法：

$P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})=\mathrm{min}\left\{1.0, \frac{P(\mathrm{n}) \alpha(\mathrm{n}\to\mathrm{o})}{P(\mathrm{o}) \alpha(\mathrm{o}\to\mathrm{n})} \right\}$
這里， $P(\mathrm{o})/P(\mathrm{n})$ 是已知的，對(duì)于NVT系綜，
$\frac{P(\mathrm{o})}{P(\mathrm{n})}=e^{-\beta\bigl\{U(\mathrm{n})-U(\mathrm{o})\bigr\}}$ 。而 $\alpha(\mathrm{o}\to\mathrm{n})$ ( $\alpha(\mathrm{n}\to\mathrm{o})$ )取決于我們?cè)O(shè)計(jì)的Monte Carlo運(yùn)動(dòng)(move)是如何由構(gòu)象 $\mathrm{o}$ ( $\mathrm{n}$ )到構(gòu)象 $\mathrm{n}$ ( $\mathrm{o}$ )的。例如對(duì)于簡(jiǎn)單的粒子體系，我們現(xiàn)在設(shè)計(jì)隨機(jī)移動(dòng)一個(gè)粒子的位置，偽代碼為

int ichosen = int(Ntot * mt19937.getReal3()); // Ntot: 粒子的數(shù)目; mt19937.getReal3()返回(0, 1)之間的隨機(jī)數(shù)
double stepsize = 0.5 * sigma_LJ; //假設(shè)LJ 粒子體系,  sigma_LJ為L(zhǎng)J直徑. 在實(shí)際模擬中一般需要調(diào)節(jié)stepsize的大小以獲得一定的嘗試運(yùn)動(dòng)成功率, 例如30%.
for (int j = 0; j < 3; ++j)
  xn[ichosen][j] = xo[ichosen][j] + stepsize * (mt19937.getReal3()-0.5); //(2.0*mt19937.getReal3()-1.0)返回(-0.5,0.5)的隨機(jī)值.

這段代碼的for循環(huán)實(shí)現(xiàn)的是以隨機(jī)選取的粒子的原有位置xo[ichosen]為中心所形成的邊長(zhǎng)為 stepsize的正方體內(nèi)隨機(jī)產(chǎn)生一個(gè)新的位置xn[ichosen]，那么， $\alpha(\mathrm{o}\to\mathrm{n}) \propto \mathrm{stepsize}^{-3}$ ，因?yàn)閟tepsize * (mt19937.getReal3()-0.5) 保證了由舊位置到新位置的嘗試運(yùn)動(dòng)幾率等于由新位置到舊位置的嘗試運(yùn)動(dòng)幾率，所以這里 $\alpha(\mathrm{n}\to\mathrm{o})=\alpha(\mathrm{o}\to\mathrm{n})$ 。有了這樣設(shè)計(jì)的嘗試運(yùn)動(dòng)再加上利用 Metropolis 算法 $P_\mathrm{acc}(\mathrm{o}\to\mathrm{n})=\mathrm{min}\left\{1.0, P(\mathrm{n})/ P(\mathrm{o})\right\}$ 來判斷是否接受這樣的嘗試運(yùn)動(dòng)，達(dá)到平衡后的構(gòu)象分布就已經(jīng)滿足我們想要抽樣的統(tǒng)計(jì)概率 $P(\Gamma)$ 分布了。

簡(jiǎn)單一句話， 通過Detailed balance設(shè)計(jì)MC嘗試運(yùn)動(dòng)，保證我們能達(dá)到平衡而且平衡后的構(gòu)象滿足我們想要的概率分布。

仍然以NVT系綜為例，任意可測(cè)量量 $\mathcal{O}$ 的統(tǒng)計(jì)定義為
$\langle \mathcal{O} \rangle_\mathrm{NVT} = \frac{\int\mathrmopq8riv\Gamma\mathcal{O}(\Gamma) e^{-\beta U(\Gamma)}}{\int\mathrm6gl8ieh\Gamma e^{-\beta U(\Gamma)}} =\sum\limits_{i=1}^{N_\mathrm{sample}} P_\mathrm{NVT}(\Gamma_i)\mathcal{O}(\Gamma_i)$
需要注意的是：我們已經(jīng)知道  達(dá)到平衡后的構(gòu)象分布就已經(jīng)滿足我們想要抽樣的統(tǒng)計(jì)概率 $P(\Gamma)$ 分布了，這就表明在MC中，
$\langle \mathcal{O} \rangle_\mathrm{NVT} = \sum\limits_{i=1}^{N_\mathrm{sample}} \mathcal{O}(\Gamma_i)$

因?yàn)镸C所產(chǎn)生的平衡構(gòu)象 $\Gamma_i$ 已經(jīng)滿足平衡分布 $P_\mathrm{NVT}(\Gamma_i)$ 了。

[ Last edited by yahoohoo on 2009-12-22 at 04:21 ]

贊一下(20人)

回復(fù)此樓

4樓2009-12-21 22:28:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求收留 +8	1943443204 2026-03-28	8/400	2026-03-31 15:00 by -迷了路啊路
[考研] 一志愿中海洋材料357 +4	麥恩莉. 2026-03-30	4/200	2026-03-31 14:35 by 記事本2026
[考研] 284求調(diào)劑 +7	小熊～～ 2026-03-31	7/350	2026-03-31 13:39 by anny19840123
[考研] 070300化學(xué)354求調(diào)劑 +13	101次希望 2026-03-28	13/650	2026-03-31 13:38 by anny19840123
[考研] 土木304求調(diào)劑 +4	兔突突突， 2026-03-31	5/250	2026-03-31 11:29 by 北風(fēng)之神.
[考研] 335求調(diào)劑 +3	321* 2026-03-31	3/150	2026-03-31 09:00 by 功夫瘋狂
[考研] 08工科求調(diào)劑286 +5	tgs_001 2026-03-28	5/250	2026-03-31 08:18 by 一只好果子?
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5	小張zxy張 2026-03-26	10/500	2026-03-31 00:29 by jp9609
[考研] 281求調(diào)劑 +5	亞克西good 2026-03-26	7/350	2026-03-30 20:42 by dophin1985
[考研] 293求調(diào)劑 +3	末未mm 2026-03-30	5/250	2026-03-30 17:23 by 王保杰33
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7	Shern—- 2026-03-28	7/350	2026-03-30 16:31 by nothing投稿中
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13	丹青奶蓋 2026-03-26	15/750	2026-03-30 12:35 by fangnagu
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7	前行必有光 2026-03-28	8/400	2026-03-28 23:22 by 小木蟲tim
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�
[考研] 085405 考的11408求各位老師帶走 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 09:19 by 樂呵呵的追夢(mèng)人
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 305求調(diào)劑 +5	哇盧卡庫 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 085602化學(xué)工程求調(diào)劑。 +4	平樂樂樂 2026-03-26	4/200	2026-03-26 17:57 by fmesaito
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5	超贊的 2026-03-24	5/250	2026-03-24 21:07 by 星空星月