来源论文: https://arxiv.org/abs/2603.15560v1 生成时间: Mar 20, 2026 22:49

执行摘要

伪能隙（Pseudogap, PG）态是强关联电子物理中最具挑战性的谜团之一，广泛存在于铜氧化物、镍酸盐及魔角石墨烯等体系中。传统的视角往往将其归结为莫特（Mott）物理或某种竞争序的波动。然而，Andreas Gleis 和 Gabriel Kotliar 在其最新论文《Pseudogapped Fermi liquids from emergent quasiparticles》中提出了一种截然不同的范式。

该研究的核心贡献在于提出了一个在任意空间维度下均可精确求解的交互模型。通过引入关联跃迁（Correlated Hopping, CH）相互作用，作者展示了单电子如何“碎裂”为三个涌现的准粒子（称为 $q$-准粒子）。这种机制在不诉诸莫特绝缘态的前提下，自然地产生了一个具有伪能隙特征的单粒子谱函数，并伴随着半填充处准粒子权重 $Z$ 的消失。此外，该工作利用高斯路径积分表象，建立了一个分析 PG-FL（伪能隙费米液体）到 Landau FL（朗道费米液体）量子相变的严谨框架，并深入探讨了 Luttinger 面（LS）在统计和动力学上的表现。这一发现不仅为理解高温超导体的正常态提供了新的理论工具，也为工程化量子系统的模拟提供了明确的目标。

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

1.1 核心科学问题：伪能隙的非莫特起源

在凝聚态物理中，费米液体理论是描述金属态的基石。然而，PG 金属挑战了这一基石：它们在费米能级附近表现出电子态密度的压制，且在光谱学实验中呈现出“费米弧”（Fermi arcs）而非闭合的费米面（FS）。长期以来，物理学界争论的焦点在于：这种态密度的缺失究竟是由于配对涨落、电荷密度波，还是由于电子靠近莫特绝缘点时的局域化效应？

本工作的核心科学问题是：是否存在一种纯粹的费米子机制，能在不依赖莫特物理的情况下，通过粒子间的动力学关联产生稳定的伪能隙态？

1.2 理论基础：关联跃迁与算符变换

作者考虑了一个具有关联跃迁项的自旋费米子哈密顿量 $H_q$。其形式不同于标准的 Hubbard 模型（局域排斥），而是侧重于电子跃迁过程受邻近位置占据态调制的效应：

$$H_q = \sum_{ij\sigma} (t^q_{ij\sigma} - \mu\delta_{ij}) c^\dagger_{i\sigma} c_{j\sigma} + \sum_{ij\sigma} t^q_{ij\sigma} (4n_{i\bar{\sigma}} n_{j\bar{\sigma}} - 2n_{i\bar{\sigma}} - 2n_{j\bar{\sigma}}) c^\dagger_{i\sigma} c_{j\sigma}$$

这里的关键在于引入了一个非线性的算符映射，将物理电子算符 $c_{i\sigma}$ 转换为所谓的 $q$-算符：

$$q_{i\sigma} = c_{i\sigma}(2n_{i\bar{\sigma}} - 1)$$

这个变换的精妙之处在于，$q_{i\sigma}$ 依然满足正则费米子反对称关系。在 $q$-算符的表象下，原本复杂的相互作用哈密顿量 $H_q$ 竟然坍缩为一个自由费米子哈密顿量：

$$H_q = \sum_{ij\sigma} (t^q_{ij\sigma} - \mu\delta_{ij}) q^\dagger_{i\sigma} q_{j\sigma}$$

这意味着，虽然 $q$-准粒子是自由运动的，但原始电子算符 $c_{i\sigma}$ 与 $q$-准粒子之间存在高度非线性的关联。电子被看作是由三个 $q$-准粒子复合而成的实体（两个粒子，一个空穴的复合过程），这种“电子碎裂化”是产生伪能隙的直接诱因。

1.3 技术难点：高斯路径积分的推广

虽然 $H_q$ 在 $q$ 表象下是平庸的，但实际材料中总会存在破坏这种精确可解性的项，如标准的单粒子跃迁 $H_c$ 或局域库仑力 $H_2$。技术上的难点在于如何建立一个扰动理论框架，既能处理 $H_q$ 的强关联本质，又能引入其他微扰。

