来源论文: https://arxiv.org/abs/2605.11823v1 生成时间: May 13, 2026 05:22

拓扑绝缘体遇见多体相互作用：基于绝热量子模拟的 SSHH 模型深度解析

0. 执行摘要

在量子计算迈向实用化的进程中，模拟强关联多体系统被认为是最具潜力的“杀手级应用”之一。传统的变分量子特征求解器（VQE）虽然在近期的噪声中等规模量子（NISQ）设备上表现出色，但在处理具有全局对称保护或复杂拓扑特征的态时，往往受限于 Ansatz 的表达能力，难以精确捕捉波函数的拓扑本质。相比之下，绝热量子模拟（AQS）通过缓慢驱动系统跨越哈密顿量空间，为获取基态提供了一条更为物理且鲁棒的路径。

本文深入探讨了一项里程碑式的工作：如何利用门阵列量子计算机对一维费米子 Su-Schrieffer-Heeger-Hubbard (SSHH) 模型进行绝热模拟。SSHH 模型不仅包含了经典的 SSH 模型所描述的单粒子拓扑相变，还引入了 Hubbard 项来刻画电子间的库仑排斥。研究证明，拓扑特征（如多体 Berry 相和子格极化）在弱相互作用下具有显著的稳定性，但随着手性对称性破缺项的增强，这些特征最终会发生瓦解。该方案在量子比特需求、门复杂度及测量成本上均实现了多项式级扩展，为未来在大规模量子硬件上探测相互作用多体系统的拓扑性质奠定了理论与算法基础。

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

1.1 核心科学问题：多体相互作用对拓扑稳定性的挑战

拓扑绝缘体的核心特征通常在单粒子图景下通过能带结构（如 Zak 相或绕数）来理解。然而，真实材料中电子之间的库仑相互作用是不可忽略的。当我们将 SSH 模型扩展为 SSHH 模型时，系统从单粒子动力学转变为真正的多体问题。科学界长期关注的一个核心问题是：在强关联环境下，受对称性保护的拓扑序是否依然稳健？ 特别是当相互作用项（如不平衡的 Hubbard U 项）显式破坏了系统的手性对称性时，拓扑态会如何演化并最终失效？

1.2 理论基础：从 SSH 到 SSHH

SSH 模型是描述聚乙炔等一维二聚化链的范式模型。其哈密顿量定义在具有 A、B 子格的单元格上，通过胞内跳跃 $v$ 和胞间跳跃 $w$ 的竞争产生拓扑非平凡相（$|w| > |v|$）和拓扑平凡相（$|v| > |w|$）。

SSHH 模型通过引入 Hubbard 项将其推向多体领域：

$$H_{Hubbard} = \sum_{j=1}^{N} (U_A n_{Aj\uparrow}n_{Aj\downarrow} + U_B n_{Bj\uparrow}n_{Bj\downarrow})$$

其中 $U_A$ 和 $U_B$ 代表不同子格上的在位排斥能。当 $U_A = U_B$ 时，系统保留手性对称性；当 $U_A \neq U_B$ 时，手性对称性被显式破坏，这为研究拓扑特征的鲁棒性提供了完美的试验场。

1.3 技术难点：希尔伯特空间的指数爆炸与能隙闭合

  1. 态空间规模：对于 $N$ 个单元格，费米子态空间维度随粒子数呈指数增长，传统的精确对角化（ED）难以处理超过 20 个轨道的系统。
  2. 绝热条件的维持：绝热演化的成功依赖于系统始终处于瞬时基态。根据绝热定理，演化速度必须远慢于能隙的平方反比。在多体系统中，相互作用可能导致能隙变窄甚至闭合（相变点），如何选择合适的演化路径 $T$ 是巨大的挑战。
  3. 费米子到量子比特的映射：利用 Jordan-Wigner 变换将费米子算符转换为 Pauli 算符时，会引入复杂的串行相位（parity strings），导致量子电路深度增加。

1.4 方法细节：绝热演化协议

