磁铁吸引金属时发生了什么？（四）：10²³ 个自旋如何达成一致，以及通往希格斯的桥

系列大结局。第（一）篇证明了经典物理禁止磁铁；第（二）篇从狄拉克方程里拽出了电子的内禀磁矩（ $$g=2$$ ）；第（三）篇把自旋一路追溯到时空本身的对称性。于是我们已经讲清了一个电子的磁性。可一块冰箱贴里有大约 $10^{23}$ 个电子，剩下的那个问题才最古怪：它们凭什么会达成一致？是什么把它们排齐的——而且，我们马上就会看到，那个负责排齐的力，根本就不是磁力。

动手之前，我先把我们正要踩进去的那个陷阱说在前头。忙活了三篇，我们好不容易挣到一个行为像小磁针的电子。这看上去像是胜利了。其实不是。把 $10^{23}$ 根这样的小磁针抓来，塞进一根铁棒，默认的结局是：热骚动把它们指得四面八方，铁棒死气沉沉——没有净磁场，吸不住你的冰箱门。要得到一块磁铁，必须有什么东西伸进这堆自旋里，强迫它们达成一致。第（四）篇讲的就是：到底是谁在强迫。而那个谜底，我现在就剧透——因为它实在太漂亮，藏着可惜——是泡利不相容原理，披着一身静电学的外衣。

8. Heisenberg Model

我们现在知道每个电子都是一个小磁针（ $$g=2$$ ）。但是，如果你把一堆小磁针放在一起，室温下的热骚动（Thermal Agitation）足以将它们的方向完全打乱，使得宏观磁矩为零（顺磁性）。要形成铁磁性，自旋之间必须存在一种极强的“耦合力”，强制它们整齐排列。经典磁偶极相互作用（Dipole-Dipole Interaction）太弱了，大约只有热能的万分之一。所以那个最显眼的嫌疑人当场就能排除：小磁针之间并不是靠磁力互相拉扯排齐的——那个力输给室温，差了整整四个数量级。真正的力量来自于我们之前提到的全同粒子统计性（Spin-Statistics）与库仑相互作用的结合。这被称为交换相互作用（Exchange Interaction）。

为了展示本质，我们考虑最简单的模型：双电子系统（例如 He 原子或者是两个相邻的 Fe 原子上的电子）。假设有两个电子 1 和 2，以及两个空间轨道 $\psi_a(\vec{r})$ 和 $\psi_b(\vec{r})$ ， $\psi_a(\vec{r})$ 局域在原子 A 附近， $\psi_b(\vec{r})$ 局域在原子 B 附近，并假设这两个轨道是正交归一的： $\langle \psi_a | \psi_b \rangle = 0$ ，注意这里的正交是我们接下来将要得到的 Heisenberg model 的前提假设，即使不正交归一我们也可以通过基底变换造出两个正交归一的新轨道，认为二者的重叠积分很小，这会使得双电子系统的计算大大简化，但是这并不意味着两个电子之间就没有交互了，因为交互作用中会伴随着交互作用项积分，一般不为零。

这个系统的总哈密顿量为： $H = H_0 + H_{int}$ ，其中 $$H_0$$ 是单电子部分（动能+核势能）， $H_{int}$ 是两个电子之间的库仑相互作用： $H_{int} = \frac{e^2}{|\vec{r}_1 - \vec{r}_2|}$ ，注意：这里完全没有磁相互作用项，只有纯粹的静电斥力。请盯住这个“没有”——这正是全部要害。最后无论是什么东西把自旋排齐了，它都得从一个压根没出现自旋的哈密顿量里冒出来。

根据量子力学全同粒子费米统计假设，电子的总波函数 $\Psi(1,2)$ 必须在粒子交换算符 $P_{12}$ 作用下变号（反对称）：

P_{12} \Psi(1,2) = -\Psi(1,2)

