 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1203 | 回復(fù): 4 | ||||
[求助]
蒙特卡洛模擬中，組態(tài)躍遷幾率的兩種表達(dá)式是怎么來的呢 已有2人參與
|
|
最近在做用蒙特卡洛ising模型的模擬，在文獻(xiàn)上看到有兩種關(guān)于躍遷幾率的表達(dá)式 一種是比較直觀的bolzman因子形式，即：p=exp(-E/KT)，當(dāng)E<0時(shí)，p=1 其中E代表躍遷前后的能量差 另一種則是我的疑惑所在了，其表達(dá)式為：P=p/(1+p) 其中小寫p=exp(-E/KT)。 書上的原話是：在統(tǒng)計(jì)力學(xué)的細(xì)致平衡原理中，沒有唯一決定躍遷幾率，因此采用metropolis抽樣法時(shí)還可以用其他適當(dāng)?shù)臍w一化幾率分布，這不會(huì)改變平衡態(tài)結(jié)果，但可以加速高溫下的收斂速度。 請(qǐng)教各位前輩們，在此謝過~ [ Last edited by 月只藍(lán) on 2014-7-30 at 13:41 ] |
回復(fù)此樓» 猜你喜歡
材料與化工求調(diào)劑
已經(jīng)有7人回復(fù)
-
304求調(diào)劑
已經(jīng)有3人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有12人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有12人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
324分 085600材料化工求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有6人回復(fù)
-
307求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有3人回復(fù)
2樓2014-07-30 11:55:32
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲號(hào): 556471
- 注冊(cè): 2008-05-11
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
潯荊(月只藍(lán)代發(fā)): 金幣+3, 感謝指導(dǎo)！ 2014-07-30 16:53:00
潯荊: 金幣+7, ★★★很有幫助 2014-07-30 22:19:54
感謝參與，應(yīng)助指數(shù) +1
潯荊(月只藍(lán)代發(fā)): 金幣+3, 感謝指導(dǎo)！ 2014-07-30 16:53:00
潯荊: 金幣+7, ★★★很有幫助 2014-07-30 22:19:54
|
建議你把Understanding Molecular Simulation里關(guān)于Monte Carlo的介紹（3.1節(jié)）仔細(xì)閱讀一下，在3.1.2節(jié)中講到了Detailed balance principle，detailed balance 無非就是 P_i * T(i->j) = P_j *T(j->i) T(i->j) 和T(j->i)在這里就是transition probability，在正則系綜中,P_i=exp(-E_i/kT)，正如書中所說，在Monte Carlo模擬中，transition probability的選擇可以有多種形式，只要它們滿足由detailed balance principle所決定的比例關(guān)系(可能還有些額外要求，具體細(xì)節(jié)我不太了解) T(i->j) /T(j->i)=P_j/P_i = exp(-(E_j - E_i)/kT)， 那么你最后得到的概率分布就會(huì)收斂到玻爾茲曼分布。采用不同的transition probability會(huì)導(dǎo)致不同的收斂速度，這個(gè)就是個(gè)很技術(shù)的問題了。 所以你列出的兩個(gè)transition probability都滿足detailed balance principle的比例關(guān)系，所以都是可以的。 |
3樓2014-07-30 13:14:25
4樓2014-07-30 22:11:33
5樓2014-07-30 22:19:39
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 潯荊 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +3 | Ncdx123456 2026-03-19 | 3/150 |
|
|
[考研] 288求調(diào)劑 +6 | 于海海海海 2026-03-19 | 6/300 |
|
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
|
[考研]
 0854，計(jì)算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|