本研究采用了以下算法流程：

  1. 初始态制备：首先在量子计算机上制备非相互作用 SSH 模型的基态。由于该态是 Slater 行列式，可以通过解析解构建高效的单粒子变换电路。利用 OpenFermion 提供的矩阵分解技术，将基态变换转化为 $O(N^2)$ 数量级的 Givens 旋转门。
  2. 时间演化算符的分解：采用一阶 Trotter-Suzuki 分解。将总演化时间 $T$ 分为 $L$ 个步长 $\delta t = T/L$。在每个步长内，演化算符被近似为： $$U_\ell \approx \exp(-i \frac{\delta t}{\hbar} H_1(t_\ell)) \exp(-i \frac{\delta t}{\hbar} H_0)$$ 其中 $H_0$ 是 SSH 轨道项，$H_1$ 是随时间线性增加的 Hubbard 相互作用项。
  3. 多体 Berry 相的测量：这是本工作的技术亮点。研究者利用位置算符 $\hat{X}$ 在计算基下的对角特性，通过单基测量估算算符 $z_N = \langle \Psi | \exp(i 2\pi \hat{X} / N) | \Psi \rangle$，从而提取 Berry 相 $\gamma_N = \text{Im} \ln z_N$。这种方法规避了复杂的辅助量子比特干涉测量。

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

2.1 Benchmark 体系设置

研究团队选择了一个具有 $N=6$ 单元格的体系，对应 24 个量子比特。模拟涵盖了两种边界条件：

  • 周期性边界条件 (PBC)：用于提取多体 Berry 相（反映体态拓扑性）。
  • 开放边界条件 (OBC)：用于观察边缘态的子格极化（反映体-边对应关系）。

2.2 Trotter 误差与保真度分析

在验证阶段，研究者通过对比经典精确解，评估了 Trotter 分解的收敛性（见论文图 2）：

  • 非相互作用情形 (U=0)：当 $T=1, 5$ 时，保真度随步数 $L$ 增加迅速饱和至 1。但在极慢演化（$T=80$）下，由于累积误差，需要极大的 $L$ 才能维持保真度。
  • 相互作用情形 (U=1)：保真度展现出类似的趋势，证明了绝热协议在处理多体项时的稳定性。最终选择 $T=1, L=40$ 作为计算 Berry 相的平衡点。

2.3 关键计算数据：Berry 相的跃迁

在体态模拟中，研究者固定 $v=1$，调节 $w$ 从 0 到 2（跨越拓扑相变点 $w=1$）：

  • 对称性平衡（$U_A = U_B = 0.01$）：Berry 相在 $w=1$ 处表现出完美的阶梯状跳变（从 $\pi$ 到 $0$），这与单粒子 SSH 模型一致，证明了弱相互作用下拓扑序的稳健性。
  • 对称性破缺（增加 $\Delta U = U_B - U_A$）：随着 $\Delta U$ 增加至 0.3，原本陡峭的阶梯变得平滑，Berry 相偏离了量子化数值（见图 3）。当 $\Delta U \gtrsim 0.1$ 时，对称性保护失效，系统不再被视为标准的拓扑绝缘体。

2.4 子格极化的空间分布

在 OBC 下，研究者观察了 $N=2N+2$（半填充加两个额外电子）情形下的子格极化 $P_j^e$：

  • 在弱 $\Delta U$ 限制下，极化强度集中在链的两端（单元格 1 和 6），而体态极化几乎为零。这是拓扑边缘态存在的直接证据。
  • 随着 $\Delta U$ 增大，边缘极化逐渐向体态渗透，揭示了强手性破缺导致的边缘态“解域化”。

3.1 软件栈选择

该研究主要基于 IBM 的 Qiskit 框架进行量子电路的构建与经典仿真，并结合 OpenFermion 进行费米子代数的预处理。

  • Qiskit (v1.0+): 用于定义量子寄存器、实现 $R_{i,j}(\theta)$（XY 相互作用）、$G_{i,j}(\theta)$（受控旋转）以及 $CP_{i,j}(\theta)$（受控相位门）等基础门集。
  • OpenFermion: 关键在于其 molecular_data 和 circuits 模块。研究者使用了 OpenFermion 中的 Givens 旋转分解算法来制备 Slater 行列式初始态。

3.2 关键电路实现逻辑

  1. Jordan-Wigner 变换： 代码中需定义映射函数将 $(A/B, j, \uparrow/\downarrow)$ 映射到 $0 \dots 4N-1$ 的整数索引。注意 parity string 的处理：在 $H_0$ 的实现中，胞内和胞间跳跃涉及到相邻轨道，JW 变换后的 Pauli 串较短，有利于减少门数量。
  2. 绝热演化循环：
    # 伪代码逻辑示例