因为总波函数由空间部分 $\phi(\vec{r}_1, \vec{r}_2)$ 和自旋部分 $\chi(s_1, s_2)$ 组成： $\Psi = \phi \otimes \chi$ 。我们首先推导自旋部分，然后根据总体的反对称性得到空间部分。我们有两个自旋为 $$1/2$$ 的粒子（比如两个电子），总共有 $2 \times 2 = 4$ 种可能的直积状态（未耦合基底），我们要把它们的自旋加在一起，看看总自旋 $\vec{S}_{tot} = \vec{S}_1 + \vec{S}_2$ 会是什么样子，找到总自旋算符 $\hat{S}^2$ 和总磁量子数 $\hat{S}_z$ 的共同本征态 $|S, M\rangle$ 。根据角动量加法规则，两个 $$1/2$$ 自旋合成的总自旋 $$S$$ 可以取： $$S = 1/2 + 1/2 = 1$$ （三重态，有 3 个 $$M$$ 值： $$+1, 0, -1$$ ）和 $$S = 1/2 - 1/2 = 0$$ （单态，有 1 个 $$M$$ 值： $$0$$ ）。

三重态对应三个分量。总磁量子数 $$M$$ 是两个粒子磁量子数之和： $$M = m_1 + m_2$$ 。要让 $$M=1$$ ，唯一的可能是两个电子都向上： $$1/2 + 1/2 = 1$$ 。所以，三重态的第一个成员是确定的： $|1, 1\rangle = |\uparrow\uparrow\rangle$ ，然后使用降算符 ( $$S_-$$ ) 得到中间态 ( $$M=0$$ ) 从 $|1, 1\rangle$ 推导出 $|1, 0\rangle$ 。利用总降算符 $\hat{S}_- = \hat{S}_{1-} + \hat{S}_{2-}$ 作用在态 $|j, m\rangle$ 上： $J_- |j, m\rangle = \hbar \sqrt{j(j+1) - m(m-1)} |j, m-1\rangle$ ，对左边的耦合态作用得到：

\hat{S}_- |1, 1\rangle = \sqrt{1(1+1) - 1(1-1)} |1, 0\rangle = \sqrt{2} |1, 0\rangle

对右边的直积态作用得到：

\begin{aligned} (\hat{S}_{1-} + \hat{S}_{2-}) |\uparrow\uparrow\rangle &= (\hat{S}_{1-} |\uparrow\rangle_1)|\uparrow\rangle_2 + |\uparrow\rangle_1 (\hat{S}_{2-} |\uparrow\rangle_2) \ = |\downarrow\rangle_1|\uparrow\rangle_2 + |\uparrow\rangle_1|\downarrow\rangle_2 \ = |\downarrow\uparrow\rangle + |\uparrow\downarrow\rangle \end{aligned}

于是得到：

|1, 0\rangle = \frac{1}{\sqrt{2}} (|\uparrow\downarrow\rangle + |\downarrow\uparrow\rangle)

最低权态 ( $$M=-1$$ ) 也很简单，只有两个都向下才能得到 $$-1$$ ： $|1, -1\rangle = |\downarrow\downarrow\rangle$ 。于是三重态自旋部分是对称的，交换都不变号，那么空间部分必须是反对称的。

|0, 0\rangle = \frac{1}{\sqrt{2}} (|\uparrow\downarrow\rangle - |\downarrow\uparrow\rangle)

单态（Singlet）的自旋部分是反对称的，对应的空间部分必须是对称的。总结为：

空间对称 $\otimes$ 自旋反对称（Singlet, 单态）能量：$E_S$
- 自旋部分（反对称）： $$S=0$$ 的单态 $\chi_S = \frac{1}{\sqrt{2}}(\uparrow\downarrow - \downarrow\uparrow)$ 。
- 空间部分（对称）： $\phi_S(\vec{r}_1, \vec{r}_2) = \frac{1}{\sqrt{2}} [\psi_a(\vec{r}_1)\psi_b(\vec{r}_2) + \psi_a(\vec{r}_2)\psi_b(\vec{r}_1)]$ 。
空间反对称 $\otimes$ 自旋对称（Triplet, 三重态）能量：$E_T$
- 自旋部分（对称）： $$S=1$$ 的三重态 $\chi_T = \{\uparrow\uparrow, \frac{1}{\sqrt{2}}(\uparrow\downarrow + \downarrow\uparrow), \downarrow\downarrow\}$ 。
- 空间部分（反对称）： $\phi_A(\vec{r}_1, \vec{r}_2) = \frac{1}{\sqrt{2}} [\psi_a(\vec{r}_1)\psi_b(\vec{r}_2) - \psi_a(\vec{r}_2)\psi_b(\vec{r}_1)]$ 。

