来源论文: https://arxiv.org/abs/2604.24887v1 生成时间: May 01, 2026 18:07

三角格点 Majorana-Hubbard 阶梯中的对称保护拓扑相：深度数值与理论解析

0. 执行摘要

本文基于 Will Holdhusen 等人的最新研究，对四条腿（four-leg）三角格点 Majorana-Hubbard (MH) 阶梯模型进行了系统的数值调研。该研究利用密度矩阵重整化群（DMRG）和变分均匀矩阵乘积态（VUMPS）算法，纠正了早期研究中关于强耦合区“能隙消失”的误判。研究发现，该模型在不同的相互作用强度下展现出极其丰富的相空间，包括多个对称保护拓扑（SPT）相（$G_1, G_3, G_4$）以及一个具有动力学临界指数 $z=2$ 的量子 Lifshitz 临界相（$GL_4$）。特别地，$G_4$ 相展示了拓扑有序与自发对称性破缺（SSB）的独特性结合。本文不仅绘制了精确的相图，还深入探讨了纠缠谱退化、费米子平移对称性和边界条件在稳定这些拓扑相中的核心作用，为在超导体表面涡旋阵列中实现 Majorana 晶格提供了关键的理论依据。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题

Majorana 零能模（MZMs）由于其非阿贝尔统计特性，被认为是容错拓扑量子计算的核心候选者。然而，大多数研究集中在非相互作用或弱相互作用的 Majorana 系统。Majorana-Hubbard 模型引入了四费米子相互作用，探讨了在强关联背景下，拓扑序是如何被修改、破坏或增强的。本文的核心问题在于：在三角格点这一具有高度几何受挫性的结构中，Majorana 费米子的相互作用会诱导哪些新奇的拓扑相？尤其是之前研究中观察到的“无能隙相”究竟是物理本质还是数值赝像？

1.2 理论基础：MH 模型与几何结构

模型哈密顿量由两部分组成：

$$H = H_0 + H_I$$

其中，$H_0$ 描述了 Majorana 费米子在三角格点上的近邻跃迁（hopping）：

$$H_0 = it \sum_{\langle p,q \rangle} \eta_{p,q} \gamma_p \gamma_q$$

$\eta_{p,q}$ 是格点上的规范场，反映了超导体中的 $\pi$ 通量。相互作用项 $H_I$ 则作用于三角形或四边形斑块（plaquettes）：

$$H_I = g (P_1 + P_2 + P_3)$$

每个 $P_\mu = \gamma_1 \gamma_2 \gamma_3 \gamma_4$ 是四格点 Majorana 算符的乘积。对于三角格点，这种四体相互作用是描述强相互作用机制的最小模型。

1.3 技术难点：几何受挫与有限尺寸效应

在 2D 三角格点上直接进行数值模拟存在巨大的计算复杂度。因此，作者采用了“四条腿阶梯”模型作为折中。然而，阶梯模型的拓扑特性高度依赖于横向（y方向）和纵向（x方向）的边界条件（OBC, PBC, APBC）。

  • 受挫感诱导的无能隙模：当纵向长度 $L_x$ 为奇数时，系统的几何结构会强行引入一个畴壁（domain wall），导致在本来有能隙的相中出现伪临界行为。这是早期研究得出错误结论的主因。
  • 对称性群的表示：识别 SPT 相需要分析哈密顿量在特定对称性群（如 $\mathbb{Z}_2 \times \mathbb{Z}_2$）下的投影表示，这在 Majorana 语言下比在自旋语言下更为微妙。

