Ising模型简述.doc

1、WORD格式Ising模型简述Lenz 曾向他的学生Ising提出一个研究铁磁性的简单模型，而Ising于 1925年发表了他对此模型求解的结果，所以这个模型被称为Ising模型。当时Ising只做出了该模型一维下的严格解，在一维情况下并没有自发磁化的发生。另外他还由此错误地推断出在更高维的情况下，这个模型也不存在自发磁化。这个推断在后来被证明是错误的。1936 年 Peierls论证了二维或三维的Ising模型存在着自发磁化，虽然当时他并没有能够给出模型的严格解。1944 年，当 Onsager 给出了二维 Ising模型的严格解之后，Ising模型开始引起人们广泛的关注。这次求解是相变理论

2、发展上的一个重要进展，它第一次清楚地证明了从没有奇异性的哈密顿量体系出发，在热力学极限下能导致热力学函数在临界点附近的奇异行为，而Onsager 本人也因此获得了诺贝尔奖。在此之后很多人又相继发表Ising模型的各种不同解法， Baxter甚至有篇论文叫Ising模型的第 399 种解法。但至今没有被学术界公认的三维Ising模型精确解。甚至有人发表论文证明无法解出三维Ising模型的精确解，因为三维Ising模型存在拓扑学的结构问题。人们通常用分子场理论及其改进理论、高温级数展开、低温级数展开、重整化群理论、蒙特- 卡罗模拟等近似计算三维Ising模型的居里温度和临界指数，而其中Wilson

3、于 1971 年发展的重整化群理论能以较高精度计算三维Ising模型的近似结果18-20。我国科学家张志东提出三维“Ising模型”精确解猜想。张志东的出发点就是拓扑学中的一个常识：低维空间的扭曲和纽结可以被高一维空间的旋转打开。通过引入第四卷曲起来的维与本征矢量上的权重这两个猜想作为处理三维Ising模型拓扑学问题的边界条件，并应用这些猜想用自旋分析法评估了三维简单正交晶格Ising模型的配分函数。当系统的对称性越高，居里温度也越高。他猜测三维系专业资料整理WORD格式统具有最高对称性的简单立方Ising 模型具有最高的居里温度黄金解，在二维系统具有最高对称性的正方Ising模型具有最高的居

4、里温度白银解。获得的结果具有一定的对称性和美学价值，并可部分返回到二维和一维的结果。当然，推定的精确解正确性取决于猜想的正确性，而且其与学术界通常接受的评价标准尚不完全吻合，有待于对相关的物理本质作进一步探讨。因此，这一工作目前还只是停留在猜想阶段。今天的 Ising模型根本不再是Ising博士论文中的模样。每年差不多有6000篇左右的论文研究这一模型。除了铁磁性之外，该模型还应用于很多方面，如合金中的有序 - 无序转变、液氦到超流态的转变、液体的冻结和蒸发、晶格气体、玻璃物质的性质，甚至于神经网络蛋白质折叠、生物膜场论甚至社会现象等广泛的领域。通过上述介绍，我们知道三维Ising模型尚未得到

5、严格解，而一维和二维情况下的解法确是多种多样的。在这里，我们将给出Ising模型的严格解，采用的是 1941 年 Kramers 和 Wannier 提出的转移矩阵方法(TransferMatrixMethod)。然后简要地说明二维 Ising模型严格解的主要结果，并且同平均场理论所得的结果进行对比。图 1.2 一维 Ising模型示意图。专业资料整理WORD格式对于如图1.2 所示的 Ising模型，自旋只能取向上或向下两个分量，它可以看作是 Heisenberg 模型的一种简化。当只考虑最近邻的交换相互作用，并认为这种相互作用在不同磁矩间是相同的，用常数 J 表示。和Heisenberg

6、模型相同，当J0 时，代表铁磁的交换相互作用，它使得近邻自旋有着同方向排列的趋向；当J0 代表铁磁相互作用，B 为 Bohr磁矩， h 是外磁场。对于一维情况，每个自旋只有两个近邻。现在采用周期性边界条件，即sN+1=s1，N 为晶格中的自旋数目。现将一维晶格弯成一个环，当N时，边界效应将不会影响到体系的热力学性质。根据如上的条件，可将哈密顿量(1-14) 写为：BhHJs i si1sisi1 ，(1-15)i2i其相应的配分函数为：N1BhQT,hexpJs i si1si si1。(1-16)s11sN1i1kBT2在这里我们引入矩阵P，其矩阵元定义为：专业资料整理WORD格式1Bhsi