先停下来，看清这套记账到底悄悄办成了什么——因为它就是一切磁性背后的秘密引擎。两个电子想自旋平行还是反平行，唯一逼它们做选择的，是总波函数的反对称性。如果自旋平行（三重态），空间部分就必须反对称，于是 $\phi_A$ 在 $\vec{r}_1 = \vec{r}_2$ 时为零——两个电子被禁止叠在一起，只好彼此保持距离，少付了库仑斥力；如果自旋反平行（单态），空间部分对称，两个电子乐得挤在一处，于是多攒了库仑能。哈密顿量里从没提过的自旋取向，就这样借着自旋被迫采取的几何位形，偷偷溜进了能量里。整个戏法就在这儿。下面我们把它写成公式。

我们的目标是：寻找一个只包含自旋算符 $\vec{S}_i$ 和 $\vec{S}_j$ 的数学表达式 $\hat{H}_{eff}$ ，使得它作用在单态和三重态上时，能够自动给出对应的能量 $$E_S$$ 和 $$E_T$$ 。要构建这个哈密顿量，最自然的积木块是两个自旋的点积 $\vec{S}_i \cdot \vec{S}_j$ 。我们需要先算出这个算符在单态和三重态下的特征值分别是多少。定义双电子系统的总自旋算符： $\vec{S}_{tot} = \vec{S}_i + \vec{S}_j$ ，对总自旋算符取平方反解出点积项： $\vec{S}_i \cdot \vec{S}_j = \frac{1}{2} \left( \vec{S}_{tot}^2 - \vec{S}_i^2 - \vec{S}_j^2 \right)$ 。利用量子力学中角动量平方算符的特征值公式 $\hat{S}^2 |s\rangle = s(s+1) |s\rangle$ （此处为了简洁暂时省略 $\hbar^2$ ，最后结果不受影响，或者视为无量纲自旋）：

对于单个电子 ( $$s=1/2$$ )： $\vec{S}_i^2 = \vec{S}_j^2 = \frac{1}{2}\left(\frac{1}{2} + 1\right) = \frac{3}{4}$
对于单态 ( $S_{tot}=0$ )： $\vec{S}_{tot}^2 |S\rangle = 0(0+1) |S\rangle = 0$
对于三重态 ( $S_{tot}=1$ )： $\vec{S}_{tot}^2 |T\rangle = 1(1+1) |T\rangle = 2 |T\rangle$ 现在，我们把这些值代回点积公式，计算特征值：
单态下的点积值： $(\vec{S}_i \cdot \vec{S}_j) |S\rangle = \frac{1}{2} \left( 0 - \frac{3}{4} - \frac{3}{4} \right) |S\rangle = -\frac{3}{4} |S\rangle$
三重态下的点积值： $(\vec{S}_i \cdot \vec{S}_j) |T\rangle = \frac{1}{2} \left( 2 - \frac{3}{4} - \frac{3}{4} \right) |T\rangle = \frac{1}{2} \left( \frac{1}{2} \right) |T\rangle = \frac{1}{4} |T\rangle$

我们假设等效哈密顿量 $\hat{H}_{ij}$ 具有如下线性形式（这是最普遍的旋转对称形式）： $\hat{H}_{ij} = C_0 + C_1 (\vec{S}_i \cdot \vec{S}_j)$ ， $$C_0$$ 只是一个与自旋构型无关的常数能量平移（Energy Shift），在研究相变和自旋动力学时，我们可以将其丢弃（或者重新定义能量零点），然后将单态和三重态的数据代入，有 $E_S=C_1(-\frac{3}{4}),E_T=C_1(\frac{1}{4})$ ，定义常数 $J \equiv E_S - E_T$ ，就得到了双粒子哈密顿量：