1.4 方法细节：DMRG, VUMPS 与 Jordan-Wigner 变换

  • 数值工具：研究使用了基于 ITensor 库的 DMRG 算法。为了消除有限尺寸效应，作者引入了 VUMPS 算法来直接处理热力学极限下的无穷矩阵乘积态（iMPS）。这种方法能够精确提取纠缠熵和纠缠谱，从而识别 SPT 相的退化特征。
  • 算符映射：通过 Jordan-Wigner (JW) 变换将 Majorana 系统映射到一维自旋链： $$\sigma_m^z = i \gamma_m^r \gamma_m^b = 2c_m^\dagger c_m - 1$$ 这种映射使得我们可以利用自旋系统的成熟理论（如 Ising 模型、Potts 模型）来解读 Majorana 相位。
  • 平均场理论：作者开发了一个 12 参数的广义平均场方案，通过自共轭循环（self-consistency loop）计算基态能量。虽然平均场在高能区表现良好，但在强相互作用 $g < g_1$ 区间，它显著高估了能量，显示了多体关联的重要性。

2. 关键 Benchmark 体系，计算所得数据与性能数据

2.1 相图分析（见图 1）

研究确定了四个主要区域：

  1. $G_1$ 相 ($g > -0.62$)：这是传统的 Majorana 链拓扑相（类似 Kitaev 链）。其特征是纠缠谱的两倍退化，且在 OBC 条件下具有 Majorana 零能边缘模。此相在弱吸引到强排斥区间内非常稳定。
  2. $G_3$ 相 ($-0.74 < g < -0.62$)：一个狭窄的有能隙相。在环面上展示出三倍基态简并（在柱面上降为两倍），具有长程自旋串序（string order）。
  3. $GL_4$ 相 ($-1.15 < g < -1.09$)：一个量子 Lifshitz 相。其能隙 $\Delta \propto 1/L_x^2$，动力学临界指数 $z=2$。这是本研究的一大发现，表明在该参数区间系统处于临界状态，且受连续变换对称性的保护。
  4. $G_4$ 相 ($g < -1.29$)：强吸引区。这是一个独特的“退化 SPT 相”，它不仅是对称保护的，还伴随着 $\mathbb{Z}_2 \times \mathbb{Z}_2$ 对称性的自发破缺。在环面上具有 4 倍简并度，在柱面上为 8 倍。

2.2 性能数据与精度

  • 键维数 (Bond Dimension)：DMRG 计算中 $D$ 最高达到 2500。在有能隙相（如 $G_1, G_4$），截断误差（truncation error）保持在 $10^{-8}$ 以下。在临界点附近，由于关联长度增加，误差维持在 $10^{-5}$ 左右。
  • 能量灵敏度 ($\partial^2 e_0$)：通过对基态能量 $e_0$ 取相互作用强度 $g$ 的二阶导数，作者成功捕捉到了相变点的奇异性（见图 3 顶图）。
  • 纠缠谱 (Entanglement Spectrum)：在 $G_4$ 相中观察到了清晰的能级对退化（Alternating large/small circles），这是识别 SPT 相最为确凿的数值证据（见图 3 中图）。

2.3 边界条件对比

  • Toroidal (x-PBC/APBC)：用于观察基态简并度和费米子宇称切换（parity switch）。
  • Cylindrical (x-OBC)：用于检测 Majorana 边缘模的物理存在。研究表明，$G_1$ 和 $G_4$ 在 OBC 下均展现出稳定的边缘零能模。

3. 代码实现细节，复现指南与开源链接

3.1 代码实现架构

本工作主要基于 ITensor（可能是 Julia 版本，因其在现代张量网络计算中的普及性）。

  • Hamiltonian 构建： 核心难点在于实现三角格点的四体算符。代码通常会将 2D 坐标 $(i, j)$ 线性化为 1D 索引 $m = 2(i-1) + j$。对于每一格点，定义 $P_\mu$ 涉及 4 个位置的 Majorana 算符，通过 JW 变换转化为复杂的 Pauli 算子乘积串。
  • VUMPS 流程：
    1. 初始化一个小尺寸的 MPS 基态。
    2. 使用 ITensor 提供的 VUMPS 迭代器，通过交替优化中心张量和环境张量来逼近热力学极限。
    3. 监控算符期望值的收敛性，尤其是 $\langle \sigma^x \sigma^x \rangle$ 长程关联。

