4. 哈伯德模型与强关联物理

在凝聚态物理中，许多最引人入胜的现象——莫特绝缘体、高温超导、量子磁性——都无法用单电子近似解释。它们的共同根源是电子之间的强库仑相互作用。哈伯德模型（Hubbard Model）正是描述这类强关联电子系统最简洁、最核心的格点模型。本节从二次量子化出发推导哈伯德哈密顿量，分析其相互作用项，讨论莫特绝缘体与磁性序的物理图像，并给出相图与平均场近似的基本方法。

一、从库仑相互作用到哈伯德哈密顿量

考虑晶体中电子在周期性势场中运动，完整的多体哈密顿量包含动能项与电子-电子库仑相互作用项。在二次量子化表象下，以格点Wannier轨道为基，哈密顿量写作：

H = Σ_{ij,σ} t_{ij} c†_{iσ} c_{jσ} + ½ Σ_{ijkl,σσ'} V_{ijkl} c†_{iσ} c†_{jσ'} c_{lσ'} c_{kσ}

其中 c†_{iσ} 在格点 i 上产生一个自旋为 σ 的电子，t_{ij} 是跃迁积分。哈伯德模型的核心近似是：只保留同一格点上的库仑排斥，即令 V_{iiii} = U，其余全部忽略。同时只保留最近邻跃迁 t。于是得到单带哈伯德哈密顿量：

H = -t Σ_{⟨ij⟩,σ} (c†_{iσ} c_{jσ} + h.c.) + U Σ_i n_{i↑} n_{i↓} - μ Σ_{iσ} n_{iσ}

这里 ⟨ij⟩ 表示最近邻对，n_{iσ} = c†_{iσ} c_{iσ} 是格点 i 上自旋 σ 的粒子数算符，μ 为化学势。第一项描述电子在格点间的跳跃，倾向于使电子离域形成能带；第二项 U n_{i↑} n_{i↓} 是同一格点上两个自旋相反电子之间的排斥能，倾向于使电子局域。t 与 U 的竞争，正是强关联物理的全部戏剧所在。

二、相互作用项的物理含义

相互作用项 U Σ_i n_{i↑} n_{i↓} 看似简单，却蕴含深刻的物理。当 U = 0 时，模型退化为紧束缚能带，电子填充到半满时系统是金属。当 U 远大于 t 时，双占据一个格点需要付出能量 U，代价极高，电子倾向于各自占据一个格点，形成局域磁矩。此时系统虽然半满，却因强排斥而无法导电，这就是莫特绝缘体。

需要强调的是，莫特绝缘体与能带绝缘体有本质区别。能带绝缘体的能隙来自晶格周期性导致的布里渊区折叠，而莫特绝缘体的能隙来自电子关联：一个电子要移动，必须与另一个电子交换位置，过程中产生双占据，需要克服能量 U。因此莫特绝缘体的能隙尺度由 U 与 t 的比值决定，而非晶格势。

三、莫特绝缘体与磁性序

半满填充（每个格点平均一个电子）是哈伯德模型最富饶的区域。在 U ≫ t 极限下，系统可以用微扰论处理。将跳跃项视为微扰，二阶微扰给出有效的自旋相互作用。通过正则变换或Schrieffer-Wolff变换，可以消去双占据态，得到有效哈密顿量：

H_eff = J Σ_{⟨ij⟩} (S_i · S_j - 1/4)

其中 J = 4t²/U > 0，S_i 是格点 i 上的自旋算符。这正是反铁磁海森堡模型。也就是说，半满强关联哈伯德模型在低能下等价于反铁磁自旋模型。这一结果解释了为什么许多过渡金属氧化物在低温下表现出反铁磁序：超交换作用通过虚跳跃过程产生，相邻格点的自旋倾向于反平行排列以降低能量。

当系统偏离半满，即掺杂少量空穴或电子时，反铁磁序被破坏，系统可能进入超导、奇异金属或自旋玻璃等相。铜氧化物高温超导体的母体正是反铁磁莫特绝缘体，掺杂后出现d波超导，这是哈伯德模型研究最重要的实验驱动力之一。

四、平均场近似与相图

严格求解哈伯德模型在二维及以上维度是NP难的，因此平均场近似是理解相图的重要工具。最常用的方法是引入磁序参量 m = ⟨n_{i↑} - n_{i↓}⟩，将相互作用项线性化：

U n_{i↑} n_{i↓} ≈ U[⟨n_{i↑}⟩ n_{i↓} + n_{i↑} ⟨n_{i↓}⟩ - ⟨n_{i↑}⟩⟨n_{i↓}⟩]

在动量空间，自旋向上和向下电子的色散关系发生劈裂：

ε_{kσ} = ε_k - σ U m / 2

其中 ε_k = -2t(cos k_x + cos k_y) 是二维方晶格的紧束缚色散，σ = +1 对应自旋向上，σ = -1 对应自旋向下。自洽条件要求 m 满足：

m = (1/N) Σ_k [f(ε_{k↑}) - f(ε_{k↓})]

其中 f 是费米分布函数。通过数值迭代求解该自洽方程，可以得到不同 U/t 和填充 n 下的磁矩大小。当 U 超过临界值 U_c 时，非零解出现，系统进入反铁磁相。平均场给出的相图定性上正确：半满时 U_c ≈ 0（因为费米面嵌套完美），偏离半满后 U_c 增大，磁性区域缩小。

五、超越平均场

平均场近似忽略了自旋涨落和空间关联，在低维和临界区域误差较大。更精确的方法包括：动力学平均场理论（DMFT），将格点问题映射为杂质模型，能正确描述莫特转变；量子蒙特卡洛（QMC），在有限温度下对哈伯德模型进行无偏模拟，但受费米子符号问题限制；张量网络方法（DMRG、PEPS），适合一维和二维小系统。这些方法共同构成了强关联物理的现代工具箱。

理解哈伯德模型的关键在于抓住 t 与 U 的竞争：t 使电子离域，U 使电子局域。半满时这一竞争产生莫特绝缘体和反铁磁序；掺杂时则孕育出高温超导等奇异量子态。掌握平均场自洽计算，是进入强关联数值模拟的第一步。