Toffoli/CCZ 门｜中性原子量子计算知识库

Toffoli/CCZ 是”逻辑与”门：两个控制都为 $∣1 ⟩$ 时才翻转/移相靶比特。里德堡阻塞让这一逻辑不必拆成 5–6 个两比特门——三维阵列排布最优几何、或时间最优控制直接驱动阻塞链，一次脉冲序列即得原生 CCZ，误差地板由里德堡衰减决定。

问题背景

在标准线路模型中，三比特 Toffoli（C $^{2}$ NOT）门至少需要 6 个 CNOT，而 CCZ 也需要分解成 5 个两比特门 Isenhower 2011 Levine 2019。对量子纠错与编译，这是实打实的开销：表面码稳定子测量、可逆运算、Grover 算法的 oracle 都大量消耗多控门（见三比特门容错稳定子测量与Geyser 原生多比特门编译）。里德堡阻塞提供了一条捷径：一个靶原子的里德堡跃迁可以被同时处于里德堡态的多个控制原子共同移出共振——“全部控制都激发”这一合取条件恰好写进了物理。原生 CCZ 的幺正为

U_{CCZ} = diag (1, 1, 1, 1, 1, 1, 1, - 1), U_{C_{k} Z} = 2 ∣0 ⟩^{\otimes k + 1} ⟨ 0 ∣^{\otimes k + 1} - I,

对靶加局部 Hadamard 即把 CCZ 换成 Toffoli。

核心难题在几何：一维链上只能实现”两控制夹一靶”（控制—靶强阻塞、两控制彼此弱作用），控制数一多，远距控制—靶相互作用按 $R^{- 6}$ 衰减，不对称阻塞条件随之崩塌。于是有两条互补的工程路线：排布最优的三维几何（Harvard 演示一维链；华东师大提出球面优化）与时间最优的脉冲设计（Jandura–Pupillo）。

理论模型：不对称阻塞的几何约束

