伊辛模型

描述物質相變的隨機過程模型

伊辛模型(Ising model)是一類描述物質相變的隨機過程(stochastic process)模型。物質經過相變,要出現新的結構和物性。發生相變的系統一般是在分子之間有較強相互作用的系統,又稱合作系統。

伊辛模型所研究的系統由多維周期性點陣組成,點陣的幾何結構可以是立方的或六角形的,每個陣點上都賦予一個取值表示自旋變數,即自旋向上或自旋向下。伊辛模型假設只有最近鄰的自旋之間有相互作用,點陣的位形用一組自旋變數來確定。常見的二維伊辛模型示意圖使用箭頭方向表示自旋方向

歷史


伊辛模型由德國物理學家威廉·楞次(Wilhelm Lenz)在1920年提出以描述鐵磁性物質的內部的原子自旋狀態及其與宏觀磁矩的關係。1924年,楞次的學生Ernst Ising求解了不包含相變的一維伊辛模型。20世紀30-40年代,勞倫斯·布拉格(Lawrence Bragg)、E. J. Williams、漢斯·貝特(Hans Bethe)、Rudolf Peierls等學者使用平均場近似理論(mean-field theory)對二維伊辛點陣模型(two-dimensional square-lattice Ising model)進行了研究。1944年美國物理學家拉斯·昂薩格(Lars Onsager)得到了二維伊辛模型在沒有外磁場時的解析解,即Onsager解。

理論


模型

這裡以二維伊辛點陣模型(two-dimensional square-lattice Ising model)為例對伊辛模型進行說明。二維伊辛點陣模型是一個空間隨機場,其中任意點的狀態 可有 兩個取值,並僅受到與其鄰接的點的影響:
式中 為與 鄰接的點的集合,求和項 被稱為外部場(external field),是模型參數,當 時,隨機場所有點的符號(正負)相同、當 時,隨機場所有點的正負相反、當 時,隨機場間的點與其鄰接點之間沒有相互作用,按均等概率取正或負值。
圖1 伊辛模型
圖1 伊辛模型
伊辛模型通常被用於模擬鐵磁性物質(鐵、鈷、鎳)的結構並對其在鐵磁性狀態和非鐵磁性狀態之間的相變(phase transition)進行理論描述。當鐵磁性物質的溫度低於居里溫度(Curie temperature)時,其內部的原子會按特定方式自旋從而產生宏觀磁矩。對應伊辛模型中 或 的情形,前者表現為鐵磁性(ferromagnetic),後者表現為反鐵磁性(anti-ferromagtic)。當溫度高於居里溫度時,原子自旋的取向非常紊亂,不產生凈磁矩,對應 的情形。

平均場理論

Weiss分子場理論(Weiss molecular field theory)
Weiss分子場理論將伊辛模型近似為點 的概率分佈 的乘積:對上式右側取對數,並考慮到概率分佈 需要表示 所有鄰接點對其自身的影響,可得如下展開:
即概率分佈可表示為該點鄰接點的狀態的期望疊加一常數。帶入二維伊辛點陣模型的表達式后可得:
由於點的概率分佈 僅有正/負兩個取值: ,由此得到:
Bragg-Williams平均場理論(Bragg-Williams mean-field theory)
Bragg-Williams平均場理論認為,某一陣點上的自旋取某一方向的幾率同近鄰陣點上的自旋取向無關,只同自旋在該方向的數目成正比。每個陣點上有一平均磁場,自旋在陣點上的取向只同該磁場有關。用這種方法可求得下列公式:式中 是每個自旋的磁矩,是每一陣點的最近鄰數,是外磁場強度,是熱力學溫度,模型參數 在式中表示自旋同向的最近鄰對之間的相互作用能,是玻耳茲曼常數。每個自旋上的磁化強度可表示為:
式中 分別代表自旋向上和向下的總陣點數, 。
將上式應用於鐵磁性物質的性質,可以得到如下結論:取 為臨界溫度(居里溫度),當溫度高於臨界溫度且沒有外磁場時, ,物質不具有鐵磁性;當溫度低於居里溫度時, ,磁化強度取正值或負值,鐵磁性物質存在相變。這個結論對一維、二維、三維點陣都應成立,
研究表明,二維、三維伊辛模型在臨界溫度以上仍有相變,即平均場理論在鐵磁性物質處於臨界溫度以上時有局限性。

