Twisting the Hubbard model into the Momentum-Mixing Hatsugai-Kohmoto Model

本文介绍了动量混合 Hatsugai-Kohmoto (MMHK) 模型，这是一个可连续变形的框架，它系统地将动量混合恢复到精确可解的 Hatsugai-Kohmoto 模型中，从而准确地重现了哈伯德模型（Hubbard model）的强关联物理特性，且与标准的有限簇技术相比具有更优的收敛速率。

要点

想象一下，你正试图理解一个房间里人群的行为。在物理学世界中，这些“人”是电子，而“房间”是晶格。描述这些电子如何相互作用的最著名规则书被称为哈伯德模型（Hubbard Model）。它是理解像铜氧化物超导体（可以无电阻导电的材料）这类材料的金标准。

然而，这里有一个难点：哈伯德模型极难求解。这就像试图预测一个“冲撞舞池”（mosh pit）中每个人的精确路径，而他们所有人都在互相碰撞。其中的数学计算变得极其复杂，以至于即使是最强大的超级计算机也难以得到完美答案，尤其是在处理二维材料（如原子薄片）时。

另一方面，存在一个更简单的、“作弊码”式的模型，叫做鸠谷-小本模型（Hatsugai-Kohmoto, HK model）。它很容易求解，但它有点“假”。它假设电子只有在处于完全相同的“座位”（动量态）时才会关心彼此，忽略了在现实世界中，电子是基于物理位置进行相互作用的。这就像是说，房间里的人只有在戴着完全相同的帽子时才会发生碰撞，而忽略了他们可能会撞到站在旁边的人。

核心思想：扭转“作弊码”

本文作者提出了一个聪明的疑问：我们能否将这个简单的“作弊码”模型进行缓慢的“扭转”，使其变成那个真实的、困难的模型，同时又不失去我们求解它的能力？

他们说：“可以。”他们创建了一个新模型，称为动量混合鸠谷-小本（MMHK）模型。

他们使用了如下类比：

• 旧方法（HK 模型）： 想象你有一个拥有 100 个座位的房间。在 HK 模型中，你根据人们的“帽子颜色”（动量）进行分组。如果两个人的帽子相同，他们就会互相排斥。但戴着不同帽子的人永远不会发生相互作用。这太简单了。
• 新方法（MMHK 模型）： 作者说，“让我们把情况搞复杂一点。”他们选取一小组座位（例如 2、4 或 10 个座位），并强迫坐在其中的人交换位置并产生相互作用。他们称之为“混合动量”。
  • 如果混合 2 个座位，你会得到一个稍好一点的近似值。
  • 如果混合 4 个座位，它会变得更好。
  • 如果混合 10 个座位，它将变得极其精确。

神奇的结果：速度与精度

他们发现的最令人惊讶的部分是这种方法运行得有多快。

通常，当科学家试图通过增加更多的组成部分（比如增加混合座位的数量）来近似一个复杂系统时，精度的提升是非常缓慢的，就像在爬一座平缓的小丘。如果你将座位的数量增加一倍，你只能稍微接近真相一点点。

作者发现，他们的 MMHK 模型就像是一枚火箭。

• 当他们将混合座位的数量从 1 增加到 10 时，该模型不仅是变好了一点，而是达到了与真实哈伯德模型 99% 的准确度。
• 他们称之为“平方律”式的提升。这意味着，如果你付出两倍的努力（混合两倍的动量），你就能获得四倍的精度。这比目前使用的标准方法要快得多。

他们证明了什么？

他们在两种情景下测试了这个新模型：

1. 一维情况（原子链）： 他们将自己的结果与唯一的已知完美解（贝特拟阵，Bethe Ansatz）进行了对比。仅使用 10 个混合动量，他们的模型就达到了与完美答案 1% 以内的误差。而标准方法需要数千个原子才能达到如此接近的程度。
2. 二维情况（平面薄片）： 这是哈伯德模型通常无法求解的“困难模式”。他们将该模型应用于方格网格。即使使用较少的混合动量（如 4 或 16 个），他们的模型也成功重现了真实材料的所有“特性”，例如：
   • 莫特转变（Mott Transition）： 材料如何突然停止导电并变为绝缘体。
   • 反铁磁性（Antiferromagnetism）： 电子自旋如何呈现出棋盘格状排列。
   • 伪能隙（Pseudogaps）： 一种神秘的状态，材料在此状态下表现得既像金属又像绝缘体。
   • 热容（Heat Capacity）： 材料如何储存热量，展现出区分电荷和自旋行为的明显峰值。

