ising model

伊辛模型
Yixin moxing
伊辛模型
Ising model
描述物质相变的一种模型。

物质经过相变，要出现新的结构和物性。

发生相变的系统一般是在分子之间有较强相互作用的系统，又称合作系统。

在铁和镍这类金属中，当温度低于居里温度（见铁磁性）时，原子的自旋自发地倾向某个方向，而产生宏观磁矩。

温度高于居里温度时,自旋的取向非常紊乱,因而不产生净磁矩。

当温度从大于或小于两边趋于居里温度时，金属的比热容趋于无限大。

这是物质在铁磁性状态和非铁磁性状态之间的相变，它并不包含在P.厄任费斯脱所分类的相变中。

伊辛模型就是模拟铁磁性物质的结构，解释这类相变现象的一种粗略的模型。

它的优点在于，用统计物理方法，对二维情形求得了数学上严格的解。

这就使得铁磁性物质相变的大致特征，获得了理论上的描述。

这个模型所研究的系统是由个阵点排列成维周期性点阵，这里=1,2,3点阵的几何结构可以是立方的或六角形的,每个阵点上都赋予一个取值+1或－1的自旋变数,如果＝+1,即第个阵点的自旋向上；如＝-1，即第个阵点的自旋向下并且认为只是最近邻的自旋之间有相互作用。

点阵的位形用一组自旋变数{}（＝1，2，…,）来确定。

图1[二维伊辛点阵模型]是一个二维伊辛模型的示意图，图中表示自旋向上,[kg2]表示自旋向下。

处理方法20世纪30年代初，不少科学家如W.L.布格 E.J.威廉斯H.A.贝特、R.E.佩尔斯等人就已从有序－无序转变问题及点阵气体等模型出发，采用平均场近似法处理伊辛模型。

布格－威廉斯平均场近似法认为，某一阵点上的自旋取某一方向的几率同近邻阵点上的自旋取向无关，只同自旋在该方向的数目成正比。

每个阵点上有一平均磁场，自旋在阵点上的取向只同该磁场有关。

用这种方法可求得下列公式
[1148-01]式中是每个自旋的磁矩,是每一阵点的最近邻数，是外磁场强度，是热
力学温度，是自旋同向的最近邻对之间的相互作用能（铁磁性物质0，非铁磁性物质0），是玻耳兹曼常数, [kg2]是每个自旋上的磁化强度，可表示为[1148-02],其中是总阵点数,[1110-17]、[1110-18]分别代表自旋向上和向下的总阵点数,显然[1110-17]+[1110-18]=。

由此研究铁磁性物质的性质，得到如下结论：存在一临界温度[1148-03]，当而＝0时，物质不磁
化，没有相变；当时，尽管仍有=0，但磁化强度可不为零（可取正值或负值），铁磁性物质存在相变。

这个结论对一维、二维、三维点阵都应成立，但严格的证明指出，二维、三维伊辛模型在临界温度以上仍有相变。

这反映了平均场近似法的简单、粗糙。

当然，用它处理临界温度以下铁磁性物质的相变，仍是一种有意义的方法。

20世纪40年代L.昂萨格对伊辛模型采用解析法得到了严格解，作出了突出的成就。

这种方法的基本点，是设每个阵点的自旋变数可取+1和－1两个值,考虑阵点上自旋的某个位形，计算每个自旋同最近邻自旋的相互作用能量以及同外磁场相互作用能量，再对全部可能的位形求和，用矩阵的方法求出配分函数，从而得到各个热力学函数。

一维情况考虑具有个自旋的直线链（图2[ 个阵]
个阵" class=image>[点的直线链]所示），每个自旋仅同它的两个最近邻自旋及外磁场相互作用。

相互作用的总能量即由{，，，…,-1}所确定的位形能量是
[1148-04],式中[1148-05]是对最近邻自旋对求和，[1148-11]
表示对所有自旋坐标求和。

配分函数可写作
[1148-06]定义一个2×2矩阵，它的矩阵元是
[1148-07]，因为=±1，=±1，故矩阵可表示为
[1148-08]如果采用周期性边界条件＝，或设想将直线链弯成闭合的
圆链，并将初端与尾端相接(图3[ 圆链])，配分函数I可表示成矩阵形式[1148-12]，式中tr表示矩阵的迹,它是[1148-13]的对角元之和。

