奇偶相位门与原生多体门

奇偶相位门把”基矢串的 Z 乘积”本身当作门目标：一个全局里德堡脉冲按比特串的奇偶性施加 $\mp \theta$ 相位，把三体、四体乃至 $N$ 体纠缠相位从”分解成一串两比特门”降为”单个脉冲”——这正是奇偶架构量子优化与快速稳定子读出所需要的原生资源。

物理图像

里德堡阻塞给出的硬约束是”阻塞半径内只能有一个里德堡激发”。传统的阻塞 CZ 门把这个约束当作误差来源小心规避，而奇偶相位门（parity phase gate）换了一个思路：直接让脉冲在”偶宇称”与”奇宇称”两个子空间之间累积相对相位。由于里德堡哈密顿量只耦合 $|1⟩$ 与 $|r⟩$、对 $|0⟩$ 无作用，演化算符在 $2^N$ 个计算基矢上是块对角的，每个基矢各拾取一个相位——于是”一个全局脉冲 = 一组对角相位”成为中性原子的天然能力 Kazemi 2025 Locher 2026。

这类似于多比特阻塞门与Toffoli/CCZ 门所利用的物理，但目标不同：多比特阻塞门天然实现”对称的”多控相位，而奇偶门家族通过几何排布与脉冲整形，可以选择性地实现纯 $N$ 体相互作用、单控多靶（fan-out）与多控（$C_{n}Z$ 型）等非对称门型，直接对接奇偶架构（parity architecture）优化问题与量子纠错中的 ${\hat{Z}}^{\otimes N}$ 型稳定子测量 Dlaska 2022 Stein 2026。

理论模型

奇偶相位门的定义

$N$ 比特奇偶相位门定义为 Kazemi 2025

$$Z_N(\theta )=\exp \left(-i\theta {\sigma }_z^{(1)}\otimes {\sigma }_z^{(2)}\otimes \cdots \otimes {\sigma }_z^{(N)}\right),$$

它对含偶数个 $|1⟩$ 的基矢施加 $e^{-i\theta }$、奇数个施加 $e^{+i\theta }$。当 $\theta =n\pi /2$ 时它退化为局域 $Z$ 旋转的乘积（不产生纠缠），而最大纠缠能力出现在 $\theta =(2n+1)\pi /4$；$\theta =\pi /4$ 是最常用也最具代表性的工作点。注意 $Z_N(\theta )$ 的生成元是纯 $N$ 体相互作用，区别于多控相位门

$$C_{N-1}P(\theta )=\exp (-i\theta \underset{j=1}{\overset{N}{⨂}}\frac{\mathbb{I}-{\sigma }_z^{(j)}}{2}),$$

后者的生成元按二项式展开混合了 1 体到 $N$ 体的所有项 Kazemi 2025。

CZZ/CZZZ、CnZ 与 CCZ 的关系

单控多靶门 $C(Z^{\otimes N})$ 在控制比特为 $|1⟩$ 时对 $N$ 个靶比特施加 $Z^{\otimes N}$ 相位，其生成元 $|1⟩⟨1{|}_c\otimes Z^{\otimes N}=\frac{1}{2}(\mathbb{I}-Z_c)\otimes Z^{\otimes N}$ 同样混合了从局域到 $N+1$ 体的成分；当 $N=2,3$ 时即为 CZZ、CZZZ 门 Stein 2026。多控相位门则常记作（$n$ 控 1 靶）

$$C_{n}Z(\alpha )=\mathrm{diag}(\underset{2^n}{\underbrace{1,\dots ,1}},e^{i\alpha }),\dim C_{n}Z=2^{n+1},$$

这是全局里德堡脉冲在三个及以上原子上能实现的所有门参数化的基本构件 Locher 2026。三者的桥梁是相位 Gadget/Z2 傅里叶展开：投影算符 $|1\cdots 1⟩⟨1\cdots 1|=2^{-N}{\sum }_{S\subseteq [N]}{\prod }_{i\in S}{Z}_i$，因此

$$CCZ=e^{i\pi |111⟩⟨111|}=\exp \left[\frac{i\pi }{8}\right(\mathbb{I}+\sum _i{Z}_i+\sum _{i<j}{Z}_i{Z}_j+Z_1{Z}_2{Z}_3\left)\right],$$

即 CCZ 在相差局域 $Z$ 旋转与全局相位意义下等于三个 $ZZ$ 奇偶门与一个 $ZZZ$ 奇偶门（各取相位 $\pi /8$）的乘积。奇偶门是对角幺正相位多项式中的单项，多控门是它们的乘积；见对角线路综合与多控相位门综合。