7、 Psi1 expJs i si1si si1，(1-17)kBT2因为 si 与 si+1 都能取1 两个值，所以P 是 22 的矩阵：si1Psi11si1Psi11Psi1Psi11si1Psi11。(1-18)JBhkBTJkBTeeJk BTJBhkBTee于是配分函数 (1-16) 可以重新写成：QT,hs1Ps2s2Ps3sN1PsNsNPs2s11sN1。(1-19)s1 PNs1TrPNs11将 P 矩阵对角化得，0P，(1-20)0+和- 即为矩阵P 的本征值，由下面的久期方程决定,JBhkBTJk BTee0 ，(1-21)JkBTJBhkBTee其解为：Jk Thcos

8、h2BhJk T 2J，(1-22)e BcoshB2e Bsinhk Tk Tk TBBB专业资料整理WORD格式要注意的一点就是+ -。专业资料整理WORD格式现在将等式 (1-20) 代入 (1-19) ，配分函数可以表达为：NQT,hNNN，(1-23)1所以，当 N时，我们得到：lim 1lnQT,hlnNN，(1-24)JlncoshBhcosh 2Bh2eJkBTsinh2JkBTkBTkBTkBT即配分函数有P 矩阵较大的本征值决定。体系的自由能和总极化强度分别为：FT,hQT,hBhBJk Tcoshk TBNNkBT，(1-25)cosh 2Bh2eJkBTsinh2JkB

9、TkBTBhsinhM1FNk TB，(1-26)NBBhTe4JkBTsinh2BhkBT其它的热力学函数也可同样由自由能求出。如图 1.3 所示，在计算中我们选取交换相互作用常数J=1kBK，对于一切T0都有 M(T,0)=0 ，也就是说Ising模型在一维的情况下不存在自发磁化，不会发生顺磁 - 铁磁转变。从物理上看，任何温度下自旋的平均取向由两个对抗的因素相互竞争决定，即能量趋向最小而熵趋向最大，使得自由能达到最小值。在一维情况下，由于近邻数低，使得自旋排列在同方向的倾向不足以对抗使熵极大的倾向，专业资料整理WORD格式结果在任何有限温度下都不能形成自发磁化。图 1.3 一维 Isin

10、g模型在不同温度下，磁化强度M随外场 h 的变化曲线。图 1.4 一维 Ising模型在有限温度下长程序被破坏的示意图。如同上文所说，当自旋排列在同方向的倾向不足以对抗使熵极大的倾向时，自旋往往会在一个较小的尺度内保持着同方向的排列，形成所谓短程序(ShortRangeOrder) ，而在较大的尺度内失去这种有序的状态，也就是破坏了所谓长程序 (LongRangeOrder) 。当我们使用蒙特卡洛方法 (MonteCarloMethod) 来计算一维 Ising模型，常常得到如图1.4所示的结果，在某个小范围内，如从格点1 到格点 4( 或者从格点5 到格点 10) 体系可以看作存在自发磁化，

11、而在整体上看( 从格点 1 到格点 N)向上的自旋和向下的自旋数目在统计上看是相等的。专业资料整理WORD格式对于二维Ising模型，我们考虑正方晶格，每个自旋有4 个最近邻。在零磁场下，系统的自由能可以表达为：FT,01lnQT,0NN。(1-27)1 1ln2sinh2J1d022其中 cosh()=cosh2cosh2-coshsinh2sinh2，而 =J， =tan -1 e-2J 。系统内能则可以写为：uT,01FJcoth2J12mK1m，(1-28)N这里的 K1(m) 是第一类椭圆积分，2dK1m，(1-29)0221msin其中 m=sinh2J/cosh 22J，m=2tanh22J-1 。而临界点由下式确定：kBTC2.269J 。(1-30)所得热容量为：CB T,0T0.4945ln1const ，(1-31)NkBTC这样的热容量在临界点处具有对数发散的奇异性。计算自发磁化的时候，我们采用杨振宁的方法。他计算了在弱磁场h 下，系专业资料整理WORD格式统的自由能，最后令h0，得到磁化强度的表达式：MT,00TTC(1-32)418。NB1 sinh2JTTC而对于平均场近似(MeanFieldApproximation，简称 MFA)所得的磁化强度可以表达为：qJM hMT,hMtanh1B，(1-33)1NBkBT其