Symmetry-enforced agreement of Kohn--Sham and many-body Berry phases in the SSH--Hubbard chain

本文通过一维 SSH-Hubbard 链的密度矩阵重整化群计算发现，尽管相互作用强度变化会显著改变多体波函数的几何响应，但在反演对称的绝缘体区域，对称性强制使得 Kohn-Sham 描述与多体 Berry 相位在弱至强耦合范围内始终保持一致，这种一致性源于$\mathbb{Z}_2$对称性导致的扇区匹配而非密度对 Berry 连接的编码。

要点

这篇论文探讨了一个非常深奥的量子物理问题，但我们可以用一些生活中的比喻来把它讲得通俗易懂。

核心故事：两个“双胞胎”的旅行

想象一下，你有一群非常调皮的电子（就像一群在操场上乱跑的孩子），它们之间会互相推搡（这就是相互作用，论文里用 $U$ 表示）。

科学家想研究这群孩子在一个环形跑道上跑一圈时，会留下什么样的“足迹”或“记忆”。这个“记忆”在物理学里叫贝里相位（Berry Phase）。它不是孩子跑了多远，而是他们在这个过程中“转了多大的弯”或者“心态发生了怎样的微妙变化”。

这篇论文主要做了两件事：

1. 真实世界（多体系统）： 科学家真的去数这群互相推搡的孩子，看他们跑一圈后的真实“记忆”。
2. 简化模型（Kohn-Sham 描述）： 科学家试图用一个更简单的模型来模拟这群孩子。在这个简化模型里，孩子们互不推搡（没有相互作用），只是乖乖地跑。但是，这个简化模型被强制要求：孩子们跑完一圈后，留在跑道上的“人数分布”（密度）必须和真实世界一模一样。

论文发现了什么惊人的巧合？

通常大家认为，如果孩子们互相推搡（强相互作用），他们的“记忆”（贝里相位）会变得非常复杂，而那个“互不推搡”的简化模型肯定算不准。

但是，这篇论文发现了一个惊人的现象：

在这个特定的“苏 - 施里弗 - 海格 - 哈伯德（SSH-Hubbard）”模型中，无论孩子们推搡得有多厉害（从完全不推搡到疯狂推搡），那个简化模型算出来的“记忆”（贝里相位），竟然和真实世界完全一致！

这就好比你让一群乱跑的孩子和一群排好队走路的孩子，虽然过程完全不同，但最后他们都在墙上留下了完全一样的签名。

为什么会有这个巧合？（关键揭秘）

科学家一开始以为，这是因为“人数分布”（密度）里藏满了所有关于“记忆”的秘密。也就是说，只要简化模型能模仿出真实的人数分布，它就能自动算出正确的“记忆”。

但论文打碎了这种幻想！

通过精密的计算，科学家发现：

1. 人数分布是“死”的： 无论孩子们怎么推搡，也无论那个“环形跑道”怎么旋转（引入磁通量 $\theta$），孩子们在每个位置的人数分布完全不变！就像一群人在操场上，不管怎么乱跑，每个角落站的人数看起来都一样。
2. 真正的“记忆”藏在别处： 虽然人数没变，但孩子们内心的状态（波函数）发生了巨大的变化。在推搡厉害的时候，孩子们变得非常“僵硬”，不敢乱动（电荷涨落被冻结），这导致他们在这个过程中的几何敏感度（量子度规）发生了剧烈变化。

结论来了：
既然人数分布根本没变，那简化模型（Kohn-Sham）根本不需要费脑子去模仿复杂的相互作用，它只需要保持最简单的状态（就像 $U=0$ 时那样）就能满足“人数分布”的要求。

那么，为什么简化模型算出的“记忆”还是对的呢？
不是因为“人数分布”告诉它答案，而是因为“对称性”强迫它必须选那个答案。

一个生动的比喻：对称的迷宫

想象一个完全对称的迷宫（具有反演对称性）。

• 真实世界（强相互作用）： 迷宫里充满了陷阱和机关，你走起来很困难，路径很曲折。
• 简化模型（无相互作用）： 迷宫里空空荡荡，你走起来很轻松，路径很直。

但是，这个迷宫有一个铁律：无论你怎么走，只要不走出迷宫的边界（能隙不关闭），你最终从起点绕一圈回到原点时，你要么转了 0 度，要么转了 180 度（$\pi$）。 这就是所谓的 $Z_2$ 拓扑保护。

• 因为迷宫是对称的，你不可能转 90 度或 120 度。
• 所以，不管你是走“困难模式”（真实世界）还是“简单模式”（简化模型），只要你还在这个迷宫里，你最后转的角度只能是 0 或 180 度。