Dlaska 等面向 LHZ/奇偶架构提出的四体门写作 $U(\gamma )$：对 4 比特 plaquette 的奇宇称态施加 $e^{+i\gamma }$、偶宇称态施加 $e^{-i\gamma }$，通过两次绝热脉冲的动力学相位满足 ${\varphi }_1+{\varphi }_3=\gamma$、${\varphi }_0+{\varphi }_4=2{\varphi }_2=-\gamma$ Dlaska 2022。类似地，基于星形图与辅助原子链的协议把输出态的宇称 ${\nu }_q$ 变成几何相位 ${\phi }_g={\nu }_q\pi$，可转换为 $C^{k}Z$ 门（见远距原子门）Delakouras 2026。

方法与实现

目前有四条互补的原生实现路线：

1. 绝热 plaquette 脉冲（四体门）。四个原子排成四面体使所有原子对相互作用 $V$ 相同；第一段绝热脉冲只耦合 $|\downarrow ⟩$ 到里德堡态、第二段只耦合 $|\uparrow ⟩$，利用 $n$ 体激发能谱的依赖关系产生宇称相关相位。采用”两次暂停”（two-pause）协议：对任意 $\gamma$ 存在解析的暂停时长 $t_{A,B}$，QAOA 的变分参数因此简化为”激光脉冲保持时长”。参数取 $V=2\pi \times 40\mathrm{MHz}$（${}^{87}\mathrm{Rb}$ $68S$、间距 $5\mu \mathrm{m}$），门误差为绝热误差与里德堡耗散的折中，亚微秒量级门时间对应百分之几的误差即可让 $p=3$ 层的奇偶 QAOA 在 20 个奇偶比特上稳健收敛；通过对被照明 plaquette 之间隔一行不照明原子的方式，整个多体相位层可在与系统尺寸无关的 9 轮照明内完成 Dlaska 2022。

2. 噪声感知 GRAPE 相位整形（$Z_N$ 奇偶门）。以 ${}^{88}\mathrm{Sr}$ 时钟态 $|1⟩={}^3{P}_0$、里德堡态 $|r⟩=5s61s{}^3{S}_1$（$C_6/h\approx -150\mathrm{GHz}\mu {\mathrm{m}}^6$）为例，${\Omega }_0/2\pi =10\mathrm{MHz}$、最小间距 $3\mu \mathrm{m}$ 给出 $V/\hslash {\Omega }_{\max }\approx 21$。只优化激光相位 $\phi (t)$（Rabi 幅度恒定）即可在线形、等边三角形、直角三角形、四面体、方形等构型上实现 $Z_2$ 到 $Z_4$；门时长约为 $N$ 个 Rabi 周期，$R_{\min }=2\mu \mathrm{m}$ 时总误差 $\sim 10^{-3}$ 量级、接近 $F=0.999$。原生 $Z_4$ 的误差谱显著偏向 $Z$ 型 Pauli 误差，而分解为单/双比特门的实现误差分布弥散——这种噪声偏置正是纠错码所需要的 Kazemi 2025 Kazemi 2025。

3. 空间阻塞工程的单控多靶门（CZZ/CZZZ）。把控制原子放在圆心、$N$ 个靶原子均布在圆周上：控制-靶距离短（强阻塞 $|V_{CT}|\gg \Omega$）、靶-靶距离长（弱阻塞 $|V_{TT}|\lesssim \Omega$），仅凭几何排布即打破置换对称性，无需额外激光或微波。GPU 加速的脉冲合成给出：CZZ 保真度最高 $99.55\%$（$575\mathrm{ns}$，半径 $3.5\mu \mathrm{m}$）、CZZZ 最高 $99.24\%$（$575\mathrm{ns}$，半径 $4.0\mu \mathrm{m}$），且已计入原子位置涨落与里德堡衰减；表面码稳定子 $Z^{\otimes 4}$ 与三角色码 $Z^{\otimes 6}$ 是直接动机 Stein 2026 Stein 2026。

4. 最优控制与几何相位 fan-out 门。基于非对称里德堡相互作用、时间最优控制与非常规几何相位（NHQC）的方案可在单控多靶结构上对靶比特执行任意单比特操作，三比特几何 FOT 门保真度 $99.51\%$、几何 FOX 门 $98.73\%$，并可用 FOX 门把 GHZ 态制备推进到 5 原子 Wei 2022。另一条面向纠错的全局脉冲路线用少数参数解析脉冲优化 $C_{n}Z$ 家族（含四比特 $CCCZ$），服务 Steane 码与 Floquet 码的稳定子读出 Locher 2026 Dlaska 2022。