求解


一維情形

考慮具有N個自旋的直線鏈,每個自旋僅同它的兩個最近鄰自旋及外磁場相互作用。相互作用的總能量即由 所確定的位形能量是:
式中 是對最近鄰自旋對求和,表示對所有自旋坐標求和。由上式可得配分函數的形式為:
定義一個的矩陣,其矩陣元具有形式: 。因為點只有正/負兩個可能取值,該矩陣可表示為:
如果採用周期性邊界條件,即,或設想將直線鏈彎成閉合的圓鏈,並將初端與尾端相接,則配分函數 有如下表示:
式中 表示矩陣的跡。由此計算可以求出配分函數。另對於長鏈,有: 。即當外磁場為0時,對於磁化強度有。即一維伊辛模型沒有自發磁化。在一維伊辛模型中,不論順磁性或逆磁性,都不會實現有序的狀態。如對逆磁性的情形,在絕對零度時,所有自旋取向相同時處於能量最低的狀態。但如果熱力學溫度不等於零,是有限的,則平均位形由兩種相反的、相互競爭的趨向所決定。一個是各個自旋的取向完全一致,使能量最低;另一個是各個自旋的取向為隨機的,使熵最大。由於一維伊辛模型中每個自旋沒有足夠多的最近鄰自旋,因而不可能出現所有自旋取向完全相同的情況。

二維情形

Onsager解
20世紀40年代Lars Onsager對伊辛模型採用解析法得到了嚴格解,作出了突出的成就。這種方法的基本點,是設每個陣點的自旋變數可取+1和-1兩個值,考慮陣點上自旋的某個位形,計算每個自旋同最近鄰自旋的相互作用能量以及同外磁場相互作用能量,再對全部可能的位形求和,用矩陣的方法求出配分函數,從而得到各個熱力學函數。
二維伊辛點陣的陣點數為。處理二維空間問題的方法與一維的類似,只需將一維的每個陣點當作一列,並逐列相加求和即可。
以Sl表示第l行的所有自旋坐標的集合 上標l()代錶行,下標代表列。邊界條件為。即要求每一行的第列的位形與第1列的相同,每一,Sl有2n個值。整個點陣的位形由確定。考慮最近鄰自旋對以及自旋同外磁場的相互作用,則配分函數可寫成 為將上式表示成矩陣的形式,引入三個矩陣,它們的矩陣元分別定義為 第一式反映不同行最近鄰自旋對的相互作用能量,它有個;第二式反映同一行最近鄰自旋對的相互作用能量,它有2n個;第三式反映同一行各個自旋與外磁場的相互作用能量,它也有2n個。為了計算方便,在補上一些“0”元素后,可把、擴大成矩陣的對角矩陣。可以證明 。當時,求矩陣的本徵值問題就變成求解矩陣的本徵值問題。H.A.克喇末和G.H.萬尼爾等人曾用數字解計算過有限的幾項,他們計算到,發現當n為有限的情形下,沒有相變。昂薩格在求解時,設外磁場強度,因而。計算結果表明:高溫時,(臨界溫度),矩陣只有一個最大本徵值υ+;低溫時,,矩陣有兩個本徵值,當時,配分函數為 並得出平均每個自旋的自由能f為 若,則上式右邊第二項被積函數θ(v)滿足。用數值計算,通過上述二式可算出時的各個熱力學量,得到以下具體結果,
其中S為每個自旋的熵。式中的臨界溫度Tc滿足方程 或。在Tc附近,每個自旋的比熱容可表示為
可見在時,自由能、熵以及內能是連續的,這意味著在時,發生的相變不包含潛熱。但是時,作為上述熱力學函數的導數,比熱容是對數發散的,無論從高溫端還是低溫端趨於Tc(即或),比熱容с的值是相同的。
為弄清處相變的細節,還需進一步考慮自發磁化(即計算自由能對磁場強度H的導數,再讓)。楊振寧於1952年採用微擾法得到了很好的結果。他證明自發磁化強度可表為 式對應於Tc的值。

其它

存在外磁場的二維伊辛模型和更高維情形的伊辛模型使用平均場近似理論進行研究,其解析解沒有被發現,或被認為沒有解析解。