这篇论文的结论就是：
简化模型之所以能算对，不是因为它“读懂”了复杂世界的密度，而是因为对称性的大棒强迫它只能选那个唯一的答案。这是一种**“被迫的一致”**，而不是因为密度里藏着所有秘密。

这篇论文告诉我们什么大道理？

1. 不要迷信“密度”： 在量子世界里，仅仅知道“哪里有多少人”（密度），并不足以告诉你“这群人经历了什么复杂的心理变化”（几何相位）。密度是“死”的，而波函数是“活”的。
2. 对称性是强大的： 在某些特定的对称保护下，即使物理过程天差地别（一个乱跑，一个排队），最终的宏观结果（拓扑相位）却可能完全一样。
3. 简化模型的局限性： 虽然在这个特定模型里，简化模型碰巧算对了，但这只是因为它被对称性“绑架”了。如果去掉对称性，或者改变模型，简化模型可能就会算错。这提醒我们，在计算复杂材料时，不能盲目依赖简单的密度泛函理论（DFT）来预测所有拓扑性质。

一句话总结：
这篇论文发现，在一个特定的量子迷宫里，不管里面的粒子是“乱成一团”还是“整齐划一”，只要迷宫是对称的，它们最后留下的“签名”就必然一样。但这并不是因为“整齐划一”的模型看懂了“乱成一团”的内心，而是因为对称性规定它们只能签同一个名字。

技术摘要

这是一份关于论文《Symmetry-enforced agreement of Kohn–Sham and many-body Berry phases in the SSH–Hubbard chain》（SSH-Hubbard 链中 Kohn-Sham 与多体 Berry 相的对称性强制一致性）的详细技术总结。

1. 研究背景与问题 (Problem)

密度泛函理论（DFT）中的 Kohn-Sham (KS) 形式通过将相互作用的多电子问题映射为非相互作用的 KS 有效系统，能够精确重现基态电子密度。然而，在强关联电子系统中，几何相位（如 Berry 相）和拓扑不变量是多体波函数的泛函，而不仅仅是密度的泛函。

核心问题：
在强关联区域，一个仅通过匹配基态密度（density-matching）构建的 KS 描述，能否重现多体波函数中的 Berry 相？

• 通常认为，几何相位并不唯一地由单粒子期望值或密度决定。
• 如果密度不编码几何信息，KS 描述为何有时能重现正确的拓扑相？这种一致性是源于密度本身包含了拓扑信息，还是源于其他机制（如对称性）？

2. 研究模型与方法 (Methodology)

模型系统：

• 研究了一维 SSH-Hubbard 模型（二聚化 Fermi-Hubbard 链），置于环上（周期性边界条件，PBC）。
• 系统处于半满填充（half-filling），具有反演对称性（inversion symmetry）。
• 通过引入 U(1) 规范扭曲（flux insertion, $\theta$）来定义 Berry 相。$\theta$ 从 $0$变化到$2\pi$，模拟绝热循环。

数值方法：

• 多体计算： 使用密度矩阵重整化群（DMRG）方法（基于 TenPy 库）精确求解多体基态 $|\Psi_0(\theta, U)\rangle$。计算了激发能隙、量子度规（quantum metric）和 Berry 相。
• KS 构建： 基于 DMRG 计算得到的多体基态密度 $n_i(\theta)$，构建非相互作用的 KS 有效哈密顿量。
  • 由于系统具有二聚化结构和反演对称性，密度约束极大地简化了 KS 势的求解问题，将其简化为单参数问题（仅寻找一个化学势偏移 $\Delta_{KS}$）。
  • 在数值精度范围内，发现多体密度 $n_i$ 对 $\theta$ 和相互作用强度 $U$ 均不敏感（保持常数）。

关键指标：

• Berry 相 ($\gamma$)： 通过 Wilson loop（相邻态重叠的乘积）的相位计算。
• 量子度规 ($g_{\theta\theta}$)： 通过相邻态重叠的模长偏差估算，用于衡量波函数在参数空间中的局部几何响应。

3. 主要结果 (Key Results)

A. 密度的“平坦性” (Density Flatness)

• 在研究的整个 $U$ 范围（从非相互作用到强耦合）和 $\theta$ 循环中，多体基态的局域密度 $n_i$ 在数值精度内保持恒定（$\langle n_i \rangle \approx 0.5$）。
• 这意味着密度对通量 $\theta$ 和相互作用 $U$ 没有依赖性。
• 推论： 密度约束无法提供关于相互作用依赖性的任何非平凡信息。因此，对称性保持的 KS 势被强制简化为 $\Delta_{KS} = 0$ 的 SSH 型二次型（自由费米子）代表。KS 哈密顿量在形式上与 $U$ 无关。

