Fock 空间碎裂下的精确刚度与动力学响应

论文信息

标题: Exact Stiffness and Dynamical Responses from Fock-Space Fragmentation

作者: Jonah Herzog-Arbeitman, Eslam Khalaf, Zhaoyu Han

发布日期: 2026-08-27

arXiv ID: 2608.27446v1

PDF 链接: 下载 PDF

3 分钟速览

• 研究问题：在量子几何嵌套（QGN）模型中，自发破缺连续对称性的刚度是否严格等于其高斯/BCS 变分值；更一般地，动力学响应能否在强关联平带体系中精确求解。
• 核心方法：利用 QGN 模型的 Fock 空间碎片化：费米子双线性生成元使激发从精确基态出发时按粒子–空穴组成 $(p,h)$ 封闭，Krylov 子空间不变。
• 关键结果：证明 QGN 超导体的超流刚度严格等于 BCS 变分上限 $D_s=(N_f/2)\nu(1-\nu)m_{\rm pair}^{-1}$（论文式 (6)）；光学电导率、结构因子和单粒子格林函数均可由少体谱精确确定。
• 主要局限：结果依赖 QGN 条件、严格平带投影和厄米生成元；对非厄米 QGN 手性平带超导体，Bootstrap 下界仍低于 BCS 变分上界。
• 适合读者：从事凝聚态理论、平带超导/铁磁、量子几何、可解多体模型和量子多体 Bootstrap 研究的读者。

论文背景和研究动机

在强关联量子多体物理中，可解模型稀缺且通常只提供部分谱信息。QGN 模型是一类特殊的平带费米子模型，具有 frustration-free 基态，同时其生成元 $S_{\mathbf{q},I}$ 与配对算符 $\eta^\dagger$ 对易，从而给出严格的 BCS 型基态。此类模型之所以引人关注，是因为在严格平带极限下，传统动能贡献消失，超流刚度完全由相互作用和量子几何产生。此前的一项量子多体 Bootstrap 工作发现，BCS 变分上界与 Bootstrap 下界数值上饱和，并猜测真正的多体刚度为 $D_s=D_s^{\rm BCS}$（见论文第 I 节）。本文正是证明了这一猜想。

刚度问题之所以困难，是因为 Kohn 公式中的顺磁项需要对全部激发谱求和。通常平均场计算只给出变分估计；若激发态混合超出变分流形，就会产生额外多体 dressing。本文的关键观察是：QGN 的 Fock 空间碎片化可以严格限制这些激发扇区，使问题回到少体谱。

核心方法和技术细节

QGN 哈密顿量写成

$$H=\sum_{\alpha} S_\alpha^\dagger S_\alpha,\qquad S_\alpha^\dagger|{\rm GS}\rangle=S_\alpha|{\rm GS}\rangle=0.$$

其中每个 $S_\alpha$ 是费米子双线性。对任意激发算符 $\mathcal O$，作用在基态上可化简为双对易子：

$$H\mathcal O|{\rm GS}\rangle =\sum_\alpha [S_\alpha^\dagger,[S_\alpha,\mathcal O]]|{\rm GS}\rangle.$$

由于与费米子双线性的对易保持粒子和空穴算符数量，激发谱按粒子–空穴组成 $(p,h)$ 分块。定义

$$\mathcal K_{p,h}= \operatorname{span}\{c^\dagger_{a_1}\cdots c^\dagger_{a_p} c_{b_1}\cdots c_{b_h}|{\rm GS}\rangle\},$$

则 $H\mathcal K_{p,h}\subseteq\mathcal K_{p,h}$。这就是本文所称的 Fock 空间碎片化（见论文第 II.2 节）。

对超流刚度的证明如下。施加恒定矢势 $\mathbf A$，展开

$$H(\mathbf A)=H+A_iH_i+\frac12 A_iA_jH_{ij}+O(A^3).$$

由于 $S_{\mathbf q,J}|0_N\rangle=0$，一阶能量修正为零；二阶 Kohn 公式为

$$\partial_i\partial_j E(\mathbf 0) =\langle0_N|H_{ij}|0_N\rangle -2\Re\sum_{E_M>0}\frac{\langle0_N|H_i|M\rangle\langle M|H_j|0_N\rangle}{E_M}.$$

