多体相显微镜：从相干性、d波超导到格林函数，基于傅里叶空间操控的冷原子新范式

来源论文: https://arxiv.org/abs/2602.12142v1 生成时间: Feb 20, 2026 05:01

执行摘要

传统的量子气体显微镜在观测量子多体系统时，主要侧重于实空间的密度和自旋关联测量，对于揭示系统深层相干性和非局域关联信息则显得力不从心。这篇题为《Protocols for a many-body phase microscope: From coherences and d-wave superconductivity to Green’s functions》的突破性研究，提出了一系列巧妙的实验协议，旨在利用物质波显微镜（matter-wave microscope）的独特傅里叶空间操控能力，实现一个前所未有的“多体相显微镜”。

这项工作核心在于构建一个多体干涉仪，通过精确控制拉曼脉冲在傅里叶空间引入动量转移，并借助物质波透镜将其转化为实空间位移，从而使系统不同部分发生干涉。通过分析自旋分辨测量中获得的干涉条纹，可以直接提取原本难以触及的非对角关联函数和相信息。文章详细阐述了三个关键应用方向：首先，测量费米子 d 波超导序参量，这对于理解高温超导机制至关重要；其次，探测非等时格林函数，进而获取材料的谱函数，揭示粒子激发动力学；最后，用于揭示分数量子霍尔（Fractional Quantum Hall, FQH）绝缘体中复合玻色子的隐藏非对角长程序，为拓扑物理研究开辟新途径。这些协议的提出，不仅克服了现有技术的局限，更预示着冷原子量子模拟在探索复杂量子现象方面将迈入一个全新的、更深层次的阶段，有望直接观测到过去只能通过理论推导或间接推断的物理量，极大地加速我们对多体量子态的理解。

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

1.1 核心科学问题：超越密度关联，直探量子相干

量子气体显微镜（Quantum Gas Microscope）自问世以来，彻底革新了我们对量子多体系统研究的方式。它提供了单原子分辨率的实空间成像能力，使得我们能够直接观测到晶格中的原子占据数，进而测量各种密度-密度关联、自旋-自旋关联，甚至是更复杂的如弦序（string order）关联等。然而，这些测量本质上都是“对角”的，即只关注粒子在特定位置或自旋态上的概率分布。对于量子系统至关重要的“非对角”信息——即不同位置或动量态之间波函数叠​​加的相干性，尤其是其相位信息——却始终难以直接获取。这是当前量子气体显微镜领域面临的核心挑战之一。

相干性是量子现象的基石，它编码了粒子波函数之间的相对相位，直接决定了超流、超导等宏观量子现象的出现。例如，超流序参量、超导序参量、格林函数等都本质上是非对角关联函数，包含了丰富的相位信息。理解这些序参量的空间分布和动力学演化，对于揭示物质的奇异量子相和激发谱至关重要。传统的手段如布拉格谱学可以提供动量空间信息，但缺乏实空间分辨率；而噪声关联（noise correlations）虽能提供一些配对信息，但通常需要复杂的反演过程且缺乏直接的相位相干测量。

本文提出的“多体相显微镜”正是为了弥补这一关键空白。它通过一系列精密的傅里叶空间操控协议，使得冷原子实验能够直接测量长程非对角关联函数，从而揭示多体系统中的相干性和相位信息。这不仅仅是对现有技术的简单升级，更是开启了一个全新的实验范式，使得许多过去仅限于理论研究或间接推断的量子现象，现在有望在实验中被直接观测和量化。

1.2 理论基础：物质波干涉与海森堡绘景

本文的核心理论基础可以概括为物质波干涉原理和海森堡绘景下算符的演化。

物质波透镜（Matter-wave Lens）与傅里叶空间: 协议的基石是物质波显微镜中独特的物质波透镜技术。一个四分之一周期（T/4）的谐振子势能演化，能够实现粒子的实空间（x）和动量空间（p）的互换，即所谓的傅里叶变换。具体来说，当粒子在谐振势中演化T/4时，其波函数在空间域上的分布将变为其初始波函数在动量域上的分布（经过一个相位因子）。这意味着，在一个T/4脉冲之后，系统从实空间图像平面转换到了傅里叶平面。随后，如果再施加一个T/4脉冲，系统则会从傅里叶平面转换回实空间图像平面，但会伴随着放大或缩小，并且可能发生反转。这种x-p互换的能力是进行傅里叶空间操控的前提。

拉曼脉冲（Raman Pulse）进行动量转移和自旋操控: 拉曼脉冲在此方案中扮演了关键角色。它能够诱导原子在两个内部自旋态之间进行跃迁，同时伴随着一个受控的动量踢（momentum kick）。通过选择拉曼光束的波矢和夹角，可以精确控制动量踢的大小和方向。在傅里叶空间中施加拉曼π/2脉冲，可以将一部分原子从物理自旋态（如↑）转移到辅助自旋态（↓），并在此过程中赋予一个特定的动量 q。这个动量 q 在后续的物质波透镜作用下，会转化为实空间中的位移 d = ħq/(mω₁)，其中 m 是原子质量，ω₁ 是物质波透镜的陷阱频率。这使得我们能够“选择性地移动”系统的一部分，而另一部分则保持不变，为构建干涉仪提供了可能。

多体干涉仪方案： 如图1所示，核心方案是一个基于马赫-曾德尔干涉仪原理的多体版本。整个协议分为以下几个步骤：

1. 初始态准备： 将原子制备在光晶格中，处于一个物理自旋态（例如↑）。
2. 傅里叶变换： 关闭晶格和相互作用，施加一个T/4物质波脉冲，将系统从实空间转换到傅里叶空间。
3. 动量和自旋操控： 在傅里叶空间中，施加一个拉曼π/2脉冲。这个脉冲将一部分原子从物理自旋态↑转移到辅助自旋态↓，并赋予其一个精确的动量踢 q。这意味着，现在系统处于一个由两部分组成的叠加态：一部分是未被转移的↑态原子，另一部分是转移到↓态并带有动量 q 的原子。
4. 实空间位移： 施加第二个T/4物质波脉冲，将傅里叶空间的动量 q 转换回实空间的位移 d。此时，带有辅助自旋态↓的原子相对于↑态原子发生了位移。
5. 干涉闭合： 施加第二个拉曼π/2脉冲，但这次不引入动量踢，只是简单地进行自旋态之间的相干叠加。这个脉冲的关键之处在于其相位 φ。
6. 自旋分辨测量： 最后，通过单原子分辨成像技术，测量每个格点上的原子自旋态（↑或↓）。

