 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1206 | 回復: 4 | |||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||
[求助]
蒙特卡洛模擬中，組態(tài)躍遷幾率的兩種表達式是怎么來的呢 已有2人參與
|
|||||
|
最近在做用蒙特卡洛ising模型的模擬，在文獻上看到有兩種關(guān)于躍遷幾率的表達式 一種是比較直觀的bolzman因子形式，即：p=exp(-E/KT)，當E<0時，p=1 其中E代表躍遷前后的能量差 另一種則是我的疑惑所在了，其表達式為：P=p/(1+p) 其中小寫p=exp(-E/KT)。 書上的原話是：在統(tǒng)計力學的細致平衡原理中，沒有唯一決定躍遷幾率，因此采用metropolis抽樣法時還可以用其他適當?shù)臍w一化幾率分布，這不會改變平衡態(tài)結(jié)果，但可以加速高溫下的收斂速度。 請教各位前輩們，在此謝過~ [ Last edited by 月只藍 on 2014-7-30 at 13:41 ] |
回復此樓» 猜你喜歡
求助
已經(jīng)有4人回復
-
0805材料320求調(diào)劑
已經(jīng)有3人回復
-
297求調(diào)劑
已經(jīng)有11人回復
-
268求調(diào)劑
已經(jīng)有9人回復
-
一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有10人回復
-
278求調(diào)劑
已經(jīng)有7人回復
-
26考研一志愿中國石油大學(華東)305分求調(diào)劑
已經(jīng)有5人回復
-
280求調(diào)劑
已經(jīng)有9人回復
-
279求調(diào)劑
已經(jīng)有5人回復
-
286求調(diào)劑
已經(jīng)有5人回復
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗: +49 - 模擬EPI: 2
- 應助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時
- 蟲號: 556471
- 注冊: 2008-05-11
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 分子模擬
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
潯荊(月只藍代發(fā)): 金幣+3, 感謝指導！ 2014-07-30 16:53:00
潯荊: 金幣+7, ★★★很有幫助 2014-07-30 22:19:54
感謝參與，應助指數(shù) +1
潯荊(月只藍代發(fā)): 金幣+3, 感謝指導！ 2014-07-30 16:53:00
潯荊: 金幣+7, ★★★很有幫助 2014-07-30 22:19:54
|
建議你把Understanding Molecular Simulation里關(guān)于Monte Carlo的介紹（3.1節(jié)）仔細閱讀一下，在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)系(可能還有些額外要求，具體細節(jié)我不太了解) T(i->j) /T(j->i)=P_j/P_i = exp(-(E_j - E_i)/kT)， 那么你最后得到的概率分布就會收斂到玻爾茲曼分布。采用不同的transition probability會導致不同的收斂速度，這個就是個很技術(shù)的問題了。 所以你列出的兩個transition probability都滿足detailed balance principle的比例關(guān)系，所以都是可以的。 |
3樓2014-07-30 13:14:25
2樓2014-07-30 11:55:32
4樓2014-07-30 22:11:33
5樓2014-07-30 22:19:39
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 336求調(diào)劑 +3 | rmc8866 2026-03-21 | 3/150 |
|
|
[考研] 材料學學碩080502 337求調(diào)劑-一志愿華中科技大學 +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 材料 271求調(diào)劑 +3 | 展信悅_ 2026-03-21 | 3/150 |
|
|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +6 | 梨花珞晚風 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿985，本科211，0817化學工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|