多比特阻塞门与猫态

阻塞半径内”只能有一个里德堡激发”的硬约束，把 $N$ 个原子变成一个共享单激发的”超原子”：一个受激控制原子能同时挡住半径内所有靶原子的激发，于是 CkNOT、CkZ 这类一控多靶门与 GHZ/猫态可以在一次全局脉冲或极少数脉冲内完成——门数不再随比特数指数增长。

物理图像

里德堡阻塞给出的最紧凑资源是 Dicke 流形：$N$ 个处于阻塞半径内的原子，可及对称态空间压缩为 $\{|gg\cdots g⟩,|W_N⟩=\frac{1}{\sqrt{N}}{\sum }_i|g\cdots r_i\cdots g⟩\}$。对称驱动下基态 $|G⟩$ 与 $|W_N⟩$ 的耦合被集体增强为

$${\Omega }_N=\sqrt{N}\Omega ,$$

即”超原子”以 $\sqrt{N}$ 倍速率吸收光子。这枚硬币有两面：集体增强让介观系综的门可以跑得很快，但也意味着逐靶操作的门必须在每个 $k$ 下重新优化脉冲（见EIT 里德堡门对这一问题的规避）；而”单激发即可全局阻塞”则开辟了另一条路——用激发的开关结构本身承载多比特逻辑 Isenhower 2011。

把二比特 CZ 的”阻塞即相位”逻辑直接铺开，就得到一族原生多比特门：一个控制原子的里德堡激发可以同时约束 $k$ 个靶原子（CkNOT/CkZ），两个控制原子可以同时约束一个靶原子（Toffoli/CCZ，见Toffoli/CCZ 门），甚至可以对整团原子按激发数奇偶性施加多体相位（见奇偶相位门）。与按奇偶性”读取”激发数不同，本条聚焦选择性阻塞：只有”全部控制都处于 $|1⟩$“（或”控制处于 $|1⟩$“）的分支才被允许/禁止翻转。

理论模型：一控多靶的门结构与误差标度

集体增强的来源

记 $|G⟩=|g{⟩}^{\otimes N}$，单激发对称态为 $|W_N⟩=\frac{1}{\sqrt{N}}{\sum }_{i=1}^N|g\cdots r_i\cdots g⟩$。全局驱动 $\Omega {\sum }_i(|r_i⟩⟨g_i|+h.c.)/2$ 只把 $|G⟩$ 与 $|W_N⟩$ 耦合起来（其余激发态非对称、不耦合），矩阵元为

$$⟨W_N|\frac{\Omega }{2}\sum _{i=1}^N{\sigma }_i|G⟩=\frac{\Omega }{2}N\cdot \frac{1}{\sqrt{N}}=\frac{\sqrt{N}\Omega }{2},$$

即有效哈密顿量 $H_{\mathrm{eff}}=\sqrt{N}\Omega (|W_N⟩⟨G|+h.c.)/2$——这就是”超原子”以 ${\Omega }_N=\sqrt{N}\Omega$ 吸收光子的原因。多比特门设计中每一个 $\sqrt{n}$ 因子都必须被脉形或序列吸收（ARP 门用双脉冲处理，EIT 里德堡门用暗态绕开）。

CkNOT 的两种序列

Isenhower、Saffman 与 Mølmer 给出了 CkNOT 门的两种实现 Isenhower 2011：

• 顺序寻址：对 $k$ 个控制原子依次施加 $|0⟩\to |s⟩$ 的里德堡 $\pi$ 脉冲（共 $k$ 个），再对靶原子执行受控交换脉冲序列，最后按逆序把控制原子送回基态（$k$ 个 $\pi$ 脉冲）——任何控制处于 $|0⟩$ 都会被先期激发的其他控制原子阻塞。共需 $2k+3$ 个里德堡 $\pi$ 脉冲；
• 同时寻址：$k$ 个控制原子被同一个脉冲同时激发到 $|s⟩$，靶原子操作后同时退激发，全程只需 5 个里德堡 $\pi$ 脉冲，与 $k$ 无关。

作为对照，电路分解极其昂贵：最优已知线路实现 CkNOT 在 $k\ge 5$ 时需约 $32k$ 个单/双比特门，且 CkNOT 至少需要 $2k+2$ 个 CNOT（三比特 Toffoli 恰需 6 个 CNOT）Isenhower 2011。

误差有两类并给出清晰的标度律 Isenhower 2011：