通过改变第二个拉曼脉冲的相位 φ，我们可以观察到不同格点上自旋态占据数的干涉条纹。这些条纹的对比度和相位偏移直接编码了系统在位移 d 上的非对角关联函数。例如，对于单粒子格林函数 g^(1)(d) = ⟨a†{x+d,↑}a{x,↑}⟩，测量到的密度 ⟨n_{x,↑}(4)⟩ 将表现出 cos(φ) 依赖的干涉项，其振幅正是 g^(1)(d) 的幅值。如果关联函数是复数（例如在存在规范场或破缺时间反演对称性时），相位偏移将给出其复数相位。

海森堡绘景下的算符演化： 为了精确地理解测量结果与关联函数之间的关系，论文采用了海森堡绘景。在这个绘景中，物质波透镜和拉曼脉冲都被视为作用在产生/湮灭算符上的单粒子幺正变换。例如，一个T/4物质波透镜将实空间产生算符 a†{x,σ} 转换为傅里叶空间算符 a†{k,σ}。拉曼脉冲则将不同自旋态的算符耦合起来，并引入动量相移。通过将所有这些幺正变换按时间顺序复合，可以得到最终测量到的密度算符 ⟨n_{x,↑}(4)⟩ 与初始态的各种关联函数之间的精确数学关系。这种方法能够清晰地揭示实验序列如何“构造”出目标关联函数，从而为实验设计和数据分析提供坚实的理论基础。

1.3 技术难点与方法细节

实现“多体相显微镜”面临着诸多技术挑战，但论文也提供了详细的方法细节来解决这些问题。

1.3.1 精确的物质波透镜操控

• 谐振子势能： 物质波透镜的核心是精确构建一个谐振子势能，其频率 ω₁ 和演化时间 T/4 必须严格控制。任何对谐振子势的偏差都可能导致傅里叶变换的不完美，从而引入不需要的相移或失真。
• 放大率： 物质波透镜还涉及放大率问题。虽然论文中假设单位放大率的协议，但实际中可能需要调整放大率以适应不同的实验需求。高倍放大可以提高实空间分辨率，但也可能引入更多像差。
• 晶格势与相互作用的开关： 在T/4脉冲期间，原子应在自由空间或纯谐振势中演化，这意味着晶格势和原子间的相互作用需要被精确、瞬时地关闭。对于费米子，Feshbach共振技术可以用来调控相互作用强度，使其趋近于零。然而，瞬时关闭高强度相互作用对实验控制提出极高要求，任何残余相互作用都可能导致退相干。

1.3.2 拉曼脉冲的精密控制

• 动量踢的精确性： 拉曼脉冲诱导的动量踢 q 必须精确控制，以确保最终实空间位移 d 是晶格间距的整数倍，从而探测到特定的长程关联。这要求拉曼光束的波矢、夹角和频率差都必须极其稳定和精确。例如，对于 ⁶Li 原子，实现10个晶格位移仅需5°的夹角，这意味着对角度的精度要求达到0.05°。
• 相位稳定性： 构成干涉仪的两个拉曼脉冲以及物质波透镜本身，其相对相位必须在整个实验序列中保持极高的稳定性。任何相位漂移都会导致干涉条纹的模糊或错误偏移，从而破坏对关联函数复相位的提取。论文建议将拉曼光束的拍频信号与拉曼脉冲的开始时刻关联起来，以保证相位参考的稳定性。
• 脉冲持续时间： 拉曼脉冲的持续时间必须远短于原子在晶格中的动力学时间尺度和物质波透镜的演化时间（T/4），以确保其操作的“瞬时性”，避免在脉冲期间引入额外的动力学或位移。

1.3.3 d波超导序参量的测量细节 (图2)

• 自旋态扩展： 为了测量d波超导序参量 Cμ,ν(d) = ⟨a†{i+eμ,↑} a†{j+eν,↓} a_{j,↓} a_{i,↑}⟩，方案扩展到使用四个自旋态：两个物理自旋态（↑，↓）和两个辅助自旋态（介，弥）。
• 独立位移和相位： 方案允许对物理自旋↑和↓分别施加独立的动量踢，从而产生独立的实空间位移 d₁ = i - j 和 d₂ = i - j + eμ - eν，以及独立的拉曼相位 φ₁ 和 φ₂。这种多维度控制是提取四点关联函数的关键。
• 四次测量组合： 论文推导指出，通过四次不同拉曼相位组合（φ₁ = φ₂ = 0, π, π/2, -π/2）的密度-密度关联测量，可以精确提取出所需的四点配对关联函数 Cμ,ν(i-j)。这一步骤通过巧妙的相位选择，抵消了不相关项，突出了目标关联函数。
• 对称性编码： d波超导的特征在于其配对波函数具有空间对称性。通过选择不同的 eμ 和 eν，可以区分 s波和d波配对。例如，d波超导中，Cx,y(d) = -Cx,x(d)，而s波超导中 Cx,y(d) = +Cx,x(d)，这可以通过方案的位移选择直接探究。

1.3.4 非等时关联函数与谱函数测量细节 (图3)

