来源论文: https://arxiv.org/abs/2605.11830v1 生成时间: May 13, 2026 13:26

0. 执行摘要

在强关联电子体系的研究中，莫特绝缘体（Mott Insulators）的掺杂演变一直是凝聚态物理的核心课题之一。传统的莫特物理关注电荷能隙的开启以及掺杂导致的超导电性。然而，一个悬而未决的关键问题是：在莫特能隙中出现的“能隙内态”（in-gap states）究竟承载了怎样的物理信息？它们是否能够作为底层自旋动态（spin dynamics）的直接探针？

本研究由南京大学物理学院及协同创新中心的科研团队完成，发表于 2026 年。文章通过量子多体数值方法——团簇摄动理论（Cluster Perturbation Theory, CPT），系统地研究了一维（1D）及准一维（quasi-1D）Kitaev-Hubbard 模型。研究发现，通过调节 Kitaev 型跳跃项（$t'$）来改变自旋各向异性，能隙内态的动力学色散表现出与底层自旋激发谱高度一致的对应关系。在 Z 链中，能隙内态随系统从 Heisenberg 模型向 Ising 模型转变而演化；在 XY 链中，能隙内态完美镜像了 Jordan-Wigner 费米子谱；而在双腿梯子（two-leg ladder）结构中，能隙内态呈现出连续谱特征，反映了自旋激发的分数化（fractionalization）以及拓扑 $Z_2$ vison 的存在。这一发现为利用 ARPES 或 STM 等电荷探测手段研究量子自旋液体提供了全新的理论路径。

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

1.1 核心科学问题

莫特绝缘体的父相通常具有长程反铁磁序，掺杂会抑制磁序并诱导高温超导等奇特物态。在光谱学上，莫特能隙中涌现的能隙内态（in-gap states）伴随着光谱权重的转移。早期的研究暗示这些态可能与自旋极化子（spin-polaron）或空穴子-背离（holon-backflow）相关，但其微观起源及与自旋激发的系统性定量关系仍不清晰。本文旨在回答：在具有各向异性交换相互作用的 Kitaev 物理背景下，电荷激发（能隙内态）与自旋激发（磁子或分数化激发）之间是否存在普适的映射关系？

1.2 理论基础：Kitaev-Hubbard 模型

研究的对象是 Kitaev-Hubbard 模型，其哈密顿量定义如下：