1. 里德堡态自发辐射：$k$ 个控制原子在门期间”值班”的里德堡占据时间随时钟变长，误差近似按 $k^2$ 增长（$E_{se}\propto k^2\pi {\Omega }_c/4B^2\tau$ 量级）；
2. 有限阻塞的旋转误差：阻塞能移 $B$ 有限使本应禁戒的激发以 $\sim (\Omega /B{)}^2$ 概率发生，多控制原子同时”值班”时该项随 $k$ 代数增长。

以 Cs 参数、偶极—偶极控制—靶相互作用 $B_{ct}\sim n^4/R^3$ 计算的数值优化表明：$k=35$ 时门误差仍可低于 $10\%$；顺序方案本征误差约为同时方案的三分之一，但同时方案脉冲数少约 15 倍、门时间 $T=3\pi /{\Omega }_t+2\pi /{\Omega }_c=1.1\mu s$（顺序方案约 $9\mu s$，$k=35$），对技术噪声更不敏感 Isenhower 2011。

多体纠缠制备是同一机制的副产品：对处于阻塞体积内的 $N$ 个原子，先用 Hadamard 把其中一个原子置于叠加态，再用阻塞 CNOT 门做扇出（每个 CNOT 只是嵌入本地旋转的原生 CZ），即可制备

$$|GHZ{⟩}_N=\frac{1}{\sqrt{2}}(|00\cdots 0{⟩}_N+e^{i\varphi }|11\cdots 1{⟩}_N),$$

即薛定谔猫态；由于阻塞约束，扇出分支上的集体激发以 $\sqrt{N}\Omega$ 增强，制备时间可以远短于里德堡寿命。按激发数奇偶性直接施加多体相位的原生途径见奇偶相位门。

方法与实现

并行多比特门的首次高保真演示（哈佛，${}^{87}$Rb） Levine 2019：比特编码于 $|0⟩=5S_{1/2},F=1$、$|1⟩=5S_{1/2},F=2$，里德堡态 $70S_{1/2}$。CZ 用两个全局脉冲实现（总时长 $0.4\mu s$，期间关闭光镊以避免里德堡态反陷阱；$\Omega \approx 2\pi \times 3.5\mathrm{MHz}$，最近邻阻塞 $2\pi \times 24\mathrm{MHz}$）：两个脉冲之间的相对相位 $\varphi$ 使 $|01⟩,|10⟩$ 各累积 $\varphi$、$|11⟩$ 累积 $2\varphi -\pi$，与 CZ 的典则形式只差可被单比特 $Z$ 吸收的单粒子相位（见门相位与局域等价）。具体序列把 CZ 嵌入回波结构以消除里德堡激光光移。五个原子对上并行制备 Bell 态：$F\ge 95.0(2)\%$（原始下界；未修正 $95.9(2)\%$），SPAM 修正后纠缠门保真度 $F_c\ge 97.4(3)\%$，CNOT 真值表 $F_c\ge 96.5(3)\%$。同一系统把三个原子排成最近邻阻塞链（边缘原子对相互作用仅 $2\pi \times 0.4\mathrm{MHz}$），用 RedCRAB 优化的幅频调制全局脉冲直接实现 CCZ，转成 Toffoli 后真值表保真度 $F_{\mathrm{Toff}}^{\mathrm{c}}\ge 87.0(4)\%$（受限层析 $F_{\mathrm{LT}}^{\mathrm{c}}\ge 86.2(6)\%$），且四个三比特门并行执行；主要残余误差来自有限原子温度与里德堡动力学期间的非共振激光散射。

多比特纠缠与算法运行（ColdQuanta/Saffman 组，${}^{133}$Cs） Graham 2022：钟态比特（微波全局 $R_{\varphi }(\theta )$ + 459 nm 光移局域 $R_Z(\theta )$）经 459/1040 nm 双光子激发到 $75S_{1/2}$，CZ 采用失谐双脉冲协议，单原子 Rabi $1.7\mathrm{MHz}$、阻塞能移 $3\mathrm{MHz}$（弱阻塞区），另加两个局域 $R_Z$ 脉冲把受控相位整定成典则 CZ。优化对上 Bell 态保真度 $92.7(1.3)\%$（SPAM 修正约 $95.5\%$）。在此基础上制备了 $N=6$ 的 GHZ/猫态（Hadamard + 五条 CNOT 边 1-2、1-3、1-4、1-5、5-6），所有 $N\le 6$ 的 GHZ 态经 SPAM 修正后保真度超过 $50\%$，相干时间超过 $600\mu s$——得益于钟态编码的长 $T_2$；GHZ 电路（含相位估计与 QAOA，编译到原生门集）运行时长 $40-180\mu s$。这是首次在门模型平台上测量 $N>2$ 的 GHZ 态保真度。

