中性原子基准测试｜中性原子量子计算知识库

中性原子基准测试的核心问题是：当一台里德堡阵列处理器拥有"全连通里德堡相互作用 + 原子搬运重构连接 + 原生多比特门"这些异于超导/离子阱的能力时，用什么指标、什么噪声模型、什么电路族才能公平地回答"这台机器跑真实算法到底行不行"。

问题背景

逐门指标（单比特/双比特门保真度、 $T_{1} / T_{2}$ 、读出保真度）刻画的是”零部件”，而用户关心的是”整车”：给定一个算法电路，处理器以多大成功率输出正确结果。随机基准（randomized benchmarking）只预测随机电路的表现，对真实应用电路的保真度预测很差 Wagner 2024。量子体积（quantum volume, QV）虽被广泛采用，但只测”宽度 = 深度”的方型随机 SU(4) 电路，无法体现中性原子平台最独特的两项资产：通过原子搬运获得的动态全连通（见FPQA 计算模型），以及以多比特阻塞门、Toffoli/CCZ 门为代表的原生多比特门。

因此，中性原子的平台级基准需要同时回答三个方法学问题：(1) 选哪套应用级基准电路族；(2) 用什么噪声模型把逐门参数映射为算法成功率；(3) 如何把”搬运重构连接”与”固定连接”的差异变成可量化的对比。目前语料中的两条代表性路线是：基于 QED-C 基准套件的全连接–近邻对比仿真 Wagner 2024，与基于量子体积 + Bernstein–Vazirani + Grover 的近期处理器算法性能评估 McInroy 2024。

理论模型

输出分布保真度。 应用级基准通常比较实验（或仿真）输出概率分布 $P_{out}$ 与理想分布 $P_{ideal}$ 。QED-C 套件采用的经典保真度定义为逐位串概率的乘积归一化形式，可写为

F_{n} (P_{ideal}, P_{out}) = (x \sum P_{ideal} (x) P_{out} (x))^{2}, F = max (F_{n}, 0),

即截断到非负的态重叠（经典限制下的纯态保真度）Wagner 2024。SupermarQ 等套件则用 Hellinger 保真度，见通用量子基准套件。

门保真度。 噪声模型中的逐门误差用平均门保真度刻画，对理想幺正 $U$ 与噪声信道 $E$ ，

F_{G} (ρ, σ) = [Tr ρ σ ρ]^{2},

其中 $ρ, σ$ 分别为理想与仿真末态密度矩阵，对 Haar 随机输入态取平均 Wagner 2024。

量子体积。 QV 定义为”能可靠实现的最大方型随机 SU(4) 电路”的规模：对宽度 $N_{Q}$ 、深度 $d = N_{Q}$ 的随机电路，统计”重输出”（heavy output，理想概率高于中位数的位串）概率 $\overset{ˉ}{h}_{p}$ ，通过判据 $\overset{ˉ}{h}_{p} > \frac{1 + ϵ}{2}$ （单边置信区间）后取

lo g_{2} V_{Q} = N_{Q} ar g max min (N_{Q}, d) .

McInroy 2024 在虚拟设备上对每个深度用 200 条随机电路评估该判据。QV 的形式化局限（只测方型电路、不含测量调度）正是 QED-C 基准把电路分为”浅层 oracle 算法 / 量子子程序 / 真实问题”三族、同时扫描宽度与深度的动机 Wagner 2024。

体积化基准（volumetric benchmarking）。 比单一 QV 分数更进一步的做法是在”宽度 × 深度”平面上逐格报告保真度，把设备的可用区域画成一张等高线图：中性原子平台因全连通而在”宽而浅”区域占优，这是方型电路指标看不到的 McInroy 2024。随机基准（如镜像电路 mirror circuits）虽也扫描更多电路，但使用随机生成的”不实用”电路，已被证明不能很好预测实际电路的保真度 Wagner 2024。

方法与实现