\hat{H}_{ij} = -J (\vec{S}_i \cdot \vec{S}_j)

如果 $$J>0$$ ( $$E_S > E_T$$ )：系数是负的，点积越大（平行，+1/4），能量越低。这是铁磁性。
如果 $$J<0$$ ( $$E_S < E_T$$ )：系数是正的，点积越小（反平行，-3/4），能量越低。这是反铁磁性。

换句话说：一个纯粹静电的相互作用，自己重组成了一个只谈自旋的项。整个计算里从没进来过磁力。 $$J$$ 的能量量级由库仑斥力定下——是电子伏特，而不是磁偶极耦合那种微电子伏特——这恰恰是它能压过室温、把 $10^{23}$ 个自旋摁在队形里的原因。这就回答了开篇那个问题：自旋之所以对齐，是因为泡利原理把它们的取向和它们要付的静电能绑在了一起，而静电，很猛。

现在我们将这个双粒子相互作用推广到整个晶格。假设每个电子 $$i$$ 只和它的最近邻（Neighbors）发生相互作用。我们需要对晶格中所有的原子求和。为了修正这种重复计数：

H_{exchange} = -\frac{J}{2} \sum_{i,j \text{ neighbor}} \vec{S}_i \cdot \vec{S}_j

最后，我们还要考虑每个电子自旋与外部均匀磁场 $\vec{B}$ 的相互作用。这是单体相互作用，不需要考虑邻居。回顾我们之前从 Dirac 方程导出的结论，电子具有自旋磁矩：

\vec{\mu}_S = -g \frac{e}{2m} \vec{S}

为了让公式更加简洁和通用，物理学家定义了一个自然常数组合，称为玻尔磁子（Bohr Magneton）。玻尔磁子是原子物理中磁矩的自然单位。它的定义为：

\mu_B \equiv \frac{e\hbar}{2m_e}

这个物理量包含了三个基本常数：基本电荷 $$e$$ 、普朗克常数 $\hbar$ 和电子质量 $$m_e$$ 。它代表了一个经典电子在氢原子基态轨道上运动所产生的轨道磁矩大小，这是国际标准公制，而在高斯单位制下，还包含光速 $\mu_B \equiv \frac{e\hbar}{2m_ec}$ 。如果我们将自旋算符 $\vec{S}$ 视为无量纲算符（即特征值为 $$1/2$$ 而不是 $\hbar/2$ ），那么真实的物理角动量是 $\hbar \vec{S}$ 。将这个 $\hbar$ 提取出来与前面的常数结合：

\begin{aligned} \vec{\mu}_S &= -g \frac{e}{2m_e} (\hbar \vec{S}_{\text{dimensionless}}) = -g \left( \frac{e\hbar}{2m_e} \right) \vec{S}= -g \mu_B \vec{S} \end{aligned}

当电子处于外部磁场 $\vec{B}$ 中时，其势能（Zeeman Energy）由经典电磁学公式 $U = -\vec{\mu} \cdot \vec{B}$ 给出。将上述磁矩表达式代入：

\begin{aligned} U_{Zeeman} &= - (-\vec{\mu}_S) \cdot \vec{B} \quad (\text{注意符号：通常势能定义为 } -\vec{\mu}\cdot\vec{B}) \\ &= - (-g \mu_B \vec{S}) \cdot \vec{B} \\ &= g \mu_B \vec{S} \cdot \vec{B} \end{aligned}

关于符号的物理约定： 在凝聚态物理中，我们通常希望哈密顿量反映能量极小值。电子带负电，磁矩 $\vec{\mu}$ 与自旋 $\vec{S}$ 反向。能量最低态是磁矩 $\vec{\mu}$ 与磁场 $\vec{B}$ 平行。这意味着自旋 $\vec{S}$ 与磁场 $\vec{B}$ 反平行。为了避免处理繁琐的负号，或者为了让自旋看起来“顺着”磁场（定义 $\vec{S}$ 指向磁矩方向而非角动量方向），有时文献会调整定义。但最标准的推导（保留电子电荷负号）给出的塞曼项通常写为 $+ g\mu_B \vec{S} \cdot \vec{B}$ 或 $- \vec{\mu} \cdot \vec{B}$ 。然而，在海森堡模型的习惯写法中，为了数学上的对称性及方便讨论（例如假设 $$g$$ 为负或重新定义自旋方向），通常将塞曼项写为负号形式，表示自旋倾向于沿着场排列（这是一种唯象的处理）：