• 粒子提取与隔离： 关键步骤是在傅里叶空间中，利用聚焦拉曼光束（其波矢对应于目标动量 ko）将一个粒子从系统的一个特定动量模式 ko 中部分提取到辅助自旋态中。这个被提取的粒子随后在实空间被位移，从而与剩余的多体系统在空间上隔离。
• 时间演化： 随后，被提取的粒子（通常只有一原子）在隔离状态下，其相互作用可以忽略，但会经历一个由单粒子能带结构决定的动态相位演化。同时，剩余的多体系统（“挖掉”一个动量模式的系统）在晶格中重新加载相互作用，并进行一段预设的动力学演化时间 t。
• 干涉闭合与谱函数： 在演化时间 t 结束后，被提取的粒子再次通过物质波透镜转换回傅里叶空间，并通过第二次拉曼脉冲与多体系统进行干涉。测量到的干涉条纹的对比度和相位，将直接给出非等时关联函数 G(ko, t) = -iθ(t)⟨a_{ko,↑}(t)a†_{ko,↑}(0)⟩。通过对时间 t 的傅里叶变换，即可得到谱函数 A(k,ω)。
• 带映射的精细处理： 在非等时关联函数协议中，晶格势的开关和粒子重加载需要通过“带映射”（band mapping）技术精心处理。当瞬时关闭晶格时，格子的准动量态会映射到连续空间的动量态。为了保持相干性，需要补偿带映射过程中引入的动态相位。这要求在傅里叶平面注入一个合适的相位。或者，也可以在关闭晶格后将粒子加载到辅助谐振阱中，避免复杂的带映射过程。

1.3.5 隐藏非对角序的测量细节 (图4)

• 局部干涉： 与图1的全局干涉不同，图4的协议中，第二个拉曼脉冲仅通过光镊聚焦到单个或少数格点上进行局部作用。这意味着只有这些局部格点会与位移后的系统发生干涉。
• 组合测量： 为了探测分数量子霍尔态中复合玻色子的隐藏序 ğ^(1)(i,j) = ⟨Π_{n≠i,j} (z_i-z_n / z_j-z_n)^{-1/ν} a†_i a_j⟩，需要将位移 d 处的单粒子相干性测量 ⟨a†_i a_j⟩ 与所有其他格点 n≠i,j 的局部占据数 {n_n} 测量结果结合起来。这通过数据后处理实现。
• 数据处理： 未被局部拉曼脉冲作用的格点，其辅助自旋态原子在数据分析时通过空间位移反向移动，叠加回原始位置，以重构总密度。只有被局部拉曼脉冲作用的格点（红色框区域）才用于提取相干信息。这使得我们能够在一个实验序列中同时收集相干性和密度信息，从而构建出复杂的非局域关联函数。

综上所述，这项研究在理论上提出了创新的多体干涉仪方案，并对每个协议的关键技术细节和挑战进行了深入分析。通过精密的物质波操控和拉曼脉冲序列设计，它有望为量子多体物理提供前所未有的实验观测能力，直接探究量子系统的深层相干性和非对角关联，为理解最前沿的量子物质状态铺平道路。

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

该研究提出的“多体相显微镜”协议并非仅仅停留在理论构想，而是基于冷原子实验中已有的高精度操控技术，并对实验参数和预期性能进行了详细的基准测试和评估。论文中列举的实验参数和性能指标，是协议在实际系统中可实现性的重要支撑。以下将围绕关键 benchmark 体系、计算所得数据和性能数据进行深入解析。

2.1 实验体系与基本参数选择

论文考虑了两种主要的冷原子体系，以适应不同的物理场景：

• 玻色子体系： 针对图1（一般相干性测量）和图4（分数量子霍尔态隐藏序测量）的协议，作者假设使用玻色子 133Cs 原子。133Cs 原子具有较大的质量，且可以通过Feshbach共振精确调控相互作用，这在物质波透镜期间关闭相互作用至关重要。
• 费米子体系： 针对图2（d波超导序参量测量）和图3（非等时格林函数与谱函数测量）的协议，作者假设使用费米子 6Li 原子。6Li 原子是研究Hubbard模型和高温超导的理想平台，其轻质量特性也为快速动力学实验提供了便利。

在两种体系中，共同的晶格间距（alat）均设定为 500 nm。这是一个在光晶格实验中常见的间距，既能保证足够的原子分辨率，又能实现相对强的相互作用。

2.2 物质波透镜与拉曼操控参数

物质波透镜的性能直接影响傅里叶变换的质量和实空间位移的精确性。拉曼脉冲的参数则决定了动量转移和相位控制的精度。

• 陷阱频率（ω₁）：
  • 对于 133Cs 原子：设定 ω₁ = 2π × 50 Hz。这意味着四分之一周期（T/4）的演化时间约为 5 ms。
  • 对于 6Li 原子：设定 ω₁ = 2π × 250 Hz。这意味着 T/4 演化时间缩短到 1 ms。更高的陷阱频率允许更快的物质波透镜操作，这对于需要长时间演化的非等时关联函数测量尤为关键。
• 拉曼波长（λ）：
  • 对于 133Cs 原子：拉曼波长选择在 870 nm 附近。
  • 对于 6Li 原子：拉曼波长选择在 671 nm 附近。
  • 拉曼波长决定了动量踢的基准大小，以及实现特定动量转移所需的光束夹角。
• 动量转移与实空间位移（d）：
  • 拉曼光束之间的夹角 θ 控制动量转移的模量 q = (2π/λ)√(2(1-cos(θ)))。这个动量 q 在物质波透镜作用下，会转化为实空间位移 d = ħq/(mω₁)。
  • 133Cs 体系的位移能力： 实现 10 个晶格间距的位移（n=10），对于 133Cs 原子只需约 30° 的拉曼光束夹角。这表明在中等角度下就能实现较大的位移。
  • 6Li 体系的位移能力： 实现 10 个晶格间距的位移，对于 6Li 原子需要非常小的拉曼光束夹角，仅约 5°。这对于光学设置提出了更高的精度要求，因为 0.05° 的角度精度对应于一个晶格间距的一部分位移。但这同时也表明，即便使用相对较小的夹角，也能在轻原子体系中实现显著的位移。
• 拉曼脉冲持续时间： 脉冲持续时间必须短于多体动力学时间尺度和物质波透镜的演化时间。论文指出，对于 ω₁ = 2π × 250 Hz（T/4 = 1 ms），可实现的拉曼脉冲持续时间短至 200 ns，这满足了瞬时操作的要求，确保在脉冲期间不会引入显著的粒子位移或动力学。