实验校准的噪声模型（NA181）。 Wagner–Saffman 等以 Nature 604, 457 (2022) 的实验数据为基准构建小规模里德堡处理器仿真：CZ 门（阻塞 CZ 门）是最易出错的操作，平均保真度 $F ≃ 0.95$ ，建模为相位翻转 + $∣1 ⟩$ 衰减 + 向亮/暗”丢失态”的损耗信道；微波旋转 $F_{R_{ϕ}} ≃ 0.9995$ ，局域 $R_{z}$ 相位门 $F_{R_{Z}} ≃ 0.995$ ； $T_{1}$ 约为门时间的 $1 0^{6}$ 倍（可忽略）， $T_{2}$ 约为门时间的 $\sim 100$ 倍（对较长电路显著）；制备与测量误差（SPAM）各约 $0.5%$ Wagner 2024。为追踪泄漏与原子丢失（见原子损失），仿真把每个原子当作四级（ququart）系统：泄漏态按读出时的亮/暗归入 $∣ l ⟩_{0}, ∣ l ⟩_{1}$ ，不参与后续操作。

模型校准的判据是与真实实验对照：3 条 4 比特 QAOA-MaxCut 电路的布居谱平均经典保真度 $98.6%$ （最差 $97.5%$ ），Bell 态保真度 $91.3%$ 与实验一致；GHZ 电路的仿真保真度平均高估 $8.0%$ ，归因于噪声参数随时空漂移而被模型假设为常数 Wagner 2024。

全连通 vs 近邻对比。 在同一噪声模型下分别模拟”全对全”（all-to-all，任意两原子可直接 CZ）与”最坏情形”（二维方格最近邻）两种连接，电路宽度至 11 比特。结论：全连通使仿真电路保真度整体提高 $10% - 15%$ ；在 7–11 比特区段，Deutsch–Jozsa 平均提升 $2.6%$ （最大 $4.3%$ ），Bernstein–Vazirani 平均提升 $3.1%$ （最大 $3.7%$ ）；Hidden Shift 电路因每次只纠缠 2 比特而无差别 Wagner 2024。

量子子程序层面的对比更鲜明：仅近邻连接时，量子傅里叶变换（QFT）基准的方法 2（制备→修改→解码）最高在 2 比特达到 $85%$ 、3 比特跌至 $67%$ ，方法 1（两次 QFT）从 $77%$ 跌至 $48%$ ；3 比特相位估计仿真保真度 $89%$ ；而所有振幅估计电路的保真度都不足 $20%$ 。这些数字给出了”哪些子程序已可用、哪些还差一个量级”的精细坐标 Wagner 2024。

近期处理器的算法性能仿真（NA182）。 McInroy 等构建了一个”存储区 + 门区”分区虚拟设备（见分区架构），参数取 $T_{1} = 4 s$ 、 $T_{2} = 1.49 s$ 、初始化误差 $0.003$ ，双比特门用 ARP 方案 CZARP（理论相位矩阵中 $∣11 ⟩$ 分量幅度 $0.9986$ ）与 LP 方案两种门、连接方式用”静态 SWAP 重排”与”相干搬运重构”两种模态；在 QuEST 库上做全密度矩阵演化，把丢失建模为非幺正过程，再用非破坏读出（见中电路测量）实现”丢失后选择”（loss correction）McInroy 2024。三项测试为：量子体积、Bernstein–Vazirani（利用单原子搬运逐对耦合 ancilla）与 Grover 搜索（利用原生 CCZ 门一次完成三比特 oracle）。

建模细节上有两点值得注意。其一，随机 SU(4) 层被分解为 3 个 CZ 与至多 15 个单比特旋转，静态模态在层间需沿行插入 SWAP 链完成配对重排（层深与门数都增加），而重构模态只需一层搬运即可整体换对——门数差异直接决定了静态模态的无后选择量子体积更低（ARP 门下 $V_{Q} = 2^{7}$ ）。其二，搬运时长被建模为固定时长 $τ$ （多原子协同搬运用 $τ_{moveA}$ 、单原子搬运用 $τ_{moveB}$ ），使误差估计不依赖于具体路由算法的优劣 McInroy 2024。

