来源论文: https://arxiv.org/abs/2607.02326v1 生成时间: Jul 03, 2026 15:52

强关联体系中的拓扑失效：单粒子格林函数零点表征多体拓扑的根本局限性深度解析

0. 执行摘要

拓扑能带理论在独立粒子体系（非相互作用体系）中取得了辉煌的成功。通过 Berry 相、Chern 数和缠绕数（Winding Number）等能带几何性质，研究者能够完美预测无能隙边界态的出现，这即是经典的“体-边界对应关系”（Bulk-Boundary Correspondence）。然而，在强关联多体系统中，由于多体基态不再能简单简化为单粒子 Slater 行列式，传统基于单粒子波函数的拓扑不变量诊断方法面临失效。

为了将拓扑分类延伸至强关联体系，近年来物理学界提出了一类基于单粒子格林函数（Green’s Function, GF）$G(\mathbf{k}, \omega)$ 的拓扑不变量。特别是在 Mott 绝缘体等强关联体系中，格林函数的“零点”（Zeroes，即满足 $\det G(\mathbf{k}, \omega_z) = 0$ 的频率点）被认为扮演着与“极点”（Poles）同等重要的拓扑角色，甚至被赋予了“涌现自旋子”（Emergent Spinons）等分数化准粒子的物理图像。人们普遍期望，格林函数零点的拓扑不变量能够完美捕捉强关联多体系统中的拓扑相变与边界态性质。

然而，最新由 Emile Pangburn 和 Olivier Gingras 撰写的研究成果《Fundamental limitations of single-particle Green’s-function zeroes as probes of many-body topology》给这一乐观期望泼了一盆冷水。该工作通过严格的解析推导与高精度数值计算（结合精确对角化 ED、密度矩阵重整化群 DMRG 以及辅助场量子蒙特卡洛 AFQMC），系统性地证明了：基于单粒子格林函数零点构建的拓扑不变量无法可靠地诊断相互作用多体系统的拓扑性质。在线性耦合的强关联多体模型（双通道 SSH 链与 2D BBH 模型）中，多体相互作用（如自旋-自旋 Heisenberg 耦合）能够在使多体基态彻底平庸化（Trivialized）的同时，不改变格林函数零点的拓扑不变量。这种失效的深层根源在于：单粒子格林函数在本质上只能探测粒子数改变（即不同 Fock 空间扇区之间）的电荷激发，而无法有效探测那些主导多体拓扑简并及自旋关联的电荷中性低能集体激发（如自旋翻转）。因此，格林函数零点并不对应物理上的拓扑准粒子，亦不能普适地刻画相互作用多体拓扑相。

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

1.1 核心科学问题

在强相互作用下，单粒子格林函数 $G$ 会因为谱权重的重新分配而出现零点。格林函数零点最初是作为拯救 Mott 绝缘体半满充填下 Luttinger 定理的一种数学机制被引入的。随后，人们尝试赋予这些零点以准粒子图像，认为它们对应着无电荷的自旋子激发（Spinons），可对热导率有贡献，并在拓扑 Mott 绝缘体中实现体-边界对应。本项工作核心要回答的问题是： 格林函数零点（以及由其定义的拓扑不变量）在强关联多体系统受到纯多体相互作用干扰时，是否依然能够忠实地反映多体基态的拓扑状态？

1.2 理论基础

1.2.1 格林函数的 Lehmann 表示与零点的产生机制

在费米子体系中，单粒子缓回格林函数在实频域的 Lehmann 表示为：

其中 $|\Psi_0\rangle$ 为 $N$ 粒子基态，能谱中的极点对应于添加或移除一个电子的多体激发能量（准粒子电荷激发）。对于格林函数的对角元 $G_{ii}(\omega)$，其单调性由其关于 $\omega$ 的导数决定：