2.3 相互作用控制与带映射

相互作用关闭： 在物质波透镜操作期间，原子之间的相互作用必须关闭。这是协议成立的关键，因为它依赖于单粒子物理。

• 133Cs 原子： 论文建议将磁场淬灭到两个零交叉点之间，以大幅降低相互作用强度。此外，可以快速关闭垂直限制势，减少密度，并在垂直方向也引入物质波条件。
• 6Li 原子： 对于 6Li，快速磁场淬灭到零可能不切实际。替代方案是利用快速拉曼自旋翻转将原子转移到弱相互作用的超精细态。这种方法利用不同自旋态之间Feshbach共振的偏移特性。另外，在物质波序列期间，可以将磁场缓慢调整到接近零，以进一步降低相互作用强度。

带映射（Band Mapping）： 在非等时关联函数测量中，当粒子从隔离态重新加载到晶格时，需要进行带映射。

• 时间尺度： 晶格势的缓变上升（或下降）时间应慢于单粒子带隙，以确保绝热转换。对于 6Li 原子，这个时间尺度约为 100 µs。这个时间需要比物质波透镜的持续时间（T/4 = 1 ms）快得多，以避免不必要的延迟。
• 相位补偿： 带映射过程中会引入动态相位 ϕ(κ)。为了保持相干性，这个相位需要通过在傅里叶平面注入一个合适的相位来补偿。这要求精确校准带映射过程。

2.4 动量模式选择与空间隔离

• 聚焦拉曼光束选择动量模式： 在非等时关联函数协议中，通过在傅里叶空间使用聚焦拉曼光束，可以选择一个特定的动量模式 ko。光束腰半径 w₀ 决定了选择模式的宽度。
• 宽度参数： 论文建议对于 6Li 原子，系统长度为 20 个晶格间距时，选择对应于倒晶格矢量 1/20 的动量模式，所需的 w₀ 约为 4 μm。这个光束腰大小可以确保选择一个足够窄的动量模式，同时又足够宽以包含一个原子。
• 空间隔离距离： 隔离的动量模式被移位到 n_{iso} 个晶格间距之外。论文提到，对于 30° 的拉曼光束夹角，可以实现 n_{iso} = 68 个晶格间距的位移。这个距离足够大，可以确保被提取的粒子在多体演化期间与剩余系统保持空间隔离，从而相互作用可以忽略。

2.5 性能预期与统计量需求

• 信噪比与光谱分辨率： 对于非等时关联函数测量，由于采用了实时方法而非线性响应谱学，预期能获得良好的信噪比和高光谱分辨率。这对于解析精细的激发谱特征至关重要。
• d波超导关联函数的统计量： 论文引用了数值模拟结果，指出在存在对角隧穿项 t’ 的情况下，d波超导关联 Cμ,ν(d) 在零温度下会在几个晶格间距内衰减到 10⁻⁴ 左右的常数值。对于有限温度下的代数衰减，检测如此微弱的信号需要非常高的统计量。作者估计，为了分辨普遍部分的长程关联，可能需要 10⁴ 次测量。即使对于短程、信号更强的非普遍部分，也能通过这些测量揭示配对对称性。
• 隐藏序的统计量： 对于分数量子霍尔态的隐藏非对角长程序测量，基于数值研究，作者预期只需几千次快照即可实现。这比d波超导的测量要求略低，表明该协议对于探测此类拓扑序具有较高的可行性。

这些详细的参数和性能评估表明，本文提出的多体相显微镜协议在当前冷原子实验技术的范畴内是可行的。虽然需要极高的实验精度和稳定性，但其潜在的物理发现价值是巨大的。通过精确控制物质波透镜、拉曼脉冲和相互作用，并结合高效的数据采集和分析，我们有望首次直接探究过去难以触及的量子相干性和非对角关联，为理解最前沿的量子多体现象提供强有力的新工具。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

尽管提供的论文主要聚焦于实验协议和理论分析，并未直接提供具体的代码实现或开源软件库链接，但作为面向量子化学科研人员的深度解析，我们可以根据论文描述的实验流程和数据处理需求，推断出实现这些协议所需的软件、硬件接口和数据分析策略。本节将从这个角度，探讨可能的代码实现细节、复现指南，并建议相关的软件包和潜在的开源方向。

3.1 推断的代码实现细节

一个完整的“多体相显微镜”实验控制与数据分析系统，将是一个复杂的软硬件集成项目，涵盖从低层硬件控制到高层物理量提取的多个层次。

3.1.1 硬件控制层

这一层主要负责与实验硬件进行实时、精确的交互，生成和同步各种实验脉冲。

• 激光系统控制：
  • 拉曼光束生成： 需要通过任意波形发生器（AWG）控制声光调制器（AOM）或电光调制器（EOM），精确生成和整形拉曼脉冲的频率、强度、持续时间和相对相位。例如，为了实现 200 ns 的拉曼脉冲，AWG 的采样率需要达到 GS/s 级别。
  • 晶格势控制： 晶格激光的强度和持续时间也需要 AOM 精确控制，实现快速的开关和缓变（带映射）。
  • 光镊控制： 如果协议需要局部拉曼作用（如隐藏序测量），则光镊激光的强度和位置也需精确控制。
• 磁场控制：
  • Feshbach共振线圈： 需要高精度电流源驱动Feshbach线圈，实现磁场的快速淬灭或缓变，以调控原子相互作用。
• 物质波透镜控制：
  • 谐振陷阱： 快速切换到谐振陷阱，并维持其频率 ω₁ 和演化时间 T/4 的稳定性。这可能涉及可编程数字模拟转换器（DAC）控制电极或额外的激光势。
• 触发与同步：
  • 所有这些硬件的动作都必须通过一个中央时序控制器（例如基于FPGA或PXI系统）进行精确同步，确保物质波透镜、拉曼脉冲、磁场变化的相对时序和相位关系严格符合协议要求。例如，拉曼脉冲的相位需要与拍频信号的某一特定点同步。

