来源论文: https://arxiv.org/abs/2605.14909v1 生成时间: May 16, 2026 00:05

揭秘费米-哈伯德系统中的隐藏关联：通过相互作用斜坡探测电荷密度波与赝能隙

0. 执行摘要

在强关联电子系统中，费米子配对（Pairing）与电荷有序（Charge Order）的交织是凝聚态物理的核心课题，尤其是在高温超导体（如铜氧化物）和过渡金属硫族化合物中。然而，在强相互作用状态下，对（Pair）的非局域性（Non-local character）往往掩盖了其空间排列的真实图景。本文解析了 Botond Oreg 等人于 2026 年发表的研究，该工作利用冷原子量子模拟平台，引入了一种创新的“相互作用斜坡”（Interaction Ramp）或“相互作用淬火”（Interaction Quench）技术。通过快速提升吸引强度，研究者将原本跨越多个格点的非局域费米子对“投影”到单格点的双占据态（Doublons）上，从而显著增强了电荷密度波（CDW）关联的可观测性。这一技术不仅揭示了在强关联交叉区隐藏的有序，还为区分无配对的费米液体与存在预成对（Preformed pairs）的赝能隙相（Pseudogap phase）提供了明确的实验依据。

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

1.1 核心科学问题：隐藏的有序

吸引费米-哈伯德模型（Attractive Hubbard Model）是研究超流性与电荷密度波竞争的经典模型。在半填充（Half-filling）状态下，该模型的超流序与 CDW 序在数学上是简并的。尽管冷原子实验已能实现该模型，但在强相互作用交叉区（$4 \lesssim U/t \lesssim 8$），由于费米子对的大小与粒子间距相当，配对呈现出强烈的非局域特征。在这种情况下，传统的原位成像只能看到由于波动导致的原子分布，而无法直接识别出成对排列的空间结构。科学界的难题在于：如何“看见”那些由于空间扩展而变得模糊的物理关联？

1.2 理论基础：吸引哈伯德模型与对关联

哈伯德哈密顿量定义为：

$$H = -t \sum_{\langle ij \rangle, \sigma} c_{i, \sigma}^\dagger c_{j, \sigma} - U \sum_i n_{i, \uparrow} n_{i, \downarrow}$$

其中 $t$ 为跃迁强度，$U$ 为吸引相互作用强度。在 $t=0$ 的极限下，所有粒子形成局域的双占据（Doublons）。当 $t > 0$ 时，粒子对通过跃迁获得动能，形成跨越相邻格点的非局域分量。这种非局域性在降低系统能量（通过 $J = 4t^2/U$ 的对跃迁尺度）的同时，也使得单点成像难以捕捉到完整的对电荷分布。特别是在交叉区，这种非局域性不仅驱动了 CDW 的形成（通过有效的一阶排斥），同时也由于对的电荷分布在两个格点间“平滑化”，掩盖了交替的双占据-空穴图案。

1.3 技术难点：动态调控与时间尺度的博弈

实施“相互作用斜坡”技术的最大难点在于时间尺度的控制。斜坡必须满足：

  1. 足够慢：使得非局域对的两个组成粒子有时间收缩到同一个格点形成双占据（受限于 $h/t$ 尺度）。
  2. 足够快：防止系统发生绝热演化产生新的配对或破坏原有的全局对排列（受限于 $h/J$ 尺度）。

在 $U \gg t$ 的强吸引极限下，$h/t \ll h/J$，窗口期较宽。但在物理意义最丰富的交叉区（$J \approx t$），这两个尺度趋于一致，实验必须精确校准。

1.4 方法细节：相互作用淬火协议

实验采用 $^{40}\text{K}$ 原子的自旋平衡费米气体。核心步骤如下：

  1. 平衡态制备：在特定的 $U/t$ 下制备平衡系统（密度 $n \sim 0.8$，温度 $T/t \sim 0.5$）。
  2. 相互作用斜坡：在 $1 \text{ms}$（约 $0.3 h/t$）内，利用 Feshbach 共振将磁场扫描至共振点，瞬间将吸引强度 $U$ 提升至极大（$U \to \infty$）。
  3. 冻结分布：在 $0.1 \text{ms}$ 内迅速将晶格深度提升至 $70 E_r$，锁定原子空间分布。
  4. 高分辨率成像：利用单原子分辨显微镜进行成像，测量双占据密度及自旋电荷关联。