for l in range(L):
    # 演化 H0 部分 (SSH 跳跃)
    circuit.append(ssh_step_circuit(delta_t))
    # 演化 H1 部分 (Hubbard 相互作用)
    # 利用 Controlled-Phase 门实现 n_up * n_down 项
    current_u = (l / L) * total_U
    circuit.append(hubbard_step_circuit(delta_t, current_u))
  3. 测量与后处理： 由于位置算子 $\hat{X}$ 是对角的，复现时仅需在 $Z$ 基下进行采样。通过对采样得到的位串 $b^{(m)}$ 计算 $\exp(i 2\pi \hat{X}(b^{(m)}) / N)$ 的均值，即可得到 Berry 相估值。

3.3 开源资源链接

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献

  1. SSH 原型：Su, W. P., Schrieffer, J. R., & Heeger, A. J. (1979). Solitons in polyacetylene. Phys. Rev. Lett. [28]。该文奠定了 SSH 模型的物理基础。
  2. SSHH 扩展：Kivelson, S., & Heim, D. E. (1982). Hubbard versus peierls and the su-schrieffer-heeger model of polyacetylene. Phys. Rev. B [27]。引入相互作用的先驱工作。
  3. 多体 Berry 相定义：Resta, R., & Sorella, S. (1999). Electron localization in the insulating state. Phys. Rev. Lett. [36]。本文测量协议的理论支撑。
  4. 电路制备算法：Jiang, Z., et al. (2018). Quantum algorithms to simulate many-body physics of correlated fermions. Phys. Rev. Applied [40]。提供了 $O(N^2)$ 的态制备方案。

4.2 工作局限性评论

尽管本工作在理论和概念验证上非常出色，但仍存在以下局限：

  1. NISQ 硬件的可扩展性限制：研究采用了 $L=40$ 的 Trotter 步数。在当前的真实量子硬件上，如此深的电路（涉及数千个 CNOT 门）会因退相干和门噪声导致信号完全淹没。实验展示目前仍依赖于无噪声的经典仿真。
  2. 绝热时间与相变能隙：在接近热力学极限时，多体能隙可能以指数级减小。这意味着对于大规模系统，所需的绝热演化时间 $T$ 可能会超出量子比特的相干寿命，形成“绝热模拟的瓶颈”。
  3. Trotter 误差阶数：一阶 Trotter 分解虽然简单，但在处理非对易算符时误差较大。未来需要引入高阶分解或变分快速绝热演化（Variational Fast Adiabatic Evolution）来优化电路深度。
  4. 单一参数路径：目前仅探索了相互作用强度 $U$ 的线性增长路径。在复杂的相空间中，可能存在更优的非线性绝热路径以避开能隙闭合点。

5. 其他补充内容：量子化学视角下的 SSHH 模型

5.1 为什么量子化学家应该关注 SSHH？

SSHH 模型虽然起源于凝聚态物理，但它本质上是一个描述共轭 $\pi$ 电子系统的简化的半经验模型。对于研究导电聚合物、有机半导体以及分子导线的量子化学家来说，该模型捕捉了单电子跳跃（离域化）与库仑排斥（局域化）之间的本质竞争。

5.2 从平均场到强关联

本文在第 V 节提到，当相互作用较弱时，Hubbard 项可以用平均场近似（Mean-field）来处理。但在强关联极限下，平均场会失效，系统可能进入莫特绝缘体（Mott Insulator）相。本项工作展示的绝热量子模拟方案，其优越性恰恰在于它不依赖于平均场假设，能够直接探测真实的多体波函数。这为研究复杂的非绝热耦合、激子动力学等量子化学前沿问题提供了新的计算工具。

5.3 展望：通往容错量子计算的桥梁

随着纠错量子计算（FTQC）的到来，本文提出的协议可以无缝切换到更高阶的算法（如量子信号处理 QSP）。该研究不仅是一个“证明”，更是一份路线图，指导我们如何从测量简单的能量值，转向测量复杂的拓扑不变量。这种从“量”到“质”的模拟转变，正是未来量子化学模拟的核心方向。

总结： 通过将绝热算法与拓扑物理结合，这项研究不仅深化了我们对 SSHH 模型的理解，更展示了门阵列量子计算在处理复杂对称性及其破缺时的巨大潜力。对于科研工作者而言，掌握这一套从初始态制备到拓扑量提取的完整流水线，是进行多体量子模拟研究的必备技能。