3.1.2 数据采集层

这一层主要负责从实验中获取原始数据，即单原子分辨率的图像。

• EMCCD/sCMOS相机控制： 相机需要精确触发，以在实验序列的末尾获取原子图像。这包括控制曝光时间、增益等参数。
• 图像传输与存储： 高速传输采集到的图像到计算机进行初步处理和存储。由于需要大量快照（例如 10⁴ 次），数据吞吐量是重要考量。

3.1.3 数据处理与分析层

这是将原始图像数据转化为物理量的核心。这一层通常在高性能计算机上执行。

• 图像预处理：
  • 降噪： 应用图像滤波算法（如高斯滤波）去除相机噪声。
  • 背景扣除： 从图像中减去背景光信号。
  • 原子识别： 基于图像中原子发出的荧光信号，开发鲁棒的原子识别算法。这通常涉及阈值分割、连通域分析或机器学习方法。每个原子会显示为一个亮点，需要精确识别其在晶格中的位置。
• 原子占据数提取：
  • 将识别出的原子位置映射回光晶格的格点位置。
  • 对于每个格点，统计是否有原子存在（二值化），以及其自旋态。这通常需要自旋分辨成像技术，例如通过梯度磁场将不同自旋态的原子在空间上分离，或通过拉曼侧带冷却选择性地移除一个自旋态的原子，然后进行差异成像。
• 关联函数计算：
  • 平均密度： 计算每个格点上原子占据数的平均值 ⟨n_x⟩。
  • 相干性提取： 根据论文中的公式，通过拟合干涉条纹（作为拉曼相位 φ 的函数）来提取关联函数 g^(1)(d) 或 Cμ,ν(i-j) 的幅度和相位。这通常涉及正弦拟合：A cos(φ - φ₀) + B，其中 A 是关联函数的幅值，φ₀ 是复数相位。
  • 高阶关联： 对于d波超导序参量和隐藏序，需要计算四点关联函数或涉及乘积的复杂关联函数。这要求将多个格点上的占据数和自旋信息组合起来。
• 谱函数计算（非等时关联）：
  • 对于非等时关联函数 G(ko, t)，需要在不同演化时间 t 下重复实验，并提取 G(ko, t)。
  • 然后对 G(ko, t) 进行傅里叶变换，得到谱函数 A(ko, ω)。这可能需要数值稳定化技术，如最大熵方法，因为时间数据通常是离散且有噪声的。

3.2 复现指南

复现这些协议需要一个高度集成且精密的冷原子实验平台，并结合专门的软件开发。以下是复现的关键步骤和建议：

1. 平台搭建与校准：

   • 冷原子气体准备： 搭建 133Cs 或 6Li 冷原子实验平台，实现原子冷却、囚禁、加载到光晶格。
   • 物质波显微镜校准： 精确校准物质波透镜，包括谐振阱频率 ω₁、T/4 演化时间、成像放大率和失真。
   • 拉曼系统校准： 标定拉曼激光器的频率、强度、光束夹角和相对相位控制。精确测量给定夹角下的动量踢 q，以及其转化为实空间位移 d 的关系。
   • 相互作用调控： 校准Feshbach线圈或拉曼自旋翻转序列，确保在物质波透镜期间能够有效关闭相互作用。
   • 单原子分辨率成像： 优化成像系统，达到单原子、单格点和自旋分辨的成像能力。

2. 协议模块化实现：

   • 基本干涉仪（图1）: 首先实现最简单的单粒子格林函数 g^(1)(d) 测量，验证物质波透镜和拉曼干涉的基本功能和相位稳定性。这是后续所有复杂协议的基础。
   • d波超导（图2）: 在基本干涉仪的基础上，引入两个辅助自旋态和独立的拉曼操控通道，实现对 d₁ 和 d₂ 的独立控制，并根据论文建议的四种相位组合进行测量。
   • 非等时关联（图3）: 实现粒子提取、空间隔离、多体时间演化和干涉闭合。重点在于带映射的相位补偿和时间演化期间多体系统的稳定性。需要对不同的演化时间 t 进行多次测量。
   • 隐藏序（图4）: 引入局部拉曼脉冲，并开发复杂的数据分析管道，将局部相干性测量与全局密度测量结合起来。

3. 数据分析与验证：

   • 统计量积累： 根据预期信号强度，进行足够次数的快照测量。例如，d波超导可能需要 10⁴ 次测量。
   • 误差分析： 对每次测量结果进行统计分析，评估系统误差和随机误差，确保结果的可靠性。
   • 理论比较： 将实验测量得到的关联函数与理论计算或数值模拟结果进行比较，验证协议的有效性。

4. 稳定性与自动化：

   • 考虑到实验的复杂性和所需的大量快照，自动化实验流程和提高系统长期稳定性是成功的关键。这包括温度、激光功率、磁场、真空环境等所有参数的稳定。

3.3 所用的软件包及开源 Repo Link (推断与建议)

由于论文未提供具体代码，以下列出的软件包是基于冷原子实验控制和数据分析的常见实践进行推断的。鼓励未来的实现者开源其代码，以促进科学发展。

3.3.1 实验控制与自动化

• LabVIEW (National Instruments): 在许多冷原子实验室中，LabVIEW 仍然是用于仪器控制、数据采集和实验序列编程的常用图形化开发环境。它与 NI 的 PXI/CompactRIO 硬件结合，可以实现高精度定时和多通道控制。
• Python (with pyvisa, qcodes, labscript etc.): 随着科学软件生态的发展，Python 已经成为越来越受欢迎的实验控制语言。pyvisa 用于与各种仪器（示波器、信号发生器等）进行通信。qcodes 或 labscript 等框架提供了一个结构化的方式来编写复杂的实验序列和数据管理。
• Julia (for real-time control where performance is critical): 对于极度要求实时性能的控制任务，如快速反馈循环，Julia 语言因其高性能和易用性，也可能成为一个有吸引力的选择。