通过这种“淬火投影”，原本分布在 $i$ 和 $j$ 格点的粒子对被压缩到其中一个格点，将“对的中心”显影为双占据点，从而放大了 CDW 的棋盘图案关联。

2. 关键 benchmark 体系，计算所得数据，性能数据

2.1 CDW 关联的显著增强

研究者对比了无斜坡（Normal imaging）和有斜坡（With ramp）的数据。在 $U/t = 5.8(3)$ 的关键点上：

  • 密度-密度相关器 $\langle n_i n_{i+\delta} \rangle_c$：无斜坡时信号微弱；施加斜坡后，呈现出清晰的交替正负关联，展示了明显的棋盘状有序。
  • 结构因子 $\chi_n(\vec{k})$：在 $(\pi, \pi)$ 点，斜坡后的峰值强度提升了近一倍。这证明了之前的实验低估了系统中实际存在的 CDW 强度。

2.2 区分相图的不同区域

利用斜坡后的数据，实验成功绘制了随 $U/t$ 变化的三个截然不同的物理区间：

  1. 费米液体区 ($U/t \lesssim 4$)：斜坡对 CDW 信号几乎没有提升。该区域主要受泡利不相容原理支配，不存在预成对，系统呈现出费米液体的关联特征。
  2. 交叉/赝能隙区 ($4 \lesssim U/t \lesssim 8$)：这是斜坡技术效果最显著的区域。斜坡揭示了大量隐藏的非局域对，证实了系统处于预成对阶段，且这些对具有明显的空间排斥力，导致了强烈的隐藏 CDW 关联。此区域对应于超流温度 $T_c$ 最高的区域。
  3. 局部对/强耦合区 ($U/t \gtrsim 12$)：在此极限下，对本身就是局域的（双占据为主），斜坡前后的成像结果趋于一致。斜坡不再显著改变关联函数，证实了配对已经高度局域化。

2.3 性能指标与校准数据

  • 单占据率（Singlon fraction）：在 $1 \text{ms}$ 斜坡下，非相互作用气体的单占据率保持不变，而强相互作用气体的单占据率降至背景水平，证明了斜坡有效地将非局域对转化为了双占据，且没有在非配对系统中凭空制造双占据。
  • 磁场精度：Feshbach 共振点位于 $202.1 \text{G}$，斜坡终点磁场确保了 $U/t \gg 100$，保证了投影的有效性。
  • 温度稳定性：通过测量总密度波动和压缩率，确认斜坡过程没有显著的热效应注入，全局温度保持在 $T/t \sim 0.5$ 左右。

注：作为实验物理工作，其“代码”主要体现在数据后处理脚本和控制程序上。实验所用的原始数据处理通常基于自定义的 Python 或 MATLAB 框架。

3.1 关联函数计算 (Connected Correlators)

复现该研究的核心在于计算连接关联函数（Connected Correlator）：

$$\langle AB \rangle_c = \langle AB \rangle - \langle A \rangle \langle B \rangle$$

对于密度关联，需计算 $\langle n_i n_{j} \rangle_c$。在实验图像中，这涉及到对上千张二值化（是否存在原子）图像的统计平均。建议使用 Python 的 numpy 和 scipy.ndimage 进行高效的矩阵卷积运算。

3.2 结构因子分析

通过对空间关联图进行离散傅里叶变换（DFT）获得结构因子：

$$\chi_n(\mathbf{k}) = \frac{1}{N} \sum_{i,j} e^{i \mathbf{k} \cdot (\mathbf{r}_i - \mathbf{r}_j)} \langle n_i n_j \rangle_c$$

在复现时，需要注意棋盘归一化，即计算整流后的关联 $\sum_{\delta} (-1)^{\delta_x + \delta_y} \langle n_i n_{i+\delta} \rangle_c$。

3.3 开源资源与类似工具

虽然本工作的具体控制代码未开源，但冷原子社区广泛使用以下工具进行此类模拟和分析：

  • QuTiP (Quantum Toolbox in Python)：用于模拟哈伯德模型的短时动力学，验证 1ms 斜坡的保真度。 Link
  • Monte Carlo 模拟：可以使用 ALPS (Algorithms and Libraries for Physics Simulations) 或基于 Julia 的 MonteCarlo.jl 进行行列式量子蒙特卡洛（DQMC）模拟，产生平衡态基准数据进行对比。
  • 图像处理：对于单原子分辨成像的数据提取，可参考哈佛 Greiner 组或 MIT Zwierlein 组早期公开的图像识别逻辑。

