二维ising模型概率c语言_二维零场 Ising model 严格解

好久没更新_(:3」∠)_

这次专栏文章展示二维零场 Ising model 的严格解，主要遵循Kaufman在1949年相关文章中的逻辑，他对Onsager的解法做出了简化。我们首先将计算配分函数的问题转化为对角化转移矩阵的问题，把转移矩阵表示为局域泡利矩阵的直积，进一步用Majorana fermion表示。在把转移矩阵按宇称在两个子空间上分解后，通过一般的观察将

维大矩阵对角化的问题约化为

维小矩阵的对角化问题，从而求出两个子空间上转移矩阵的本征值，也就求出系统的配分函数。

---

一、二维Ising model的转移矩阵

我们考察的对象是经典的二维Ising model，spin放置在m行n列的方形格子上，只有向上与向下两个取向；每个格点有四个近邻，行间铁磁相互作用强度为

，而行内为

；两个方向均采取周期性边界条件。记第

行spin的构型为

，它给出这一行所有spin的朝向

，那么逆温

下系统的配分函数可以形式化地表示为对所有不同构型玻尔兹曼因子的求和：

其中

表示行间相互作用的能量，而

表示行内相互作用的能量。

转移矩阵（transfer matrix）的想法是，将一行的所有（

）构型视为一组基，那么单行的玻尔兹曼因子，最多只与两个构型（本行与下一行）有关，它“吃”进去两个构型（两个态）就“吐”出来一个数，因此就可以看做构型空间上的矩阵。具体来说，行间：

而行内的玻尔兹曼因子只与一个构型有关，所以是单位阵：

其中的

与

是无量纲化的相互作用强度：

在这种表示下，给定整个格子的构型

， 玻尔兹曼权重是

其中每一个

，都联系了邻近两行的自旋