3.3.2 数据采集与图像处理

• Andor SDK / Basler Pylon SDK 等相机供应商提供的 SDK： 用于直接控制 EMCCD/sCMOS 相机，优化图像采集参数。
• Python (with OpenCV, scikit-image, numpy, scipy):
  • numpy 和 scipy：进行高效的数值计算和科学数据处理。
  • OpenCV (Open Source Computer Vision Library)：进行图像处理，如滤波、阈值分割、连通域分析等，以识别原子。
  • scikit-image：提供了更多的图像处理算法，特别是用于科学图像分析的功能。

3.3.3 数据分析与可视化

• Python (with matplotlib, seaborn, pandas):
  • matplotlib 和 seaborn：生成高质量的图表，可视化干涉条纹、关联函数、谱函数等结果。
  • pandas：进行数据管理和结构化存储，便于处理大量的实验数据。
  • lmfit 或 scipy.optimize：用于对干涉条纹进行非线性拟合，提取振幅和相位。
• Julia / MATLAB / R： 同样可以用于数据分析和可视化，具体取决于实验室的偏好和已有工具链。

3.3.4 数值模拟与理论计算（辅助）

• Python (with QuTiP, TenPy, Julia for high-performance computing):
  • 为了验证实验结果和指导实验设计，通常需要进行数值模拟，例如通过精确对角化（Exact Diagonalization）、密度矩阵重整化群（DMRG）或蒙特卡洛（Monte Carlo）方法模拟 Hubbard 模型。
  • QuTiP：一个用于量子信息和开放量子系统演化的Python库。
  • TenPy：一个基于张量网络（Tensor Network）的Python库，适用于模拟一维和准二维强关联系统。

3.3.5 开源 Repo Link (建议)

论文中未提供任何开源代码库链接。 这是当前冷原子领域的一个常见现象，但随着“可复现性危机”的凸显，越来越多的研究者开始重视开源。如果未来有团队成功复现了这些协议或基于此开发了相关工具，强烈建议他们将代码开源到GitHub或类似的平台上，例如：

• github.com/YourLab/ManyBodyPhaseMicroscope-Control (用于硬件控制和实验序列)
• github.com/YourLab/ManyBodyPhaseMicroscope-Analysis (用于图像处理和数据分析)
• github.com/YourLab/ManyBodyPhaseMicroscope-Simulations (用于数值模拟验证)

通过开源代码，可以极大地降低其他团队的复现门槛，加速科学发现的进程。这将有助于将这些前沿的协议转化为社区共享的强大工具。

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

这项工作站在了冷原子量子模拟领域的巨人之肩上，广泛引用了该领域最前沿和最基础的研究成果。同时，任何开创性的研究都不可避免地存在其自身的局限性，这些局限性既是当前技术的瓶颈，也是未来研究的方向。

4.1 关键引用文献解析

论文的引用列表反映了其跨越量子气体显微镜、超导物理、谱学和拓扑物理多个前沿领域的特性。以下是几个对理解这项工作至关重要的引用类别：

1. 量子气体显微镜的奠基性工作 ([1-3]):

   • [1] W. S. Bakr et al., Science 329, 547 (2010): 这篇论文首次展示了基于光学晶格的单原子分辨率成像，开辟了量子气体显微镜的新时代。
   • [2] J. F. Sherson et al., Nature 467, 68-72 (2010): 另一篇同期发表的开创性工作，展示了对Hubbard模型的单格点分辨测量。
   • [3] C. Gross and W. S. Bakr, Nature Physics 17, 1316 (2021): 对量子气体显微镜近年发展和应用的综述或里程碑式工作。 这些工作为本文提出的“多体相显微镜”提供了最基本的实验平台和技术支撑。

2. 物质波透镜与傅里叶空间操控 ([16, 17, 18]):

   • [16] L. Asteria et al., Nature 599, 571 (2021): 这篇论文首次实现了物质波透镜的傅里叶空间成像能力，是本文所有协议的核心技术基础。它证明了通过物质波透镜可以实现实空间和动量空间的转换。
   • [17] S. Brandstetter et al., Phys. Rev. Lett. 135, 103401 (2025): 可能进一步发展了物质波透镜技术，例如高放大率或与其他技术的结合。
   • [18] L. Asteria, PhD thesis (2022): 作者之一的博士论文，很可能包含了物质波透镜的详细理论和实验细节。 这些工作直接证明了本文提出的傅里叶空间操控的可行性。

3. 局部电流与相位的早期探索 ([13, 14, 15]):

   • [13] A. Impertro et al., Phys. Rev. Lett. 133, 063401 (2024): 可能展示了局域电流算符的测量。
   • [14] A. Impertro et al., Nature Physics 21, 895 (2025): 进一步发展了局域电流测量方法。
   • [15] J. C. Brüggenjürgen et al., arXiv:2410.10611 (2024): 提出了一个用于量子气体的相位显微镜，与本文思想有共通之处，但可能侧重于局域相位涨落。 这些工作是本文之前，对超越密度测量、获取相位信息的初步尝试，为本文的更宏大目标提供了铺垫。

4. d波超导与Fermi-Hubbard模型 ([19, 20, 21, 28]):

   • [19] C. Roth et al., arXiv:2511.07566 (2025): 关于神经量子态在Fermi-Hubbard模型中识别d波超导。
   • [20] M. Xu et al., Nature 642, 909 (2025): 可能在Fermi-Hubbard体系中实现了新的冷却技术，达到了低温区域。
   • [21] T. Chalopin et al., arXiv:2412.17801 (2024): 探索了赝能隙机制。
   • [28] H. Xu et al., Science 384, eadh7691 (2024): 强调了对角隧穿项 t’ 对d波超导稳定的重要性，这是本文强调传统映射失效的原因。 这些论文代表了当前在Fermi-Hubbard模型和d波超导研究上的前沿进展，为本文提出的测量协议提供了重要的物理背景和验证动机。

