 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1205 | 回復(fù): 4 | ||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
蒙特卡洛模擬中，組態(tài)躍遷幾率的兩種表達(dá)式是怎么來(lái)的呢 已有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) 其中小寫(xiě)p=exp(-E/KT)。 書(shū)上的原話是：在統(tǒng)計(jì)力學(xué)的細(xì)致平衡原理中，沒(méi)有唯一決定躍遷幾率，因此采用metropolis抽樣法時(shí)還可以用其他適當(dāng)?shù)臍w一化幾率分布，這不會(huì)改變平衡態(tài)結(jié)果，但可以加速高溫下的收斂速度。 請(qǐng)教各位前輩們，在此謝過(guò)~ [ Last edited by 月只藍(lán) on 2014-7-30 at 13:41 ] |
回復(fù)此樓» 猜你喜歡
336求調(diào)劑
已經(jīng)有3人回復(fù)
-
286求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿華中科技大學(xué)071000，求調(diào)劑
已經(jīng)有3人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有14人回復(fù)
-
332求調(diào)劑
已經(jīng)有3人回復(fù)
-
0856材料專(zhuān)碩353求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué)
已經(jīng)有5人回復(fù)
-
085601調(diào)劑 358分
已經(jīng)有4人回復(fù)
-
0856 271求調(diào)劑
已經(jīng)有3人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有15人回復(fù)
4樓2014-07-30 22:11:33
2樓2014-07-30 11:55:32
lsloneil
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲(chóng)號(hào): 556471
- 注冊(cè): 2008-05-11
- 性別: GG
- 專(zhuān)業(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 無(wú)非就是 P_i * T(i->j) = P_j *T(j->i) T(i->j) 和T(j->i)在這里就是transition probability，在正則系綜中,P_i=exp(-E_i/kT)，正如書(shū)中所說(shuō)，在Monte Carlo模擬中，transition probability的選擇可以有多種形式，只要它們滿(mǎn)足由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ù)的問(wèn)題了。 所以你列出的兩個(gè)transition probability都滿(mǎn)足detailed balance principle的比例關(guān)系，所以都是可以的。 |
3樓2014-07-30 13:14:25
5樓2014-07-30 22:19:39
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +3 | Faune 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 0856 271求調(diào)劑 +3 | 展信悅_ 2026-03-21 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚(yú) 2026-03-18 | 7/350 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 | ￣^￣゜汗 2026-03-19 | 4/200 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！ +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 085601專(zhuān)碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|