来源论文: https://arxiv.org/abs/2602.21199v1 生成时间: Feb 26, 2026 11:10

拓扑 Floquet 格林函数零点：理论、方法与量子模拟实现

0. 执行摘要

传统的拓扑物态分类通常基于单粒子 Hamiltonian 的能带结构（极点，Poles）。然而，在强关联电子系统中，相互作用可以导致格林函数（Green’s Function）出现零点（Zeros），这些零点在拓扑性质的定义中起着与极点同等重要甚至更核心的作用。本文解析的论文《Topological Floquet Green’s function zeros》将这一前沿概念扩展到了非平衡态的 Floquet 系统。作者证明了在 Floquet-BDI 对称类中，格林函数零点不仅存在，而且即使在无相互作用的情况下也会因驱动特性而出现。通过引入基于格林函数的拓扑不变量 $N_1$ 和 $ ilde{N}_1$，作者建立了 Floquet 系统的体-边对应关系，并详细阐述了对称质量产生（Symmetric Mass Generation, SMG）如何使拓扑边界态消失，转而代之以拓扑边界零点。最后，文章给出了在 NISQ 量子硬件上利用 Jordan-Wigner 变换和受限 Toffoli 门实现该物理过程的路线图。

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

1.1 核心科学问题：超越极点的拓扑观

在平衡态强关联系统中，格林函数零点通常与 Mott 物理、Luttinger 体积丢失以及非费米液体行为相关。最近的研究发现，这些零点可以形成“零点带”，并具有拓扑荷。本研究的核心问题是：在周期性驱动的非平衡 Floquet 系统中，格林函数零点如何演化？它们是否能像在平衡态那样通过对称质量产生（SMG）来消除极点并维持拓扑稳定性？

1.2 理论基础：Floquet-BDI 对称性与格林函数定义

研究集中在 1D Majorana 链（Kitaev-like 链），其具有费米子宇称对称性 $\hat{P}$ 和时间反演对称性 $\hat{T}$，且 $\hat{T}^2 = 1$，属于 BDI 对称类。

离散时间格林函数定义： 不同于连续时间，Floquet 系统在离散时间 $n \in \mathbb{Z}$ 下定义格林函数：

$$G_{R/A}^{AB}(n) = \mp i heta(\pm n) \langle \{ \gamma_A(n), \gamma_B(0) \} angle$$

其中 $\gamma_A(n) = \hat{U}_F^{-n} \gamma_A \hat{U}_F^n$。关键在于作者定义了两个不同的 Floquet 演化算符 $\hat{U}_F$ 和 $\hat{ ilde{U}}_F$，分别对应从 $t=0$ 和 $t=T_F/2$ 开始的微运动（micromotion），这对于捕捉 Floquet 拓扑的完整信息（包含 $0$ 和 $\pi$ 模式）至关重要。

1.3 技术难点：Floquet 特有的零点

在平衡态自由费米子系统中，格林函数行列式永远不为零（只有极点）。但在 Floquet 系统中，由于能量（准能量）的周期性，在扩展布里渊区方案中，极点之间必然夹杂着零点。作者指出，这意味着在 Floquet 拓扑中，零点是内禀的，必须将其与相互作用诱导的零点区分开来。

1.4 方法细节：拓扑不变量的推广

作者推广了平衡态的拓扑绕数（Winding Number）公式，引入了两个基于格林函数的不变量：

$$N_1 = \int_{-\pi}^{\pi} rac{dp}{4\pi i} ext{tr}[\Gamma G_R^{-1} \partial_p G_R]_{\Omega=0}$$

$$ ilde{N}_1 = \int_{-\pi}^{\pi} rac{dp}{4\pi i} ext{tr}[\Gamma ilde{G}_R^{-1} \partial_p ilde{G}_R]_{\Omega=0}$$

其中 $\Gamma$ 是手征对称性算符。这两个不变量共同决定了 Majorana 零能模（MZMs）和 $\pi$ 能模（M$\pi$Ms）的数量：$ u_0 = (N_1 + ilde{N}_1)/2, u_\pi = (N_1 - ilde{N}_1)/2$。

2. 关键 Benchmark 体系与数据分析

2.1 自由费米子相图 (Benchmark)

作者首先建立了 $N_f=1$ 的无相互作用 Floquet-Kitaev 链作为基准。通过调节化学势参数 $g$ 和耦合参数 $J$，系统展现出四个不同的相：$(0,0), (1,0), (0,-1), (1,-1)$。图 2 清晰地展示了极点和零点的能带结构。在相变点，不仅极点带会闭合（Gap closing of poles），零点带也会闭合（Gap closing of zeros），这说明零点带的变化同样可以触发拓扑相变。

2.2 相互作用诱导的对称质量产生 (SMG)

作者引入了 Fidkowski-Kitaev (FK) 型相互作用。FK 相互作用是一种特殊的 8 费米子项，它能在不破缺对称性的情况下，将 $N_f=8$ 的拓扑边模转化为平凡相。在 Floquet 设定下，作者通过微扰论计算了 $N_f=8$（对于 $0$ 或 $\pi$ 模）和 $N_f=4$（同时具有 $0$ 和 $\pi$ 模）的情况。

