Plane-Wave Photon-Fock Cavity QED-DFT: Chiral-Cavity-Induced Topology in Graphene

技术摘要：石墨烯中的平面波光子-福克腔量子电动力学密度泛函理论 (Plane-Wave Photon-Fock Cavity QED-DFT)

问题陈述
扩展周期性系统的光与物质相互作用的理论描述，需要一个能够将电子与光子置于平等地位的框架。虽然量子电动力学密度泛函理论 (QEDFT) 已成功应用于使用基于波函数的计算方法（如耦合簇、构型相互作用）和局域基组的分离分子系统，但将其能力扩展到周期性固体仍面临挑战。现有的用于固体的平面波 (PW) 实现通常通过有效密度泛函数相互作用或实空间表示来处理光子耦合，缺乏能够精确携带极化子轨道的显式福克空间 (Fock-space) 拟设。这一差距限制了对周期性晶格中真空驱动物质进行系统、第一性原理计算的能力，特别是在研究由腔场诱导的拓扑和对称性破缺方面。

方法论
作者开发了一种针对周期性系统的量化速度规范 Pauli–Fierz 哈密顿量的显式平面波 (PW) $\times$ 光子-福克 (Photon-Fock) 方法。关键方法论特征包括：

• 哈密顿量公式化： 该方法利用长波（偶极）近似下的速度规范 Pauli–Fierz 哈密顿量。由于量化的矢量势 $\hat{A}$ 在空间上是均匀的，它与晶格平移对易，从而保持晶体动量 $k$ 作为一个良好的量子数。
• 算符位移： 光-物质耦合简化为晶体动量的算符值位移：$\mathbf{K} \to \mathbf{K} + \hat{A}/c$。这使得现有的标准固态 DFT 平面波机制可以在极小修改的情况下被复用。
• 基组： 极化子 Kohn–Sham 轨道在平面波（针对物质）与光子数（针对光）的张量积基组中展开：
  $$\Psi_{nk}(\mathbf{r}) = \sum_{\mathbf{G}} \sum_{m=0}^{M} C_{nk}(\mathbf{G}, m) \frac{1}{\sqrt{\Omega}} e^{i(\mathbf{k}+\mathbf{G})\cdot\mathbf{r}} |m\rangle$$
  其中 $M$ 是福克空间截断值。
• 实现细节：
  • 动能项和局部势项在光子指数上是对角化的，而非局部伪势则需要一致的速度规范位移（$\beta(\mathbf{K}) \to \beta(\mathbf{K} + \hat{A}/c)$）。
  • 计算成本扩展效率高，仅需 $M+1$ 次快速傅里叶变换 (FFT) 处理局部势以及处理光子耦合的带状运算。
  • 该方法采用单一腔模。虽然曾考虑使用投影算符 $\hat{Q}$ 来消除因子化拟设中不物理的自位移通道，但生产计算采用了稳定的收敛、未投影处理，以避免在强耦合下的反馈不稳定性。
  • 忽略了平均场水平以上的电子-光子交换相关效应；光子在每个极化子轨道层面上被精确处理，但泛函仍保持仅含物质的特性。

核心贡献

1. 用于腔 QED 的第一性原理 PW 框架： 本工作将此前用于分子的电子-光子福克空间构建扩展到了周期性、第一性原理平面波伪势设置中。它将周期性固体腔 QED 置于与已建立的分子福克空间多体方法相同的概念基础上。
2. 对称性互补驱动： 该方法应用于单层石墨烯在两种线性及手性（圆偏振）腔场中的情况，展示了极化矢量 $\varepsilon$ 如何从定性上改变系统的拓扑结构和对称性。
3. 拓扑不变量计算： 作者对占据流形进行了完整的布里渊区计算，以贝里曲率 (Berry curvature) 和陈数 (Chern number)，直接证实了腔诱导能隙的拓扑性质。

结果

• 手性腔对比线性腔：
  • 手性腔： 打破了时间反演对称性，在狄拉克点处打开了一个哈尔德能隙 (Haldane gap)（在 $A_0/c = 0.04$ a.u. 时，$\Delta_{\text{Hal}} = 129$ meV）。该能隙通过虚耦合在物质部门产生，能隙形成态主要保持物质特性（$w_0 \gtrsim 0.93$）。
  • 线性腔： 保留了时间反演对称性，但打破了晶格的三重旋转对称性。它使狄拉克点保持无能隙状态（在数值精度范围内），但诱导了光学各向异性。
• 陈数阶梯： 占据流形携带一个量化的陈数 $C$，该数随耦合强度增加呈现出非单调序列：$C = 1 \to 3 \to -1 \to 1 \to 2 \to -1 \to 1$。
  • 该序列包含了在耦合强度增加过程中出现的极窄、符号反转的窗口（例如 $C=-1$），其宽度约为 $1 \times 10^{-4}$。
  • 这些转变对应于远离狄拉克点的局部能隙关闭与重新开启事件，已通过贝里曲率图和 Wilson 环计算得到证实。
• 光学特征：
  • 手性腔： 产生强烈的圆二色性（$\varepsilon^2_+ - \varepsilon^2_-$）和非零的非对角介电张量分量（$\varepsilon_{xy}$），作为打破时间反演对称性和哈尔德质量的光学指纹。
  • 线性腔： 产生强烈的线二色性，但圆二色性微乎其微。
• 态密度与实空间密度： 总态密度呈现出光子副本的梳状结构。物质投影密度恢复了预期的价电子计数。实空间密度的变化显示出电荷向 C-C 键的重新分布，且在线性驱动与手性驱动下具有截然不同的模式。

意义与主张
本文声称建立了一个稳健、系统的路径，用于周期性固体的第一性原理腔 QED 计算，并充分利用了现有的平面波基础设施。其主要意义在于证明了扩展的、真空驱动的物质可以表现出拓扑性质（量化的陈数）以及在有限分子系统中不存在的非单调相位序列。

作者强调，手性腔诱导了一个真正的拓扑相变，其特征是具有极化选择性的哈尔德能隙以及复杂的、符号反转的陈数阶梯。这为“真空驱动”拓扑提供了一个理论基准，有别于摄动或经典修正。这项工作表明，该框架可以扩展到其他周期性系统，包括多层/莫尔结构以及拓扑材料（如 HgTe），并且潜在地可以扩展到多模腔，尽管目前的实现仅限于单模、平均场处理。