查看“易辛模型”的源代码

{{noteTA
|T=zh-hans:伊辛模型; zh-hant:易辛模型;
|1=zh-hans:蒙特卡洛; zh-hant:蒙地卡羅;
|2=zh-hans:马尔可夫; zh-hant:馬可夫;
|3=zh-hans:伊辛; zh-hant:易辛;
}}
'''易辛模型'''（{{lang-en|'''Ising model'''}}，{{IPAc-en|ˈ|aɪ|s|ɪ|ŋ}}，{{IPA-de|ˈiːzɪŋ|lang}}），是一個以物理學家[[恩斯特·伊辛]]為名的[[數學模型]]，用於描述物質的[[鐵磁性]]。該模型中包含了可以用來描述單個原子[[磁矩]]的參數<math>\sigma_i</math> ，其值只能為+1或-1，分別代表[[自旋]]向上或向下，這些磁矩通常會按照某種規則排列，形成[[晶格]]，並且在模型中會引入特定[[基本交互作用|交互作用]]的參數，使得相鄰的自旋互相影響。雖然該模型相對於物理現實是一個相當簡化的模型，但它卻和鐵磁性物質一樣會產生[[相變]]。事實上，一個二維的{{le|方晶格易辛模型|square-lattice Ising model}}是已知最簡單而會產生相變的物理系統。<ref>See {{harvtxt|Gallavotti|1999}}, Chapters VI-VII.</ref>