$$H = -\sum_{\langle ij \rangle_\alpha} c_i^\dagger \left( \frac{t + t' \sigma_\alpha}{2} \right) c_j + H.c. + U \sum_i n_{i\uparrow} n_{i\downarrow} - \mu \sum_{i,\sigma} n_{i\sigma}$$

其中，$t$ 是传统的各项同性跳跃常数，$t'$ 是 Kitaev 型跳跃项，$\sigma_\alpha$ 是 Pauli 矩阵（$\alpha = x, y, z$ 取决于键的取向）。在强耦合极限（$U \gg t, t'$）下，该模型可以映射到著名的 Kitaev-Heisenberg 模型：

$$H = \sum_{\langle ij \rangle_\alpha} K S_i^\alpha S_j^\alpha + J S_i \cdot S_j$$

其中有效交换相互作用为 $K = 2t'^2/U$ 和 $J = (1-t'^2)/U$。这种映射为连接电子态与自旋态提供了天然的理论框架。

1.3 技术难点：强关联与多体激发的计算

计算莫特绝缘体中的单粒子谱函数 $A(\mathbf{k}, \omega)$ 面临多重挑战：

  1. 非微扰特性：$U$ 项使得单粒子图像失效，无法直接应用能带理论。
  2. 热力学极限与有限尺寸效应：精确对角化（ED）受限于团簇大小，而动力学平均场理论（DMFT）在低维体系中忽略了动量空间的重要关联。
  3. 掺杂态的复杂性：掺杂一个空穴后的基态可能处于完全不同的希尔伯特子空间。

1.4 方法细节：团簇摄动理论 (CPT)

为了克服上述难点，作者采用了团簇摄动理论。CPT 的核心思想是将无限格点划分为完全相同的有限尺寸团簇（本文采用 12 格点团簇）。

  • 步骤 1：在团簇内部利用 ED 计算精确的格林函数 $G_c(z)$。$z$ 是复频率。
  • 步骤 2：将团簇间的跳跃项 $V(\mathbf{k})$ 作为摄动，通过 Dyson 方程重构全空间的格林函数： $$G_{cpt}^{-1}(\mathbf{Q}, z) = G_c^{-1}(z) - V(\mathbf{Q})$$
  • 步骤 3：通过周期化（Periodization）恢复平移对称性，得到动量相关的格林函数 $G(\mathbf{k}, z)$。
  • 步骤 4：取虚部得到谱函数 $A(\mathbf{k}, \omega) = -\frac{1}{\pi} \text{Im} G(\mathbf{k}, \omega + i\eta)$。

本文采用了 $\eta/t = 0.1$ 的展宽因子，并在 $t' = 0.99$ 处模拟 Kitaev 极限，以避免数值奇异性。

2. 关键 Benchmark 体系与计算数据分析

2.1 Z 链：从 Heisenberg 到 Ising 的演化

Z 链只包含 $z$-型键（$\sigma_\alpha = \sigma_z$）。

  • $t' = 0$ 极限：对应各项同性 Heisenberg 模型，自旋激发是无隙的。计算得到的谱函数显示，能隙内态在 $k = \pi/2$ 附近表现出无隙特征。
  • $t' \to 1$ 极限：对应反铁磁 Ising 模型，自旋激发具有有限能隙。计算结果发现，随着 $t'$ 增加，能隙内态逐渐变平且开启能隙。
  • 定量拟合：图 3 展示了电荷能隙 $\Delta \omega$ 随 $1/U$ 的变化。数据明确显示，$\Delta \omega$ 完美收敛于 $2t'^2/U$，这正是 Ising 模型的自旋翻转激发能。这证明了在 Z 链中，能隙内态的能量尺度直接受控于自旋交换能。

2.2 XY 链：费米子色散的完美镜像

XY 链具有交替的 $x$ 型和 $y$ 型键。在纯 Kitaev 极限（$t'=1$），自旋模型可通过 Jordan-Wigner 变换映射为无相互作用的费米子。

  • 能带结构：计算结果（图 4d, 4e）显示，能隙内态分裂为一个色散支（dispersive branch）和一个近乎平坦的支（flat branch）。
  • 理论验证：通过 Jordan-Wigner 变换得到的解析谱为 $E(k) = \pm \sqrt{K_x^2 + K_y^2 + 2K_x K_y \cos(k)}$。作者发现，数值计算出的电荷谱色散曲线与该解析公式在定性及半定量上高度吻合。特别是 $k=0$ 处的激发峰，精确对应了 Majorana 费米子的激发边界。

2.3 二腿梯子：连续谱与 Vison 激发

这是本文最具突破性的部分。Kitaev 梯子是二维 Kitaev 蜂窝格子的准一维简化，支持分数化激发。

  • 连续谱的涌现：与链状结构不同，梯子体系的能隙内态不再是孤立的相干准粒子能带，而是呈现出弥散的连续谱（图 5d, 5e）。
  • 物理机制：这种连续谱源于局部自旋翻转被分解为巡游的 Majorana 费米子和静态的 $Z_2$ 规范场通量（vison）。
  • Vison 隙的观测：谱函数中存在一个约 $0.4t$ 的能隙，将能隙内态与 LHB 分开。通过 Appendix B 的解析计算对比，确认该能隙直接对应于产生一对 vison 所需的最小能量。这是首次在电荷谱中明确识别出拓扑 vison 激发的贡献。

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

3.1 核心算法：CPT + ED

复现本工作的核心在于构建一个高性能的 CPT 求解器。主要的步骤如下：

  1. 团簇构建：对于 12 格点体系，希尔伯特空间维度在半填充附近约为 $\binom{12}{6} \times \binom{12}{6} \approx 8.5 \times 10^5$。虽然可以在个人工作站完成，但建议使用并行化的 Lanczos 算法。
  2. ED 求解核心：
    • 使用 Lanczos 迭代计算基态 $|GS\rangle$。
    • 计算连分数形式的格林函数 $G_c(z) = \langle GS | c \frac{1}{z - (H - E_0)} c^\dagger | GS \rangle$。
    • 由于涉及 $t'$，哈密顿量矩阵是复数矩阵（因为 $\sigma_y$ 项）。需要使用 std::complex<double> 进行运算。
  3. CPT 周期化：这涉及到对 $V(\mathbf{k})$ 矩阵的逆运算。对于 $N_c = 12$ 的团簇，每一步 $\mathbf{k}$ 点需要求逆一个 $24 \times 24$ 的复数矩阵。可以使用 LAPACK 或 Eigen 库。

3.2 软件包建议

虽然作者未提供开源库链接，但以下开源项目可以用于构建类似的计算流：

  • QuSpin (Python)：非常适合构建 Kitaev-Hubbard 哈密顿量并进行小尺寸 ED 计算。
  • TRIQS (C++/Python)：具有强大的多体格林函数处理框架，其 CPT 插件可以扩展用于此类计算。
  • DMRG (ITensor)：如果需要研究更长的链或更复杂的梯子，DMRG 是比 ED 更好的选择，尽管计算单粒子谱函数的动力学过程更为昂贵。

3.3 复现参数表

参数数值备注
团簇尺寸 ($N_c$)12周期性边界条件或开边界条件
跳跃常数 $t$1.0能量单位
Kitaev 项 $t'$0.0 - 0.99调控磁各向异性
库伦斥力 $U/t$3.0 - 5.0确保处于莫特相
展宽 $\eta/t$0.1用于获取平滑的光谱线型
掺杂浓度 $\delta$1/12单空穴掺杂

4. 关键引用文献与局限性评论

4.1 关键文献回顾

  1. Kitaev (2006): 奠定了 Kitaev 模型的理论基石，定义了 Majorana 费米子和 $Z_2$ 规范场。
  2. Jackeli & Khaliullin (2009): 提出了在铱氧化物（iridates）中实现 Kitaev 相互作用的物理机制，是本文研究 Kitaev-Hubbard 模型的实验背景。
  3. Senechal et al. (2000, 2002): CPT 方法的奠基性工作，本文的数值框架即基于此。
  4. Kohno (2015): 讨论了莫特绝缘体中能隙内态的起源，本文在此基础上将其扩展到了具有各向异性交换的体系。

4.2 局限性评论

尽管本文结论显著，但仍存在以下局限：

  • 一维性的约束：研究集中在 1D 和准 1D。虽然作者讨论了 vison，但在二维系统中，vison 是真正局域的拓扑缺陷，电荷激发与它们的耦合机制可能更加复杂。
  • 有限尺寸效应：12 格点的团簇虽然捕捉了局部短程关联，但对于具有长程关联的相（如某些自旋液体相），CPT 的精度可能会受限。图中 $k=0$ 附近的一些震荡特征可能源于团簇边界。
  • 单空穴近似：研究仅考虑了极低浓度的掺杂。在有限浓度下，空穴间的相互作用以及它们对背景磁序的反馈（如形成极化子晶格）未被考虑。
  • 计算复杂性：Kitaev-Hubbard 模型的非对角项导致希尔伯特空间不能简单地通过 $S_{total}^z$ 划分（除了 Z 链），这限制了更大规模团簇的模拟。

5. 补充内容：从实验室到理论的桥梁

5.1 实验观测的可能性

本文最具吸引力的观点是：“能隙内态可以作为自旋液体的探针”。对于实验物理学家而言，这意味着：

  • ARPES (角分辨光电子能谱)：如果我们在掺杂的 $Na_2IrO_3$ 或 $\alpha-RuCl_3$ 中观察到能隙内态分裂为明显的 dispersive 和 flat 两支，这可以直接证明底层 Kitaev 物理的存在。
  • STM (扫描隧道显微镜)：局部态密度（LDOS）的能量分布应反映 vison 隙。特别是在梯子体系中，连续谱的起始位置直接锁定了 vison 的能量尺度。

5.2 Jordan-Wigner 变换的直观理解

为什么电荷谱会镜像费米子自旋谱？这可以从“自旋-电荷复合”图像理解。在莫特背景下，移动的空穴（holon）会留下一个改变了磁背景的“磁子”或“分数化费米子”。如果磁激发是费米性质的（如 Kitaev 极限下），那么空穴必须与其耦合以维持体系的量子数平衡。这种耦合导致空穴的动能被磁激发的色散关系调制。因此，观察到的电荷色散本质上是“空穴 + 分数化自旋激发”的复合体，而由于空穴动能在某些路径上被抑制，自旋激发的特征便凸显出来。

5.3 结论与展望

这项工作成功地在电子结构与奇特自旋动力学之间架起了桥梁。它不仅加深了我们对莫特绝缘体谱权转移的理解，更关键的是，它为寻找“自旋液体”这一凝聚态物理的圣杯提供了新的搜索策略。未来的研究可以进一步扩展到非平衡态动力学，观察超快激光激发如何调控这些能隙内态，从而实现超快自旋操控。