複製鏈接
請複製以下鏈接發送給好友

量子蒙特卡羅方法

鎖定
在Hartree-Fock理論中,很重要的一項誤差就是在Hartree-Fock理論中,電子之間的關聯項被全部忽略了。這其中的原因是因為我們將體系的基態波函數近似為單行列式的形式,電子所感受到的是平均場的作用,這使得我們可以將多體的波函數分解成為單電子波函數簡化計算。如果我們用無窮多個行列式的組合作為基函數,則理論上也可以得到精確的基態能量。而量子蒙特卡羅(Quantum Monte Carlo)則從一個全新的角度考慮問題。在量子蒙特卡羅中,體系的基態波函數顯式地寫成關聯的波函數,也就是説波函數是電子--電子之間距離的顯式函數。 [1] 
中文名
量子蒙特卡羅方法
所屬學科
電子學
性    質
電子學術語

目錄

量子蒙特卡羅方法實現

對於波函數的期望值的計算則通過蒙特卡羅方法實現。試想我們有N個電子的體系,則座標的自由度為3N,要在這麼大的維數的空間裏進行普通的積分是不可行的。在變分蒙特卡羅中,計算的精度完全由探試波函數的精度所決定。而在擴散蒙特卡羅中,則在計算過程中會通過投影方式改善試探波函數,從而提高計算精度。
我們知道,由於多體波函數中含有3N個電子座標的自由度,針對這麼多維空間的積分,普通的方法,如辛普森方法不再適用。因為這類積分方法的誤差對步長M通常按照收斂,其中d是積分的維數。所以當積分的維數非常大時(比如多體波函數的維數),一般的積分方法不再適用。量子蒙特卡羅,顧名思義,計算的積分方法用的是蒙特卡羅方法,其最重要的理論基礎就是中心極限定理(central limit theorem)。下面首先對中心極限定理做一介紹。

量子蒙特卡羅方法定理

假設我們用表示3N維空間中的一個矢量,在R點找到電子的概率用P(R)表示。顯見P(R)滿足如下的條件
我們接着定義一個依賴於電子座標R的函數f(R),則其平均值和方差應由如下公式給出
在實際計算中,由於dR涉及到針對很高維數的積分,所以通過定義(6.1.2)求解並不現實。為了近似,我們隨機選取一組互不關聯的滿足P(R)分佈的R,不妨記作R1,R2,R3,...,RM,並且計算函數f(R)針對這組{Rm: m=1,M}的平均值。
中心極限定理告訴我們,只要M取得足夠大,Z會以的速度趨向於,而和P(R)的具體形式無關,也和積分的維數無關。
那麼如何將此中心極限定理具體運用到求算符對多體基態波函數的期望值呢?比如説哈密頓的期望值就是基態能量。假設我們需要求哈密頓量對於基態波函數的期望值(注意,下式中的是不能分離變量的多體波函數)
其中叫做局域能量。而概率分佈則滿足方程(6.1.1)的限制條件。所以,我們可以套用中心極限定理,隨機產生一組滿足P(R)分佈的R,並且在每個R點上計算局域能量,則在取足夠多點的情況下,哈密頓量的期望值可以由局域能量的平均值近似。

量子蒙特卡羅方法其他信息

其他算符對於基態多體波函數的期望值也可以通過類似的方法計算。最後,我們指出,量子蒙特卡羅中的基態多體波函數通常是以Hartree-Fock或者密度泛函的波函數為基礎,加上一定的關聯修正得到。具體的如何產生基態多體波函數我們將在下一節介紹。通常,計算能力可以處理大到含有1000個電子的蒙特卡羅計算。
參考資料