H_{Zeeman} = -g \mu_B \sum_i \vec{S}_i \cdot \vec{B}

（这意味着能量最低时， $\vec{S}_i$ 与 $\vec{B}$ 同向。这蕴含了我们重新定义了自旋方向，或者 $$g$$ 取负值。在唯象模型中，我们只关心：磁场倾向于把自旋拉向哪个方向。）

现在，我们将两部分拼图合二为一：内部相互作用（由泡利原理和库仑力产生的交换能）与外部相互作用（由相对论量子效应产生的磁矩与外场耦合），将它们相加，我们终于得到了描述固体磁性最核心的哈密顿量——海森堡模型（Heisenberg Model）：

\boxed{H = -\frac{J}{2} \sum_{\langle i,j \rangle} \vec{S}_i \cdot \vec{S}_j - g \mu_B \sum_i \vec{S}_i \cdot \vec{B}}

这个公式是现代磁学的基石。第一项（ $$J$$ ）解释了为什么磁铁会有磁性（自发磁化，自旋整齐排列）。第二项（ $$B$$ ）解释了磁铁如何被外界控制（磁化过程，磁滞回线）。 $\mu_B$ 和 $$g$$ 将微观的量子常数（ $\hbar, e, m_e$ ）与宏观的可观测磁场联系在了一起。请留意这两项是怎么悄悄把整个系列收口的：第一项是第（一）篇全同粒子统计性的遗产，第二项是第（二）篇 $$g=2$$ 的遗产。那块冰箱贴，就坐在这两者的交汇点上。

9. Ising Model

我们已经完成了微观机制的构建（Dirac $\to$ Spin $\to$ Exchange Interaction $\to$ Heisenberg Hamiltonian）。现在，我们要从微观走向宏观。我们手里有了一个想让自旋对齐的哈密顿量——可”想”和”是一块磁铁”不是一回事。一根烧到十亿度的铁棒，交换耦合一点没少，它却不是磁铁。在”自旋偏好对齐”和”自旋真的对齐了”之间，横着一个温度；而这两种状态之间的过渡，来得突然、来得锋利，是物理里最美的东西之一。为了看见它，我们需要处理 Heisenberg 模型。但在二维或三维中精确求解 Heisenberg 模型极度困难（因为它包含非对易算符）。因此，我们需要引入 Ising 模型作为近似，并使用平均场论（Mean-Field Theory）来展示对称性是如何破缺的。

我们对 Heisenberg model 进行无量纲化处理，并分析其结构。假设外磁场沿着 $$z$$ 方向 $\vec{B} = (0, 0, B)$ 。我们将自旋算符 $\vec{S}$ 的点积展开为纵向（ $$z$$ ）和横向（ $$x, y$$ ）分量：

\vec{S}_i \cdot \vec{S}_j = S_i^z S_j^z + (S_i^x S_j^x + S_i^y S_j^y)

为了看清横向部分的物理意义，我们引入升降算符（Ladder Operators）：

S_i^+ = S_i^x + i S_i^y, \quad S_i^- = S_i^x - i S_i^y

于是，海森堡模型可以重写为两部分：

H = \underbrace{\left[ -\frac{J}{2} \sum_{\langle i,j \rangle} S_i^z S_j^z - h \sum_i S_i^z \right]}_{\text{Ising Part}} + \underbrace{\left[ -\frac{J}{4} \sum_{\langle i,j \rangle} (S_i^+ S_j^- + S_i^- S_j^+) \right]}_{\text{Flip Part}}

（其中 $h = g\mu_B B$ ）。这两部分有着截然不同的物理含义：