5. 非等时关联函数与谱函数 ([29, 34, 35, 39, 40-42]):

   • [29] P. T. Brown et al., Nature Physics 16, 26 (2020): 使用量子气体显微镜测量谱函数，但可能面临统计量挑战。
   • [34] A. Damascelli et al., Rev. Mod. Phys. 75, 473 (2003) & [35] J. A. Sobota et al., Rev. Mod. Phys. 93, 025006 (2021): ARPES的综述，提供了固态系统中谱函数测量的标准。
   • [40-42] M. Cetina et al., Science 354, 96 (2016); L. A. Pena Ardila et al., Phys. Rev. A 99, 063607 (2019); M. G. Skou et al., Nature Physics 17, 731 (2021): 使用拉姆齐干涉法测量量子杂质或极化子的格林函数。 这些工作展示了在冷原子系统中测量谱函数和格林函数的不同尝试和挑战，凸显了本文方法在直接性、全动量分辨和实时性方面的优势。

6. 隐藏序与分数量子霍尔态 ([43, 44, 46, 47]):

   • [46] N. Read, Phys. Rev. Lett. 62, 86 (1989): Read的经典工作，将分数量子霍尔态的拓扑序与复合玻色子的非对角长程序联系起来。
   • [47] F. Pauw et al., Phys. Rev. Res. 6, 023180 (2024): 一篇关于探测分数量子霍尔态隐藏序的最新理论工作，可能与本文提出的局部干涉协议紧密相关。 这些论文为本文对拓扑相中隐藏序的探测提供了理论框架和实验灵感。

4.2 对这项工作的局限性评论

尽管本文提出的协议具有巨大的潜力，但其实现和应用也面临一系列不可忽视的局限和挑战：

1. 极高的实验复杂性与精度要求：

   • 多步骤序列： 所有协议都涉及多个物质波透镜和拉曼脉冲的序列，每个步骤都必须精确控制其时间、强度、频率和相位。任何一个环节的微小偏差都可能导致整个干涉仪的失败。
   • 相位稳定性： 物质波干涉对相位漂移极其敏感。在整个实验序列（可能持续数毫秒到数十毫秒）中，所有激光器和磁场的相对相位必须保持在远低于 π 的水平，这是一个巨大的技术挑战。
   • 瞬时开关： 晶格势和相互作用的“瞬时”开关并非易事。快速变化的磁场可能引入涡流或加热，而激光脉冲整形也存在带宽限制。这些非理想因素可能导致原子偏离理论预测的单粒子演化轨迹，从而引入退相干。

2. 相互作用控制的非理想性：

   • Feshbach共振的限制： 尽管Feshbach共振可以调控相互作用，但在极短时间内将相互作用精确归零，且不引入额外损失或加热，仍然是挑战。尤其是在多自旋态系统中，需要仔细选择磁场点，以最小化所有自旋通道的相互作用。
   • “隔离”的近似： 在非等时关联函数协议中，被提取的粒子被假设与剩余系统“隔离”且相互作用可忽略。然而，完全的隔离可能难以实现，特别是在有限大小的阱中。残余相互作用可能导致被提取粒子的非平凡动力学。

3. 带映射的挑战：

   • 绝热性要求： 带映射需要缓变晶格势，以确保准动量态的绝热演化。但协议要求在物质波透镜（毫秒级）与带映射（微秒级）之间进行切换，这在时间尺度上存在冲突。虽然论文建议通过在傅里叶平面注入相位来补偿，但这要求对带映射过程的精确表征和控制。
   • 相位相干性： 在带映射过程中，需要保持原子波函数的相位相干性，这意味着晶格势的相位也必须是稳定的。

4. 统计量与测量时间：

   • 微弱信号： 特别是对于长程d波超导关联和某些拓扑序的隐藏关联，信号强度可能非常微弱（例如 10⁻⁴）。为了从噪声中提取这些信号，需要进行数千甚至上万次的快照测量。这会使单个实验序列的总时间非常长，受限于原子寿命、系统漂移和实验稳定性。
   • 单次快照信息量： 尽管论文提到可以通过重叠区域的密度差从单个快照中提取相干性，但对于非局域和低信噪比的关联，仍然需要大量的平均。

5. 理论近似与模型限制：

   • 理想哈密顿量： 协议的理论推导假设了一个理想的哈密顿量和幺正演化。实际系统中，存在各种非均匀性、无序和损耗，这些都可能偏离理想情况。
   • 有限尺寸效应： 实验在有限尺寸的晶格上进行，而理论有时假设无限大系统。有限尺寸效应可能改变关联函数的衰减行为和临界性质。

6. 推广到更复杂系统的挑战：

   • 尽管协议原则上通用，但将其推广到更复杂的多组分、多轨道或具有长期非局域相互作用的系统时，可能需要更复杂的操控和额外的修正。

综上，这项研究在概念上非常强大和创新，但其实验实现是一个巨大的工程挑战，需要冷原子实验技术的进一步发展和完善。未来的工作需要专注于克服这些具体的实验瓶颈，从而充分发挥“多体相显微镜”的巨大潜力。

5. 其他你认为必要的补充

5.1 开启量子模拟的新篇章：从现象到机制的桥梁

“多体相显微镜”的提出，无疑为冷原子量子模拟领域开启了一个全新的篇章。传统量子气体显微镜的优势在于其对实空间原子占据数和自旋的直接观测，这对于理解密度波、自旋序等“对角”物理现象至关重要。然而，许多更深层次的量子物态，例如超流、超导、拓扑序等，其本质特征在于非对角长程关联和量子相干性。这些现象的微观机制，往往隐藏在粒子波函数之间的相对相位中，是仅凭密度测量无法完全捕捉的。

本文提出的协议，正是架起了从宏观量子现象到微观量子机制之间的桥梁。通过直接测量格林函数、配对序参量和复合玻色子关联函数，实验者将能够：