由于 $\Re G_{ii}(\omega)$ 在任意两个相邻极点 $\omega_1$ 和 $\omega_2$ 之间是严格单调递减的，且在极点两侧其符号分别趋于 $+\infty$ 和 $-\infty$，因此根据介值定理，$\Re G_{ii}(\omega)$ 在区间 $](\omega_1, \omega_2)[$ 内必然且仅会穿过一次零点（参见原文中的 Fig. S1 示意图）。格林函数零点的存在是强关联体系中谱权重分裂（如 Mott 分裂形成上下 Hubbard 带）的直接结果。

1.2.2 零点哈密顿量 $\tilde{H}_0(\mathbf{k})$ 的引入

利用 Dyson 方程，格林函数可以写为自能 $\Sigma(\mathbf{k}, \omega)$ 的函数：

由此可见，格林函数 $G$ 的零点对应于自能 $\Sigma$ 的极点。在 Mott 绝缘极限下，自能可以近似展开为：

此处，$\tilde{H}_0(\mathbf{k})$ 即为表征零点色散的有效“零点哈密顿量”。它保留了原始单粒子哈密顿量 $H_0(\mathbf{k})$ 的所有对称性（如手征对称性、晶格对称性），但其参数经过了多体关联效应的重整化。

1.2.3 手征对称性保护下的零点拓扑缠绕数

对于具有子格点（手征）对称性的 1D 系统（如 SSH 链），满足 $\hat{S} H \hat{S}^{-1} = H$，其中子格点对称算符使得不同子格点上的产生消灭算符带上相反的相位。在这种对称性约束下，非相互作用的 Bloch 哈密顿量可以 spectrally flattened 并写为纯反角块形式：

由此定义拓扑缠绕数（Winding Number）：

零点哈密顿量 $\tilde{H}_0(\mathbf{k})$ 在手征对称性下也具有完全相同的反角块形式，因而可以定义零点的拓扑不变量 $\nu[\tilde{H}_0]$。在开放边界条件（OBC）下，若 $\nu[\tilde{H}_0] \neq 0$，则在零能处 $\omega = 0$ 会出现受拓扑保护的格林函数零点边界态（Zero-frequency boundary zeroes）。

1.3 技术难点与失效机制的物理本质

该研究的核心难点在于揭示为何零点拓扑不变量与实际的多体基态拓扑会发生脱钩。其深层物理机制在于Fock 空间扇区的不匹配性：

1. 单粒子格林函数的局限：格林函数 $G_{ij}(\omega)$ 的算符形式为 $c_i$ 和 $c_j^\dagger$。这意味着它测量的所有激发谱都是连接 $N$ 粒子基态与 $N \pm 1$ 粒子激发态的通道（即电荷激发扇区）。
2. 强关联体系的边界自由度：在 Mott 极限下（$U \gg t$），电荷波动被强烈抑制，系统边界处的低能自由度退化为纯自旋-$1/2$ 自由度（自旋态）。这些边界态之间的多体拓扑简并及低能激发属于电荷中性激发（例如自旋翻转激发 $S^+ = c^\dagger_\uparrow c_\downarrow$），其将 $N$ 粒子扇区内的状态连接到同一 $N$ 粒子扇区内的其他激发态。
3. 多体相互作用的干涉：当引入纯多体相互作用（如 Heisenberg 自旋耦合 $J_{12} \mathbf{S}_1 \cdot \mathbf{S}_2$）来耦合两个拓扑边界自旋时，这些自旋会两两配对形成一维自旋单态（Singlets）。这使多体基态彻底平庸化（退简并，打开物理自旋能隙 $\Delta_s > 0$）。然而，由于自旋-自旋相互作用不涉及电荷转移，它在格林函数所探测的 $N \pm 1$ 粒子电荷激发扇区中几乎不产生一阶修正。这就导致格林函数零点边界态对这种多体平庸化过程“毫无察觉”，依旧在零能处保持无能隙（Gapless）状态。这种脱钩现象直击了单粒子格林函数零点理论的盲区。

1.4 方法细节：哈伯德算符（Hubbard Operator, HO）方法与 Roth 解耦

为了精确求解强关联系统中的自能极点和零点哈密顿量 $\tilde{H}_0(\mathbf{k})$，论文采用了哈伯德算符方法（HO Method）结合 Roth 解耦方案。该方法的基本思路是将局部物理空间中的状态完全基底化，定义局域哈伯德算符 $\Psi_i^{\mathrm{loc}}$，使得运动方程在局域项下是闭合的：

而对于非局域的跃迁项 $H^t_{ij}$，运动方程会产生非闭合的高阶算符。HO 方法的核心近似在于将这些高阶算符线性投影回原算符基底：

通过构造自共轭的重叠矩阵（Overlap Matrix） $I_{ij,nm} = \langle \{ \Psi_{i,n}^{\mathrm{loc}}, \Psi_{j,m}^{\mathrm{loc}\dagger} \} \rangle$ 和动力学矩阵 $M_{ij,nm} = \langle \{ [\Psi_{i,n}^{\mathrm{loc}}, H], \Psi_{j,m}^{\mathrm{loc}\dagger} \} \rangle$，可以自洽地求解格林函数色散。通过 Roth 解耦方案，可以精确计算出包含两点及四点关联函数（如 $\langle S_i^- S_j^+ \rangle$ 和双占有率 $\langle n_{i\uparrow} n_{i\downarrow} \rangle$）的自洽动力学。对于半满充填的单轨道 Hubbard 模型，利用 HO 方法可以直接求得零点哈密顿量为：

其中 $p = \langle n_{i\sigma} n_{j\sigma} \rangle + \langle S_i^- S_j^+ \rangle - \langle \Delta_i \Delta_j^\dagger \rangle$ 描述了相邻键上的关联强度。该理论方法为零点的解析分析提供了强大的支撑。

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

论文中设计了两个具有里程碑意义的强关联多体 Benchmark 体系，用以无可辩驳地证明零点拓扑不变量的失效。

2.1 Benchmark 体系 1：双通道强关联耦合 SSH 链 (Two coupled SSH chains)

2.1.1 体系哈密顿量

考虑两条平行的、具有 Hubbard $U$ 相互作用的 Su-Schrieffer-Heeger (SSH) 链（标记为 $l = 1, 2$），单链哈密顿量为：

其中跃迁项呈现交替模式：$t_i = t + (-1)^i \delta t$。当 $\delta t > 0$ 时，系统处于拓扑平庸相（图 1a）；当 $\delta t < 0$ 时，系统处于拓扑非平庸相（图 1b），此时在开放边界（OBC）下，其两端会各留出一个未配对的格点，相当于两个 Hubbard 孤立点，在强关联下形成自旋-$1/2$ 边界态。

接下来，将这两条处于拓扑相的链进行耦合（图 1c），耦合项包含单粒子通道层间跃迁 $t_{12}$ 和纯多体通道自旋-自旋相互作用 $J_{12}$：

2.1.2 关键计算数据分析

研究者利用 DMRG 计算了长度为 $L=24$ 的耦合链体系在不同耦合强度下的物理自旋能隙 $\Delta_s$；同时利用精确对角化（ED, $L=6$）和自洽哈伯德算符方法（HO）计算了格林函数零点的空间展宽量 $Z[G](\omega) = [\min |\lambda_G(\omega)|]^{-1}$。计算采用的物理参数为：$t = 0.75$，$\delta t = 0.25$，$U = 8$。

• 物理可观测量的变化（DMRG 结果，图 2a）：
  • 当单独引入层间跃迁 $t_{12}$ 或自旋相互作用 $J_{12}$ 时，自旋能隙 $\Delta_s$ 均迅速打开（对于 $J_{12} = 0.25$，$\Delta_s \approx 0.12$；对于 $t_{12} = 0.25$，$\Delta_s \approx 0.28$）。自旋能隙的打开意味着边界处的非平庸双重简并自旋单态配对成了自旋单态，系统在多体基态层面上已经彻底平庸化。
• 格林函数零点响应（ED 与 HO 结果，图 2b）：
  • 情况 A：单粒子跃迁 $t_{12} = 0.25, J_{12} = 0$（绿色虚线）。由于 $t_{12}$ 直接耦合单粒子电荷算符，它打破了子格点（手征）对称性。格林函数零点不变量成功捕捉到了这一变化，零点能隙迅速打开（$\omega = 0$ 处 $Z[G]$ 降为 0），拓扑不变量失效，零点与多体拓扑同步走向平庸。
  • 情况 B：多体相互作用 $J_{12} = 0.25, t_{12} = 0$（橘色点线）。尽管在物理上边界自旋已经配对、基态已完全平庸化，但在 $\omega = 0$ 处，格林函数零点依然坚挺地存在，无能隙特征完好保留，且其拓扑缠绕数 $\nu[\tilde{H}_0] = 2$ 保持恒定。这明确展示了格林函数零点在面对纯多体相互作用干扰时给出了完全错误的拓扑诊断结果。

2.2 Benchmark 体系 2：二维高阶拓扑 BBH 绝缘体中的零点失效

为了证明该局限性在更高维度下的普适性，研究者将研究拓展至二维 Benalcazar-Bernevig-Hughes (BBH) 四极矩拓扑绝缘体模型。

2.2.1 体系模型与相互作用

二维 BBH 模型在非相互作用极限下是一个典型的高阶拓扑绝缘体（HOTI），在开放边界下其四个角落会各出现一个零维的零能电荷模式（角落态，Corner States）。引入原位 Hubbard 相互作用 $U$ 后，在半满 Mott 绝缘体状态下，角落的电荷自由度被冻结，退化为四个局域自旋-$1/2$ 自由度。论文考虑将两层 BBH 体系通过层间自旋-自旋相互作用 $J_{12}$ 进行耦合。

2.2.2 高精度数值计算结果 (AFQMC vs. HO)

由于二维量子多体计算的复杂性，研究者使用了辅助场量子蒙特卡洛（AFQMC）来求解高精度的多体基态，并结合自洽 HO 方法求解零点行为（计算尺寸 $18 \times 18$，参数：$\gamma = 0.5, \lambda = 1.0$）：

• AFQMC 能量与能隙数据（图 S12）：随着 $U$ 的增加，拓扑相的电荷能隙 $\Delta_c$ 立即打开（当 $U=4$ 时，$\Delta_c \approx 2.5$），而自旋能隙 $\Delta_s$ 保持为 0（因角落自旋尚未耦合）。当引入层间自旋相互作用 $J_{12}$ 后，四个角落的自旋由于两两配对形成自旋单态，其自旋激发能隙 $\Delta_s$ 立即打开，系统整体的多体拓扑特性被彻底平庸化。
• 格林函数特征：与 1D 链的情况完全平行，层间多体自旋相互作用 $J_{12}$ 无法破坏零能处的格林函数零点。零点不变量继续错误地指示系统处于高阶拓扑状态。这表明格林函数零点拓扑的失效性与系统维度无关，具有普遍的物理本质。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

为了让量子化学与凝聚态物理计算领域的科研人员快速复现论文中的核心结论，以下提供基于开源社区最主流工具链的复现指南。

3.1 推荐开源软件包与计算工具

1. ITensor (Julia 语言版本)：用于高精度 DMRG 计算。Julia 版 ITensor 拥有极佳的易用性与卓越的张量网络收缩性能。
   • GitHub 链接：https://github.com/ITensor/ITensors.jl
2. SAFIRE：用于辅助场量子蒙特卡洛（AFQMC）的高效并行计算框架。
   • 官方链接：https://safire.flatironinstitute.org/

3.2 1D 耦合 SSH 链复现指南 (基于 ITensor Julia)

第一步：环境配置与依赖包导入

首先确保安装了 Julia（建议 >= 1.9），并在 REPL 中执行：

using Pkg
Pkg.add("ITensors")

第二步：定义哈密顿量与 DMRG 求解

以下是复现双通道耦合 SSH 链基态能及自旋能隙的核心 Julia 脚本：

using ITensors

function run_coupled_ssh_dmrg(L::Int, U::Float64, J12::Float64, t12::Float64; t=0.75, dt=0.25)
    # 总格点数：2条链，每条链L个格点，共2L个格点
    # 格点编号：1 到 L 属于链1，L+1 到 2L 属于链2
    sites = siteinds("Electron", 2*L; conserve_qns=true)
    
    # 构建 MPO 形式的 Hamiltonian
    ampo = OpSum()
    
    # 1. 链 1 内部跃迁与 Hubbard U
    for i in 1:(L-1)
        ti = t + (-1)^i * dt
        # 跃迁项
        ampo += -ti, "Cdagup", i, "Cup", i+1
        ampo += -ti, "Cdagup", i+1, "Cup", i
        ampo += -ti, "Cdagdn", i, "Cdn", i+1
        ampo += -ti, "Cdagdn", i+1, "Cdn", i
    end
    for i in 1:L
        ampo += U, "Nupdn", i
    end
    
    # 2. 链 2 内部跃迁与 Hubbard U
    for i in 1:(L-1)
        ti = t + (-1)^i * dt
        idx = i + L
        # 跃迁项
        ampo += -ti, "Cdagup", idx, "Cup", idx+1
        ampo += -ti, "Cdagup", idx+1, "Cup", idx
        ampo += -ti, "Cdagdn", idx, "Cdn", idx+1
        ampo += -ti, "Cdagdn", idx+1, "Cdn", idx
    end
    for i in 1:L
        ampo += U, "Nupdn", i+L
    end
    
    # 3. 两链之间的耦合项
    for i in 1:L
        # 单粒子杂化项 t12
        if t12 != 0.0
            ampo += -t12, "Cdagup", i, "Cup", i+L
            ampo += -t12, "Cdagup", i+L, "Cup", i
            ampo += -t12, "Cdagdn", i, "Cdn", i+L
            ampo += -t12, "Cdagdn", i+L, "Cdn", i
        end
        
        # 多体自旋耦合项 J12
        if J12 != 0.0
            # S.S = SzSz + 0.5*(S+S- + S-S+)
            ampo += J12, "Sz", i, "Sz", i+L
            ampo += 0.5*J12, "S+", i, "S-", i+L
            ampo += 0.5*J12, "S-", i, "S+", i+L
        end
    end
    
    H = MPO(ampo, sites)
    
    # 定义半满填充状态下的基态（N_up = L, N_dn = L）
    state = [isodd(n) ? "Up" : "Dn" for n in 1:(2*L)]
    psi0 = randomMPS(sites, state; linkdims=10)
    
    # 设置 DMRG 扫描参数
    sweeps = Sweeps(6)
    setmaxdim!(sweeps, 50, 100, 200, 400, 800, 1024)
    setcutoff!(sweeps, 1e-10)
    
    energy, psi = dmrg(H, psi0, sweeps; silent=true)
    return energy, psi, sites
end

# 运行基态能量计算
L = 24
U = 8.0
energy, psi, sites = run_coupled_ssh_dmrg(L, U, 0.25, 0.0)
println("System Ground State Energy: ", energy)

第三步：计算自旋能隙 $\Delta_s$

自旋能隙定义为 $\Delta_s = E_0(N_\uparrow+1, N_\downarrow-1) - E_0(N_\uparrow, N_\downarrow)$。可以通过调整基态搜索时的总自旋量子数来实现。在上述代码的基础上，构建一个处于激发自旋扇区（例如总自旋为 1）的随机 MPS 作为初始状态进行 DMRG，所得能量差值即为自旋能隙。

3.3 零点的 Lanczos 精确对角化（ED）求解步骤

由于需要计算零频率处的 $\det G(\omega=0) = 0$，对于 $L=6$ 的小体系，可以采用 Lanczos 方法直接计算 Lehmann 表示中的格林函数。其核心步骤如下：

1. 求解基态：用 Lanczos 求解 $N$ 粒子基态 $|\Psi_0\rangle$ 与能量 $E_0^N$。
2. 应用消灭算符：构造初态 $|\phi_{i,\sigma}\rangle = c_{i,\sigma} |\Psi_0\rangle$，在 $N-1$ 粒子扇区使用 Lanczos 算法产生克里洛夫（Krylov）子空间，进而求出 $G_{ii}(\omega)$ 的极点和留数。
3. 应用产生算符：构造初态 $|\psi_{i,\sigma}\rangle = c^{\dagger}_{i,\sigma} |\Psi_0\rangle$，在 $N+1$ 粒子扇区进行同样的操作。
4. 零点定位：通过对角元 $G_{ii}(\omega)$ 的零位交叉点确认零点频率，并验证行列式 $\det G(\omega=0) = 0$ 的零特异值行为。

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

4.1 关键引用文献

本项研究建立在多体物理、能带拓扑和强关联自洽方法等多维领域的基石之上。以下为其核心引用文献：

1. W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. Lett. 42, 1698 (1979)：著名的 SSH 模型，奠定了一维拓扑绝缘体几何与电荷孤立子特性的基石。
2. W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, Science 357, 61 (2017)：定义了高阶四极矩拓扑绝缘体（BBH模型），将体-边界对应关系推广到多极矩矩边界态。
3. I. Dzyaloshinskii, Phys. Rev. B 68, 085113 (2003)：阐述了非费米液体与 Mott 绝缘体中，由于关联效应产生的 Luttinger 费米面重构与格林函数零点存在的必然性。
4. L. Fidkowski and A. Kitaev, Phys. Rev. B 81, 134509 (2010)：指出相互作用会将一维拓扑超导和费米子拓扑分类从 $\mathbb{Z}$ 简化至 $\mathbb{Z}_8$（而在带有手征对称的一维费米子体系中，多体作用将拓扑分类重整为 $\mathbb{Z}_4$ ）。

4.2 对本项工作局限性的客观评论

尽管该工作对强关联多体物理中的拓扑不变量应用给出了极具震撼力的修正，但作为严谨的科学讨论，我们依然应当看到其局限性所在：

1. 分析模型的高度对称性（粒子-空穴对称性）限制： 为了避免先前研究所指出的“格林函数不变量随化学势改变而产生的非物理依赖性”（参见文献 J. Zhao et al., PRL 131, 106601），作者刻意将模型限制在严格的粒子-空穴对称性（Particle-Hole Symmetry）下。这使得结论极具数学纯洁性。然而，真实的物理材料（如过渡金属氧化物、重费米子材料等）往往普遍伴随有强烈的粒子-空穴非对称性（例如次近邻跃迁项 $t'$ 的存在）。在缺乏这种人工对称性约束的现实体系中，零点不变量的失效方式可能比本文阐述的还要复杂，这也削弱了本结论直接指导具体多轨道材料计算的便利度。

2. 哈伯德算符（HO）方法的非变分性缺陷： 在对二维 BBH 体系进行求解时，由于量子蒙特卡洛在远离半满时面临严重的负符号问题（Sign Problem），作者只能大量依赖 HO 方法辅以 Roth 自洽解耦。需要指出的是，HO 方法在本质上是非变分的，其求得的基态能量可能低于精确值。在强关联区（$U/t \gg 1$），HO 方法中的顺磁近似完全忽略了可能涌现的反铁磁长程磁序（Antiferromagnetic Ordering）。这种物理相的遗漏可能导致对高维系统实际多体相图界面的微小偏差。

3. 缺乏替代性不变量的统一计算框架： 文章在结尾处提出，要正确诊断强关联 Mott 绝缘体中的多体拓扑，应该引入多粒子格林函数（如动态自旋易受率 $\chi(\mathbf{q}, \omega)$ ）。然而，多粒子格林函数的数值计算成本呈指数级增长，目前在计算物理界缺乏类似单粒子格林函数的成熟且通用的计算工具箱。文章指出了一条正确的路线，但未能给研究者提供一个在实际复杂材料计算中低成本、高效率的替代方案。

5. 其他必要的补充：格林函数零点的历史脉络与凝聚态物理启示

为了帮助读者更全面地理解这一研究的来龙去脉，我们对格林函数零点在凝聚态物理中的角色进行深度的历史补充。

5.1 从 Luttinger 定理说起

在费米液体理论中，Luttinger 定理是一个极具分量的结论：它指出，费米海所包围的动量空间体积在相互作用下保持不变，严格正比于系统的电子密度。其背后的数学构造是 Luttinger-Ward 泛函 $Y[G]$。然而，在以铜氧化物高温超导体和 Mott 绝缘体为代表的强关联体系中，实验观测到的费米面体积似乎违背了这一经典定理。 为了“拯救” Luttinger 定理，Dzyaloshinskii 证明了，不仅格林函数的极点（代表物理费米子）可以贡献费米面体积，格林函数的零点（满足 $G=0$ 的亚表面，称为 Luttinger 表面）同样可以在数学上等效地计入 Luttinger 体积的积分中。这就为 Mott 绝缘体提供了一个优雅的数学描述：一个半满的 Mott 绝缘体之所以不导电，是因为其费米面被一圈“零点表面”所占据，电荷激发在该表面上完全被自能的发散所压制。

5.2 铜氧化物超导体中的“假想自旋子”

在高温超导体的伪能隙（Pseudogap）研究中，科学家们观测到了所谓的“费米弧”（Fermi Arcs）现象。许多理论学家（如 M. Fabrizio, S. Sachdev 等）提出，格林函数零点的色散 $\tilde{H}_0(\mathbf{k})$ 实际上代表了量子自旋液体中涌现的分数化准粒子——自旋子（Spinons）。由于自旋子仅携带自旋而不携带电荷，它们无法贡献电导，但这在格林函数零点的数学结构上表现得极为协调。这也是格林函数零点拓扑学说之所以广受追捧的原因：它承诺了一条从单粒子能带结构直接通往分数化多体激发的“捷径”。

5.3 本项工作的终极启示

然而，这篇论文彻底厘清了这一科学幻觉：格林函数零点仅仅是电子作为基本激发单元在强相互作用下发生崩溃（Breakdown）的数学信号，而不代表真正的低能分数化物理准粒子。

格林函数零点哈密顿量虽然在数学上形式优美，且在非相互作用极限下完美契合拓扑性质，但一旦系统在强关联下发生电荷与自旋的彻底分离（Spin-Charge Separation），多体系统的拓扑保真度便由电荷中性的自旋相干性主导。而由于单粒子格林函数天生带有电荷特征（以电子消灭和产生为探针），它被死死地禁锢在电荷激发 Fock 扇区中，从而对自旋多体配对产生的拓扑平庸化完全“失明”。

这一结论对未来量子化学计算、强关联电子材料设计具有深远的指导意义：在探索新型拓扑材料（如拓扑 Mott 绝缘体、近藤拓扑绝缘体）时，我们不能再盲目套用基于单粒子自能或单粒子格林函数零点计算得到的拓扑不变量（如缠绕数、Chern 数等）。在关联电子计算中，发展计算成本可行、可直接反映多体关联相关性的多体本征泛函和多粒子响应函数，才是走向下一个强关联拓扑材料时代的必然选择。