Ising Part (纵向项)：这一项只涉及 $$S^z$$ 。由于不同格点的 $$S_i^z$$ 是对易的（ $$[S_i^z, S_j^z]=0$$ ），它们表现得像经典的标量变量。这描述了自旋在 $$z$$ 轴上的静态排列。
Flip Part (横向项/翻转项)：这一项涉及 $$S_i^+ S_j^-$$ 。它的作用是：将第 $$j$$ 个格点的自旋向下翻（ $$S^-$$ ），同时将第 $$i$$ 个格点的自旋向上翻（ $$S^+$$ ）。 $S_i^+ S_j^- |\downarrow_i \uparrow_j\rangle = |\uparrow_i \downarrow_j\rangle$ 。这物理上代表了自旋激发的移动。就像一个自旋反转的状态在晶格中跳跃一样，这对应于自旋波（Spin Wave）或磁振子（Magnon）。这是一种量子涨落，它赋予了系统“动能”，倾向于破坏有序排列。

在许多真实的磁性材料中，由于晶格结构的对称性，存在磁各向异性（Magnetic Anisotropy）。这意味着自旋在某些方向（如 $$z$$ 轴，易轴）上的能量比在 $$x,y$$ 平面内更低。如果各向异性足够强，或者我们只关心经典极限下的相变行为，我们可以忽略 Flip Part（量子涨落项），只保留纵向项。这便是著名的 Ising Model。

此时，我们将算符 $$S_i^z$$ 替换为经典变量 $\sigma_i = \pm 1$ （吸收系数到 $$J$$ 中）：

H_{Ising} = - \frac{J}{2} \sum_{\langle i,j \rangle} \sigma_i \sigma_j - h \sum_i \sigma_i

这是一个巨大的简化：我们将一个非对易的量子矩阵问题，变成了一个经典的统计排列组合问题。2D Heisenberg model 至今是没有办法精确解的，我们可以计算方法求解。不过其中的 Ising model 相对简单一些，1D Ising model 的解很简单，在 1D 没有相变（Ising, 1924），但是 2D 也是很难解，直到有一个奇人 Lars Onsager 发了个很简单的论文说这个他解了，但是没有给出怎么解的，但是他给出了那个 critical temperature，且指出低温时铁磁，高温是顺磁，这是第一个可以精确解出的铁磁相变的 model，他一直没有公布解答过程，直到杨振宁看到这篇论文，给出了解答，不过一样很困难。为了直观理解物理图像，我们这里采用平均场近似。

我们面临的问题是多体耦合： $\sigma_i$ 的状态取决于 $\sigma_j$ ，而 $\sigma_j$ 又取决于 $\sigma_k$ … 这种连锁反应使得配分函数难以计算。平均场的思路偷懒得可爱，而它的力量恰恰在这份偷懒里：当我们关注第 $$i$$ 个原子时，我们不关心邻居 $$j$$ 到底是在 $$+1$$ 还是 $$-1$$ 之间跳变，我们只关心邻居的平均影响。我们将 $\sigma_j$ 写为平均值加涨落： $\sigma_j = \langle \sigma \rangle + \delta \sigma_j$ 。忽略二阶涨落项 $\delta \sigma_i \delta \sigma_j \approx 0$ ，哈密顿量可以线性化为单体形式：

H_{MFA} = - \sum_i \sigma_i \underbrace{\left( J \sum_{j \in \text{neigh}} \langle \sigma \rangle + h \right)}_{h_{eff}}

这就定义了有效分子场（Molecular Field） $h_{eff}$ ：

h_{eff} = J z \langle \sigma \rangle + h

其中 $$z$$ 是配位数（每个原子的邻居数量）。现在，问题变成了一个自旋在“外场” $h_{eff}$ 中的统计分布。根据玻尔兹曼分布，该自旋向上的概率正比于 $e^{\beta h_{eff}}$ ，向下的概率正比于 $e^{-\beta h_{eff}}$ （ $\beta = 1/k_B T$ ）。于是，该自旋的热力学平均值 $\langle \sigma_i \rangle$ 为：

\langle \sigma_i \rangle = \frac{(+1)e^{\beta h_{eff}} + (-1)e^{-\beta h_{eff}}}{e^{\beta h_{eff}} + e^{-\beta h_{eff}}} = \tanh(\beta h_{eff})

