自旋模型简述

自旋模型简述1、自旋的基本概念与表述自旋是电子的基本性质之一，是电子内禀运动量子数的简称。电子自旋的概 念是由Uhlenbeck和Goudsmit为了解释碱金属原子光谱的精细结构以及反常 Zeeman效应而提出的。他们认为电子的运动与地球绕太阳运动相似，电子一方 面绕原子核运动，从而产生了相应的轨道角动量；而另一方面它又有着自转，其 自转的角动量为 /2，并且它在空间任何方向的投影都只能取两个值，即/2（也就是自旋向上和向下两个状态），与自旋相对应的磁矩则是e /2mc。当然，这样带有机械性质的概念是不正确的，而自旋作为电子的内禀 属性，是标志电子等各种粒子（如质子、中子等）的一个重要的物理量。对于自旋这个自由度，我们一般用算符表示（这里的记号表示算符，在下文中为了简便我们将略去这一记号）。因为自旋角动量与轨道角动量有着相同 的特征，所以一般也认为它们具有相同的对易关系，即S s=i S。在这里我们 引入泡利算符S=/2。由于S沿任何表象的投影都只能取/2两个值，即沿任何方向的投影只能取1这两个值，所以泡利算符 的每个分量都可以用 2 2的矩阵来表示。我们一般采用分量对角化的表象，得到其矩阵表示：Z0 10 i10（x 10 5 y i 0 5 z 01 51 1这样的表示就是著名的Pauli矩阵。2、自旋模型的形式物质的磁性与自旋模型 由于原子核的磁矩很小，物质的磁矩可以看成其轨道磁矩和自旋磁矩之和。 电子的总磁矩（轨道磁矩+自旋磁矩），直接体现为物质的宏观磁性。而对于过渡 金属的原子或离子，因为轨道角动量的冻结，其磁性主要来源于未配对电子的自 旋磁矩。对于物质的磁性，很早以来就有着广泛的研究，比如 Langevin 的顺磁理论， Wiess的分子场理论，Bloch的自旋波理论。这些理论中，原子（离子）都具有磁 矩，而磁矩之间存在着一定的相互作用。在绝对零度以上，每个原子都在做热振 动，磁矩的方向也在作同样的振动，而磁矩间的相互作用又使得每个磁矩趋向于 某种有序的排列，这就是物质宏观磁性的来源。磁矩之间的相互作用有很多种：1、经典的磁偶极子之间相互作用。2、交换相互作用（也称直接交换相互作用）。氢分子模型、海森堡交换模型 就是采用这一类的相互作用。交换相互作用没有相应的经典对应，它来自于电子 间的库伦作用以及量子力学的全同粒子系特性。下面所说的其它种类的交换相互 作用也是基于这样的原理。3、超交换相互作用（也称间接交换相互作用）。这种相互作用由 Kramers 于 1934 首先提出，用于解释反铁磁性的自发磁化的起源。它是阳离子的电子以氧 离子的 p 电子为媒介进行间接的相互作用。4、RKKY相互作用。这种相互作用由Ruderman、Kitt el、Kasuya以及Yosida 提出的，是一种以巡游电子为媒介，使得磁性原子（或离子）中的局域电子自旋与 其邻近的磁性原子（或离子）中的局域电子自旋产生的交换相互作用。5、双交换相互作用。以氧离子为媒介，两个不同价态的过渡族粒子间之交 换相互作用。在锰氧化物中，这种相互作用就起到了十分重要的作用。以上的几种相互作用中，除了偶极间的相互作用是一种经典的相互作用，而 其余的几种交换相互作用却是基于体系的量子特性，即全同粒子的特征。这样的 相互作用，在我们研究物质的磁性以及其它以磁性相关的性质，或者以磁性变化 为主导的相变时，起着至关重要的作用。对于这类的磁性原子体系，我们认为它们位于某种晶格格点位置上，通过磁 矩进行相互作用，我们可以建立一种自旋模型来进行描述，其最基本的形式可以 写成如下的哈密顿量：HJ SxSxSySyJzSzSz ,i, j i ji ji ji ji ji j(1-2)这里的s.是自旋算符，上标x, y以及z为s.在三个方向的投影，可以分别对应1 1泡利矩阵的三个分量。J和J z代表着位于格点i和格点j上原子磁矩之间的 1jij相互作用，这样的相互作用可以是我们之前所说的几种交换相互作用中的任意一种。对于自旋算符S.，我们可以做如下的变换：s =sx isy,这里s+算符将自旋1向上的状态转换为自旋向下的状态，而S-算符则是将自旋向下的状态转换为自旋向上的状态。哈密顿量(1-2)可以重新表达为如下的形式：JS S S S2 i j i j i jJz SzSz oi, j i j1 j(1-3)对于哈密顿量(1-2)和(1-3),当我们取J二Jz时，就成了 Heisenberg模ij ij型；取J =0时，这个哈密顿量是对角化的，自旋只有向上和向下两种取向，即 ij/2，这样的模型就是著名的Ising模型；而如果Jz =0时，得到的则是XYij模型。Heisenberg模型简述(1-4)在Weiss提出分子场假说20年后，Heisenberg提出了电子间的交换相互作 用导致了自发磁化的产生，并且按这一模型，也就是Heisenberg模型计算了自 发磁化随温度变化的性质，为铁磁性量子理论的发展奠定了基础1。假设系统 是均匀的，只考虑最近邻相互作用，Heisenberg模型的哈密顿量可以写成如下 形式：H J s s h s ,i j Bi其中si表示位于格点i处的自旋，是矢量(Sx, Sy, Sz), J是交换相互作用常数,1这里我们取J0代表铁磁相互作用，而J0代表反铁磁相互作用，为Bohr磁B矩，h是外磁场。图氢分子电子云分布示意图。我们知道交换相互作用和任何一种经典的相互作用都没有对应，它是一种量 子行为，是由粒子的全同性产生的相互作用。为了理解这种特殊的相互作用，通 常以氢分子这种最简单的模型为例来说明。如图所示，一个氢分子的系统，由两个原子核a b以及两个电子1, 2组 成。r代表粒子m与粒子n之间的距离，其中m，n=a，b，1，2，且m n。氢 mn分子的哈密顿量可以写为：2H22m 1e2ra1e2ra2e2rb1e2rb2e2r12e2rab(1-5)H1a22me221ra1H22.2e2(1-6)b2m2r b2e2e2e2V1,2rrr12a2b1其中H (1)与Hb(2)是两个孤立氢原子的哈密顿量，V(l, 2)是两原子间的相互作 b用，而哈密顿量(1-5)中的最后一项是常数，对本征态没有影响。我们用(1)与b(2)分别代表H (1)与Hb(2)两个本征态，则(1) b(2) babab是下列Schr dinger方程的本征态：H 1 H 2 E ，(1-7)ab000有粒子的全同性可知(2)(1)也是Schr dinger方程(1-7)的本征态(相ab当于将电子1与电子2互换)。考虑到电子自旋的波函数，以及对Fermi子波函 数反对称的要求，我们可以写出氢分子基态波函数的近似形式：1 2 1 2 1,2I abbaA(_8)1 2 1 21,2 II abbaS这里(1,2)与(1,2)分别是电子波函数的单态和三重态。所以， 态两电子ASI反平行，而n态两电子平行。从哈密顿量(1-5)出发，这两个波函数所对应的 能量为：E 2E 竺 UJI 0 r 12ab，(1-9)E 2E 竺 U_JII 0 r 12ab其中U是库伦排斥能，而J是交换相互作用能，为a，b原子波函数的重叠积分：这里1 2V 1,22d12 V 1,2b2 1d db12(1-10)d为电子1，22的全空间积分，而 显然是在0和1之间的。因此对于能量E和E我们，当J0时，E 0IIII II时，E E，两电子趋向平行排列。这和我们之前对于Heisenberg模型的描述是II I致的。从以上对于氢分子模型的计算，交换相互作用能是量子力学的结果，虽然它 和库伦势同样和电子间的库伦相互作用有关，但我们并不能找到其经典的对应。 从直观上看，交换相互作用能来源于电子波函数的交叠，在交叠区域电子是不可 识别的，如果我们仍然按照处理宏观粒子的方式对其标识，则不可避免的会引入 交换相互作用能。从Heisenberg模型我们知道了 Weiss的分子场实际上是交换相互作用 的平均场。而Heisenberg模型作为一种量子磁性最基本的相互作用模型已 经被运用于极为广泛的领域。在某一些物质中，部分未配对电子会被局域在离子实周围。当这部分局域电 子数目较多（一般来说大于3）的时候，我们可以在一定程度上忽略其量子效应， 即电子的自旋s可以不再被看成是算符而是一个经典的磁矩。比如在锰氧化物 中，Mn离子的3d轨道上的电子可分为两类，一种叫e电子，另一种叫t电子。g2g其中e电子是可以通过双交换相互作用在Mn离子之间移动的，而t电子则是局g2g域在Mn离子周围的。这一体系的哈密顿量可以简单的写成2:Hta d d J s S J SS,(1-11)i i aHi i AFi i aiaiia其中s.代表e电子，用Pauli算符；而s.代表t电子，由于有三个t电子局igi2g2g域在Mn离子周围，所以S可以看作是自旋为3/2的磁矩。哈密顿量的第一项代i表双交换相互作用；第二项是Mn离子内部的Hund相互作用，也就是e电子(s ) gi与t电子(S )的相互作用；第三项不同格点t电子间的Heisenberg相互作用。 2gi2gHeisenberg 模型的许多拓展也在物理学各个领域中得到运用。如果我们在 Heisenberg 模型引入了各向异性能，就得到了各向异性 Heisenberg 模型 (Anisotropic Heisenberg Model):H J s s h s Dsx 2 ，(1-12)i j B iii,jii在此哈密顿量中，我们引入了 x方向的单轴各向异性晶体场D。在Heisenberg模型中还可以加入了 Dzyaloshinskii-Moriya相互作用(DM 相互作用)3,4:H J s s D a r s s h s，(1-13)i ji i j B ii,ji,ji其中r是位于格点i的原子(离子)的的位移，Da(r )是DM矢量。这一模型被认 ii为是多铁性材料5-7以及磁性金属材料8-13中螺旋自旋序产生的根源。关于 螺旋自旋序及其相关问题我们将在第二章内做详细地讨论。Ising模型简述Lenz曾向他的学生Ising提出一个研究铁磁性的简单模型，而Ising于1925 年发表了他对此模型求解的结果14，所以这个模型被称为 Ising 模型。当时 Ising只做出了该模型一维下的严格解，在一维情况下并没有自发磁化的发生。 另外他还由此错误地推断出在更高维的情况下，这个模型也不存在自发磁化。这 个推断在后来被证明是错误的。1936年Peierls论证了二维或三维的Ising模 型存在着自发磁化，虽然当时他并没有能够给出模型的严格解15。 1944 年， 当Onsager给出了二维Ising模型的严格解之后，Ising模型开始引起人们广泛 的关注16,17。这次求解是相变理论发展上的一个重要进展，它第一次清楚地 证明了从没有奇异性的哈密顿量体系出发，在热力学极限下能导致热力学函数在 临界点附近的奇异行为，而 Onsager 本人也因此获得了诺贝尔奖。在此之后很多 人又相继发表Ising模型的各种不同解法，Bax ter甚至有篇论文叫Using模型 的第399种解法但至今没有被学术界公认的三维Ising模型精确解。甚至有 人发表论文证明无法解出三维Ising模型的精确解，因为三维Ising模型存在拓 扑学的结构问题。人们通常用分子场理论及其改进理论、高温级数展开、低温级 数展开、重整化群理论、蒙特-卡罗模拟等近似计算三维I sing模型的居里温度 和临界指数，而其中Wilson于1971年发展的重整化群理论能以较高精度计算三 维Ising模型的近似结果18-20。我国科学家张志东提出三维“Ising模型” 精确解猜想21。张志东的出发点就是拓扑学中的一个常识：低维空间的扭曲和 纽结可以被高一维空间的旋转打开。通过引入第四卷曲起来的维与本征矢量上的 权重这两个猜想作为处理三维Ising模型拓扑学问题的边界条件，并应用这些猜 想用自旋分析法评估了三维简单正交晶格Ising模型的配分函数。当系统的对称 性越高，居里温度也越高。他猜测三维系统具有最高对称性的简单立方 Ising 模型具有最高的居里温度黄金解,在二维系统具有最高对称性的正方Ising模型 具有最高的居里温度白银解。获得的结果具有一定的对称性和美学价值，并可部 分返回到二维和一维的结果。当然，推定的精确解正确性取决于猜想的正确性， 而且其与学术界通常接受的评价标准尚不完全吻合，有待于对相关的物理本质作 进一步探讨。因此，这一工作目前还只是停留在猜想阶段。今天的Ising模型根本不再是Ising博士论文中的模样。每年差不多有6000 篇左右的论文研究这一模型。除了铁磁性之外，该模型还应用于很多方面，如合 金中的有序-无序转变、液氦到超流态的转变、液体的冻结和蒸发、晶格气体、 玻璃物质的性质，甚至于神经网络蛋白质折叠、生物膜场论甚至社会现象等广泛 的领域。通过上述介绍，我们知道三维Ising模型尚未得到严格解，而一维和二维情 况下的解法确是多种多样的。在这里，我们将给出Ising模型的严格解，采用的 是1941年Kramers和Wannier提出的转移矩阵方法(Transfer Matrix Method)22-24。然后简要地说明二维Ising模型严格解的主要结果，并且同平均场理论所得的结果进行对比。I I I SII1“图一维Ising模型示意图。对于如图所示的Ising模型，自旋只能取向上或向下两个分量，它可以看作 是Heisenberg模型的一种简化。当只考虑最近邻的交换相互作用，并认为这种 相互作用在不同磁矩间是相同的，用常数J表示。和Heisenberg模型相同，当 J0时，代表铁磁的交换相互作用，它使得近邻自旋有着同方向排列的趋向；当 Jo代表铁磁相互作用，B为BohrB磁矩，h是外磁场。对于一维情况，每个自旋只有两个近邻。现在采用周期性边 界条件，即s =s, N为晶格中的自旋数目。现将一维晶格弯成一个环，当NN+11时，边界效应将不会影响到体系的热力学性质。根据如上的条件，可将哈密顿量(1-14)写为:J ssi i 1ih-Bs s2i i 1i(1-15)其相应的配分函数为:Q T,hsi1exp1 i 11kTBJssi i 1(1-16)在这里我们引入矩阵P,其矩阵元定义为:silP!si1-1exp k TBJssi i 1(1-17)因为s与s都能取1两个值，所以P是2 2的矩阵:ii+1P汽1咔i11咔i11汽1咔i11叩i11e jh k TBBeJkBTeJ k TJ Be jBh kBT(1-18)于是配分函数(1T6)可以重新写成:Q T,h(s1lPls2Xs2lPls3/s1 sN 1,(s |P n |Tr P ns1 1|P|s Xs |P|s)N * N i(1-19)将P矩阵对角化得,(1-20)0P 0 ，和 即为矩阵P的本征值，由下面的久期方程决定,+ -(1-21)eJ Bh；kBTe JkBT0e JkBTeJ Bh /kBT其解为:(1-22)eJ kBT cosh 丄cosh? 丄2e J kBT sinh-2k Tk Tk TBBB要注意的一点就是。+ -现在将等式(1-20)代入(1-19)，配分函数可以表达为:(1-23)Q T,h N N N 1所以，当N 时，我们得到:lim1 lnQ T,h ln n N(1-24)J ln cosh 丄 k Tk TBB:cosh2 丄 2e J kBr sinhk TB k T1BB即配分函数有P矩阵较大的本征值决定。体系的自由能和总极化强度分别为:FT,h kTQ T,hB NJ k T coshBhBk TB(1-25):cosh2hBk TBsinhk TBsinh(1-26)其它的热力学函数也可同样由自由能求出。如图所示，在计算中我们选取交换相互作用常数J=1k$,对于一切T0都有 M(T, 0)=0,也就是说I sing模型在一维的情况下不存在自发磁化，不会发生顺 磁-铁磁转变。从物理上看，任何温度下自旋的平均取向由两个对抗的因素相互竞争决定，即能量趋向最小而熵趋向最大，使得自由能达到最小值。在一维情况 下，由于近邻数低，使得自旋排列在同方向的倾向不足以对抗使熵极大的倾向， 结果在任何有限温度下都不能形成自发磁化。1.00.5c-0.50510r = ioK图一维Ising模型在不同温度下，磁化强M随外场h的变化曲线。图一维Ising模型在有限温度下长程序被破坏的示意图。如同上文所说，当自旋排列在同方向的倾向不足以对抗使熵极大的倾向时， 自旋往往会在一个较小的尺度内保持着同方向的排列，形成所谓短程序 (Short Range Order)，而在较大的尺度内失去这种有序的状态，也就是破坏了所谓长程 序(Long Range Order)。当我们使用蒙特卡洛方法(Mon te Carlo Met hod)来计算 维Ising模型，常常得到如图所示的结果，在某个小范围内，如从格点1到格 点4(或者从格点5到格点10)体系可以看作存在自发磁化，而在整体上看(从 格点1到格点N)向上的自旋和向下的自旋数目在统计上看是相等的。对于二维Ising模型，我们考虑正方晶格，每个自旋有4个最近邻。在零磁场下，系统的自由能可以表达为：lnQ T,0F T,0N(1-27)其中 cosh ( )=cosh2cosh2 -cosh sinh2 sinh2二t an-1e-2 J。系统内能则可以写为:u T ,0J coth2 J 1(1-28)这里的K1 (m)是第一类椭圆积分,0 J1 m 2 sir?(1-29)(1-30)其中m=sinh2 J/cosh22 J, m =2tanh22 J-1。而临界点由下式确定:k T 2.269J。B CC T,0 BNkB所得热容量为：(1-31)这样的热容量在临界点处具有对数发散的奇异性。计算自发磁化的时候，我们采用杨振宁的方法25。他计算了在弱磁场h 下，系统的自由能，最后令h 0,得到磁化强度的表达式:(1-32)M T,00T TcN1 sin!2 J 4 18T TBC而对于平均场近似(Mean Field Approximation,简称MFA)所得的磁化强度可以表达为:M T,hB十 hqJMtanh1k TB(1-33)其中q是最近邻自旋的数目，对于二维正方晶格来说q=4。HK)图严格解与平均场近似所得二维sing模型磁化强度的比较。在图中，与平均场近似所得的解(Tc=4K)相比，严格解(TC=有着更低的临界温 度而且磁化强度在T TjO有着更陡的温度变化率。平均场近似忽略了系统 的涨落，而涨落是倾向于破坏有序的，所以在平均场近似下所得的-是高于实 际体系的，磁化强度的变化也反映的这一特点。在这里我们需要注意，在平均场 近似所得的结果中，Ising模型在一维的情况下存在着自发的磁化，这个结果是 错误的。在前文我们提过，在一维系统中由于近邻数目低，系统的涨落完全抑制 了有序，而平均场近似忽略了系统的涨落才得到了有序相。Ising模型自Ising提出后，有了很多的发展，不仅是解法的多样化，其具 体形式也发生了不少变化。Ising模型的一些重要拓展成为描述相变(比仅仅是 磁性)等问题重要工具。比如说横场伊辛模型(Transverse-field Ising model， 简称TIM)，它是于1963年由Gennes提出26。TIM是在Ising模型的基础上考 虑了横向外场的作用，所谓横向外场，是指外场的方向垂直于Ising模型中 自旋投影的方向。在这里横向外场可以看作是晶体内部横向遂穿效应的一种等 效，从而可以应用于零温的量子相变。TIM的哈密顿量可以写成：H J szszsx，(1-34)i jii, ji这里 是横向场，也可以看作是遂穿积分，决定了从一个势能极小态到另一个势 能极小态的遂穿几率。TIM被广泛运用到多种体系，比如量子自旋玻璃，量子弛 豫，量子磁滞，量子铁电等，详细的内容我们将在第三章中做具体的阐述。我们再说Ising模型的另一种拓展形式弹性伊辛模型(Elastic IsingModel，简称EIM)。我们知道，晶体中的原子(离子)受到其周围的原子(离子)作 用，被束缚在晶格格点，也就是平衡位置周围做微小的振动，可以用劲度系数 (stiffness factor)k来描述。在EIM中，交换相互作用J不再是常数，而是同 离子间的距离存在某种类似于弹性能的关系，新的弹性作用使得原子产生位移， 出现新的平衡位置。当自旋以某种有序状态排列时，弹性交换相互作用使得晶格 中的某些原子产生同向的位移。这样的特性对于解释某些自旋共线排列的多铁性 材料有着重要的意义，我们将在第四章中作具体的介绍。XY模型简述XY模型是I sing模型的延伸，允许自旋在(xy)平面内转动。当仅考虑最近领交互作用时，体系的哈密顿量可以表示为:H J S S J cos ， (1-35) i j i jij ij这里J表示交换作用常数，S为位于格点i的自旋，求和包括所有的最近邻晶格i点，&为第i个自旋相对于x轴的夹角。对于XY模型，最为重要的就是它的 iKT相变。早在1966年，Mermin和Wagner就严格证明了对于二维或二维以下具有连 续对称性和短程相互作用的自旋体系，在有限温度下由于热扰动或量子波动不可 能形成长程序。所以在二维XY模型中，不可能存在长程的磁有序。1973年，J. M. Kosterlitz和D. J. Thouless提出了拓扑序和拓扑性相变的概念，建立了 二维XY模型的KT相变理论。他们在研究二维XY模型时发现：低温下局域自旋 会形成正反涡旋两两配对，当温度升高到相变的临界温度时，这种配对被热运动 所拆散，出现独立运动的正反涡旋。这种新的相变被称为Kos terli tz-Thouless 相变(KT相变)。它的特点是：尽管体系在相变点两侧的序参量总为零，但是温 度的变化会引起关联函数和拓扑性质的本质性改变，是一种拓扑相变。具体来说， 体系存在两个相：低温下的准长程有序，关联函数以幂指数形式衰减，所以也称 为代数磁有序；高温的无序相，关联函数以e指数形式衰减。前人的工作表明二 维四方晶格铁磁XY模型的KT相变温度严。自旋模型的研究方法自旋模型看似简单，但严格的求解却十分困难。就连其最简单的形式Ising 模型的求解也经历了重重困难，而直到现在Ising模型在三维情况下的解依然是 一个难以攻克的难关。对于实际的体系，模型将更加复杂，而严格求解也更为困 难，因此一系列近似求解的方法便发展起来。近似的方法可以简化计算，但免不 了产生这样那样的误差，有着各自的优缺点，因此对于方法的选用就须要格外注 意。下面我们将介绍几种近似的处理方法：1) 平均场理论(Mean Field Approximation，简称MFA)。平均场近似的方法可以说是最早被提出求解自旋模型的方法，比如说Weiss 的分子场理论就是一种平均场理论。这个方法忽略了各种格点自旋之间的关联， 把近邻自旋的作用用一种集体的平均作用来进行代替，从而忽略了涨落。这样就 可以把多粒子相互作用的系统转化为近独立粒子系统，大大简化了计算。但是由 于这种近似忽略了涨落，在涨落较强的体系中，比如一维体系或准一维体系中， 它往往难以得到令人满意的结果。我们前文也说了，用平均场近似处理一维 Ising模型时可以产生自发磁化，这是与事实完全相背离的。所以在遇到类似的 情况时，我们要尽量的避免使用平均场近似，或者对其形式做出某种修正。2) 相关有效场理论(Correlated Effective Field Theory)。这种方法在平均场理论的基础上考虑了静态自旋关联的多体效应，通过微分 算符技术和退耦近似等数学处理来对自旋模型的计算进行简化27，所以在很大 程度上比平均场近似更能反映自旋模型所描述的物理现象的微观本质。3) 级数展开法(Expansion Met hod)。这种方法就是利用适当变量的幂级数展开来计算配分函数，只要展开项足够 多，就可以很好地反映临界点附近的热力学性质。对于某些较为理想的模型仅仅 计算几项就足够了，对于大多数模型则十分繁琐，比如说计算三维Ising模型时， 需要40多项的展开式才能达到令人满意的精度28,29。4) 重整化群方法 (Renormalization Group Technique)。重整化群理论建立在标度理论和普适性的基础上，提供了从微观上计算临界 指数的系统方法。重整化群的基本思想是：在临界点关联长度趋于无穷大，因此 体系应该具有尺度变换下的不变性，由此不需要直接计算配分函数，而是找尺度 变化下的不变性，从而确定临界点并计算出临界指数18-20。重整化群方法对 于自旋模型的计算可以达到一个很高的精度。5) 蒙特卡罗模拟(Monte Carlo Simulation)。蒙特卡罗模拟方法是建立一个统计模型，通过随机抽样来进行数字模拟实 验，从而可以得出这个系统的各种热力学性质。Metropolis等人于1953年在研 究限制在一个盒子中的液体和浓稠气体的行为时，提出了一种有效的蒙特卡罗抽 样方法(后来人们称之为Metropolis抽样方法)30。蒙特卡罗方法的精度抽样 样本的大小和计算时间有关，当抽样样本足够大，计算时间相当长的时候可以得 到很高的精度。而在不需要高精度的时候也可以将它们相应地减小，所以这是一 种很灵活的方法。蒙特卡罗方法的应用非常广泛，不但用于自旋模型的计算，还 被推广到统计物理与凝聚态物理研究的很多领域中31-35。参考文献：1. W. Heisenberg, Z. Physik, 49, 31 (1928).2. E. Dagotto, T. Hotta, and A. Moreo, Phys. Rep. 344, 1(2001).3. I. Dzyaloshinskii, Sov. Phys. JETP 19, 960 (1964).4. T. Moriya, Phys. Rev. 120, 91 (1960).5. I. A. Sergienko and E. Dagotto, Phys. Rev. B 73,094434 (2006).6. I.A.Sergienko,C.en,andE.Dagotto,Phys.Rev.Lett.97,227206(2006).7. Q. C. Li, S. Dong, and J. -M. Liu, Phys. Rev. B 77, 054442 (2008).8. C. Pfleiderer, S. R. Julian, andG. G. Lonzarich, Nature (London) 414, 427 (2001).9. N. Doiron-Leyraud, I. R. Walker, L. Taillefer, M. J. Steiner, S. R. Julian, and G. G. Lonzarich, Nature (London) 425, 595 (2003).10. C.Pfleiderer,D.Reznik,L.Pintschouvius,H.v.L hneysen,M.Garst, and A. Rosch, Nature (London) 427, 227 (2004).11. Y. Ishikawa, K. Tajima, D. Bloch, and M. Roth, Solid State Commun.19, 525 (1976).12.13.14.15.16.17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32.G. Shirane, R. Cowley, and C. Majkrzak, Phys. Rev. B 28, 6251 (1983).M. Uchida, Y. Onose, Y. Matsui, and T. Tokura, Science 311, 359 (2006). E. Ising, Z. Phys. 31, 253 (1925).R. E. Peierls, Proc. Camb. Phil. Soc. 32, 477 (1936)L. Onsager, Phys. Rev. 65, 117 (1944).B. Kaufman and L. Onsager, Phys. Rev. 76, 1244 (1949).K. G. Wilson, Phys. Rev. B 4, 3174 (1971).K. G. Wilson and J. Kogut, Phys. Reports 12C, 75 (1974).K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975).Z-D Zhang, Philosophical Magazine, 87(34), 5309 - 5419 (2007).H. A. Kramers and G. H. Wannier, Phys. Rev. B 60, 252 (1941).R. J. Baxter, ExactlySolvedModelsinStatisticalMechanics(Academic Press, London, 1982).北京大学物理系量子统计物理学编写组，量子统计物理学(北京大学出 版社，1987).C. N. Yang, Phys. Rev. 85, 809 (1952). de Gennes, Solid State Commum. 1, 132 (1961).T. Kaneyoshi, I. P. Fittipaldi, R. Honmura, and T. Manabe, Phys. Rev B 24, 481 (1981).M.F.Sykes,J.W.Essam,andD.S.Gaunt,J.Math.Phys.6,283(1965).M.F. Sykes, D. S. Gaunt, J. W. Essam, and C. J. Elliott, J. Phys. A 6, 1507 (1973).N.Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, J. Chem. phys. 21, 1087 (1953).K.Binder,MonteCarlomethodinStatisticalPhysics(Springer-Verlag, 1978).K. Binder,Application of the Monte Carlo Simulation in StatisticalPhysics (Springer, 1987).33. K. Binder and D. M. Heermann, Monte Carlo Simulation in StatisticalPhysics (Springer, 2002).34. D. P. Landou and K. Binder, A Guide to Monte Carlo Simulation inStatistical Physics(Cambridge, 2000).35. J. S. Liu, Monte Carlo Strategiesin ScientificComputing(Springer,2001).