多控阻塞门要求不对称相互作用：每个控制—靶对的阻塞强度 $U_{c_{j} t}$ 必须远大于控制—控制之间的相互作用 $U_{c_{j} c_{j^{'}}}$ （否则控制原子互相挡路，“合取”条件失效）：

U_{c_{j} t} ≫ U_{c_{j} c_{j^{'}}} ⟺ χ_{j} = \frac{U _{c_{j} t}}{U _{c_{j} c_{j^{'}}}} ≫ 1.

Yu 等人把这一约束变成一个可优化问题：单个靶原子置于球心、 $n$ 个控制原子放在回转椭球面上，用演化算法在球面上搜索控制原子的最优分布，最大化 $min_{j} χ_{j}$ Yu 2022。以 $^{87} Rb$ 、偶极型 $C_{3} /2 π = 4.194 GHz \cdot μ m^{3}$ 与范德瓦尔斯 $C_{6} /2 π = - 12 GHz \cdot μ m^{6}$ 、控制—靶距离 $R_{c t} = 5.0 μ m$ 为例，优化后的 C $_{6}$ NOT 构型（4 个控制位于赤道面、2 个位于两极）给出

(min U_{c_{j} t}, max U_{c_{j} c_{j^{'}}}) /2 π = (33.552, 0.096) MHz, j min χ_{j} = 349.5 ≫ 100,

不对称阻塞误差被压到 $1 0^{- 5}$ 以下。数值模拟给出 C $_{6}$ NOT 门保真度 $\overset{ˉ}{F}_{6} = 0.9920$ ——剩余误差几乎全部来自里德堡态自发衰减，几何与阻塞误差已可忽略；拉伸到 $R_{c t} = 7 μ m$ 时 C $_{12}$ NOT 仍有 $\overset{ˉ}{F} = 0.9537$ 。对照组十分直观：同样的门放在二维蜂窝几何（6+1）上，因不对称条件被破坏，保真度跌至 $0.2821$ 。椭球构型还对原子位置涨落表现出极强的鲁棒性（位置误差可忽略），这是把相互作用各向异性 $C_{3} (θ) \propto (1 - 3 cos^{2} θ)$ 用足的结果。

理论模型：时间最优脉冲

第二条路线问的是：给定阻塞强度 $B$ 、最大 Rabi 频率 $Ω_{m a x}$ ，实现 CCZ 的最短全局脉冲是什么？Jandura 与 Pupillo 用 GRAPE 数值搜索 + Pontryagin 极大值原理（PMP）半解析证明，时间最优解是光滑的相位调制脉冲（振幅固定在 $Ω_{m a x}$ 、只连续地调相位 $φ (t)$ ），并可用少数参数解析刻画 Jandura 2022。

在无穷阻塞 $B = \infty$ 极限下：

CZ：时间最优全局脉冲 $T Ω_{m a x} = 7.61$ ，门误差 $1 - F = 3.1 \times 1 0^{- 10}$ ；它比传统需要单点寻址的脉冲还快约 $10%$ ，且放宽到允许单点寻址也不再提速——单比特分辨能力对两比特门并非必需；
C $^{2}$ Z（等边三角形排布，三对原子阻塞相同）：找到两类定性不同的最优脉冲， $T Ω_{m a x} = 16.43$ 与 $16.53$ ，误差 $3.1 \times 1 0^{- 7}$ 与 $2.8 \times 1 0^{- 6}$ 。

有限阻塞 $B$ 下的结构与代价：

把 $B$ 当作已知参数，可以对特定 $B$ 做出精确门： $B = 10 Ω_{m a x}$ 时 CZ 与 C $^{2}$ Z 都有 $1 - F < 1 0^{- 9}$ 的脉冲（CZ $T Ω_{m a x} = 7.58$ ；C $^{2}$ Z $T Ω_{m a x} = 16.2 - 17.1$ ）；
反过来，有限 $B$ 也可用来加速：允许双里德堡激发态（反阻塞区）参与时， $B = 10 Ω_{m a x}$ 的 C $^{2}$ Z 可缩短到 $T Ω_{m a x} = 10$ ，但对 $B$ 极其敏感—— $B$ 减小 $10%$ 误差即恶化到 $1 - F = 0.59$ ；
若要求脉冲只依赖 $B = \infty$ 的解并对 $B$ 涨落鲁棒，则一阶误差可写成

(1 - F) B^{2} ≃ \frac{1}{4} ⟨ 11∣ ψ_{11}^{(1)} (T)⟩^{2},

即误差按 $B^{- 2}$ 标度；以 $min (1 - F) B^{2} T^{2}$ 为目标优化出的脉冲把有限 $B$ 的效应完全归约为可补偿的 AC Stark 移。当 $B ≫ Ω_{m a x}$ （实验已演示 $B /2 π = 180 MHz$ ，理论上限 $B /2 π \sim 3 GHz$ ）时该方案给出与容错阈值相容的保真度。若三原子排成直线而非三角形，外—外距离加倍使 $B^{'} = B /64$ （ $R^{- 6}$ ），需要重新设计脉冲。

时间最优 CCZ 脉冲的激光相位 φ(t) 与三角/直线两种阻塞几何（图源：Jandura 2022 Fig. 2）

方法与实现

一维链上的原生 CCZ/Toffoli 演示 Levine 2019：三个 $^{87}$ Rb 原子排成最近邻阻塞链，两个外层控制原子各以 $2 π \times 24 MHz$ 阻塞中间靶原子，外层彼此仅 $2 π \times 0.4 MHz$ 。用 RedCRAB 远程优化算法数值优化幅频调制的全局激光脉冲直接实现 CCZ，嵌入回波序列消除光移；转成 Toffoli 后对全部 8 个计算基态测量真值表，保真度 $F_{Toff} \geq 83.7 (3) %$ 、SPAM 修正后 $F_{Toff}^{c} \geq 87.0 (4) %$ ，受限层析给出 $F_{LT}^{c} \geq 86.2 (6) %$ ；四个三比特门在四个原子簇上并行执行（见并行纠缠门）。

实现质量的上游约束是基元里德堡激发的相干性。de Léséleuc 等对单原子光镊中 $∣ g ⟩ \leftrightarrow ∣ r ⟩$ Rabi 振荡的系统分析表明 de Léséleuc 2018：对比度损失来自制备/读出误差，而振荡衰减是多普勒退相（ $30 μ K$ 下双光子失谐起伏 $\sim 2 π \times 120 kHz$ ）、中间态自发散射与激光相位噪声各自百分之几贡献的叠加—— $2 π \times 2 MHz$ Rabi 频率下 $1/ e$ 衰减时间仅约 $5 μ s$ ，使单原子 $π$ 脉冲保真度约 $95%$ ；里德堡寿命则按 $n^{3}$ 标度、所用态超过 $100 μ s$ 。这份”误差清单”直接决定了多比特门能达到的保真度上限，也指明改进方向：更低温、更大中间态失谐与更窄线宽激光（与误差预算与真实建模的框架一致）。

几何升级：三维可重构阵列（见原子装配无缺陷阵列）使”球心靶 + 球面控制”的椭球构型成为现实选项，把多控门从一维链的 3 比特扩展到 6–12 比特规模 Yu 2022。

参数与量级

量	典型值	来源
电路分解代价	Toffoli $\geq 6$ CNOT；CCZ 分解需 5 个两比特门	Isenhower 2011、Levine 2019
一维链阻塞比	控制—靶 $2 π \times 24 MHz$ vs 控制—控制 $2 π \times 0.4 MHz$	Levine 2019
演示 Toffoli 保真度	$F_{Toff}^{c} \geq 87.0 (4) %$ （四个三比特门并行）	Levine 2019
球面最优几何不对称度	$min χ_{j} = 349.5$ （ $U_{c_{j} t} /2 π = 33.552 MHz$ ， $U_{c_{j} c_{j^{'}}} /2 π = 0.096 MHz$ ）	Yu 2022
原生 C $_{6}$ NOT / C $_{12}$ NOT	$\overset{ˉ}{F}_{6} = 0.9920$ （衰减极限）； $\overset{ˉ}{F}_{12} = 0.9537$ （ $R_{c t} = 7 μ m$ ）	Yu 2022
时间最优 CZ / C $^{2}$ Z	$T Ω_{m a x} = 7.61$ （ $1 - F = 3 \times 1 0^{- 10}$ ）/ $16.43$ 、 $16.53$	Jandura 2022
有限阻塞下的精确门	$B = 10 Ω_{m a x}$ ： $1 - F < 1 0^{- 9}$ ；鲁棒脉冲误差 $\propto B^{- 2}$	Jandura 2022
单原子激发衰减	$1/ e \sim 5 μ s$ （ $Ω/2 π = 2 MHz$ ）； $π$ 脉冲约 $95%$	de Léséleuc 2018

与其他概念的关系

单控制版本的一控多靶门见多比特阻塞门与猫态；把合取条件推广为”按激发数奇偶相位”的是奇偶相位门。
一维链 CCZ 与阻塞 CZ 门共享双全局脉冲框架，其并行化讨论见并行纠缠门。
时间最优方案是量子速度极限思想在里德堡门上的落地，脉冲的光滑性与鲁棒性分析工具属于鲁棒最优控制与GRAPE 梯度上升脉冲工程。
最优几何依赖阵列的可编程排布能力，见原子装配无缺陷阵列与FPQA 计算模型；原生多比特门如何被编译器利用见Geyser 原生多比特门编译。
表面码与稳定子测量对原生 Toffoli 的需求见三比特门容错稳定子测量；衰减与位置误差的建模框架见误差预算与真实建模。
相位分配（CCZ 的 $- 1$ 放在哪个分支、局部 $Z$ 如何吸收多余相位）遵循门相位与局域等价。

参考文献

D. Yu, H. Wang, J.-M. Liu, S.-L. Su, J. Qian, W. Zhang, “Multiqubit Toffoli Gates and Optimal Geometry with Rydberg Atoms”, Phys. Rev. Applied 18, 034072 (2022). NA075 · DOI 10.1103/PhysRevApplied.18.034072 · arXiv:2203.14302
S. Jandura, G. Pupillo, “Time-Optimal Two- and Three-Qubit Gates for Rydberg Atoms”, Quantum 6, 712 (2022). NA077 · DOI 10.22331/q-2022-05-13-712 · arXiv:2202.00903
S. de Léséleuc, D. Barredo, V. Lienhard, A. Browaeys, T. Lahaye, “Analysis of imperfections in the coherent optical excitation of single atoms to Rydberg states”, Phys. Rev. A 97, 053803 (2018). r3_leseleuc_2018 · DOI 10.1103/PhysRevA.97.053803
H. Levine et al., “Parallel Implementation of High-Fidelity Multiqubit Gates with Neutral Atoms”, Phys. Rev. Lett. 123, 170503 (2019). NA063 · DOI 10.1103/PhysRevLett.123.170503
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

中性原子百科

探索

Toffoli/CCZ 门