伊辛模型

关键词#

伊辛模型; 楞次; 昂萨格; 自旋; 铁磁相变; 临界温度; 自发对称性破缺; 普适类; 临界指数; 蒙特卡洛; 平均场; 二维精确解

"非黑即白的玩具模型"——它却解开了相变的本质#

标题："一个只有上下两种状态的玩具模型"——它却是统计物理里被研究得最透彻的真理来源

伊辛模型听起来简单到近乎幼稚：在一个格点网络上，每个格点放一个只能取"向上"或"向下"两个值的自旋（写成 $+1$ 或 $-1$），相邻自旋倾向于指向同一方向。就这么点设定。很多人第一次见到它会想：这种"非黑即白"的玩具，能说明什么真实的物理？

答案是：它说明了相变这件事的本质。把温度调高，热扰动让自旋乱指，整体不显磁性（顺磁态）；把温度降到某个临界值以下，自旋会自发地大面积对齐，材料整体出现磁化——铁磁态诞生了。关键在于"自发"：在临界温度 $T_c$ 以下，即使没有任何外磁场偏向某个方向，系统也会自己"选边站"。这就是自发对称性破缺（spontaneous symmetry breaking），是从超导、超流到希格斯机制的共同数学骨架。

外场为零时，哈密顿量在全体自旋同时翻转下不变；低温相却必须挑选一个方向，磁化才不是零。对称性还写在方程里，不写在基态里。伊辛是离散的 $\mathbb{Z}_2$，希格斯机制是连续规范对称，骨架同类，细节不同。

更深刻的误解是以为伊辛模型"只关于磁铁"。其实它是一台通用的相变机器：液气共存、合金有序-无序转变、神经网络、舆论传播、蛋白质折叠的某些侧面，都能映射成伊辛型问题。它之所以重要，不是因为它逼真，而是因为它最简单却仍保留了相变的全部精髓。

把格点上的占据数写成 0 或 1，再线性换成 $\pm 1$，液气共存就可以读成伊辛模型。临界点附近起雾的那杯水，和退磁的那块铁，微观力完全不同，映射改的是变量名字，不是普适类。

映射要守住对称性：液气的粒子–空穴对称与伊辛的上下翻转对应；乱改守恒律或长程相互作用，就会滑出三维伊辛类，指数不再承诺一致。

被导师"判了死刑"，被昂萨格精确解救#

标题：被自己导师"判了死刑"的模型

模型的名字来自恩斯特·伊辛（Ernst Ising），但点子来自他的导师威廉·楞次（Wilhelm Lenz）——1920年楞次提出这个磁性模型，并把求解它当作博士课题交给了学生伊辛。

1925年，伊辛在一维情形（自旋排成一条链）求出了精确解，得到一个令人沮丧的结论：一维伊辛模型在任何有限温度下都没有相变。热涨落总能把任何已经形成的有序区域打断，链条永远无法在有限温度自发磁化。伊辛由此（错误地）推断，这个模型在更高维度也不会有相变，于是放弃了它，转行去当中学教师。这个模型差点就此被埋没。

转机来自鲁道夫·派尔斯（Rudolf Peierls）。1936年，派尔斯用一个漂亮的几何论证（后世称"派尔斯论证"）证明：在二维及以上，伊辛模型确实存在有限温度的相变。他的思路是估算"翻转一块自旋畴"所需的能量代价与熵收益的竞争——在二维，畴壁的能量代价足以在低温压住熵的诱惑，有序态得以稳定存在。

派尔斯论证的要点是能量与熵抢一块翻转畴。一维里畴壁只是一个点，能量代价固定，位置却有整条链可选，熵随长度增长，任何有限温度都会把秩序打断——这与伊辛 1925 年的精确解一致。二维里畴壁是一条线，能量随周长增加，低温时能量赢过熵，有序才能站住。平均场看不见这种几何，所以会错误地给一维也安排一个临界温度。

1941 年，克拉默斯与万尼尔（H. A. Kramers、G. H. Wannier）用对偶变换把高温展开和低温展开互相对调，在精确解出现之前就写出了正方格子的临界条件。各向同性时 $\sinh(2J/k_B T_c)=1$，与后来昂萨格的 $T_c = 2J/[k_B\ln(1+\sqrt{2})]$ 是同一句话。对偶不算出比热曲线，却把相变温度钉在一个可检验的数上。

