作者: D. J. Hoare, P. Dunniece, J. Rodewald, K. K. Voges, M. R. Tarbutt

1. 核心物理图象

• 核心物理图象：中性原子量子计算中，原子被随机装载进光镊阵列，必须实时"重排"成无缺陷的整齐阵列才能用于计算。传统方法用声光偏转器逐个搬运原子，速度随原子数线性增长；本文提出一种基于空间光调制器（SLM） 的全并行搬运方案——所有原子同时移动，重排时间只随阵列尺寸的平方根增长。

• 关键突破：作者不依赖迭代傅里叶变换算法，而是解析地算出每一帧全息图（把每个光镊所需的复光场直接叠加、取相位），并让每个光镊的光学相位随位置线性变化。这样既保证了相邻全息图之间的相位连续（避免开关瞬态导致原子丢失），又使光镊位置可以是连续变量而非网格点。

• 一句话总结：用单块GPU在0.284 ms内算出一帧全息图（快于SLM 0.696 ms的刷新周期），从而让重排速度完全由硬件而非计算决定；实验上从256个随机装载位点组装出100原子无缺陷阵列，耗时不到20 ms，单原子存活率99.31%。

---

2. 关键术语解释

1. SEH（Superposition of Elementary Holograms，基元全息图叠加）：一种全息图计算方法。每个光镊对应一个由波矢（决定横向位置）和相位偏移决定的复光场，把N个这样的场直接相加后取相位，即得到SLM相位掩模。相比迭代算法，它能显式控制每个光镊的位置和光学相位，且位置是连续参数，不受傅里叶网格限制。

2. I-SEH（Interpolated SEH，插值SEH）：本文提出的重排协议。在初始和最终全息图之间，对每个光镊的位置、权重、相位做线性插值，得到K个中间全息图（帧）。相位插值走圆周最短路径，把每步相位变化降到∆θ/K，使瞬态强度下降从最坏情况的100%降到K=10时的2.4%。

3. 线性相位插值（LPI）与相位连续性问题：SLM切换全息图时存在瞬态，两帧光场同时显示并干涉。若同一光镊在两帧间光学相位差∆φ=π，瞬态强度会降为零，原子丢失。LPI通过控制相位轨迹使∆φ尽可能小。本文数值验证相位误差标准差仅0.043 rad，比CNN方案（0.2 rad）小一个数量级。

---

3. 主要贡献 (Key Contributions)

1. 解析全息图计算实现硬件极限的并行重排：用单块NVIDIA RTX 5090在0.284 ms内算出每帧全息图，显著快于SLM 0.696 ms刷新周期。这意味着重排速度由SLM刷新率决定，而非计算速度——对最多约2200个光镊的阵列，计算不会成为瓶颈。

2. 相位连续性与强度均匀性的双重精确控制：通过线性相位插值和每步一次权重均衡迭代，将重排过程中的光镊强度非均匀性从约20%降至1.24%（两次迭代降至0.48%），且相位控制几乎不受影响（相位误差标准差从0.0432 rad仅增至0.045 rad）。这使得瞬态强度下降主要由空间重叠而非相位差决定。

3. 实验验证高效2D重排：从随机装载的16×16阵列（填充率55.26%）出发，在18.23 ms内组装出10×10无缺陷阵列，单原子重排存活率99.31(12)%，原子温度仅升高9(3) μK（远低于480 μK阱深），表明原子在搬运后仍保持冷且深束缚。

4. 首次演示3D并行重排与"跳蛙"传输：通过在全息图中加入二次相位项，实现轴向（z方向）搬运。演示了将一组原子轴向抬起、从另一组原子上方"跳"过、再落回焦平面的传输序列，静止原子存活率99.4(3)%，移动原子95.1(4)%。

---

4. 研究方法 (Methodology)

