01 它问的是什么问题(The Question It Poses)
一块铁在 770 摄氏度以下是磁铁,以上就不是了。这件事看起来该由铁的化学性质负责——毕竟不同金属的这个温度差得很远。
但 Wilhelm Lenz 在 1920 年问的是另一个问题:要让"整体突然有了方向"这件事发生,最少需要几条规则? 能不能干脆不要化学,只留下"每个单元有两种状态"和"挨着的单元倾向于一致"这两句话?
他把这个问题交给学生 Ernst Ising。答案后来变得比问题重要得多:这个几乎什么都没有的模型,不但会发生相变,而且它的临界行为跟真实的磁体、跟液气临界点、跟合金的有序化完全一样 。它成了"普适类"这个概念的原型。
02 规则本身(The Rules)
取一张方格纸,每一格放一个只有两种取值的变量:+1 或 −1 (想成"朝上/朝下",或者"赞成/反对")。
相邻两格如果取值相同,这一对贡献能量 −J (J > 0,是一个正数);如果不同,贡献 +J 。也就是说:跟邻居一致更省能量。
整个系统的能量就是所有相邻对的贡献之和。
给定一个温度 T 。随机挑一格,算出把它翻过来会让能量变化多少,记作 ΔE。
如果 ΔE ≤ 0(翻了更省),就翻。如果 ΔE > 0(翻了亏本),也仍然以概率 e−ΔE/kT 翻(k 是玻尔兹曼常数,一个把温度换算成能量的固定系数)——温度越高,越容易做亏本的事。
重复第 4、5 步,一直跑到统计性质不再变化为止。
规则只有这些。没有电子、没有量子力学、没有铁和镍的区别。整个模型里只剩两个数:J(邻居之间有多黏)和 T(随机扰动有多强),而真正决定一切的其实只有它们的比值。
第 5 条是全部戏剧性的来源。如果只允许"往省能量的方向走",系统会立刻冻死在某个格局里。允许一定概率地亏本行动,系统才会不停试探,两种倾向——"跟邻居一致"和"随便乱翻"——才有得较劲。温度就是这场较劲的天平。
规则只管一件事:和邻居一样,省能量
四个邻居都同向
能量最低 → 稳
四个邻居都反向
能量最高 → 想翻
能量差 8J
温度决定
多大概率
容忍「亏本翻」
翻一格要付的代价,只取决于它的四个邻居里有几个跟它一致。整个模型的全部物理就在这一张图里。
03 跑起来会看到什么(What You See When It Runs)
把 T 从低往高扫一遍,你会看到三种截然不同的画面。
同一套规则,只改一个温度
低温 kT/J = 1.6
几乎全体一致
临界 kT/J = 2.269
每种尺寸的团块都有
高温 kT/J = 3.4
团块只有几格大
三张都是 64×64 的格子,规则一模一样,只有温度不同。中间那张的团块从几格大到横跨全场都有——没有"典型尺寸"。
低温 :几乎全场一个颜色,只有零星的反向小点。系统选定了一个方向,而且守得住。
高温 :一片杂点,团块最多几格大。任何一致性都撑不过几个格子的距离。
中间那个特定的温度 :团块什么尺寸都有。你把图放大或缩小,看到的东西统计上长得一样。这就是关联长度(correlation length)发散的样子——相隔多远的两格都还在互相牵扯。
把"整体朝向"(所有格子的取值取平均,再取绝对值)对温度画出来,就是下面这条曲线。二维方格上,它有精确解。
昂萨格 1944 年算出来的精确曲线
温度 kT/J →
整体朝向 |m|
1
2
3
4
1
0
Tc = 2/ln(1+√2) ≈ 2.269
低温:几乎全体一致
高温:整体朝向恒为 0
贴近 Tc 处 |m| ∝ (Tc−T)^(1/8)
曲线在 Tc 处不跳,但坡度是无穷大 — 这正是「连续相变」的样子
Onsager 1944 年的精确解。临界温度是 kTc/J = 2/ln(1+√2) ≈ 2.269,曲线在这一点连续但斜率无穷大。
贴近 Tc 的那一段,整体朝向按 (Tc − T) 的 1/8 次方趋于 0。这个 1/8 是精确值,不是拟合出来的。三维的情况没有精确解,但数值方法与共形自举给出的指数是 0.3264 左右——正是真实磁体和液气临界点测到的那个数。
04 它解释了现实中的什么(What It Explains)
铁磁性的居里点。 这是它的老本行:为什么磁性不是逐渐变弱,而是在一个确定的温度上归零。
液气临界点。 把"这一格有分子/没分子"当成 +1/−1,伊辛模型就变成了"格子气"(lattice gas),这个对应是 Lee 与 Yang 在 1952 年写清楚的。于是水的临界点和磁体的居里点在数学上是同一件事。
合金的有序-无序转变。 黄铜里的铜原子和锌原子,低温下规规矩矩地间隔排列,高温下随机占位——把"这个格点上是铜还是锌"当成 ±1,又是同一个模型。
两种液体的分层。 油和水在某个温度以上互溶、以下分层,同样落在这个类里。
这四件事的共同点,不是它们像,而是它们的临界指数在实验精度内相同 。这就是"三维伊辛普适类"这个名字的全部含义。也正因为如此,它是复杂性科学里少数几个"跨领域的相同"能被精确到小数点后三位的例子——大部分跨领域类比做不到这一点。
它不能解释什么(What It Cannot Explain)
不能解释真实磁性的来源。 为什么电子自旋倾向于对齐,是量子力学(交换相互作用)的事,伊辛模型直接把这件事塞进常数 J 里当输入。而且真实的自旋可以指向三维空间的任意方向,不是只有上下两种;那种情况属于另一个普适类(海森堡类),指数不一样。
不能给出 Tc 是多少。 模型能算出"kTc/J 等于多少",但 J 是多少完全由材料决定,模型不管。所以它永远说不出铁的居里点是 770 度。它给的是形状,不是位置。
它本身没有时间。 模型定义的是平衡态的统计分布,不是一个动力学过程。第 4、5 步那套翻转规则(Metropolis 算法)是人为编出来的采样手段,不是真实物质的运动方式。所以拿它去谈"转变要花多久""怎么演化过去的",是在问模型没定义的东西。
一维版本什么都不会发生。 Ising 1924 年算出的正是这个:一条线上,任何有限温度都不会有相变。他由此断定任何维度都没有——这个结论是错的,但错得很有名。(也因此,"我的系统像一维链"是个真实的、会取消结论的边界条件。)
不能直接套到人身上。 把 ±1 读成"赞成/反对"、把 T 读成"社会噪声",念起来很顺,但推导用到的前提在人的系统里基本不成立:需要平衡态(社会不在)、单元同质(人不是)、耦合对称(我影响你多少你就影响我多少,社会里不是)、规模大到有限尺度效应可忽略(一个组织几十上百人,差得远)。更要紧的是社会网络里到处是长程连接,而射程会直接换掉普适类。作为一个"局部一致如何长出全局一致"的比喻 ,它很好用;作为一个能算数的模型,用之前得先把这四条逐条检查过。