1. 直接验证理论预测： 许多理论模型对量子相变和相图做出了精确的预测，特别是关于序参量的形式、空间衰减行为和动力学。多体相显微镜将提供直接的实验数据来验证这些预测，甚至可能发现理论未曾预料到的新现象。
2. 揭示隐藏的量子序： 对于分数量子霍尔效应等拓扑相，其序参量是非局域的，难以用传统方法测量。该协议通过巧妙的局部-全局测量组合，有望揭示这些“隐藏”的拓扑序，从而加深我们对拓扑物质本质的理解。
3. 实时探测动力学： 非等时格林函数的测量允许我们实时追踪粒子的激发和传播，从而直接获取谱函数。这对于研究量子淬火（quantum quench）、热化过程和非平衡态动力学具有重要意义，能够直接观察到例如极化子形成、多体局域化（Many-Body Localization, MBL）等复杂动力学过程。

5.2 潜在应用领域与未来展望

这项研究的深远影响不仅限于冷原子物理本身，更可能辐射到凝聚态物理、量子信息科学等多个领域。

5.2.1 凝聚态物理的前沿探索

• 高温超导机理： 理解铜氧化物高温超导的微观机理是凝聚态物理学的圣杯。Fermi-Hubbard模型被认为是描述这一现象的关键模型。本文提出的d波超导序参量测量协议，将为直接研究Hubbard模型中的d波配对提供前所未有的实验工具。通过探测配对波函数的空间分布、衰减行为和对称性，可以深入了解d波超导的起源，并探索新材料的潜力。
• 拓扑物质态： 除了分数量子霍尔效应，许多其他拓扑物质态也以其独特的非局域量子序为特征。例如，在1D拓扑系统中，Majorana零模的相干性可以作为拓扑相存在的标志。本文的协议原则上可以扩展到探测这些Majorana零模之间的非局域关联，甚至可能在多维拓扑系统中寻找更复杂的拓扑序。
• 量子磁性与自旋液体： 在量子磁性系统中，自旋液体等奇特物态也可能具有非对角关联或隐藏序。多体相显微镜有望帮助识别和表征这些复杂磁序，推动量子磁性研究。

5.2.2 量子信息与非平衡态动力学

• OTOCs与量子混沌： 论文提及的一个重要未来方向是扩展协议以测量负时间演化（通过反转哈密顿量的符号），从而访问非等时序关联函数（Out-of-Time-Order Correlators, OTOCs）。OTOCs被认为是量子混沌和信息搅动（information scrambling）的诊断工具，在量子计算和黑洞物理等领域具有重要意义。通过测量OTOCs，可以实验性地区分量子混沌系统和可积系统，探究多体局域化（MBL）和热化之间的转变。
• 纠缠熵测量： 论文最后提到通过Hong-Ou Mandel干涉，从两个系统副本中提取纠缠熵。虽然细节未在此处展开，但物质波显微镜的相干性和傅里叶空间操控能力，原则上为测量非局域的纠缠量提供了可能性，这对于量子信息科学至关重要。

5.2.3 广义量子模拟工具箱

• 传输系数提取： 结合对局域电流的测量（如参考文献[13, 14]），以及本文提出的非局域关联测量，可以构建更全面的图像，从而提取传输系数（如电导率、热导率等），这对于理解材料的输运性质至关重要。
• 连续系统与非晶格系统： 尽管协议主要针对晶格系统设计，但其核心思想——物质波干涉与傅里叶空间操控——原则上可以推广到连续或非晶格系统，例如均匀费米气体或玻色气体，从而研究其长程关联。

5.3 对实验和理论界的启示

这项研究不仅提出了新的实验协议，也对冷原子实验和理论研究产生了深远影响：

• 实验技术推动： 协议对物质波透镜的精密性、拉曼光束的相位稳定性、Feshbach共振的瞬时开关、带映射的相干控制提出了极高要求。这些挑战将推动冷原子实验技术向更高精度、更高稳定性方向发展。例如，对拉曼光束角度 0.05° 的精度要求，将促使光学系统和校准方法更加精细化。
• 新理论工具的验证： 许多关于超导、拓扑和非平衡态的理论模型，可以通过这些协议进行直接的实验验证，从而加速理论的发展。例如，张量网络态（Tensor Network States）或量子蒙特卡洛（Quantum Monte Carlo）模拟的结果，将有直接的实验对应物。
• 多学科融合： 这项工作体现了冷原子物理与凝聚态物理、量子信息科学的深度融合。通过冷原子平台解决凝聚态物理的核心问题，同时探索量子信息前沿，将激发跨学科的合作与创新。

5.4 与其他方法的比较及独特优势

与现有的冷原子测量方法相比，多体相显微镜具有独特的优势：

• 噪声关联 (Noise Correlations)： 噪声关联是测量配对关联的强大工具，但它通常提供的是动量空间中关联函数的幅值信息，而非直接的相位信息，且通常需要复杂的反演才能得到实空间信息。本文的方法直接在实空间或傅里叶空间中通过干涉测量相位和幅值。
• 布拉格谱学 (Bragg Spectroscopy)： 布拉格谱学能够探测动量空间中的激发谱，但缺乏实空间分辨率，且通常是集成式的测量。本文的非等时关联函数协议提供了动量分辨的谱函数，且潜在地可以结合实空间信息（通过选择性提取粒子）。
• 拉姆齐干涉 (Ramsey Interferometry) 测量格林函数： 之前的工作（如[40-42]）利用拉姆齐干涉测量量子杂质或极化子的格林函数。本文将这种思想推广到没有辅助杂质的体系统，且具有全动量分辨能力。
• 线性响应谱学： 传统的谱学方法通常基于线性响应理论，这在强关联系统中可能不完全适用，并可能限制信号强度。本文的实时方法避免了这些限制，有望提供更好的信噪比和光谱分辨率。

总而言之，“多体相显微镜”代表了冷原子量子模拟技术的一个重大飞跃，它将允许研究者直接深入探测量子多体系统的深层奥秘，从而加速对高温超导、拓扑物质和量子混沌等一系列最复杂、最迷人的量子现象的理解。虽然实验挑战巨大，但其潜在的科学回报无疑是巨大的，预示着冷原子领域激动人心的未来。