易辛模型最早是由物理學家[[威廉·楞次]]在1920年發明的，他把該模型當成是一個給他學生恩斯特·易辛的問題。易辛在他一篇1924年的論文中求得了一維易辛模型的[[解析解]]，並且證明它不會產生相變。<ref name="Ising">[http://www.hs-augsburg.de/~harsch/anglica/Chronology/20thC/Ising/isi_fm00.html Ernst Ising, ''Contribution to the Theory of Ferromagnetism'']</ref> 
二維方晶格易辛模型相對於一維的難出許多，因此其解析的描述在一段時間之後才在1943年由[[拉斯·昂薩格]]給出。一般來說，二維易辛模型的解析解可由{{le|傳遞矩陣法|transfer-matrix method}}求得，不過也有幾個和[[量子場論]]有關的解法。對於大於三維的易辛模型目前還沒有找到解析解，但其近似解可由諸多方法求得，例如[[平均場論]]。

==定義==
令Λ為所有晶格點的[[集合]]，其中每個晶格點都有一個所有和它相鄰的晶格點的集合(在數學上稱之為[[圖 (數學)|圖]])並使這些晶格點形成一個d維的晶格。對於每個晶格點 ''k''&nbsp;∈&nbsp;Λ 都有一個[[離散]][[變數]] σ<sub>''k''</sub> ，其中 σ<sub>''k''</sub> &nbsp;∈ {+1, −1}，代表一個晶格點的自旋。而所有變數的集合σ&nbsp;=&nbsp;(σ<sub>''k''</sub>)<sub>''k''∈Λ</sub>則稱作 '''自旋組態'''。

對於兩個相鄰的晶格點''i'',&nbsp;''j''&nbsp;∈ Λ ，我們可以引入一個交互作用參數''J''<sub>''ij''</sub>，此外，我們可以假設每個自旋''j''&nbsp;∈&nbsp;Λ都和外加的[[磁場]] h<sub>j</sub> 作用。則整個系統的[[哈密頓量]]可寫成:

:<math>H(\sigma) = - \sum_{<i~j>} J_{ij} \sigma_i \sigma_j -\mu \sum_{j} h_j\sigma_j</math>

其中 ''<ij>'' 代表晶格點 ''i'' 和晶格點 ''j'' 是相鄰的的晶格點。因此哈密頓量的第一項為對每一對相鄰晶格點的總和(每一對只算一次)，代表所有自旋之間交互作用的[[能量]]，而第二項則是磁場和自旋交互作用的能量。µ是晶格點[[磁矩]]的值，值得注意的是，[[電子]]的[[磁矩]]和他的自旋方向相反，所以哈密頓量的第二項應該要是正號比較合理，但在習慣上，還是會令第二項為負號。<ref>See {{harvtxt|Baierlein|1999}}, Chapter 16.</ref> 

該系統的的'''組態機率''' P(σ)為在[[熱平衡]]下某個特定自旋組態 σ 的[[機率]]，為[[波茲曼分布]]:

<math>P_\beta(\sigma) ={e^{-\beta H(\sigma)} \over Z_\beta},</math>

其中 β&nbsp;=&nbsp;(''[[波茲曼常數|k<sub>B</sub>]][[溫度|T]]'')<sup>−1</sup>，而:

<math> Z_\beta = \sum_\sigma e^{-\beta H(\sigma)}</math>

是該機率分布的[[歸一化常數]]，在[[統計力學]]中又稱做[[配分函數]]。對於有為自旋組態[[函數]]的[[物理量]] f(σ) ，其期望值可表示為:

<math> \langle f \rangle_\beta = \sum_\sigma f(\sigma) P_\beta(\sigma) \,</math>

===參數===
''H''(σ) 中兩項前的負號是約定俗成的。因為第一項為負號，因此參數 ''J''<sub>''ij''</sub> 的正負號決定了該系統的性質，對於每一對 ''i'' 和 ''j'' :
::<math>J_{ij} > 0 </math> ，則該系統為[[鐵磁性]]。
::<math>J_{ij} < 0 </math> ，則該系統為[[反鐵磁性]]。
::<math>J_{ij} = 0 </math> 表示自旋間無交互作用。

除此之外該系統為 '''非鐵磁性'''。

在鐵磁性的易辛模型中，相鄰[[自旋]]同方向時能量較低，因此自旋會傾向於同向排列，反之，在反鐵磁性的易辛模型中，相鄰[[自旋]]反向的能量較低，因此自旋會頃向於反向排列。

''H''(σ) 中的第二項為負號，表示自旋頃向於和外加[[磁場]]同向，因此 h<sub>j</sub> 的正負也決定自旋頃向排列的方向。對於所有的j ，如果:

::<math>h_j>0</math>, 則晶格點 ''j'' 頃向於朝向正向。
::<math>h_j<0</math>, 則晶格點 ''j'' 頃向於朝向負向。
::<math>h_j=0</math>, 表示沒有外加[[力場]]作用在[[自旋]]上 。

===簡化===
一個常見的簡化是假設沒有外加的[[磁場]]作用在易辛模型上，也就是說，對於所有的 j∈Λ，''h<sub>j</sub>''&nbsp;=&nbsp;0 。利用這項簡化，其[[哈密頓量]]可以寫成:

<math>H(\sigma) = - \sum_{<i~j>} J_{ij} \sigma_i \sigma_j. </math>

此時易辛模型在反轉所有自旋之下是[[對稱]]的:一個外加的力場會破壞這種對稱。

另一個常見的簡化是假設所有相鄰晶格點的交互作用都是相等的，因此可以設 ''J<sub>ij</sub>'' = ''J'' 對於所有相鄰的 ''i'', ''j'' ∈Λ，而其哈密頓量可以寫:
<math>H(\sigma) = -J\sum_{<i~j>}\sigma_i \sigma_j. </math>

==特性==
[[Image:Ising-tartan.png|thumb|right|一維平移不變鐵磁性在無外加場下的易辛模型的所有組態之[[機率密度]]投射在二維平面上。]]
被研究的最透徹易辛模型是在 d [[維度|維]]晶格上，平移不變、[[鐵磁性]]並且無外加場的模型。也就是 ：Λ&nbsp;=&nbsp;'''Z'''<sup>''d''</sup>, ''J''<sub>''ij''</sub>&nbsp;=&nbsp;1, ''h''&nbsp;=&nbsp;0。

易辛在他1924年的博士論文中，解決了在 d=1 時的情況 ，這個一維的模型可以想像成一排的[[自旋]] ，而每個自旋都只和它左右兩邊的自旋交互作用。這個一模型不會產生[[相變]]<ref>{{Cite web |url=http://users-phys.au.dk/fogedby/statphysII/no-PT-in-1D.pdf |title=存档副本 |access-date=2014-09-07 |archive-url=https://web.archive.org/web/20140407154639/http://users-phys.au.dk/fogedby/statphysII/no-PT-in-1D.pdf |archive-date=2014-04-07 |dead-url=yes }}</ref> 
，換句話說，對於所有正值的β，任意兩自旋的[[相關係數]] <σ<sub>''i''</sub>σ<sub>''j''</sub>> 都對 |''i''−''j''| 呈[[指數|指數衰減]]:

<math> \langle \sigma_i \sigma_j \rangle_\beta \leq C \exp(-c(\beta) |i-j|)</math>

其中c(β)是一個只和β有關的函數。

由此可知。這個系統是[[無序相|無序]]的。根據一維模型的結論，易辛錯誤的認為任何維度的易辛模型都不會有[[相變]]。但事實上，在二維或二維以上的模型中，該系統可以歷經從[[無序相]]轉變成[[有序相]]的[[相變]]。基本上在β值小(高溫)時，系統處在無序相，而β值大(低溫)時，系統處在有序相中。換句話說，當系統在有序相時：

<math> \langle \sigma_i \sigma_j \rangle_\beta \geq c^{\prime}(\beta) > 0.\,</math>

其中c'(β)也是一個只和β有關的函數。

這個性質是首先被[[魯道夫·佩爾斯]](Rudolf Peierls)在1936年證明的<ref>http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=2027260</ref>
，他的證法後來被稱為'''佩爾斯論述'''(Peierls argument)。

在零磁下的二維方晶格易辛模型的[[解析解]]後來在1943年被[[拉斯·昂薩格|昂薩格]]解出，他證明了該模型的{{le|相關函數|correlation function}}和[[自由能]]可由一個無交互作用的[[格點費米子場]](noninteracting lattice fermion)來界定。昂薩格在1949年發表了一個決定了自發磁化現象的公式，但卻沒有給出推導過程。後來是[[楊振寧]]在1951年發表了第一個正式的推導過程，其中裡面用到了包括{{le|塞格極限定理|Szegő limit theorems}}和{{le|弗雷德霍姆行列式|Fredholm determinant}}等數學工具<ref name="Montroll 1963 pages=308-309">{{harvnb|Montroll|Potts|Ward|1963|pages=308–309}}</ref>。
==數值模擬==
===蒙地卡羅方法===
易辛模型一般來說很難直接進行數值計算，因為他的自旋組態非常之多。考慮一個擁有  ''L'' 個晶格點的模型，每個晶格點 σ<sub>''j''</sub> 有 ±1 兩種可能，因此所有的自旋組態共有 ''2''<sup>''L''</sup> 種可能<ref name = "Newman">Newman MEJ, Barkema GT, "Monte Carlo Methods in Statistical Physics, Clarendon Press, 1999</ref>，這個數字會隨著 ''L'' 的增加而進行[[指數]]增長。這也是為什麼一般在做易辛模型的數值模擬時，都會採用[[蒙地卡羅方法]]（Monte Carlo method）<ref name="Newman" />

用蒙地卡羅方法來模擬易辛模型所用的[[哈密頓量]]為:
:<math>H(\sigma) = - J\sum_{<ij>}\sigma_i \sigma_j -h\sum_{j}\sigma_j.</math>
此外，可以假設外加磁場 (''h'') 為零以簡化模型，因為大部分的問題都只需要用到零磁下的模型。因此，近一步簡化的[[哈密頓量]]為

:<math>H(\sigma) = - J\sum_{<ij>}\sigma_i \sigma_j.</math>
可以從這個模型計算出的[[物理量]]包括該系統在定溫下的[[比熱]]和[[磁化強度]]等等。<ref name="Newman" />

===Metropolis 演算法===
[[Metropolis–Hastings演算法]]是在數值模擬易辛模型時最常用的一種蒙地卡羅方法。<ref name="Newman" />這種方法首先要選定一個'''選擇機率''' ''g''(μ, ν)，代表系統在狀態μ下，在所有可能狀態中選到狀態ν的機率。另外還需定義出一個'''接受機率''' ''A''(μ, ν) ，也就是說當系統在狀態μ下，接受系統跳到態ν的機率。如此的設計是為了讓系統達到[[細緻平衡]]。如果狀態ν被接受了，則整個系統就會跳到狀態ν，並且選擇和決定下一個要跳到的狀態。如果狀態ν被沒被接受，系統則留在狀態μ，一樣重新選擇下一個要跳到的狀態。這樣的步驟一直重複直到某些條件達成為止，譬如說整個易辛模型完全被[[鐵磁性|磁化]]，也就是所有的[[自旋]]都只到同一個方向。<ref name="Newman" />另外在實行這種演算法，有一點需要注意的是需要選到適當 ''g''(μ, ν) 以保證整個過程的[[遍歷性]](ergodicity)。

在[[熱平衡]]時，整個系統的[[能量]]只會有小幅度的擾動<ref name="Newman" />，這點促成了在演算時採用'''單一自旋反轉法'''進行計算，也就是說每次系統轉換其狀態時，只改變其中一個[[自旋]]的方向。對自旋數很多的一易辛模型來說，系統在不同的狀態之間跳躍時，其[[能量]]改變的幅度都很小。事實上，對於每個晶格點都和 ''c'' 個晶格點相鄰的模型來說，每次能量改變的最大幅度為 2''cJ''。此外，採用這種單一自旋反轉法可以保證演算過程的[[遍歷性]]，因為任意一個狀態都可以藉由逐次的反轉相異的自旋，而變成任意其他狀態<ref name="Newman" />。

在下面的演算過程我們會採用[[周期性边界条件]]使的每個晶格點的相鄰數 ''c'' 都相等。<ref name="Newman" />

==== 演算過程 ====
整個用於數值模擬易辛模型的演算過程，可由下列的方法建立。

因為共有 ''L'' 個晶格點，而且單一自旋反轉法是唯一系統可以從一個狀態跳到另一個狀態的方法，所對於任何一個狀態 μ，有 ''L'' 個新狀態 ν 它可以跳去。可以假設 μ 跳到這 ''L'' 個新狀態的機率是相等的，因此 ''g''(μ, ν) = 1/''L''。 為了滿足{{le|細緻平衡|detailed balance}}，所以下面的等式要成立:

:<math> \frac{P(\mu,\nu)}{P(\nu,\mu)} = \frac{g(\mu,\nu)A(\mu,\nu)}{g(\nu,\mu)A(\nu,\mu)} = \frac{A(\mu,\nu)}{A(\nu,\mu)} = \frac{P_\beta(\nu)}{P_\beta(\mu)} = \frac{\frac{1}{Z}e^{-\beta(H_\nu)}}{\frac{1}{Z}e^{-\beta(H_\mu)}} = e^{-\beta(H_\nu-H_\mu)}.</math>

因此我們希望接受機率滿足:

:<math> \frac{A(\mu,\nu)}{A(\nu,\mu)}=e^{-\beta(H_\nu-H_\mu)}.</math>

如果 ''H''<sub>ν</sub> > ''H''<sub>μ</sub>  則 ''A''(ν, μ) > ''A''(μ, ν)，因此 Metropolis 另 ''A''(μ, ν) 和 ''A''(ν, μ) 中較大的為 1。按照這樣的定可以得到 ''A''(μ, ν) 的值為: <ref name="Newman" />

<math>A(\mu,\nu)=\begin{cases}
 e^{-\beta(H_\nu-H_\mu)}, & \text{if }H_\nu-H_\mu>0 \\
 1, & \text{otherwise}.
\end{cases}</math>

最後，整個演算法最基本的形式為: 
#用 ''g''(μ, ν) 選出一個自旋，並且計算所有和其自旋相關的能量貢獻。(也就是所有和它相鄰自旋交互作用的能量貢獻。)
#反轉其自旋，再計算一次所有和它相關的能量貢獻。
#如果其能量貢獻下降，保持自旋反轉。
#如果貢獻能量上升，則令自旋有 <math> e^{-\beta(H_\nu-H_\mu)}</math> 的機率保持反轉。
#重複步驟一。

整個系統能量的改變量 ''H''<sub>ν</sub>−''H''<sub>μ</sub> 只和反轉的自旋和和它相鄰的自旋有關，所以只要相鄰的自旋數不要太多，計算的速度算是相當快的。而整個系統會逐漸的趨近於某個平衡。

===將易辛模型視為馬可夫鏈===
將易辛模型比擬為[[馬可夫鏈]]是一件很容易的事情，因為下一刻狀態 ν 的轉移機率  ''P''<sub>β</sub>(ν) 只和目前狀態 μ有關 。事實上，Metropolis 演算法就是[[馬可夫蒙地卡羅]](Markov Chain Monte Carlo)的一個版本。而且因為演算時採用的是單一自旋反轉法，每個狀態可視為和另外 ''L'' 個狀態相連，而每一個狀態轉移則是將其中一個[[自旋]]反轉。<ref>{{cite journal|last=Teif V.B. |title=General transfer matrix formalism to calculate DNA-protein-drug binding in gene regulation|journal=Nucleic Acids Res. |year=2007|volume=35|pages=e80|doi=10.1093/nar/gkm268}}</ref> 此外，因為不同狀態之間 ''H''<sub>σ</sub> 變化只和相鄰自旋的交互作用 '''J''' 有關，所以易辛模型也可視為{{le|選民模型|Voter model}}的一種。

==一維易辛模型==
在一維易辛模型系統中，假設每個带有自旋的原子分布在一維的圓圈中，且原子仅和鄰居發生交互作用，交互作用均為<math>J</math>，能量可表示為

::<math>H(\sigma) = -J \sum _{i=1}^N  \sigma_i\sigma_{i+1}+H\sum _{i=1}^N \sigma_i</math>

其中''h''为外加磁场的强度，''J''为相邻原子的耦合强度。该系统的[[自由能]]为：

:<math>Z(\beta, h)=-\lim_{N\to \infty} \frac{1}{\beta L} \ln (Z(\beta))=-\frac{1}{\beta} \ln\left(e^{\beta J} \cosh \beta h+\sqrt{e^{2\beta J}(\sinh\beta h)^2+e^{-2\beta J}}\right)
</math>
相邻自旋的{{en-link|相关函数_(统计力学)|Correlation function (statistical mechanics)|相关函数}}为：
:<math>\langle \sigma_i \sigma_j\rangle-\langle \sigma_i \rangle\langle\sigma_j\rangle=C(\beta)e^{-c(\beta)|i-j|}</math>
其中，C(β)和''c''(β)是关于温度的函数，当温度''T'' > 0时取正值。当''T'' → 0时，''c''(β)趋于零。

===证明===

易辛自己的论文中已给出了一维易辛模型[[配分函数]]精确解。<ref name="Ising"/>无外加磁场，即''h'' = 0时，计算是很简单的。

此时能量简化为：
:<math>H(\sigma)=-J(\sigma_1\sigma_2+\cdots+\sigma_{N-1}\sigma_N).</math>
对于相邻两原子自旋朝向的四种情况，之间的能量只有两种状态：同向和反向。应用变量替换，
:<math>\sigma'_j=\sigma_j\sigma_{j-1} \qquad j\ge 2.</math>
代入配分函数定义式，得：

:<math> Z(\beta) =\sum_{\sigma_1,\ldots, \sigma_N} e^{\beta J\sigma_1\sigma_2}\; e^{\beta J\sigma_2\sigma_3}\; \cdots e^{\beta J\sigma_{N-1}\sigma_N}= 2\prod_{j=2}^N \sum_{\sigma'_j}  e^{\beta J\sigma'_j} =2\left[ e^{\beta J}+e^{-\beta J}\right]^{N-1}. </math>

因此，系统的自由能为

:<math>Z(\beta,0)=-\frac{1}{\beta } \ln\left[e^{\beta J}+ e^{-\beta J}\right]. </math>

同样的变量替换法，得间隔N-1个自旋间的相关函数：

:<math>\langle \sigma_{j}\sigma_{j+N}\rangle=\left[\frac{e^{\beta J}- e^{-\beta J}}{e^{\beta J}+ e^{-\beta J}}\right]^N</math>

''T'' ≠ 0时相关性呈指数衰减。而仅仅在绝对零度时''T'' = 0下，换句话说，β → ∞ 时，系统保持长程相关性。

如果外磁场''h'' ≠ 0，配分函数的计算需要引入[[传递矩阵法]]。[[周期性边界条件]]的近似下，配分函数为
:<math>Z(\beta)=\sum_{\sigma_1,\ldots, \sigma_N} e^{\beta h \sigma_1}e^{\beta J\sigma_1\sigma_2}\; e^{\beta h \sigma_2}e^{\beta J\sigma_2\sigma_3}\; \cdots e^{\beta h \sigma_N}e^{\beta J\sigma_N\sigma_1}=  \sum_{\sigma_1,\ldots, \sigma_L} V_{\sigma_1,\sigma_2}V_{\sigma_2,\sigma_3}\cdots V_{\sigma_N,\sigma_1}. </math>
系数<math>V_{\sigma, \sigma'}</math>可视为2x2矩阵的元素。相邻两自旋各有四种可能状态，其玻尔兹曼因子分别为：
:<math> V_{\sigma, \sigma'} = e^{\frac{\beta h}{2} \sigma} e^{\beta J\sigma\sigma'} e^{\frac{\beta h}{2} \sigma'}</math>
或写作传递矩阵的四个元素：
:<math>V=\begin{bmatrix} e^{\beta(h+J)}&e^{-\beta J}\\ e^{-\beta J}&e^{-\beta(h-J)} \end{bmatrix}.</math>
配分函数的展开式恰好和传递矩阵自乘后对角元素之和（[[矩阵的迹]]）一致。而矩阵的迹可通过求解传递矩阵的特征值得出，因而：
:<math> Z(\beta)={\rm Tr} V^N= \lambda_1^N + \lambda_2^N= \lambda_1^N\left[1+ \left(\frac{\lambda_2}{\lambda_1}\right)^N\right]</math>
其中λ<sub>1</sub>是''V''最大的特征值，λ<sub>2</sub>是该矩阵另一个特征值。因而
:<math>\lambda_1=e^{\beta J} \cosh \beta h+ \sqrt{e^{2\beta J} (\sinh\beta h)^2  +e^{-2\beta J}} </math>
由于|λ<sub>2</sub>| < λ<sub>1</sub>，当''N''很大时对配分函数的贡献可忽略。

===结论===
系统最低能量为−''N''，此时所有自旋朝向相同。而其它构型相较基态的能量增量等于一维序列中自旋转向的次数''k''。例如，能量次低的状态相较基态，能量差为2''k''。由于能量与转向数之间是线性关系，则两个相邻自旋方向相反的概率符合玻尔兹曼分布：

:<math>{p \over 1-p} = e^{-2\beta J}.</math>

藉由[[統計力學]]的[[配分函數]]可以計算再給定溫度下<math>T</math> (<math>
\beta =1/ (k T)
</math>)的每個原子的磁矩期望值為
:: <math>
M (H,T)=\frac{\sinh (\beta H )}{\sqrt{\sinh
   ^2(\beta H )+e^{-4\beta J  }}}
</math>.

所以一維易辛模型並沒有[[居里温度]]、不會發生[[相變]]，即沒有{{le|自发磁化|spontaneous magnetization}}的現象。
:: <math>
M (0,T)=0
</math>.

==二維易辛模型==

设二维方晶格易辛模型横纵两方向的交互作用能分别为<math>J_1</math> and <math>J_2</math>。[[拉斯·昂萨格]]求得无外磁矩，即<math>h=0</math>时自由能的解析解：
:<math> -\beta f = \ln 2 + \frac{1}{8\pi^2}\int_0^{2\pi}d\theta_1\int_0^{2\pi}d\theta_2 \ln[\cosh(2\beta J_1)\cosh(2\beta J_2) -\sinh(2\beta J_1)\cos(\theta_1)-\sinh(2\beta J_2)\cos(\theta_2)]. </math>

从自由能的偏导数可得到各种热力学函数。

特别地，二维易辛模型有一个相变点，临界温度<math>T_c</math>满足以下方程：

: <math> \sinh\left(\frac{2J_1}{kT_c}\right)\sinh\left(\frac{2J_2}{kT_c}\right) = 1 </math>.

若<math>J_1 = J_2</math>，则<math> T_c = \frac{2 J}{k\ln(1+\sqrt{2})} </math>，或 <math> \beta = \frac{1}{kT_c} \approx \frac{0.44}{J} </math>.

== 参见 ==
* [[圖模式]]
* [[马尔可夫链]]
* [[马尔可夫逻辑网络]]

==延伸閱讀 ==
* Kerson Huang, Introduction to Statistical Physics.
* I. A. Stepanov. Exact Solutions of the One-Dimensional, Two-Dimensional, and Three-Dimensional Ising Models. - Nano Science and Nano Technology: An Indian Journal. 2012. Vol. 6. No 3. 118 - 122.(這篇文章可在該期刊的網站免費閱讀)
* Barry A. Cipra, "The Ising model is NP-complete", [[SIAM News]], Vol. 33, No. 6; [https://web.archive.org/web/20070926223950/http://www.siam.org/pdf/news/654.pdf online edition (.pdf)](一篇文章闡述為何任意的易辛模型無法有一般性的精確解，因為非平面的易辛模型是[[NP完全]]的。)

==相關連結==
* [http://scienceworld.wolfram.com/physics/IsingModel.html Science World article on the Ising Model]
* [https://web.archive.org/web/20120716232344/http://www.phy.syr.edu/courses/ijmp_c/Ising.html An Ising Applet by Syracuse University]
* [http://physics.ucsc.edu/~peter/ising/ising.html A dynamical 2D Ising Applet by UCSC]
* [https://web.archive.org/web/20140907212614/http://personal-pages.ps.ic.ac.uk/~achremos/Applet2-page.htm A nice dynamical 2D Ising Applet]
* [http://www.physics.uci.edu/~etolleru/IsingApplet/IsingApplet.html A larger/more complicated 2D Ising Applet]
* [https://web.archive.org/web/20101111022621/http://dtjohnson.net/projects/ising A nice HTML5 Ising Model simulation]
* [http://demonstrations.wolfram.com/IsingModel/ Ising Model simulation]by Enrique Zeleny，由 [[Wolfram 演示項目]] 提供。
* [http://ibiblio.org/e-notes/Perc/contents.htm Phase transitions on lattices]
* [http://www.sandia.gov/media/NewsRel/NR2000/ising.htm Three-dimensional proof for Ising Model impossible, Sandia researcher claims]
* [http://dx.doi.org/10.1016/j.cpc.2010.05.005 Multi-GPU accelerated multi-spin Monte Carlo simulations of the 2D Ising model]
* [http://www.fisica.unige.it/~cavaliere/soft/MacIsing.zip Interactive dynamical simulation for MacOs of the 2D ising model on a square lattice]
* [http://www.unige.ch/math/folks/velenik/smbook/ Equilibrium Statistical Mechanics of Classical Lattice Systems: a Concrete Introduction, Chapter 1 (The Ising model)]
* [http://isingspinwebgl.com Interactive Monte Carlo simulation of the Ising, XY and Heisenberg models with 3D graphics(requires WebGL compatible browser)]

==註解==
{{Reflist}}

{{统计力学}}
{{Stochastic processes}}

[[分類:統計力學]]