B. 多体波函数的几何响应 (Many-body Geometric Response)

• 尽管密度不变，多体波函数表现出显著的几何响应：
  • 量子度规 ($g_{\theta\theta}$)： 在中等 $U$ 区域依赖于 $\theta$；在强耦合极限（大 $U$）下，由于电荷涨落被冻结（charge fluctuation freezing），量子度规被强烈抑制并趋于零。
  • 这表明多体波函数在希尔伯特空间中经历了非平凡的几何演化，而不仅仅是全局相位旋转。

C. Berry 相的一致性 (Agreement of Berry Phases)

• 核心发现： 尽管 KS 描述是 $U$ 无关的自由费米子模型，而多体系统随 $U$ 变化，但两者计算出的 Berry 相 $\gamma$ 在整个能隙区域（gapped regime）内完全一致（均为 $0$或$\pi$，模$2\pi$）。
• 机制解释： 这种一致性并非因为密度编码了多体 Berry 联络（Berry connection）。
• 根本原因： 这是对称性强制的 $Z_2$ 扇区匹配（symmetry-enforced $Z_2$ sector matching）。
  • 在具有反演对称性的一维能隙系统中，Berry 相是 $Z_2$ 拓扑保护的（量子化为 $0$或$\pi$）。
  • 只要体带隙不闭合且对称性保持，同一绝热类中的所有哈密顿量（包括相互作用的 SSH-Hubbard 和对应的非相互作用 KS 模型）必须拥有相同的量子化 Berry 相。
  • KS 模型恰好属于同一个 $Z_2$ 对称保护拓扑（SPT）类，因此其 Berry 相被对称性“锁定”为与多体系统相同。

D. 对比与反例

• 论文指出，如果相互作用项本身显式依赖于通量（例如 $U \to U(\theta)$），则密度可能保持平坦，但多体 holonomy（整体几何相位）可能变得非平凡且不再与 KS 描述一致。这证明了密度匹配本身并不保证几何等价性。

4. 主要贡献 (Key Contributions)

1. 澄清了密度与几何相位的关系： 明确证明了在强关联绝缘体中，基态密度本身并不包含多体 Berry 相的几何信息。密度可以是“平坦”的（与 $U$ 和 $\theta$ 无关），而波函数却具有复杂的几何结构。
2. 揭示了 KS 与多体一致性的真正机制： 指出在 SSH-Hubbard 模型中观察到的 KS 与多体 Berry 相的一致性，并非 DFT 能够重构相互作用拓扑的普遍证据，而是由对称性保护的拓扑量子化（Symmetry-protected topological quantization）所强制的。
3. 区分了局部几何响应与全局拓扑不变量： 展示了量子度规（局部几何指标）随 $U$ 的剧烈变化（在强耦合下被抑制），而 Berry 相（全局拓扑不变量）却保持恒定。这揭示了拓扑保护与局部动力学响应之间的层级结构。
4. 提供了受控的理论框架： 利用一维 SSH-Hubbard 模型作为受控平台，通过 DMRG 精确分离了密度约束和波函数几何演化，为理解关联拓扑系统中的 DFT 局限性提供了清晰的物理图像。

5. 意义与启示 (Significance)

• 对 DFT 在拓扑材料计算中的应用提出警示： 在强关联体系中，仅依靠密度匹配的 KS 描述可能无法正确捕捉相互作用依赖的几何性质（如量子度规），除非该性质被对称性强制量子化。如果对称性破缺或系统处于非量子化区域，KS 描述可能会给出错误的拓扑预测。
• 深化了对拓扑相的理解： 强调了拓扑不变量（如 Berry 相）的鲁棒性源于对称性和能隙，而非具体的微观相互作用细节。只要处于同一拓扑类，非相互作用的代表模型（KS 系统）就能捕捉到正确的拓扑分类。
• 未来方向： 论文建议进一步研究没有保护对称性的区域（此时 Berry 相不再量子化），以及在其他相互作用分类（如 $Z_8$ 分类）下的 KS 与多体一致性，以验证“对称性强制匹配”是否是更普遍机制的特例。

总结：
该论文通过精确的数值模拟证明，在 SSH-Hubbard 链中，Kohn-Sham 描述与多体描述在 Berry 相上的一致性是由对称性保护的拓扑量子化强制实现的，而非因为电子密度编码了多体几何信息。这一发现揭示了密度泛函理论在处理强关联拓扑系统时的内在机制与局限性。