真正的高潮在1944年。拉斯·昂萨格（Lars Onsager）给出了二维零场伊辛模型的精确解——这是统计物理史上的丰碑级成就。他不仅证明了相变存在，还精确算出了正方格子的临界温度 $T_c = 2J / [k_B \ln(1+\sqrt{2})]$，并发现比热在 $T_c$ 处呈对数发散。一个"玩具"模型，第一次让人类在一个真实有相互作用的多体系统里，把相变算到了底。

无量纲临界温度是 $k_B T_c/J = 2/\ln(1+\sqrt{2})\approx 2.269$。二维正方格子上的蒙特卡洛只要扫温度，比热峰就应落在这附近；差一截，多半是有限尺寸或抽样出了问题。

昂萨格用的是转移矩阵：把二维格子看成一条条自旋链依次叠上去，配分函数变成巨大矩阵的本征值问题。考夫曼（Bruria Kaufman）1949 年用旋量语言把同一套代数写得更干净。

自发磁化当时还没印成论文。昂萨格 1948–1949 年在黑板上写出 $M=(1-k^2)^{1/8}$，正式推导由杨振宁于 1952 年发表，$\beta=1/8$ 由此落地。比热的对数发散意味着 $\alpha=0$，能量连续，对比热却没有平均场预言的那种有限跳跃。

外场一旦打开，二维伊辛的精确解就停了：昂萨格解的是零场。有场时对称已被外场打破，相变线上的奇异性结构不同，至今没有同等地位的闭式解。教科书把「二维零场」单独写成丰碑，原因在这里。

从配分函数到临界指数：为什么维度会改变一切#

标题：从配分函数到临界指数

伊辛模型的全部物理藏在一个量里——配分函数 $Z = \sum_{\{s\}} e^{-E(\{s\})/k_B T}$，对所有自旋构型求和，其中能量 $E = -J\sum_{\langle ij\rangle} s_i s_j$（$J>0$ 表示相邻自旋同向降低能量）。一旦知道 $Z$，磁化、磁化率、比热全都能导出。难点在于：自旋构型数随格点数指数增长，直接求和不可能。

物理学家发展了几条互补的攻击路线：

• 平均场近似（mean-field）：假设每个自旋只感受到周围自旋的"平均"作用。它给出相变的定性图像和一组临界指数，但在低维定量错误——它高估了序的稳定性（错误地预言一维也有相变）。
• 昂萨格精确解：仅二维零场可解，是检验一切近似方法的"标准答案"。
• 蒙特卡洛模拟（Monte Carlo）：用计算机随机抽样构型，是处理三维及更复杂情形的主力。Metropolis 算法（1953）与后来的集团算法（Wolff、Swendsen-Wang）让人们能数值逼近临界行为。

三维没有精确解，不是机器不够快，而是二维那套自由费米子结构不再把门打开。平均场假定每个自旋只看见均匀背景，等于先把涨落洗掉；伊辛的上临界维是四，四维及以上平均场才定量可靠，二维、三维的指数会被涨落改写。临界点处关联变长，单自旋翻转的 Metropolis 抽样会越来越慢；集团算法把整块同向自旋一起翻转，三维指数表才有了可核对的位数。任何新算法的第一道标准答案，仍是二维昂萨格解。

1953 年洛斯阿拉莫斯的 Metropolis 等人用蒙特卡洛给二维伊辛抽样，那是计算机进入统计物理的早期现场。他们要的不是再解一次昂萨格，而是证明随机抽样能复现已知的相变。三维没有闭式答案，这条路线后来变成主力。见蒙特卡洛方法。

二维精确指数是分数：$\beta=1/8$、$\gamma=7/4$、$\nu=1$。三维伊辛的 $\beta$ 大约 0.326、$\nu$ 大约 0.630，夹在二维分数与平均场之间。维度往上走，涨落变弱，指数被拉向平均场，却在四维以下到不了。

有限尺寸缩放把“格子有多大”也写进标度律：关联长度撞上边界时，比热峰会变钝并略微移位。蒙特卡洛必须外推到热力学极限，否则测到的“临界温度”只是有限格子上的有效峰位；昂萨格的无穷格子答案，正是这条外推线应当指向的终点。

