Estimating applied potentials in cold atom lattice simulators

技术摘要：冷原子晶格模拟器中应用势能的估算

问题陈述
冷原子在光晶格中的应用使其成为通用的量子模拟平台，能够通过量子气体显微镜（QGM）实现具有位点解析度读取能力的哈伯德模型（Hubbard models）。一个关键目标是实现任意的、位点依赖的光势，以支持量子模拟和信息处理的应用。然而，由于光晶格的衍射极限（其中晶格间距通常低于衍射极限），精确施加和校准这些势能仍具有挑战性。

在实际操作中，实验者在试图施加目标势 $V$ 时，往往实现的是一个略有不同的势 $W$。两者之差 $V - W$ 构成了一种系统偏差，而非波动噪声。这种差异会导致模拟误差随时间累积，从而阻碍高精度量子模拟的实现。虽然从动力学中学习哈密顿量提供了一种解决方案，但现有的提议通常需要复杂的态制备、多局部基底下的测量，或对多体动力学的经典模拟，而这在目前的超冷原子装置中难以实现，因为目前的测量通常局限于固定的粒子数基底（粒子位置）。

方法论
作者提出了一种简单且高效的实验方案，仅利用粒子数读出即可估算施加的势能（及跳跃率）。该方案依赖于两个关键的物理与计算要素：

1. 相互作用抑制：该方案利用某些原子种类可以通过费什巴赫共振（Feshbach resonance）关闭相互作用的能力，将系统简化为非相互作用（自由费米子）模型。这使得系统的随时间演化在经典计算上是高效的。
2. 数据采集：实验涉及制备一个已知的初始态（具体为电荷密度波），允许系统在未知的哈密顿量下演化若干时间 $t_j$，并使用 QGM 收集位点占据情况的快照。

作者提出了两种从随时间变化的局部密度 $D_i(t) = \langle c^\dagger_i(t)c_i(t) \rangle$ 重构势能的方法：

• 严谨解析重构法：该方法通过拟合测量数据中的多项式，通过相关矩阵在 $t=0$ 处的二阶导数来推断哈密顿量参数。虽然在理论上是成立的，但数值分析表明，该方法对统计噪声高度敏感，若要实现准确重构，需要极大量的快照（$R \sim 10^7$）。
• 高效启发式方案：为了克服噪声敏感性，作者提出通过最小化代价函数来进行梯度优化。代价函数 $K(v; C(0))$ 衡量了基于试验势 $v$ 预测的占据数与测量数据之间的平方差。该方法避免了计算高阶导数，并利用标准的梯度下降算法（如 L-BFGS）。

主要贡献与结果
论文证明了该启发式方案具有鲁棒性、可扩展性且在实验上可行。数值实验的关键发现包括：

• 最优参数：重构误差在最大演化时间 $t_{max} \approx 5$（以跳跃时间为单位）和较少的采样时间点（$S \gtrsim 5$）时趋于饱和。进一步增加这些参数带来的改进微乎其微。
• 可扩展性：平均重构误差遵循 $\epsilon_{MRE} \sim 1/\sqrt{M}$ 的比例关系（其中 $M$ 是总快照数），并且在很大程度上与系统尺寸 $N$ 无关。对于高达 100 个位点的系统，使用大约 $3 \times 10^4$ 个快照即可实现 $\le 5\%$ 的平均重构误差。
• 对误差的鲁棒性：该方案对态制备误差和跳跃率的不确定性具有鲁棒性。当跳跃率不匹配时，通过对位点势能和跳跃率进行联合优化，可以有效减轻由此产生的偏差。
• 初始化鲁棒性：即使初始猜测与真实势能存在显著差异（初始化不确定度 $\lambda \lesssim 1$），优化仍能收敛至真实的势能。虽然较大的偏差可能需要多次初始猜测，但该方法不会失效。

意义与主张
作者声称，该方案通过提供一种可扩展的哈密顿量校准工具，为实现具有任意势能的精密量子模拟铺平了道路。该方法的显著之处在于：

• 实验可行性：它不需要额外的实验修改（如局部控制），并利用标准的 QGM 读出和全局猝灭（global quenches）。
• 经典高效性：通过利用非相互作用动力学，后处理过程在计算上是高效的，避免了与多体模拟相关的“维度灾难”。
• 低样本复杂度：该方案能以可控的实验次数获得高精度，适用于基准测试多体定位（many-body localization）和随机化协议。

论文将这项工作定位为一种稳健的、与模拟兼容的影子断层扫描（shadow-tomography）方案，用于数值学习高斯自由费米子的真实反演映射。虽然目前的实现侧重于非相互作用系统和静态势能，但作者指出，将此方法扩展到相互作用系统和随时间变化的势能是未来的研究方向。他们还提到了借鉴随机测量思想的“影子校准”方法的潜在用途，尽管在目前的研究中并未实现这些扩展。