Topological Phase Transition in the Two-Leg Hubbard Model: Emergence of the Haldane Phase via Diagonal Hopping and Strong Interactions

本研究通过利用密度矩阵重整化群（DMRG）模拟证明，在双链哈伯德模型中，由对角跳跃诱导的几何受挫与强在位库仑排斥之间的相互作用，驱动了从平凡绝缘体到对称保护的哈尔丹相（Haldane phase）的拓扑相变，并由独特的磁性、边缘相关性和非零的弦序参数所证实。

要点

想象一个由“梯子”组成的微观世界，它由两条平行的轨道构成。在这个世界里，微小的粒子——电子——正试图四处奔跑。通常情况下，电子喜欢留在自己的轨道上，在相邻的横档之间跳跃。但在这一项特定的研究中，研究人员增加了一个转折：他们允许电子进行“对角线跳跃”，即以一定的角度从一条轨道跳到另一条轨道，就像斜着穿过马路而不是走人行横道一样。

研究人员想要观察，当这种“对型跳跃”与一种让电子讨厌同时出现在同一个位置的规则（强排斥作用）相结合时，会发生什么。

以下是他们发现的故事，通过简单的概念进行了拆解：

1. 设置：一个受挫的梯子

把电子想象成试图在两条车道舞池上跳舞的人。

• 规则： 他们可以在自己的车道上向前跳舞，可以跳到另一条车道，或者——这是新加入的部分——可以斜着跨越车道跳舞。
• 冲突： 电子还有一个规则，那就是他们非常不喜欢共享同一个位置（就像两个人试图坐在一把椅子上）。
• 目标： 科学家们想看看这种特定的对角线移动与“不共享”规则的结合，是否能创造出一种特殊的、隐藏的物质状态。

2. 发现：“哈尔丹”（Haldane）相

他们发现，当对角线跳跃足够强，且“不共享”规则足够严格时，电子会进入一种非常特殊的态，称为哈尔丹相。

你可以把这个相位想象成电子进行的一种秘密握手。

• 在普通状态下： 电子只是随机地跳舞，或者呈现出一种简单的模式。
• 在哈尔丹态中： 电子形成了一种隐藏的、长程的有序结构。这就像一队人，每个人都以一种特定的、复杂的模式牵着手，你无法仅仅通过观察两个相邻的人来理解这种模式。你必须观察整条队伍才能理解其中的奥秘。

3. 他们是如何知道它很特别的（证据）

研究人员并非凭空猜测；他们使用了一种强大的计算机方法（称为 DMRG）来模拟系统，并寻找这种特殊相位的“指纹”：

• “边缘”效应（漂浮的幽灵）：
  在普通的梯子模型中，如果你观察两端（边缘），不会发生任何特别的事情。但在这种哈尔丹相中，梯子的末端表现得好像有自己的小“幽灵”电子在漂浮。尽管梯子的中间部分很平静，但末端却是躁动且具有磁性的。这就像一根绳子，中间被固定住了，但两端却在自由摆动。

• “弦”序（隐形的线）：
  他们发现了一种连接电子的数学“弦”。如果你拉着一根线穿过梯子中间，电子会对这种连接做出特定反应，从而证明它们被这条隐形的线连接在一起。这是“拓扑”序的一个标志——一种基于形状的属性，非常难以被破坏。

• “能隙”（能量墙）：
  在物理学中，“能隙”就像是一堵能量墙，你需要跨越它才能改变系统的状态。

  • 在梯子中间，存在一个坚实的墙（能隙），使电子保持稳定。
  • 但在非常边缘处，这堵墙消失了，允许“幽灵”电子自由移动。这种组合（稳定的中间，自由的边缘）是这种拓扑相的特征。

• “纠缠”（双生连接）：
  当他们将梯子从中切开，观察左侧和右侧之间的连接时，他们发现了一种完美的对称性。左侧的电子与右侧的电子以一种产生完美镜像（简并）的方式“纠缠”在一起。这就像拥有两个双胞胎，无论相隔多远，他们总是保持完美的同步运动。

4. “穹顶”形状

他们地图中最有趣的部分是这个特殊区域的形状。

• 如果没有对角线跳跃，电子就处于一种普通的、乏味的态。
• 如果相互作用过强或者对角线跳跃过少，情况依然是普通的。
• 但是，在他们的地图中间有一个穹顶状区域。在这个穹顶内部，特殊的哈尔丹相存在。这就像是一个“金发姑娘区”（适中区），在这里，对角线跳跃和电子间的排斥力恰到好处地结合在一起，创造了这种神奇的、隐藏的秩序。

总结

论文表明，通过让电子在双轨梯子上进行对角线跳跃，并强制它们保持距离，你可以迫使它们进入一种拓扑相。这种相位很特别，因为它具有：

1. 隐藏的秩序（你无法通过简单的局部检查来观察到它）。
2. 边缘处的磁性“幽灵”（这在中间是不存在的）。
3. 鲁棒性，意味着即使稍微改变条件，它也能保持稳定。

研究人员通过展示边缘处的电子行为与中间处不同，并证明了连接它们的“弦”依然稳固，证实了这一点。这有助于我们理解复杂材料可能如何表现，并为科学家构建未来的量子材料提供了新的目标。

技术摘要

技术摘要：双链 Hubbard 模型中的拓扑相变

问题陈述
本研究探讨了双链 Hubbard 模型中几何挫折与强电子-电子相互作用之间的相互作用。具体而言，作者研究了在引入对角跳跃（$t_d$）以及结合最近邻跳跃（$t$）和梯级跳跃（$t_\perp$）的情况下，结合在位库仑排斥力（$U$），如何影响拓扑相的涌现。虽然先前的研究已经确立了自旋-1 海森堡链中 Haldane 相的存在，以及特定极限下的 Hubbard 梯子（例如铁磁梯级耦合或非配对边缘位点），但在存在电荷涨落和非均匀对角跳跃的情况下，Haldane 相的行为仍是一个开放性问题。核心问题在于，当系统处于偏离严格海森堡极限（即电荷涨落被抑制的极限）且对角跳跃振幅不同于梯级跳跃（即 $t_d \neq t$）时，Haldane 相是否依然存在。