精彩的自我指涉就在这里。每个自旋都在响应邻居布下的平均场——可每个自旋本身就是别人的邻居。因为晶格是均匀的， $\langle \sigma_i \rangle$ 必须等于场源本身的平均值 $m = \langle \sigma \rangle$ 。一个自旋，得跟它自己帮忙造出来的那个平均值达成一致。代入 $h_{eff}$ ，我们得到了著名的自洽方程（Self-consistent Equation）：

m = \tanh\left( \frac{J z m + h}{k_B T} \right)

让我们考虑最关键的情况：没有外部磁场 ( $$h=0$$ )。方程简化为：

m = \tanh\left( \frac{T_c}{T} m \right)

其中我们将常数打包定义为 居里温度（Curie Temperature） $$T_c = Jz/k_B$$ 。这是一个超越方程，我们可以通过作图法（寻找 $$y=m$$ 和 $y=\tanh(\frac{T_c}{T}m)$ 的交点）来分析解的行为。一切都系于一个比较—— $\tanh$ 曲线在原点的斜率 $$T_c/T$$ ，对阵直线 $$y=m$$ 的斜率，也就是 $$1$$ ：

高温相 ( $$T > T_c$$ )： $\tanh$ 曲线在原点的斜率 $$T_c/T < 1$$ 。直线与曲线只有一个交点 $$m=0$$ 。 物理意义： 热骚动剧烈，没有外场就没有磁性。这是顺磁相（Paramagnetic）。
低温相 ( $$T < T_c$$ )： $\tanh$ 曲线在原点的斜率 $$T_c/T > 1$$ 。原点 $$m=0$$ 变成不稳定解，出现了两个新的稳定非零解 $m = \pm m_0$ 。 物理意义： 即使 $$h=0$$ ，系统也会自发地产生非零的磁化强度 $$m_0$$ 。这是铁磁相（Ferromagnetic）。

这就是自发对称性破缺（Spontaneous Symmetry Breaking）： 哈密顿量在 $$h=0$$ 时具有 $\sigma \to -\sigma$ 的翻转对称性。但是，当温度低于 $$T_c$$ 时，大自然被迫在“全体向上”和“全体向下”中选择一个。这种选择不是由外力强加的，而是系统为了降低能量（由交换相互作用 $$J$$ 驱动）而自发做出的集体决策。这正是磁铁吸铁这一宏观现象背后的统计力学本质。请在这个“被迫”上多停一会儿——第（一）篇里我们见过它，相对论逼着狄拉克接受矩阵；马上我们还要在宇宙尺度上和它做最后一次照面。同一套逻辑反复回响：一个在高能下成立的对称性，不是被主动放弃的，而是因为继续保持对称变成了更贵的那个选项，才破缺掉。

10. Magnetic Domains

我们根据 Ising 模型和平均场论得出了结论：当温度低于居里温度 $$T_c$$ 时，电子自旋会自发排列，产生巨大的宏观磁化强度 $$M$$ 。可这下理论给我们出了个难堪。如果你去五金店买一颗铁钉（室温显然远低于铁的居里温度 $$1043K$$ ），它并没有磁性。它不会去吸别的东西。按我们刚证完的一切，它本该是块磁铁——它却懒洋洋地躺在料斗里。这是因为我们忽略了最后一项能量竞争。

我们之前的哈密顿量只考虑了交换能（Exchange Energy）和塞曼能（Zeeman Energy）。但在宏观尺度上，还有一个经典的静磁能（Magnetostatic Energy）。如果一块铁中的所有 $10^{23}$ 个原子都整齐朝上排列，这块磁铁会在周围空间建立一个巨大的磁场。磁场包含了能量密度 $B^2/2\mu_0$ 。把所有磁力线撑开需要耗费巨大的能量。为了降低这个静磁能，材料会自发地分裂成许多微小的区域，称为磁畴（Magnetic Domains）。虽然在磁畴内部，交换相互作用让自旋整齐排列（满足微观铁磁性）；但在整体上，各个磁畴的磁矩矢量和为零（ $\sum \vec{M}_i = 0$ ），外部没有磁感线，从而极大地降低了静磁能。所以那颗铁钉不是没被磁化，它是被完美地磁化了——只不过分成了上百万个互相抵消的小块。