3.4 实验参数复现指南

  1. 原子数与密度：维持单格点填充率 $n \approx 0.8$。
  2. 晶格电势：设置 $V_x = V_y = 4.3 E_r$，确保 $t/h \approx 340 \text{Hz}$。
  3. 斜坡速度：磁场扫描速率应设为约 $7 \text{G/ms}$，从远离共振处扫描至共振点附近，对应 $U/t$ 从初始值突变至 $> 100$。

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献

  1. [7] Hartke et al., Science 381, 82 (2023)：本项目组的前期工作，奠定了吸引哈伯德模型在冷原子平台的基础测量。
  2. [31] B. Kyung et al., Phys. Rev. B 64, 075116 (2001)：关于赝能隙理论的早期重要论述。
  3. [18] H. Shiba, Prog. Theor. Phys. 48, 2171 (1972)：建立了吸引与排斥哈伯德模型之间的对称性转换（Shiba Transformation），是理解 CDW 与超流简并的基础。
  4. [33] M. Xu et al., Nature 642, 909 (2025)：同期关于利用自旋关联探测赝能隙的研究，提供了互补的视角。

4.2 局限性评论

尽管这项工作在技术上非常巧妙，但仍存在以下局限：

  • 单格点映射的不确定性：虽然斜坡旨在将对压缩到中点，但如果费米子对跨越了 3 个或更多格点，或者在压缩过程中发生了格点跳跃，单占据态可能会被错误地分配到相邻格点。虽然文章中提到这种概率由于强 $U/t$ 而很小（$\sim 10^{-4}$），但在极强关联区，这对长程有序的精细结构测量可能仍有影响。
  • 探测损耗问题：实验观察到约 13% 的双占据原子损耗。在高精度关联测量中，损耗往往会系统性地削弱 CDW 信号的绝对强度，尽管可以通过校准修正，但引入了额外的背景噪声。
  • 温度限制：目前的实验温度 $T/t \sim 0.5$ 仍高于 BKT 转变温度。这意味着观测到的是 CDW 的强关联波动（Fluctuations）而非长程有序的基态。如何将此技术推广至更低温度以直接观测相变，是未来的重大挑战。
  • 自旋平衡限制：目前主要针对自旋平衡系统。在自旋不平衡（Spin-imbalanced）系统中，对的结构会更加复杂（如 FFLO 相），斜坡技术的投影效果是否依然清晰尚待验证。

5. 其他必要的补充

5.1 对偶性：从吸引到排斥

本研究的一个深远意义在于它对排斥哈伯德模型（Repulsive Hubbard Model）的启示。根据 Shiba 变换，吸引模型的 CDW 对应于排斥模型在半填充附近的磁有序。文章在结尾展望，类似的技术可以应用于掺杂的排斥哈伯德模型，以揭示那些被空穴稀释隐藏的“条纹序”（Stripe Order）。条纹序被认为是高温超导机制的关键，但其在冷原子实验中的直接观测一直受限于极低的温度要求和关联模糊。相互作用斜坡可能成为探测这些复杂有序的“显微镜”。

5.2 局部测温计（Local Thermometry）

该技术还意外提供了一种全新的测温手段。斜坡后的单占据率（Singlon fraction）直接反映了系统初始状态下的熵含量。相比于传统基于非局域关联对比理论模型的测温法，这种方法更加直接且对数据量的依赖较低，可作为冷原子系统的一种高效“局部温度计”。

5.3 对量子计算的意义

在量子模拟向容错量子计算过渡的过程中，这种“态投影”测量方法展示了如何通过动力学演化将难以直接观测的量子纠缠或非局域关联转化为易于读取的局部观测值。这为量子算法中的算符测量优化提供了物理层面的借鉴。

5.4 总结

Oreg 等人的这项工作证明了：在量子模拟中，主动的动态干预（相互作用斜坡）往往比静态的观测更能揭示自然界的深层规律。通过将非局域性这一“物理障碍”转化为投影测量的“物理杠杆”，我们不仅看清了 CDW 的真面目，也为理解高温超导中的预成对现象迈出了实验性的一大步。