| 5 | 1/1 | 返回列表 |
| 查看: 1214 | 回復(fù): 4 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
蒙特卡洛模擬中,組態(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 ] |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
建議你把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)系,所以都是可以的。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 4/200 |
|
|
[考研] 299求調(diào)劑 +6 | 嗯嗯嗯嗯2 2026-03-27 | 6/300 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
[考研] 尋找調(diào)劑 +5 | 倔強(qiáng)芒? 2026-03-21 | 8/400 |
|
|
[考研] 調(diào)劑310 +3 | 溫柔的晚安 2026-03-25 | 4/200 |
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|