临界点附近，物理量按幂律发散或消失，指数称为临界指数（如磁化 $\sim (T_c-T)^\beta$，磁化率 $\sim |T-T_c|^{-\gamma}$，关联长度 $\sim |T-T_c|^{-\nu}$）。昂萨格的二维解给出 $\beta = 1/8$、$\gamma = 7/4$ 等精确分数值，与平均场的 $\beta=1/2$ 明显不同——这告诉我们：维度真的会改变临界行为，平均场在低维不可信。

普适性：为什么沸水和退磁的铁在临界点"长得一样"#

标题：普适性——为什么水和磁铁在临界点"长得一样"

伊辛模型最震撼人心的遗产，是普适性（universality）概念的诞生。临界指数不再属于某一种材料，而属于一类对称性和一种维度。

实验上有个惊人的事实：液气临界点附近的临界指数，和单轴铁磁体（伊辛型）临界点附近的临界指数，在实验误差内完全相同。两种看上去毫不相干的物理系统——一杯沸腾的水和一块退磁的铁——在临界点表现出同一套数字。这绝非巧合。

肯尼斯·威尔逊（Kenneth Wilson）在1971—1972年用重整化群（renormalization group）给出了解释，并因此获得1982年诺贝尔物理学奖。核心洞见是：临界点处关联长度发散，系统在所有尺度上自相似，于是微观细节（晶格类型、相互作用的具体形式）被"洗掉"，只剩下少数几个全局特征——空间维度和序参量的对称性——决定临界指数。具有相同维度和对称性的系统属于同一个普适类（universality class），共享同一套临界指数。液气共存与单轴铁磁同属"三维伊辛普适类"，所以它们的临界指数一致。

重整化群的操作可以想成：把格子分成小块，每块用一个有效自旋代替，再看这个粗粒化系统的耦合如何随尺度跑。临界点是跑动里的不动点。微观的晶格形状、次近邻耦合，只要不改对称性和维度，都会流到同一个不动点，指数因此相同。平均场相当于在不动点处丢掉了涨落。

2014 年 Kos、Poland 与 Simmons-Duffin 用混合关联函数的共形 bootstrap，把三维伊辛临界点收成参数平面上的一座小岛：$\Delta_\sigma\approx 0.518$、$\Delta_\varepsilon\approx 1.413$，对应 $\eta\approx 0.036$、$\nu\approx 0.630$。这不是又一次蒙特卡洛拟合，而是交叉对称与幺正性把允许的三维共形场论逼进一个角落。实验上的三维伊辛类物质，指数落在同一角落里。见相变与临界现象。

普适类预言可以拿去实验室对质。氙的液气临界点、β-黄铜的有序无序、单轴铁磁体，测量到的 $\beta$ 都靠近 0.32 而不是平均场的 0.5。指数对上了，才叫同一类，不是因为它们看起来都像磁铁。

至今仍活跃的前沿包括：用共形场论（conformal field theory）和"共形 bootstrap"方法，对三维伊辛模型的临界指数做超高精度的解析-数值确定（2014年后该方法把三维伊辛指数算到了前所未有的精度）；以及把伊辛型自发对称破缺的思想，迁移到量子相变（绝对零度、由量子涨落驱动）——横场伊辛模型就是量子相变的"果蝇"，是研究纠缠与量子临界性的标准平台。

横场伊辛模型在绝对零度仍可相变：量子涨落扮演温度的角色，序在参数空间里消失。它被称作量子相变的果蝇，是因为二维经典伊辛的精确解，通过量子-经典对应，给一维量子链提供了可核对的极限。热涨落和零点运动可以讲成同一套配分函数，局部规则如何长出全局秩序，经典与量子共用这份语法。

量子临界点附近，纠缠熵常带面积律的对数修正；相变因此不只是热力学奇异性，也接到可测的纠缠结构上。同一套对称性破缺语言，实验室里既可用低温磁体，也可用离子阱或超导比特去模拟。

一个模型，连起磁性、宇宙、生命与机器#

标题：一个模型，连起磁性、宇宙、生命与机器

伊辛模型是少数几个真正"跨学科"的物理结构，它的触角伸向许多看似无关的领域：

