第6章蒙特卡羅方法(3)_第1頁
第6章蒙特卡羅方法(3)_第2頁
第6章蒙特卡羅方法(3)_第3頁
第6章蒙特卡羅方法(3)_第4頁
第6章蒙特卡羅方法(3)_第5頁
已閱讀5頁,還剩27頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領

文檔簡介

1、蒙特卡羅方法(三)6.3 Metropolis 算法算法目的:在積分變量 X 的空間(可能是多維的)內(nèi)產(chǎn)生一組按概率密度 (X ) 分布的點Metropolis 算法:假想有一個隨機行走者在 X 的空間中運動。該隨機行走過程相繼各步經(jīng)過的點產(chǎn)生出一個序列:X0, X1, ;隨著行走的路程越長,這些點的分布會逼近所要求的分布 (X ) 。注意在這里,行走不是完全無規(guī)的。每一步可能達到的位置依賴且僅依賴于上一步的位置。這種隨機行走可以用下列躍遷概率描述狀態(tài)ABCDA1/21/401/4B01/32/30C0100D1/201/20考慮大量行走者從不同的初始點出發(fā)在 X 空間獨立的隨機行走。若 Nn

2、(X ) 是 n 步后這些行走者在 X 點的密度, Nn(Y) 是 n 步后這些行走者在 Y 點的密度,那么在下一步從 X 點走到 Y 點的行走者凈數(shù)目為X 點上的一個行走者轉(zhuǎn)移到點上的一個行走者轉(zhuǎn)移到Y(jié) 的概率的概率當 時,行走者的布局不會有凈變化,系統(tǒng)將會達到平衡。這個條件稱為細致平衡條件。有多種具體方案給出躍遷概率,其中最重要的是 Metropolis方案細致平衡條件 并不能唯一的定出 躍遷概率為了使行走者的分布達到給定的概率密度分布 (X ) ,需要規(guī)定適當?shù)能S遷概率設行走者處于序列中的 Xn 點上,為了產(chǎn)生 Xn+1,行走者邁出試探性的一步到一新點 Xt 。 該新點可以用任何方便的方

3、法選取,如可以在點 Xn 周圍的一個邊長為 d 的多維立方體中均勻地隨機選取。然后按比值來決定是“接受”還是“拒絕”該試探步。Metropolis 的方案Metropolis提出的隨機行走的規(guī)則這樣產(chǎn)生出 X n+1 之后,可以再從X n+1 出發(fā)邁出一個試驗步,按照同樣的過程產(chǎn)生 X n+2 ,不斷重復以上步驟,得到整個隨機行走序列。如果 r 1,接受該步(即取 X n+1 = X t );如果 r 1,則以概率 r 接受這一步:把 r 和一個在0, 1區(qū)間上均勻分布的隨機數(shù) h 比較,若 h r 就接受這一步(即取 X n+1 = X t ) ,否則就舍棄該步(即取 X n+1 = X n

4、 )證明 Metropopis方案能導致平衡分布 (X ) 其中 T 是從 X 試探到 Y 的概率,A 是接受這一試探步的概率。根據(jù)Metropopis方案,從 X 到 Y 的轉(zhuǎn)移概率為因此,Metropopis隨機行走者的平衡分布滿足首先試探概率滿足若 (X ) (Y ),則 A (YX )=1 , 細致平衡條件并沒有定出躍遷概率的具體形式,因此除了 Metropolis 方案外,還存在其它一些選擇,例如 Barker 方法容易證明 Barker 方法確實滿足細致平衡條件。其它的選擇3. 構(gòu)成隨機行走的點 X0 , X1 , X2 , ,由于產(chǎn)生的方法,彼此不是獨立的。1. 隨機行走從何出發(fā)

