以蒙地卡罗模拟法研究临界现象.pdfVIP

以蒙地卡羅模擬法研究臨界現象 以蒙地卡羅模擬法研究臨界現象 胡進錕 摘要 1.引言 2.硬核粒子的相變與展透性變遷 3.階層蒙地卡羅模擬法 4.晶胞至晶胞階層蒙地卡羅重整群算法 5.邊界條件及晶格形狀與普適尺度函數 6.討論與結語 本文評述以蒙地卡羅模擬法研究臨界現象的重要發展。主題包括： 1. 用蒙地卡羅模擬法發現硬核粒子的相變是展透性變遷； 2. 階層蒙地卡羅模擬法主要觀念及計算程序； 3. 階層蒙地卡羅重整群法及其應用； 4. 用階層蒙地卡羅模擬法發現有限大小尺度函數和邊界條件與晶格形狀的關係以及找出普適尺度函數的程 序。 1.引言 　　楊振寧除了在粒子物理與場論很有貢獻以外，在統計物理的臨界現象理論也 頗有貢獻。臨界現象 道認為它可以做為液體和氣體之間相 變的模型 1。1952年楊振寧發表二度空間方型晶格易行相變模型 自發磁性強度的 準確解2，他說這是他所從事極艱巨的計算工作之一。他發現在臨界溫度Tc 附近 ，當 而M(T)是易形模型的秩序參數。 　　楊振寧所研究的易行模型，水平方向自旋之間的耦合常數J1與垂直方向自旋 之間的耦合常數J2相 同，即J1=J2。他建議張承修研究J1胡2的情形。1952年張承 修發表一篇短文3，他發現當　0 　　當溫 度T趨近於臨界溫度 Tc，比熱Ch、磁化率c、相關長度x等物理量會趨近 於無限大。為了表示這些物 理量趨近於無限大的行為，科學家乃定義了臨界指數 。令t = (T-Tc)/Tc，則當溫度由TTc接近 Tc 磁化率等 物理量，可以由自由能的偏微分求得。1965年，Widom 假設自由能在臨界點附近 包含一個 以蒙地卡羅模擬法研究臨界現象 　　由這些等式，自然可以得到前述的尺度關係。因此理論上，只要會算yt和yh 就可以得到各種臨界 指數。 二階相變模型可 得到頗準確的臨界點和臨界指數，但是用於一階相變模型則無法得到準確的臨界 點。 準確算法可得到準確的臨界點和臨界指數，但通常只有一些簡單的一度空間 和二度空間的模型才有準 確解。1972 年楊振寧曾說5：“我認為找到不平常的二 度空間量子統計力學模型或三度空間古典統計 力學模型的準確解，是極為遙遠的 。”從1972年到現在，二十多年過去了，楊振寧當年講的話，到現 在基本上還是 對的。但是二十多年來，電腦的硬體和數值計算的算則有極大的進步。現在我們 可以用 數值計算的方法，極準確地算出三度空間相變模型的臨界點、臨界指數、 秩序參數等物理量，也可以 藉由電腦的幫助，找出自然界的規律性。本文的目的 ，就是簡單的介紹我們在這方面研究所得到的一 些成果。我們的數值計算用到蒙 地卡羅模擬法以及重整群法。 　　本文下面各節結構如下：第 2節簡要地介紹我們用蒙地卡羅模擬法發現硬核 粒子的相變是展透性 變遷；第3 節介紹我們發展的階層蒙地卡羅模擬法並用它算 出位置展透模型的存在或然率和展透或然 率。第4 節介紹我們將階層蒙地卡羅模 擬法與重整群法結合成階層蒙地卡羅重整群法並用它算出相變 模型極準確的臨界 點、臨界指數、與秩序參數；第5 節介紹我們用階層蒙地卡羅模擬法發現展透模 型 的有限大小尺度函數和邊界條件及晶格形狀有密切關係。適當的選擇每一晶格 的寬高比及非普適參數 2.硬核粒子的相變與展透性變遷 　　在硬核粒子模型，一個晶格點被某一原子占了以後，與這晶格點相距某一範 圍內（如最鄰近，或 最鄰近與次鄰近）的晶格點，都不容許有原子占住。滿足前 述條件之下，兩最鄰近原子交互作用能量 為-u。某一原子周圍不能被別的原子占 住的範圍可用點線表示。平面方型晶格上最鄰近排斥的硬核粒 角晶格上「最鄰近排斥」的硬核粒子模型相當於硬六邊形的排列問題，因此 又叫硬六邊形模型(Hard pc，Fn 和Fu都趨近於無限大；p由pc繼續變大，Fn和Fu則又逐漸變小。 有1-p的機率不被接通，凡是被連通鍵 接在一起的點就是同一集團。由晶格一邊伸展至另一邊的集團叫 展透集團，其他 的集團叫非展透集團。非展透集團的平均大小 S，叫集團平均大小。晶格點屬於 展透 集團的比率叫展透或然率，以P表示。兩點相距r而且在同一非展透集團的或 然率G，叫連接函數，它 數。 為。 以蒙地卡羅模擬法研究臨界現象 　　若G上每一點有ps的機率被粒子占住，有1-ps 的機率空著且兩個最鄰近的粒 子視為同一集團，則 態出現的機率，除了和