关键数据点 (来自 Fig 4 & 6)：

  • 边界格林函数： 在强相互作用极限下，边界 MZMs 的格林函数 $G_R^{(0)}(\Omega)$ 在 $\Omega=0$ 处展现出清晰的零点。当相互作用强度 $w > 0$ 时，原本在 $\Omega=0$ 处的极点消失，谱权重被转移到高能区，而零点被稳固地钉扎在 $\Omega=0$ 和 $\Omega=\pi$。
  • 体格林函数： 计算得到的体格林函数 $G_R(\Omega, p)$ 在小 $w$ 极限下表现为： $$G_R(\Omega, p) \simeq rac{1}{2} rac{\sin(\Omega) + rac{\pi}{2} \cos(7\pi w) h(p)}{\cos(7\pi w) - \cos(\Omega^+)}$$ 这证明了相互作用诱导的零点带在动量空间中具有与极点相反的拓扑绕数，从而在整体上贡献了拓扑不变量。

2.3 性能数据：谱权重转移

在图 4 中，作者对比了无限温度和零温度下的格林函数行为。结果显示，零点不仅是基态属性，在有限温度下依然作为格林函数实部改变符号的特征点存在。这增强了该现象在 NISQ 设备（通常具有显著热噪声）中被观测到的可能性。

3. 代码实现细节与复现指南

3.1 物理体系映射：Jordan-Wigner 变换

为了在量子芯片上复现，必须将 Majorana 算符映射为 Pauli 算符：

  • 定义局部费米子宇称：$X_{j,\alpha} = i \gamma_{j,L}^\alpha \gamma_{j,R}^\alpha$
  • 费米子产生/湮灭算符对应的串：$Z_{j,\alpha}, Y_{j,\alpha}$ 作者推荐使用 bottom-to-top 的 JW 串排序，以优化电路的局域性。

3.2 电路构建块

Floquet 循环由三个主要酉算符组成：$\hat{U}_F = \hat{U}_{g/2} \hat{U}_{J/2} \hat{U}_w \hat{U}_{J/2} \hat{U}_{g/2}$。

  • 单比特门： $\hat{U}_g$ 对应旋转 $R_x$。
  • 双比特门： $\hat{U}_J$ 对应 $ZZ$ 耦合（可通过两个 CNOT 和一个 $R_z$ 实现）。
  • 相互作用门 $\hat{U}_w$： 这是实现 SMG 的难点。作者提出 $\hat{U}_w$ 可以分解为两部分：
    1. $\hat{U}_1$：一系列 $fSim$ 门，用于处理 Pa/Pb 项。
    2. $\hat{U}_2$：广义 Toffoli 门（CCCNOT），用于实现 4 个费米子相互作用项。在 $w=1/8$ 等特殊点，这可以简化为标准受控非门电路。

3.3 复现流程与开源工具建议

  1. 初始化： 将量子比特初始化为随机比特串（模拟无限温度平均）。
  2. 循环执行： 应用 $n$ 次 Floquet 门序列。
  3. 测量： 测量边界比特的自相关函数 $\langle Z_{1,L}(n) Z_{1,L}(0) angle$。
  4. 数据处理： 对时间序列进行离散傅里叶变换，提取 $G_R(\Omega)$，寻找符号改变点（零点）。

推荐工具链：

  • Cirq / Qiskit: 用于构建门级电路。
  • Fermionic Quantum Emulator (OpenFermion): 用于验证 JW 映射的正确性。
  • 推荐 Repo 参考： 读者可参考 Google Quantum AI 在 Nature 612 (2022) 中关于 Floquet 边模实现的电路结构进行修改。

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

4.1 关键引用

  • Fidkowski & Kitaev (2010): 定义了 $\mathbb{Z}_8$ 分类，是本文相互作用设计的基石。
  • Gurarie (2011): 奠定了格林函数拓扑不变量的理论框架。
  • Google Quantum AI (2022/2024): 提供了在 NISQ 硬件上实现类似 Floquet 链的实验先例。

4.2 局限性评论

  1. 维度限制： 本文仅讨论 1D 链。在 2D 或 3D 系统中，格林函数零点可能形成复杂的闭合面（Luttinger surfaces），其拓扑性质的 Floquet 推广将极度复杂且难以解析处理。
  2. 微扰论依赖： 文章中关于体格林函数的结果主要基于一阶 Floquet 微扰论。在强驱动强相互作用的非微扰区，格林函数零点的稳定性尚待数值验证（如使用 Floquet-DMRG）。
  3. NISQ 噪声： 虽然自相关函数在 $T o \infty$ 时也显示零点，但量子门的退相干会模糊这些精细的谱特征，特别是零点（格林函数实部为零处）对相位噪声极为敏感。

5. 补充：量子化学视角下的“零点”意义

对于量子化学背景的读者，理解本文的一个捷径是将其类比为强关联分子中的自然轨道占据数。在简单的单参考方法（如 Hartree-Fock）中，轨道要么全占（极点），要么全空。而在处理具有多参考特性的过渡金属配合物或化学键断裂时，相互作用会导致占据数偏离 0 或 1。格林函数零点在此物理意义上，标志着传统的单粒子图像（准粒子描述）彻底失效。本文的贡献在于，它证明了即使在动力的、受驱动的过程中，这些“准粒子失效”的点依然遵循严谨的拓扑规律。这为未来设计基于量子拓扑保护的动态化学模拟提供了理论支撑。例如，我们是否可以利用 Floquet 驱动下的格林函数零点来模拟特定的催化循环中由于强关联效应导致的电子转移受阻？这是一个非常值得探索的方向。