磁畴之间的边界被称为畴壁（Domain Wall）。在畴壁内，自旋不是突然反转的，而是渐进旋转的。这又是一次能量的博弈：交换能希望自旋平行，不想让它们转弯（希望畴壁越宽越好）；磁各向异性（Anisotropy）希望自旋沿着易轴，不想让它们指向中间方向（希望畴壁越窄越好）。这两者的平衡决定了畴壁的厚度（通常是几百个原子层）。磁畴的形成并不是像自旋那样来自“第一性原理”的推导，而是来自连续介质场论下的能量最小化（Energy Minimization）属于微磁学（Micromagnetics）的范畴，这里就不展开了。

现在，我们终于可以完整描述“磁铁吸铁”的宏观过程了——也就是开启这整个系列的那一声“咔哒”，如今从一个电子的时空对称性一路拼到了这里：

初始状态：铁钉内部充满了杂乱无章的磁畴，宏观磁矩为零。
外场介入：当你拿磁铁靠近铁钉时，提供了一个外磁场 $\vec{B}_{ext}$ 。
畴壁移动（Domain Wall Motion）：这就打破了平衡。那些方向与 $\vec{B}_{ext}$ 一致的磁畴拥有更低的塞曼能（ $E = -\vec{M} \cdot \vec{B}$ ）。于是，这些“顺从”的磁畴开始吞并周围“不顺从”的磁畴。畴壁发生了移动。
宏观磁化：铁钉迅速获得了一个巨大的净磁矩。
梯度力做功：这个感应出来的宏观磁矩 $\vec{m}_{total}$ 受到磁铁产生的非均匀磁场的梯度力 $\vec{F} = \nabla(\vec{m}_{total} \cdot \vec{B})$ 的牵引。
咔哒：铁钉飞向了磁铁。

Conclusion: The Deep Symmetry of the Universe

当你手中把玩两块磁铁，感受它们之间的排斥与吸引时，你感受到的不仅仅是力。你是在用手触摸量子力学的本质和宇宙演化的秘密。让我们回顾这段旅程，看看我们是如何重建物理直觉的：

经典的崩塌：我们发现洛伦兹力不做功，且经典统计物理禁止磁性（Bohr-van Leeuwen 定理）。
相对论的修正：Dirac 方程揭示了电子必须是 4 分量的旋量，并携带 $$g=2$$ 的内秉磁矩。磁性，是相对论效应在低速世界的残留。
量子统计的强力：泡利不相容原理结合库仑斥力，产生了一种等效的“交换相互作用”，迫使自旋平行排列。
对称性的破缺：Ising 模型告诉我们，当温度降低时，系统为了生存（降低能量），不得不破坏旋转对称性，选择一个方向。

最后，值得一提的是，我们在 Ising 模型中看到的自发对称性破缺（SSB），其意义远超固体物理。它是现代物理学理解宇宙的核心范式。在宇宙大爆炸之初（极高温），物理规律具有极高的对称性。所有的基本粒子都没有质量，就像高温下的铁块没有磁性一样（顺磁相）。随着宇宙冷却，当温度跌破某个临界值时，充满宇宙的**希格斯场（Higgs Field）**发生了相变。就像电子自旋突然选择朝向一个方向一样，希格斯场在真空中获得了一个非零的期望值（Vacuum Expectation Value）。我们在原点画下的那张自洽方程图，换一组坐标轴的标签，就是早期宇宙冷却时求解的同一个方程。

在铁磁体中，对称性破缺赋予了材料磁性。
在标准模型中，对称性破缺赋予了基本粒子质量。

所以，下一次当你看到磁铁吸起一枚回形针时，请意识到：你正在目睹一个微缩版的宇宙创生时刻。那个让铁钉获得磁性的机制，正是让你身体里的夸克和电子获得质量、让这个宇宙得以存在的同一机制。

磁力不做功，是时空的几何在做功；磁铁的吸力，是量子的幽灵在宏观世界的舞蹈。