• 宇宙学与粒子物理：宇宙早期的对称性破缺（如电弱相变、希格斯机制）在数学结构上与铁磁相变同源——安德森（Philip Anderson）在1962年首先点出这一类比，为电弱统一理论铺路（见《相变与临界现象》《希格斯机制》）。
• 凝聚态深化：把伊辛模型推广到带"无序"和"阻挫"的相互作用，就得到自旋玻璃——一种没有简单有序态、却有无穷多亚稳态的复杂相，是统计物理最难啃的问题之一（见《自旋玻璃与无序系统》）。
• 神经科学与机器学习：Hopfield 网络本质上是一个带可调相互作用的伊辛模型，用自旋构型存储记忆；玻尔兹曼机（Boltzmann machine）则把伊辛配分函数直接搬进了机器学习。2024年的诺贝尔物理学奖部分正是表彰这条"统计物理→人工神经网络"的脉络。
• 复杂系统：从舆论极化、流行病传播到鸟群相变，"局部相互作用涌现全局秩序"的伊辛思路被反复借用。

一个"非黑即白"的玩具，最终成了人类理解"群体如何涌现秩序"的通用语言。这正是物理学追求简单模型的胜利。

跨域连接#

• 相变与临界现象：伊辛模型是普适性概念的摇篮：液气临界点与单轴铁磁在临界点附近共享同一套指数，微观细节被"洗掉"，只剩维度与对称性。推论：改晶格结构而不改普适类，临界指数就不变；改序参量对称性，指数必变——这是重整化群最锋利的预言。
• 神经网络：霍普菲尔德网络本质上是一个耦合可调的伊辛模型，记忆被存储为能量地形上的吸引子；玻尔兹曼机则直接把配分函数搬进机器学习。推论：网络的记忆容量应有类似相变的崩溃阈值——存得超过极限，回忆质量不是渐差而是骤降，这正是统计物理预言的行为。
• 蒙特卡洛方法：三维伊辛模型没有精确解，而构型数随格点指数增长，随机抽样成了理论研究的正规军。推论：临界指数的数值精度随算法代际不断提高——一个"玩具模型"的指数表，记录的其实是计算方法的进步史。
• 自旋玻璃与无序系统：把伊辛模型中整齐的耦合换成随机符号的耦合，有序基态就碎成指数级多的亚稳态。推论：磁化行为从"有确定的转变温度"变为"依赖历史的冻结"——无序加阻挫，把最简单的磁性模型变成统计物理最难的问题之一。
• 社会科学：舆论极化、流行病传播常被映射为伊辛型模型，借的是"局部相互作用涌现全局秩序"的机制。但要写明边界：这类比只在机制层面有效，社会系统的耦合并非常数，定量的临界指数预言多数不成立——它提供问题框架，不提供可靠数值。

参考文献#

• Onsager, L. "Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition." Physical Review 65 (1944): 117–149. DOI: 10.1103/PhysRev.65.117
• Brush, S. G. "History of the Lenz-Ising Model." Reviews of Modern Physics 39 (1967): 883–893. DOI: 10.1103/RevModPhys.39.883
• Wilson, K. G. "The renormalization group and critical phenomena." Reviews of Modern Physics 55 (1983): 583–600. DOI: 10.1103/RevModPhys.55.583
• Kramers, H. A., & Wannier, G. H. "Statistics of the Two-Dimensional Ferromagnet. Part I." Physical Review 60 (1941): 252–262. DOI: 10.1103/PhysRev.60.252
• Kaufman, B. "Crystal Statistics. II. Partition Function Evaluated by Spinor Analysis." Physical Review 76 (1949): 1232–1243. DOI: 10.1103/PhysRev.76.1232
• Yang, C. N. "The Spontaneous Magnetization of a Two-Dimensional Ising Model." Physical Review 85 (1952): 808–816. DOI: 10.1103/PhysRev.85.808
• Kos, F., Poland, D., & Simmons-Duffin, D. "Bootstrapping mixed correlators in the 3D Ising model." Journal of High Energy Physics 2014 (11): 109. DOI: 10.1007/JHEP11(2014)109

延伸阅读#

• McCoy, B. M., & Wu, T. T. The Two-Dimensional Ising Model. Harvard University Press, 1973.