作者通过引入基于 $q$-相干态的高斯路径积分解决了这一问题。通过定义 $|\chi\rangle_q = \prod_{i\sigma} e^{-\chi_{i\sigma} q^\dagger_{i\sigma}} |0\rangle$，作者推导出了分区函数 $Z$ 的路径积分表示。在这种表示下，$H_q$ 对作用量的贡献是二次型的，这允许研究者使用标准的多体物理工具（如费曼图、平均场理论）来研究偏离精确解极限的行为。这在处理强关联系统时是极其罕见的成就。

1.4 方法细节：格林函数的精确求解

单电子格林函数 $G_{k\sigma}(z)$ 的计算是理解谱学特征的关键。作者展示了 $G_{k\sigma}$ 可以分解为两个部分：

  1. 相干准粒子极点：其权重为 $Z_\sigma = (2n_{\bar{\sigma}} - 1)^2$。在半填充时，$n_{\bar{\sigma}} = 0.5$，$Z_\sigma$ 恰好变为 0，意味着费米面上不再有寿命无限长的物理电子。
  2. 非相干连续谱：描述了电子碎裂为三个 $q$-准粒子的动力学过程。这一部分的谱函数 $A^{(3)}_{k\sigma}(\omega)$ 在低能下遵循 $\omega^2 + \pi^2 T^2$ 的特征，这正是伪能隙形成的数学源头。

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

2.1 二维方晶格上的谱学特征

作者选取了具有近邻跃迁 $t$ 和次近邻跃迁 $t'$ 的方晶格作为 Benchmark。在半填充（Half-filling）状态下：

  • 谱函数（Figure 1a）：在 $\omega=0$ 附近观察到明显的态密度压制（伪能隙）。与 Landau 费米液体不同，此处不存在清晰的费米极点。
  • 自能量（Figure 1b）：自能量 $\Sigma_{k\sigma}(\omega)$ 表现出极点行为，这对应于 Luttinger 面（Luttinger Surface, LS）。LS 的存在直接导致了格林函数的零点，从而压制了态密度。
  • 局部谱函数（Figure 1d）：随温度升高，伪能隙逐渐被填充。在低温下，由于相空间限制，电子无法有效地碎裂为 $q$-准粒子，导致能隙打开；而高温下热涨落打破了这种限制。

2.2 Luttinger 积分与和规则（Sum Rule）

通过计算 Luttinger 积分 $I_L$，作者发现该模型显著违反了传统的 Luttinger 和规则（Figure 2a）。在非微扰极限下，$I_L$ 从 0 偏离，反映了系统中隐藏电荷的存在。这一 Benchmark 数据有力地支持了“伪能隙金属中存在 Luttinger 面”的假说，且作者通过数值计算展示了 LS 如何随填充（Filling）演化并最终与 FS 湮灭，这为解释实验中的“费米弧”提供了几何上的直观图像。

2.3 从杂质模型到格点模型的扩展

为了验证模型的健壮性，作者使用了 NRG（数值重整化群） 对单杂质和双杂质 Anderson 模型进行了 Benchmark。结果显示（Figure 3a）：

  • 当杂质-浴耦合以 $q$-型为主时，系统进入伪能隙相。
  • 当以 $c$-型（标准电子跃迁）为主时，系统回归到具有 Kondo 峰的常规费米液体。
  • 在临界点 $V_c = V_q$，谱函数高度恰好是 Friedel 和规则值的一半，揭示了一个二通道 Kondo 型的非费米液体相。

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

3.1 核心算法：卷积结构加速

计算单电子格林函数的难点在于动量空间的双重积分（四维积分）。作者利用了动量空间的卷积结构：

$$S_{k\sigma}(t) \propto \sum_{r} e^{-ik\cdot r} [\bar{F}_{\bar{\sigma}}(r,t) F^*_{\bar{\sigma}}(r,t) \bar{F}_{\sigma}(r,t) + ...]$$

这里的 $F(r,t)$ 是费米统计权重的傅里叶变换。通过 快速傅里叶变换（FFT），算法复杂度从 $O(N^2)$ 降低到 $O(N \log N)$。对于 $1024 \times 1024$ 的格点，这是复现该工作的技术门槛。