规模上限：理论分析指出在合理的里德堡参数下（$n\sim 75$ 顺序寻址、$n\sim 60$ 同时寻址，晶格周期约 $4\mu \mathrm{m}$），CkNOT 门可推广到数十个原子的控制规模而误差保持在百分之几量级；液氮温区运行可再降约 $30\%$ 的里德堡辐射误差 Isenhower 2011。

GHZ 态 2–6 比特的 parity 振荡与保真度随 N 的衰减（图源：Graham 2022 Fig. 2）

参数与量级

量	典型值	来源
集体 Rabi 增强	${\Omega }_N=\sqrt{N}\Omega$	Pelegrí 2022
CkNOT 脉冲数	顺序 $2k+3$；同时寻址 5 个（与 $k$ 无关）	Isenhower 2011
电路分解代价	$k\ge 5$ 时约 $32k$ 个基本门；下界 $2k+2$ 个 CNOT	Isenhower 2011
CkNOT 规模与误差	$k=35$ 时误差 $<10\%$（Cs 参数）	Isenhower 2011
同时寻址门时间	$1.1\mu s$（$k=35$；顺序约 $9\mu s$）	Isenhower 2011
并行双比特门（${}^{87}$Rb）	双全局脉冲 $0.4\mu s$，$F_c\ge 97.4(3)\%$，五对并行	Levine 2019
原生 Toffoli（RedCRAB）	$F_{\mathrm{Toff}}^{\mathrm{c}}\ge 87.0(4)\%$，四个三比特门并行	Levine 2019
GHZ 猫态	$N=6$，修正后 $F>50\%$，相干 $>600\mu s$	Graham 2022

与其他概念的关系

• 集体 $\sqrt{N}$ 增强与 Dicke 流形是里德堡阻塞的直接推论；介观集体编码的”透明开关”版本见EIT 里德堡门。
• “两个控制约束一个靶”的三比特版本即Toffoli/CCZ 门；对激发数整体施加相位的推广是奇偶相位门。
• 多簇门同时执行而不互相串扰，是并行纠缠门的出发点与实验基础。
• 阻塞分支上的动力学可以用与阻塞 CZ 门、ARP 门相同的整形/最优化工具改善（解析整形脉冲、鲁棒最优控制）。
• 原生 CkZ/Toffoli 门改变编译层的门计数与调度，相关消费方见Geyser 原生多比特门编译与三比特门容错稳定子测量；免测量纠错中的多比特门需求见免测量纠错。
• GHZ 态是qLDPC 码架构与磁测量/相位估计等应用的原生资源。

参考文献

• L. Isenhower, M. Saffman, K. Mølmer, “Multibit C${}_k$NOT quantum gates via Rydberg blockade”, Quantum Inf. Process. 10, 755 (2011). NA053 · DOI 10.1007/s11128-011-0292-4 · arXiv:1104.3916
• H. Levine, A. Keesling, G. Semeghini, A. Omran, T. T. Wang, S. Ebadi, H. Bernien, M. Greiner, V. Vuletić, H. Pichler, M. D. Lukin, “Parallel Implementation of High-Fidelity Multiqubit Gates with Neutral Atoms”, Phys. Rev. Lett. 123, 170503 (2019). NA063 · DOI 10.1103/PhysRevLett.123.170503 · arXiv:1908.06101
• T. M. Graham, Y. Song, J. Scott, C. Poole, L. Phuttitarn, K. Jooya, P. Eichler, X. Jiang, T. Marra, B. Grinkemeyer, M. Kwon, M. Ebert, J. Cherek, M. T. Lichtman, M. Gillette, J. Gilbert, D. Bowman, T. Ballance, C. Campbell, E. D. Dahl, O. Crawford, N. S. Blunt, B. Rogers, T. Noel, M. Saffman, “Multi-qubit entanglement and algorithms on a neutral-atom quantum computer”, Nature (2022). NA074 · DOI 10.1038/s41586-022-04603-6