方法论
作者采用密度矩阵重整化群（DMRG）方法对双链 Hubbard 梯子的哈密顿量进行数值求解。

• 模型： 该系统由两个耦合的 Hubbard 链组成，包含内梯跳跃 $t$、梯级跳跃 $t_\perp$ 和对角跳跃 $t_d$。哈密顿量包含在位库仑排斥项 $U$。研究在半填充状态下进行，设定 $t$ 为能量单位，且 $t_\perp = t$。
• 参数： 研究覆盖了广泛的参数空间，对角跳跃比 $t_d/t$ 从 0（非挫折二部系）到 1 进行变化，相互作用强度 $U/t$ 最高达 16。
• 系统尺寸： 计算在长度范围为 $L=28$ 到 $L=112$ 个梯级的系统上进行。DMRG 模拟使用了高达 $m=2000$ 个密度矩阵特征态，截断误差在 $10^{-9}$ 到 $10^{-11}$ 之间，并守恒总粒子数和总磁化强度 $S_z$。
• 分析工具： 作者使用多种观测物理量来表征物相：自旋-自旋相关函数（梯级、梯级间及对角相关）、自旋结构因子、边缘相关、弦序参数、自旋能隙（特别是第一和第二能隙）以及纠缠谱，以及电荷能隙。应用了有限尺寸缩放分析以验证热力学极限行为。

主要贡献与结果
论文证明了在费米子双链 Hubbard 模型中，由对角跳跃引起的挫折与强相互作用共同驱动，产生了对称保护拓扑（SPT）Haldnane 相。

1. 相图： 在 $(t_d/t, U/t)$ 平面上识别出一个明显的穹顶状区域，在该区域内存在 Haldane 相。该相由拓扑平凡的莫特绝缘体所界定。
2. 磁性性质与相关性：
   • 相关性的符号改变： 转变的一个关键特征是对角自旋-自旋相关函数 $\langle S^z_{0,r} S^z_{1,r+1} \rangle$ 符号的变化。在平凡相中，这些相关性为正；而在 Haldane 相中，它们变为负值，在强耦合极限下趋于理论推导的自旋-1 海森堡链的值（$\epsilon_H/12 \approx -0.117$）。
   • 结构因子： 代表反铁磁相关性的自旋结构因子 $S_z(\pi, \pi)$ 在穹顶内显著受到抑制，而 $S_z(\pi, 0)$ 表明了与梯级三重态占优一致的磁有序转变。
3. 边缘态与相关性：
   • 边缘相关长度： 边缘相关性的相关长度（$\xi_{edge}$）在相边界处发散。在 Haldane 相内部，边缘相关性衰减缓慢并在相反边界处再次上升，这表明存在局域化的自旋-1/2 边缘态。
   • 自旋分布： 对剩余自旋分布的分析表明，在 Haldane 相中，自旋激发在梯子边缘呈指数级局域化，而在平凡相中，自旋则集中在体相中。
4. 弦序： 通过 Kennedy-Tasaki 变换定义的弦序参数在穹顶内非零，并在热力学极限下保持稳定。相反，Néel 序参数随系统尺寸增加而消失，证实了不存在传统的自发对称性破缺以及该序的拓扑性质。
5. 自旋能隙：
   • 第一能隙 ($\Delta E_1$)： 对应于从基态（$S_{tot}=0$）到第一激发三重态（$S_{tot}=1$）的跃迁，该能隙在 Haldane 相内随系统尺寸呈指数级闭合，指示了无能隙边缘激发（Kennedy 三重态）的存在。
   • 第二能隙 ($\Delta E_2$)： 在热力学极限下保持有限，代表体相 Haldane 能隙。
6. 纠缠谱： 纠缠谱在 Haldane 相内表现出特征性的偶数简并，这是对称保护（特别是 $Z_2 \times Z_2$、时间反演或反演对称性）的标志。Schmidt 谱中的最大能隙（$\Delta \lambda$）随着系统尺寸和相互作用强度的增加趋于零，这与在平凡相中观察到的有限能隙形成对比。
7. 电荷能隙： 系统在所研究的参数范围内始终保持为莫特绝缘体，其电荷能隙 ($\Delta E_c$) 在强耦合机制下随 $U$ 线性增长。

意义与主张
作者声称，他们的工作明确展示了 Haldane 相可以通过几何挫折（由对角跳跃引入）和强在位库仑排斥力的互补效应在费米子系统中产生。

• 鲁棒性： 一个关键发现是，尽管存在通常预期会破坏此类相的电荷涨落，Haldane 相依然保持稳健。该相在广泛的对角跳跃振幅范围内都是鲁棒的，超出了将系统精确映射为自旋-1 海森堡链的特定点（$t_d = t$）。
• 机制： 研究阐明了这种转变是由挫折与相互作用之间的相互作用驱动的，导致了有效实现模拟自旋-1 链的梯级三重态。
• 实验相关性： 作者指出，他们的发现为未来在强关联梯子材料和量子模拟平台中的实验探索提供了指导，特别提到了近期关于梯子几何构型中超冷费米子原子的实验。他们强调，其结果有助于定位这些系统中平凡相与拓扑相之间的相边界。

论文总结认为，对角跳跃与强相互作用的结合是产生对称保护拓扑序的一种可行机制，为理解相互作用费米子系统中的拓扑相提供了更深入的理解。