• 全息图生成：将目标光镊阵列的逆傅里叶变换写成解析形式——每个光镊在SLM平面对应一个线性相位梯度，波矢 ( k_n = 2\pi M r_n /(\lambda f) ) 精确编码其横向位置。N个光镊的复场叠加 ( U(r') = \sum_n w_n e^{i(k_n \cdot r' + \theta_n)} )，取相位即得SLM掩模。权重 (w_n) 和相位 (\theta_n) 通过类似加权Gerchberg-Saxton的迭代优化，但每步只需一次 (M \times N) 矩阵乘法（GPU高效）。

• 重排协议（I-SEH）：

1. 预先优化初始和最终阵列的全息图；
2. 荧光成像获取占据图；
3. 用Jonker-Volgenant线性分配算法将占据位点匹配到目标位点，熄灭未匹配光镊；
4. 对位置、权重、相位做K步线性插值，逐帧计算全息图；
5. 每步可选做一次权重均衡迭代（相位固定）。

• 3D扩展：在复场中加入二次相位项 ( h_n|r'|^2 )，其中 ( h_n = \pi M^2 z_n/(\lambda f^2) )，即可控制光镊的轴向位置。传播矩阵相应修改，但仍是矩阵-向量乘法，GPU效率不变。

• 计算-显示流水线：生成线程在GPU上算全息图并推入队列，写入线程从队列取出并写入SLM，两者并行。写入节奏被限制在0.700 ms（SLM刷新0.696 ms + 4 μs余量），防止液晶未完全刷新就覆盖下一帧。

---

5. 实验结果与结论 (Results and Conclusion)

• 2D重排：16×16（256位点）随机装载 → 10×10（100原子）无缺陷阵列，平均需20.78帧，总重排时间18.23 ms（从收到初始图像到最后一帧上传），单原子存活率99.31(12)%。最大步长δ≲1.0 μm时效率稳定在~99%，超过后因空间重叠减少导致瞬态阱深不足而下降。重排引起温升9(3) μK（26.6 → 35.1 μK），远小于480 μK阱深。

• 3D传输：将棋盘格子晶格轴向移动∆z=5、10 μm，原子仍被束缚（只是成像离焦）。轴向步长δz=2.0 μm时往返存活率98.7(2)%，滚降点~3.0 μm与瑞利长度3.69 μm吻合。"跳蛙"传输中，轴向距离d>10 μm时静止原子存活率99.4(3)%、移动原子95.1(4)%；d=5 μm时因光镊过于接近产生串扰和碰撞而效率降低。每步一次强度均衡可使跳蛙存活率提高10-20%。

• 数值表征：85帧、20×20→14×14重排模拟中，相位误差标准差0.0432 rad；一次强度均衡迭代将强度非均匀性从19.85%降至1.24%，两次降至0.48%，而相位误差几乎不变（0.045 rad）。

• 扩展性分析：SEH计算成本随光镊数N线性增长（O(MN)），FFT方案为O(M log M)。N≲2000时SEH可在SLM刷新周期内完成（无均衡迭代）；N≳5000且零填充因子P=4时FFT更快。但重排帧数K∝√N（目标阵列占初始阵列固定比例时），因此总重排时间∝√N，优于顺序搬运的线性增长。

• 开放问题与未来方向：

• 搬运过程中量子比特相干性是否保持？作者预计当强度非均匀性（~1%）小于 (k_B T/U_0)（~5%）时，退相干率与静态阱中相当，但尚待实验验证。
• 如何利用3D运动避免原子碰撞或光镊干涉？可将最小间距约束加入JV分配算法，并通过平滑绕行实现。
• 全息图计算还可加速：消除零权重光镊、融合GPU核函数、用CUDA图捕获每帧、多GPU并行等。
• 该方法适用于连续装载架构、分子阵列组装（分子填充率更低，并行重排收益更大）、以及双物种阵列中辅助量子比特的穿梭（对纠错和读出很重要）。

---

6. 论文标签 (Tags)

中性原子, 里德堡原子, 物理硬件, 量子信息, 编译与优化

📄 点击此处展开/折叠原文 PDF

原文链接： Fast, parallel, three-dimensional rearrangement of atom arrays using superpositions of elementary holograms