通过双对易子结构，作者证明 $H_i|0_N\rangle$ 仅由 $[\partial_iS_{-\mathbf q,I},S_{\mathbf q,J}]|0_N\rangle$ 贡献，而后者完全落在 1 玻色子（即一粒子–一空穴）片段内。因此 Lehmann 求和中只有 1-boson 态贡献。

进一步，作者证明 1-boson 激发恰好等于 BCS 变分流形的切空间。利用 BdG 等式

$$\gamma_{\mathbf k,n}\eta^{\dagger N}|0\rangle =2N\sum_m F_{nm}(\mathbf k) \gamma^\dagger_{-\mathbf k,m}\eta^{\dagger(N-1)}|0\rangle,$$

可以得到一个中性激发在凝聚体上只变形一个 Cooper 对。BCS 型试探态 $(\eta^\dagger[F+\delta F])^N|0\rangle$ 的一阶变分正是所有 1-boson 激发。因此，精确一阶形变包含在 BCS 切空间内，变分曲率穷尽了 Kohn 响应。这给出等式（见论文第 III.3 节，式 (38)）

$$[D_s]_{ij} =\frac{N_f}{2}\nu(1-\nu)[m_{\rm pair}^{-1}]_{ij}.$$

同期，论文将同一机制用于动力学关联函数。光学电导率、投影结构因子和单粒子格林函数分别表现出封闭性。最终得到

$$\sigma_{ij}(\omega) =2N_f\nu(1-\nu)\omega \sum_{\mu>0}\operatorname{Tr} \mathcal P_\mu\partial_i\mathcal P_0\partial_j\mathcal P_0 \,\delta(\omega-E_\mu),$$

其中 $\mathcal P_\mu$ 是 1-boson 本征模投影子（论文式 (44)）。单粒子格林函数精确为

$$G(\mathbf k,\omega) =\frac{1-\nu}{\omega-\Sigma(\mathbf k)+i0^+} +\frac{\nu}{\omega+\Sigma(\mathbf k)+i0^+},$$

其中 $\Sigma(\mathbf k)$ 由两点自能精确给出（论文式 (48)）。

创新点和贡献

本文的第一个主要贡献是将此前仅作为猜想的超流刚度等式变成严格证明。证明的关键不是直接计算谱求和，而是把精确一阶波函数嵌入 BCS 切空间。这解释了为什么类似 Bootstrap 数值下界会饱和。

第二个贡献是系统化了 QGN 模型的 Fock 空间碎片化，并将其从基态扩展至动力学响应。文中指出：任何保持固定 Fock 空间组成的微扰，其响应函数都只由对应的少体扇区决定。因此，光学电导率、动态结构因子和单粒子格林函数都可由一点或两点问题精确求解（见论文第 IV 节）。

第三个贡献是把碎片化推广到自旋系统。文中构造了自旋-1 局域矩的 $SU(2)$ 生成元 $K_j^a$，基态为 $|0\rangle_j$ 直积，并满足

$$K_j^a|{\rm GS}\rangle=(K_j^a)^\dagger|{\rm GS}\rangle=0.$$

此时 Krylov 空间按 $SU(2)$ 不可约表示分块，而不是粒子–空穴组成。这为其他可解响应函数的构造提供了方向（见论文第 V 节）。

局限与待解决问题

该结果高度依赖 QGN 模型的特殊结构。具体地，需要严格平带、无限远带隙、投影相互作用和 QGN 条件。论文未讨论远程带耦合、有限温度、无序或偏离平带的扰动对刚度公式的修正。

更关键的限制是厄米性。论文明确指出，对手性平带超导体中由非厄米双线性平方构造的 frustration-free 模型，生成元的共轭不湮灭基态，双对易子约化失效，微扰可以耦合到更高粒子–空穴扇区。因此这些模型中 Bootstrap 下界仍严格低于 BCS 变分上界；作者并未将其归为可解响应（见论文第 VI 节）。

此外，自旋模型的碎片化结构虽已建立，但论文未给出相应的动力学响应精确求解，只将其留作未来工作。更一般的碎片化条件、更多可精确求解响应函数的模型类别，也还没有被完整刻画。