为什么这很重要？

把 MMHK 模型看作是一个高保真模拟器。

• 旧的模拟器： 为了获得清晰的图像，你需要一台庞大且昂贵的超级计算机运行数天，而且你仍然不敢确定结果是否完美。
• MMHK 模拟器： 你可以使用一个微小、简单的设置，获得 99% 清晰度的图像。它捕捉到了复杂物理学的“灵魂”（即莫特物理学），同时保持了数学上的可解性。

作者总结道，该模型为物理学家提供了一个全新的、强大的工具。它允许研究人员以一种此前无法实现的快速且精确的方式，去研究强电子相互作用（这是理解高温超导体的关键），而这一切仅仅是通过将少量动量态混合在一起实现的。

简而言之： 他们找到了一种方法，能将一个简单的、可解的玩具模型，转化为一个高度精确的、能够模拟真实复杂电子世界的模型，而且他们做这件事的效率极高。

技术摘要

技术摘要：通过动量混合将哈伯德模型扭转为哈茨盖-科莫托模型

问题陈述
哈伯德模型（Hubbard model）是理解强关联电子系统（如铜氧化物超导体）的基础理论框架。然而，由于动能项与势能项的非对易性，该模型在 $d=2$ 及 $d>1$ 维下的精确解仍然难以求得。这种非对易性阻碍了通过双占据数对本征态进行直接排序，使得标准的解析技术难以奏效。虽然数值方法（如量子蒙特卡洛 QMC、密度矩阵重整化群 DMRG 和动力学平均场理论 DMFT）提供了深刻见解，但它们面临着符号问题、有限尺寸缩放问题或需要集群扩展以捕捉非局部关联等挑战。相反，精确可解的哈茨盖-科莫托（Hatsugai-Kohmoto, HK）模型通过动量空间相互作用捕捉了莫特物理（Mott physics），但由于它假设动量空间具有局域性（在实空间具有长程性）且缺乏动量混合，因此无法重现关键的哈伯德特征（例如反铁磁性、特定的谱权重转移）。本文所解决的核心挑战是寻找一种系统的方法，在不牺牲精确可解性的前提下，在可解的 HK 模型与未解的哈伯德模型之间进行插值。

方法论：动量混合 HK (MMHK) 模型
作者提出了动量混合哈茨盖-科莫托（MMHK）模型，它是标准 HK 模型的一种连续形变，重新引入了动量混合。

