本帖產(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í)
蟲(chóng)號(hào): 74894
注冊(cè): 2005-06-15
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

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

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

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

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

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

歡迎大家多提寶貴意見(jiàn)。

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

回復(fù)此樓

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

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

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

zhang_jaj

» 猜你喜歡

08工科275分求調(diào)劑已經(jīng)有13人回復(fù)
材料工程322分已經(jīng)有6人回復(fù)
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑已經(jīng)有7人回復(fù)
274求調(diào)劑已經(jīng)有5人回復(fù)
【求調(diào)劑】新能源材料本科，一志愿211，初試321 已經(jīng)有3人回復(fù)
材料考研調(diào)劑已經(jīng)有3人回復(fù)
材料調(diào)劑已經(jīng)有12人回復(fù)
英一數(shù)一408，總分284，二戰(zhàn)真誠(chéng)求調(diào)劑已經(jīng)有14人回復(fù)
085410 一志愿211 22408分?jǐn)?shù)359求調(diào)劑已經(jīng)有4人回復(fù)
271求調(diào)劑已經(jīng)有19人回復(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法做分子模擬，用什么編程語(yǔ)言實(shí)現(xiàn)比較好？已經(jīng)有25人回復(fù)
【資源】Introduction to Monte Carlo Method 已經(jīng)有7人回復(fù)
【討論】大家對(duì)未來(lái)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ò)有關(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的回帖

coy123

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 949.1
散金: 423
帖子: 500
在線: 155.2小時(shí)
蟲(chóng)號(hào): 309796
注冊(cè): 2006-12-23
專業(yè): 理論和計(jì)算化學(xué)

★
lei0736(金幣+1):歡迎參與討論 2010-01-30 17:01

聽(tīng)說(shuō)過(guò)好多MC的種類。格點(diǎn)模型的，非格點(diǎn)模型的。已經(jīng)摟住說(shuō)的kinetic MC，但不知道 KMC是否都是格點(diǎn)模型的

贊一下(2人)

回復(fù)此樓

8樓2009-12-23 04:30:11

已閱回復(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í)
蟲(chóng)號(hào): 62582
注冊(cè): 2005-03-17
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 量子化學(xué)

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

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

贊一下(1人)

回復(fù)此樓

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

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

qzhaosdu

金蟲(chóng) (著名寫(xiě)手)

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

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

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

贊一下(2人)

回復(fù)此樓

Anewday,anewhour,anewminute,anewpeople.

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

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

yahoohoo

鐵桿木蟲(chóng) (著名寫(xiě)手)

模擬EPI: 10
應(yīng)助: 0 (幼兒園)
貴賓: 1.55
金幣: 7632.4
散金: 251
紅花: 4
帖子: 1176
在線: 167小時(shí)
蟲(chóng)號(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:
祝賀開(kāi)貼
順便問(wèn)下細(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\mathrmqlpkuoj\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)單的語(yǔ)言描述，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\}$ 來(lái)判斷是否接受這樣的嘗試運(yùn)動(dòng)，達(dá)到平衡后的構(gòu)象分布就已經(jīng)滿足我們想要抽樣的統(tǒng)計(jì)概率 $P(\Gamma)$ 分布了。

簡(jiǎn)單一句話， 通過(guò)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\mathrmvlkosrq\Gamma\mathcal{O}(\Gamma) e^{-\beta U(\Gamma)}}{\int\mathrmzcwha9q\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ā)表

[考研] 材料調(diào)劑 +9	一樣YWY 2026-04-01	9/450	2026-04-01 22:57 by peike
[考研] 295材料工程專碩求調(diào)劑 +19	1428151015 2026-03-27	19/950	2026-04-01 22:34 by peike
[考研] 0703化學(xué)338求調(diào)劑！ +8	Zuhui0306 2026-03-26	9/450	2026-04-01 21:59 by xhai2011
[考研] 273求調(diào)劑 +19	李芷新1 2026-03-31	19/950	2026-04-01 21:49 by chyhaha
[考研] 生物與醫(yī)藥考研調(diào)劑 +5	鐵憨憨123425 2026-03-31	5/250	2026-04-01 18:01 by syh9288
[考研] 化學(xué)308分調(diào)劑 +13	你好明天你好 2026-03-30	14/700	2026-04-01 16:19 by yanflower7133
[考研] 086000生物與醫(yī)藥初試274求調(diào)劑 +5	小叮當(dāng)來(lái)了 2026-03-30	6/300	2026-04-01 14:25 by mxj1203
[考研] 302求調(diào)劑一志愿北航070300，本科鄭大化學(xué) +8	圣日耳曼條 2026-04-01	10/500	2026-04-01 14:03 by ZXlzxl0425
[考研] 考研調(diào)劑 +11	Amber00 2026-03-31	11/550	2026-04-01 11:32 by wangjy2002
[考研] 考研生物與醫(yī)藥調(diào)劑 +7	鐵憨憨123425 2026-03-31	7/350	2026-04-01 08:45 by JourneyLucky
[考研] 353求調(diào)劑 +3	江上楓_26 2026-03-28	3/150	2026-03-31 15:53 by jp9609
[考研] 277跪求調(diào)劑 +8	1915668 2026-03-27	13/650	2026-03-31 14:58 by 王亮_大連醫(yī)科大
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3	080700調(diào)劑 2026-03-30	4/200	2026-03-31 12:00 by KLMY666
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3	utopiaE 2026-03-30	3/150	2026-03-30 23:42 by ms629
[考研] 328求調(diào)劑 +8	嗯滴的基本都 2026-03-27	8/400	2026-03-30 17:20 by Wang200018
[考研] 342求調(diào)劑 +4	加油a李zs 2026-03-26	4/200	2026-03-30 16:39 by 晶體之美
[考研] 291求調(diào)劑 +5	Y-cap 2026-03-29	6/300	2026-03-29 13:18 by mumin1990
[考研] 340求調(diào)劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無(wú)際的草原
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�
[考研] 08開(kāi)頭275求調(diào)劑 +4	拉誰(shuí)不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012