四体里德堡门协议：四面体 plaquette 能谱随失谐的演化与两段绝热脉冲、宇称相关的动力学相位（图源：Dlaska 2022 Fig. 2）

非对称阻塞的几何排布：控制原子居中、靶原子圆周分布，C-T 强阻塞而 T-T 弱阻塞，实现 CZZ/CZZZ（图源：Stein 2026 Fig. 1）

参数与量级

量	典型值	来源
四体门相互作用 $V$	$2\pi \times 40\mathrm{MHz}$（Rb 68S，$5\mu \mathrm{m}$）	Dlaska 2022
四体门对应的 QAOA 误差容忍	门误差达百分之几仍稳健收敛（$p=3$，20 比特）	Dlaska 2022
Sr 奇偶门 $C_6/h$	$\approx -150\mathrm{GHz}\cdot \mu {\mathrm{m}}^6$（${}^{88}\mathrm{Sr}$ ${}^3{S}_1$）	Kazemi 2025
$Z_2$/$Z_4$ 门总误差	$8.8\times 10^{-4}$ / $4.0\times 10^{-3}$（$R_{\min }=2\mu \mathrm{m}$）	Kazemi 2025
$Z_4$ 门时长	${\Omega }_{0}T/2\pi =3.90$（约 4 个 Rabi 周期）	Kazemi 2025
CZZ 门保真度	$99.55\%$（$575\mathrm{ns}$，Rb）	Stein 2026
CZZZ 门保真度	$99.24\%$（$575\mathrm{ns}$，Rb）	Stein 2026
三比特几何 fan-out 门	FOT $99.51\%$ / FOX $98.73\%$	Wei 2022

与其他概念的关系

• 物理根基是里德堡阻塞与里德堡相互作用：阻塞约束给出块对角相位结构，$C_6/R^6$ 的可调性决定门型设计自由度。
• 对称的多控版本见多比特阻塞门与Toffoli/CCZ 门；非对称的单控多靶与奇偶门是对它们”打破置换对称”的推广。
• 奇偶门是对角幺正相位多项式的单项，与相位 Gadget、对角线路综合和多控相位门综合互为编译层面的对应物；Geyser 原生多比特门编译研究如何用原生多体门减少脉冲数。
• 面向纠错：快速 ${\hat{Z}}^{\otimes N}$ 读出直接降低三比特门容错稳定子测量的测量深度，多控门族被免测量纠错（Steane 码 $C_{2}Z$/$C_{3}Z$）与 Floquet 码读出采用，噪声偏置特性与擦除转换互补。
• 面向优化：四体奇偶门是里德堡模拟优化之外”数字式”解决组合优化的原生门资源，与qLDPC 架构同属”硬件-算法协同设计”的思路。
• 脉冲层面的噪声建模（里德堡衰减、反冲、范德瓦尔斯力致运动）见误差预算与真实建模与GRAPE 梯度上升脉冲工程；门相位约定与局域等价的处理见门相位与局域等价。

参考文献

• C. Dlaska, K. Ender, G. B. Mbeng, A. Kruckenhauser, W. Lechner, R. van Bijnen, “Quantum Optimization via Four-Body Rydberg Gates”, Phys. Rev. Lett. 128, 120503 (2022). NA076 · arXiv:2106.02663
• J. Kazemi, M. Schuler, C. Ertler, W. Lechner, “Multiqubit parity gates for Rydberg atoms in various configurations”, Phys. Rev. Research (2025). NA084 · arXiv:2506.09686
• D. F. Locher, J. Old, K. Brechtelsbauer, J. Holschbach, H. P. Büchler, S. Weber, M. Müller, “Multiqubit Rydberg Gates for Quantum Error Correction”, PRX Quantum (2026). NA086 · arXiv:2512.00843
• S. Stein, C. Liu, S. Kan, E. Crane, Y. Ding, Y. Mao, A. Schuckert, A. Li, “Multitarget Rydberg gates via spatial blockade engineering”, Phys. Rev. Research (2026). NA087 · arXiv:2504.15282
• J.-F. Wei, F.-Q. Guo, D.-Y. Wang, Y. Jia, L.-L. Yan, M. Feng, S.-L. Su, “Fast multiqubit Rydberg geometric fan-out gates with optimal control technology”, Phys. Rev. A 105, 042404 (2022). r5_fanout_optimal · DOI 10.1103/PhysRevA.105.042404