• 构建： 该方法涉及将 $n$ 个不同的动量归为一组“单元（cell）”并进行杂化。在 $n$-MMHK 模型中，每个位于缩减布里约拉区（rBZ）中的动量态 $k$ 都被 $n$ 个位点（或轨道）所装饰。
• 哈密顿量： 动能项变为一个描述单元内 $n$ 个位点间跳跃的矩阵 $g_{\alpha\alpha'}(k)$，而相互作用项在混合基组下保持对角，但作用于每个动量态对应的 $n$ 个位点上。
• 极限： 当 $n \to N$（其中 $N$ 是总格点数）时，缩减布里约拉区收缩为单个点，此时相互作用项在实空间变为纯局域的，从而恢复为标准的单位点哈伯德模型。
• 可解性： 尽管引入了混合导致了非对易性，但由于该哈密顿量可以简化为针对每个 $k$ 点求解 $n$ 位点集群的本征态问题，因此它对于适中的 $n$ 值仍是精确可解的。

主要贡献与结果

1. 一维情况下的快速收敛：

   • 在一维情况下，$n$-MMHK 模型的基态能量以极快的速度收敛至哈伯德模型的精确贝特拟阵（Bethe ansatz）结果。
   • 当 $n=10$ 时，偏差小于 1%。
   • 收敛速度遵循 $1/n^{\gamma}$ 规律，其中 $\gamma \approx 2$（具体取决于 $U$，范围在 $1.83$到$2.07$ 之间），其表现显著优于标准有限集群 DMRG 计算（其标度为反比例线性 $1/L$）。
   • 该模型正确捕捉到了当 $U \to 0$ 时莫特能隙（Mott gap）消失的现象，这是使用开边界条件下的有限尺寸计算经常遗漏的特征。

2. 二维物理的重现：

   • 将 $n$-MMHK 模型应用于方格点阵（使用 $n=4$ 和 $n=16$ 集群）重现了最先进的模拟结果。
   • 莫特转变： 该模型正确预测了在半填充条件下，对于任何非零 $U$ 都会发生向绝缘态的转变（由于费米面嵌套），并捕捉到了正确的双占据数行为——其起始值为 $1/4$（符合非相互作用哈伯德极限），而非 $1/2$（带状 HK 模型的值）。
   • 磁性关联： 不同于带状 HK 模型的铁磁不稳定性，$n$-MMHK 模型表现出反铁磁（AF）关联作为主要的自旋响应，这与哈伯德物理一致。
   • 谱函数： 使用 $n=16$（即每个 $k$ 点对应一个 $4 \times 4$ 集群）时，其谱函数 $A(k, \omega)$ 与 QMC 和集群摄动理论的结果定量匹配，捕捉到了哈伯德能带的变窄以及在 $t' \neq 0$ 时的粒子-空穴不对称性。
   • 伪能隙： 当引入次近邻跳跃（$t'$）时，该模型成功地在欠掺杂区域的态密度（DOS）中产生了伪能隙，这与最近关于扭转边界条件集群的研究一致。
   • 热力学： 热容表现出将电荷和自旋自由度分开的双峰结构，其低温行为遵循与电荷中性激发一致的幂律。

3. 动力学谱权重转移 (DSWT)：

   • $n$-MMDK 模型捕捉到了上、下莫特子带之间非平凡的动力学混合。
   • 它重现了下哈伯德带权重超过简单电子计数（$1+x$）的现象，这一特征此前仅能由复杂的数值方法捕捉，而无法由解析可解模型实现。

标度与理论依据
论文指出，$n$-MMHK 模型的效率源于通过动量混合引入的量子涨落。收敛速率在 1D 中遵循 $1/n^2$ 标度，在准 2D 梯形晶格中可能更快（达到 $1/n^3$），这表明在高维情况下，只需较少的混合动量即可实现收敛。从理论上讲，作者认为 HK 模型和哈伯德模型属于同一普适类，该类由莫特不动点处离散 $Z_2$ 对称性的破缺所支配。由于动量混合构造构成了一种保持 Luttinger 表面（格林函数零面）的局域形变，因此电荷动力学是稳健的，并能快速收敛至哈伯德极限。

意义与主张
作者声称，MMHK 模型为研究强关联量子物质提供了一个强大的替代工具。其主要意义在于提供了一个系统、可控且精确可解的框架，在可解的 HK 极限与复杂的哈伯德模型之间进行插值。

• 效率： 它通过较少的混合动量（$n=4$ 到 $16$）即可实现定量准确性，相比大型集群数值方法具有计算优势。
• 洞察力： 它架起了解析可解性与数值复杂性之间的桥梁，使得能够对 DSWT 和伪能隙等现象进行解析研究，而这些现象在可解模型中通常是无法触及的。
• 普适性： 该工作表明，莫特转变的核心电荷物理对于局域形变是稳健的，这强化了这样一个观点：电荷能隙是莫特绝缘体的普适特征，而自旋扇区（如 AF 与 FM）的差异则源于具体的相互作用细节。

论文总结道，MMHK 模型可作为莫特/哈伯德物理的高效模拟器，能够以一种简化但忠实的平台，研究关联、拓扑与非平衡动力学之间的相互作用。