3.2 软件包建议

  • NRG 部分：推荐使用开源的 MuNRG 软件包或作者提到的基于 QSpace 的张量库。这些工具能够处理复杂的非阿贝尔对称性。
  • GDMFT（广义动力学平均场理论）：复现格点计算需要构建一个自洽循环，将格点模型映射到多轨道的有效杂质模型。需要实现特定的算符变换（如 Hubbard 算符到正则费米子的变换）。

3.3 复现指南步骤

  1. 定义 $q$-表象：构建基于 $q$ 算符的自由能带结构 $\epsilon_k$。
  2. 时域积分：利用等式 (S22)-(S26) 计算时域相关函数 $S_{k\sigma}(t)$。注意 $t$ 的采样范围需要足够大以保证能量分辨率，典型取值 $t_{max} = 2048 \delta t$。
  3. Kramers-Kronig 变换：通过虚部重建实部，建议使用分段线性解析积分法（见论文附录 S-III.B）以增强数值稳定性。
  4. 求解 LS：寻找 $\text{Re} G_{k\sigma}(\omega=0, T=0)$ 的零点轨迹。

4. 关键引用文献与评论

4.1 关键引用

  1. Hirsch & Marsiglio (1989)：关联跃迁概念的鼻祖。本工作将其从超导研究引申到了伪能隙的底层构建。
  2. Dzyaloshinskii (2003)：关于 Luttinger 面及其对和规则贡献的奠基性理论工作。
  3. Fabrizio (2020/2022)：提出了涌现准粒子可能出现在 LS 上的设想。本工作对此进行了微观验证，并纠正了 Fabrizio 关于比热贡献的部分推导错误（见附录第 11 页）。

4.2 局限性评论

  • 模型特异性：该模型之所以精确可解，是因为 $H_q$ 的特殊对称性。在真实的材料（如铜酸盐）中，$H_c$ 和局域排斥 $U$ 通常与 $H_q$ 处于同一量级。虽然作者证明了 PG-FL 在小微扰下的稳定性，但当 $H_c$ 足够强时，系统必然发生一级量子相变进入朗道费米液体，这意味着该模型捕捉的是一种“理想极限”。
  • 超导配对：论文虽提及 Cooper 通道存在不寻常的配对对称性，但未给出具体的超导转变温度（$T_c$）预测。伪能隙如何与超导竞争或演化，仍需后续的超流密度计算。
  • 维度的影响：虽然理论上在任意维度成立，但在一维（1D）中，van Hove 奇异性的存在可能会导致 $q$-准粒子物理与 Luttinger 液体物理的混淆，需要更精细的有限尺寸缩放分析。

5. 补充：关联跃迁与 Hubbard 算符的深层关系

为了更直观地理解为何 $q$-算符如此有效，我们需要考察其与 Hubbard X-算符 的关系。在强关联物理中，电子被禁止双占据。而 $q_{i\sigma} = X^{0\sigma} - X^{\bar{\sigma} 2}$ 实际上是将物理电子拆分成了在两个不同占据子空间（空穴-单占据，单占据-双占据）中跃迁的组合。

本工作最深刻的启示在于：伪能隙不一定要来自某种“消失”的费米子，它也可以来自费米子“碎裂”后的重新组合。 这种组合方式受到了局部希尔伯特空间的几何限制（即那项 $2n_{i\bar{\sigma}}-1$）。

此外，对于量子化学家而言，这种基于算符变换的精确解思路为处理分子轨道中的强关联项（如多组态自洽场方法中的活性空间）提供了借鉴。如果能在大分子体系中识别出类似的“准自由”涌现算符，将极大地降低计算关联能的复杂度。

最后，论文在附录中给出的“无自由度测量”（Non-freeness）是一个非常新颖的角度，利用冯·诺依曼熵来量化基态偏离单 Slater 行列式的程度。结果显示在该模型中，非自由度接近最大值（$\ln 2$），这再次证明了即使在动量空间看起来是“伪能隙”的简单结构，其波函数的量子纠缠背景是极其复杂的。这是一个典型的“大道至简”的物理案例。