由此计算可以求出I，对于很大的链,则有[1148-09][1148-30]。

磁化强度是
[1148-10]，(,
伊辛模型
)同的关系如图4[一维伊辛模型磁化强度随磁场的变化曲线]
所示。

当→0时,对于0,有(,0)＝0。

可见，一维伊辛模型没有自发磁化即不显示铁磁性，因而不发生相变。

在一维伊辛模型中，不论铁磁性或反铁磁性，都不会实现有序的状态。

如对于0的铁磁性物质，在绝对零度时，所有自旋取向都是相同的，此时，处于能量最低的状态。

然而，如果热力学温度不等于零，是有限的，则平均位形由两种相反的、相互竞争的趋向所决定。

一个是各个自旋的取向完全一致，使能量最低；另一个是各个自旋的取向为随机的，使熵最大。

由于一维伊辛模型中每个自旋没有足够多的最近邻自旋，因而不可能出现所有自旋取向完全相同的情况,而是如图5[一维伊辛模型0时自旋的排列示意图]
0时自旋的排列示意图" class=image>所示。

二维情况（昂萨格解）二维伊辛点阵的阵点数为×=,为便于计算，画成图6[二维伊辛点阵]
所示的情形。

处理二维空间问题的方法与一维的类似，只需将一维的每个阵点当作一列，并逐列相加求和即可。

以表示第行的所有自旋坐标的集合
[1149-08]上标(=1，2，…，)代表行,下标(1，2，…，)代表列。

边界条件为[1149-09]。

即要求每一行的第+1列的位形与第1列的相同,每一[1149-12]，有2个值。

整个点阵的位形由{，，…，}确定。

考虑最近邻自旋对以及自旋同外磁场的相互作用，则配分函数可写成
[1149-01]为将上式表示成矩阵的形式,引入三个矩阵、、，它们的矩阵元分别定义为
[1149-02]第一式反映不同行最近邻自旋对的相互作用能量，它有2×2个；第二式反映同一行最近邻自旋对的相互作用能量，它有2个；第三式反
映同一行各个自旋与外磁场的相互作用能量，它也有2个。

为了计算方便，在补上一些“0”元素后，可把、扩大成2×2矩阵的对角矩阵。

可以证明
＝tr()。

当→∞时,求×矩阵的本征值问题就变成求解2×2矩阵的本征值问题。

H.A.克末和G.H.万尼尔等人曾用数字解计算过有限的几项，他们计算到=5，发现在为有限的情形下，没有相变。

昂萨格在求解时，设外磁场强度=0,因而=1。

计算结果表明：高温时,(临界温度)，矩阵＝只有一个最大本征值；低温时，，矩阵=有两个本征值，当→∞，→∞时，配分函数为
[1149-03]并得出平均每个自旋的自由能为
[1149-04]若令[1149-002]，则上式右边第二项被积函数()满足ch()＝ch2ch2-cos sh2sh2。

用数值计算,通过上述二式可算出→时的各个热力学量，得到以下具体结果，
-＝-＝(0.9296…)，
＝-＝ln(1.358)。

其中为每个自旋的熵。

式中的临界温度满足方程
[1149-06]或－=2.269185。

在附近，每个自旋的比热容可表示为
[1149-07]
可见在＝时，自由能、熵以及内能是连续的，这意味着在＝时，发生的相变不包含潜热。

但是当[1149-13]时，作为上述热力学函数的导数，比热容是对数发散的,无论从高温端还是低温端趋于(即→+0或→-0)，比热容的值是相同的。

图7[二维伊辛模型比热容]
给出了比热容随温度变化的曲线，并且同时画出了布格—威廉斯平均场近似法所得结果的曲线（图中虚线），以作比较。

为弄清=处相变的细节，还需进一步考虑自发磁化（即计算自由能对磁场强度的导数，再让=0）。

杨振宁于1952年采用微扰法得到了很好的结果。

他证明自发磁化强度(0,)可表为
[1150-01]式中[1150-02]对应于的值为[1150-03]，自发磁化强度随温度变化的曲线如图8[ 二维伊辛模型自发磁化强度曲线]
所示。

至于存在外磁场的情形,以及三维空间的解析解,虽经许多理论物理学家多年的工作，但至今还未取得令人满意的结果。

参考书目
Kerson Huang,Statistical echanics,John Wiley & Sons,New York, London, 1963.
欧阳容百。