横向坐标。 仿真的 $V_{Q} = 2^{9}$ （9 比特设备的理论上限，也是更大系统的下界）已与 IBM Prague 超导处理器报告的 $2^{9}$ 相当，而全部平台中的最好成绩来自离子阱的 $2^{16}$ ——这为”里德堡处理器是有竞争力的近期平台”提供了可比较的坐标 McInroy 2024。

方法学局限

仿真式基准的结论强度取决于噪声模型的诚实程度。NA181 的模型自检给出了三处已知偏差：实验中噪声参数随比特位置与时间漂移而模型取常数；个别低误差信道被简化以降低模拟复杂度（如把微波旋转简化为去极化信道，并不能完全刻画该门的真实行为）；由此 GHZ 电路保真度平均被高估 $8.0%$ Wagner 2024。NA182 则在两处明确声明其结果是下界： $V_{Q} = 2^{9}$ 是 9 比特设备的上限，但更大 $N_{Q}$ 的判据判定超出其计算范围；若不做丢失后选择，数字会更低 McInroy 2024。此外，两篇工作都未建模串扰与关联错误——这些恰是通用量子基准套件中”并行度”特征要探测的对象。

参数与量级

量	典型值	来源
CZ 门平均保真度（噪声模型）	$F ≃ 0.95$	Wagner 2024
微波 / $R_{z}$ 门保真度	$0.9995$ / $0.995$	Wagner 2024
SPAM 误差	制备与测量各 $\sim 0.5%$	Wagner 2024
全连通对近邻的保真度增益	$10% - 15%$ （DJ $+ 2.6%$ 、BV $+ 3.1%$ 均值，7–11 比特）	Wagner 2024
量子体积（9 比特虚拟设备）	$V_{Q} = 2^{9}$ （该规模上限；静态 SWAP+ARP 仅 $2^{7}$ ）	McInroy 2024
Bernstein–Vazirani 成功率	$> 0.95$ （9 数据 + 1 ancilla，未做丢失校正原始值）	McInroy 2024
Grover 搜索成功率	$0.97$ （丢失校正后，6 数据 + 3 ancilla，原生 CCZ）	McInroy 2024
设备相干参数（仿真）	$T_{1} = 4 s$ ， $T_{2} = 1.49 s$ ，初始化误差 $0.003$	McInroy 2024

全对全 vs 最近邻连接下 Deutsch-Jozsa/Bernstein-Vazirani/Hidden Shift 基准保真度热图（图源：Wagner 2024 Fig. 2）

量子体积仿真的重输出概率与丢失概率随电路尺寸的变化，含 ARP/LP 与静态/重构两种模态（图源：McInroy 2024 Fig. 3）

与其他概念的关系

基准电路族与特征向量指标的通用设计（QASMBench、SupermarQ、MQT Bench）见通用量子基准套件；把基准用于对比编译器而非硬件的方法学见编译器公平评测与可视化。
噪声模型中的 CZ 门即阻塞 CZ 门；Grover oracle 使用的三比特门见Toffoli/CCZ 门与多比特阻塞门，并行执行方式见并行纠缠门。
“动态全连通”这一被基准量化的能力来自原子搬运与重构，其形式化模型见FPQA 计算模型，搬运本身见原子搬运与相干输运。
丢失后选择依赖非破坏读出，见中电路测量与非破坏读出；丢失信道建模与原子损失与容损失架构直接相关。
把逐门误差预算分解到物理机制的建模方法见误差预算与真实建模；以真实硬件跑到逻辑层面的基准结果见逻辑量子处理器。

参考文献

N. Wagner, C. Poole, T. M. Graham, M. Saffman, “Benchmarking a neutral-atom quantum computer”, Int. J. Quantum Inf. (2024). NA181 · arXiv:2311.04141 · DOI 10.1142/S0219749924500011
K. McInroy, N. Pearson, J. D. Pritchard, “Benchmarking the algorithmic performance of near-term neutral atom processors”, 预印本 (2024). NA182 · arXiv:2402.02127

中性原子百科

探索

中性原子基准测试