5、,即選擇何處為 X0 。原則上,任何位置都合適,但實際中,合適的起始點是 值大的地方。Metropolis 方法要注意的地方太小了關聯(lián)太強,需要很多步才能達到一個統(tǒng)計獨立的構(gòu)型。太大了會使大多數(shù)嘗試都失敗,導致更新緩慢。經(jīng)驗的做法是使得大約的一半試探步入選。2. 步長 d 如何選取?關聯(lián)函數(shù)實際做蒙特卡羅模擬時,需要在生成的點列中進行采樣。相鄰兩個樣本點之間應有足夠大的的間隔,使得 C(l) 足夠小。關聯(lián)函數(shù)定義為自相關函數(shù)圖例其中 R=(r1, r2, , rN ) 為3N 維坐標矢量,假設粒子間為對勢應用實例:單原子氣體Metropolis 算法最早被用在經(jīng)典流體的的模擬中( Metrop

6、olis, 1953)。這里討論最簡單的經(jīng)典流體單原子氣體??紤]一定體積,溫度為T 的單原子氣體,這是一個正則系綜, 其物理量 A 的期望值為這里沒有考慮 U 對速度的依賴,因為速度的分布由麥克斯韋分布給出。 任何依賴于速度的物理量可以直接用該分布解析的計算。下面只考慮依賴于位形的物理量, 如總勢能等。 其中Ri 為相空間中依據(jù)下面的分布函數(shù)抽樣得到的樣本點用蒙特卡洛方法計算 的方法為其中 hk 為第 k 個方向的步長,、 為 0, 1 區(qū)間的隨機數(shù)很多時候,更新一次構(gòu)型的只改變一個粒子的坐標效率較高,特別是當系統(tǒng)非常接近平衡態(tài)或接近相變時。構(gòu)型的更新隨機的選擇第 i 個粒子,其坐標被更新為這

7、步更新被接受的概率為新構(gòu)型中,只有第 i 個粒子的坐標變化了,因此沒有必要為了得到概率 p 而計算整個 (Rn+1)我們將 U 表示為新舊構(gòu)型的能量差. 構(gòu)造系統(tǒng)的一個初始態(tài). 隨機選擇一個粒子 i,產(chǎn)生一個試探步 ri=ri+. 計算這一試探位移引起的能量變化 E. 如果 E 0,生成一個隨機數(shù)r, 滿足0r1. 如果 r exp(- E/KT),接受這一試探步;否則拒絕. 執(zhí)行第 (2) 步。 單原子氣體的具體算法為了有效的降低計算量,需要對位勢引入截斷, rij rc 。勢函數(shù)的截斷其截斷長度可取為 rc =3,其中 為勢函數(shù)為 0 時的距離長程相互作用, 例如庫侖勢,不能被截斷,需要額

8、外的方法來處理。例如簡單流體中典型的相互作用為 Lennard-Jones 勢 步長 h 如果太小: 粒子只能移動很小的距離, 所以需要很多步才能達到一個統(tǒng)計獨立的構(gòu)型。步長 h 如果太大:粒子雖然平均來說可以移動更大的距離,但是絕大多數(shù)移動對系統(tǒng)構(gòu)型的改變過大,以至于能量會顯著的增加。 這導致拒絕率過高,從而構(gòu)型的更新緩慢。步長的選擇一個經(jīng)驗法則是接受率平均在 0.4 和 0.6 之間。 對于硬球, 接受率應該低些,約為 0.1。 用來取平均的樣本點之間典型的應有 10-15 個間隔樣本點。經(jīng)典的二維 Ising 模型:在一個二維正方格子上,每個格點 i 上有一個自旋,可以取值 +1 或 -

9、1。相鄰自旋通過一個交換耦合能 J 相互作用,此外還存在一個外磁場 B。其中 代表對最近鄰求和。應用實例Ising 模型系統(tǒng)的哈密頓量為伊辛模型最早被用來研究磁相變,另一個有趣的應用是二元合金。零溫下,當沒有外磁場時 J 0 : 如果所有自旋都朝同一方向,系統(tǒng)能量會最低, 對應鐵磁態(tài)J0:如果相鄰兩個自旋朝向相反,系統(tǒng)能量會最低。對應反鐵磁態(tài)交換耦合能 J二維 Ising 模型示意圖周期性邊條件二維 Ising 模型不同溫度下典型的自旋構(gòu)型二維 Ising 模型的解析結(jié)果三維伊辛模型仍然缺乏精確解Lars Onsager 于1944 年得到二維 Ising 模型的精確解,并證明存在相變點,這是統(tǒng)計物理學發(fā)展過程中的里程碑。其臨界指數(shù)為磁化強度磁化率內(nèi)能比熱需要計算的物理量正則系綜下磁化強度的平均值被定義為磁化強度可以通過蒙特卡羅方法計算其中 =1,2,M 根據(jù)下面的分布函數(shù)取樣正則系綜下物理量的計算其中 為構(gòu)型 的平均自旋。注意與每一個格點相聯(lián)系的 可以儲存并在模擬中更新。另外,接受概率只有 5 種可能的值,可以事先計算并儲存起來,以避免重復進行指數(shù)運算。隨機的在所有格點生成 1 或 -1, 然后隨機的選擇一個格點來更新它, 接受率為其中j 對 i 的所有鄰居求和。.

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論