3.2 复现指南

  1. 环境配置：安装 Julia 以及 ITensors.jl 库。
  2. 格点定义：定义一个具有 XC4 嵌入的四腿阶梯（见论文图 2）。特别注意三角形跳跃算符的符号 $\eta_{p,q}$，必须符合 $\pi$-flux 配置。
  3. 相互作用项实现：
    # 示例伪代码：构建 Plaquette 相互作用
for i in 1:Lx, j in 1:2
    # 定义四个算符的位置
    sites = [idx(i,j), idx(i+1,j), idx(i+1,j+1), idx(i,j+1)]
    # 构建四体项 P = γ1*γ2*γ3*γ4
    os += g, "Majorana", sites[1], "Majorana", sites[2], "Majorana", sites[3], "Majorana", sites[4]
end
  4. 相变点扫描：在 $g \in [-3.0, 0.0]$ 区间进行步长为 0.05 的扫描，观察 entropy 和 gap 的变化。

4. 关键引用文献与局限性评论

4.1 关键引用文献

  1. A. Yu. Kitaev (2001)：Majorana 链的奠基之作，定义了 $G_1$ 相的原型。
  2. Fu & Kane (2008)：提出了在 TI/超导体界面实现 Majorana 的物理图景，是本文背景的基石。
  3. Tummuru, Nocera, & Affleck (2021)：早期关于 TLMH 模型的 DMRG 研究，本论文在其基础上进行了修正。
  4. Pollmann et al. (2010)：定义了通过纠缠谱识别 1D SPT 相的标准方法。

4.2 局限性评论

  • 维度限制：虽然 4-leg 阶梯比 1D 链更接近 2D，但它仍然无法完全捕获 2D 拓扑序（如手性 Majorana 边缘模）。阶梯系统的 SPT 相在真正 2D 极限下是否演变成更复杂的非阿贝尔拓扑有序态仍不确定。
  • 未标记区域的模糊性：在 $g \approx -0.9$ 附近的“未标记区域”，DMRG 结果表现出强烈的不一致性。这可能是由于长程有序带来的巨大原胞（unit cell），超出了当前张量网络的模拟能力。作者坦诚地将其留作未来的工作，这说明该区域可能隐藏着更复杂的空间非均匀相（如电荷密度波或条纹相）。
  • 实验可行性：虽然理论上这些相可以在涡旋晶格中实现，但调整相互作用强度 $g$ 在实际材料（如 $Bi_2Te_3/NbSe_2$）中极具挑战，需要极其精确的化学势和磁场控制。

5. 补充内容：从自旋物理到 Majorana 拓扑的深度关联

5.1 费米子宇称切换的物理意义

在 $G_1$ 到 $G_4$ 的相变过程中，一个显著的特征是系统的基态费米子宇称（Parity）会随着系统尺寸或边界条件的变化而切换。在有能隙的拓扑相中，这种宇称切换直接关联到 Majorana 零能模的存在：当边界条件从 PBC 变为 APBC 时，Majorana 模的占据态会发生改变，从而导致总宇称的变化。这一现象是拓扑序在有限尺寸系统中的“指纹”。

5.2 $G_4$ 相：拓扑与破缺的共生

$G_4$ 相最引人注目的地方在于它不是纯粹的 SPT 相。在自旋语言下，它对应于一个 8 态 Potts 模型（$g > 0$）或类似的关联态。它自发地破缺了某种平移对称性（表现为 6 格点或 4 格点的周期性，见图 5），但同时保留了纠缠谱的退化。这种“具有自发对称性破缺的 SPT 相”打破了传统的 Landau 范式，暗示了在 Majorana 晶格中，相互作用可以同时诱导凝聚态物理中最强大的两种机制：对称性破缺和拓扑保护。

5.3 对量子计算的启示

本文确定的 $GL_4$ 相（具有 $z=2$ 临界指数）具有极高的纠缠度，虽然它是无能隙的，但其拓扑分数的某种形式可能被用于非局部信息编码。此外，明确了强吸引区 $G_4$ 的有能隙本质，实际上为在该区域寻找稳定的拓扑比特提供了更广阔的参数空间，而非此前担心的临界不稳定性。