量子速度极限

量子速度极限（quantum speed limit, QSL）回答一个问题：在给定的相互作用与陷阱物理下，一个门最快能多快？它由两种时间尺度中的较大者决定——相互作用积累受控相位所需的时间，以及门结束时必须恢复的运动（振动）自由度的时间；最优控制则把这一极限变成可操作的数值判据：把门时间一压再压，直到优化算法再也找不到高保真解。

问题背景

“门越快越好”在中性原子里有直接的物理动机：门时间短，里德堡态衰变、激光相位噪声、多普勒失谐等技术误差的积累时间就短（见误差预算与真实建模）。但加速不是免费的：受控相位 $\chi$ 需要相互作用能 $V$ 在时间 $T$ 上积累（量级上 $\chi \sim VT/\hslash$），而把原子瞬间激发到相互作用强的态会通过电子–核运动耦合把振动激发打进陷阱——门结束时若不能把运动态复原，量子比特寄存器就”漏”了。

Goerz、Calarco 与 Koch 对光晶格中两个钙原子的受控相位门做了首个系统性研究：激光共振把 $|00⟩$ 激发到相互作用强的电子激发态（偶极相互作用 $1/R^3$），用最优控制计算脉冲形状，并完全计入相互作用的 $R$ 依赖。他们发现门时间由两种因素竞争决定：要么受激发态相互作用强度限制，要么受陷阱基态振动运动的恢复限制，“两者中给出较大时间的那个”说了算 Goerz 2011。这正是 QSL 的物理内容：它不是单一公式，而是一组物理约束下的最小值问题。

理论模型

从态转移界到门时间界

标准量子速度极限界给出态从 $|{\psi }_0⟩$ 演化到正交态的最短时间，由平均能量与能量不确定度表示（Mandelstam–Tamm 与 Margolus–Levitin 两类界）；幺正门目标可化为若干态到态转移的约束后套用同一思想。NA055 指出：解析求值一般只有简单模型系统才可行，数值上可以把”优化收敛的崩溃点”定义为 QSL——先把门时间设得足够长获得高保真解，再逐步压缩 $T$，直到优化算法无法再找到高保真实现为止 Goerz 2011。这一”操作型定义”把 QSL 与量子最优控制框架直接连在一起。

两类时间尺度的竞争

在 NA055 的模型中，两原子处于间距 $d$ 的谐振陷阱，激光驱动 $|0⟩\to |a⟩$（电子激发态），激发态上有偶极–偶极相互作用 $-C_3/R^3$。非局域相位定义为 $\chi ={\varphi }_{00}-{\varphi }_{01}-{\varphi }_{10}+{\varphi }_{11}$（见门相位与局域等价）。相互作用越强，相位积累越快：对 Ca 原子的 $B^1{\Sigma }_u^{+}$ 势（$C_3=16.04a.u.$），在 $d=5\mathrm{nm}$ 处相互作用能约 $4\mathrm{c}{\mathrm{m}}^{-1}$，足以在几十皮秒内积累 $\pi$ 相位；但要把 $d=200\mathrm{nm}$ 的标准 UV 晶格调到同样相互作用能，需要 $C_3\sim 10^{6}a.u.$——已接近高激发里德堡态的 $C_3\sim 3\times 10^{6}a.u.$ 量级（对应 $4\mu \mathrm{m}$ 间距约 $1.3\times 10^{-3}\mathrm{c}{\mathrm{m}}^{-1}$）Goerz 2011。

代价出现在运动自由度：共振激发使波包在激发态势能面上加速，振动激发随相互作用增强而放大。数值实验显示，把 $C_3$ 提到 $4\times 10^{8}a.u.$ 虽可获得 $\pi$ 非局域相位，但振动保真度（末态在振动目标态上的投影 $F_{00}$）跌到 $75\%$ 以下——因为 $0.5-1\mathrm{ps}$ 的门远短于陷阱振动周期 $T_v$，优化算法根本无法分辨并恢复陷阱本征态。要恢复运动态，陷阱频率须不低于 $250\mathrm{kHz}$（对应 $T_v\approx 2\mathrm{ns}$），而任何现实的陷阱频率都给出远大于皮秒的运动时间尺度 Goerz 2011。结论：对中性原子，超快里德堡门的时间下界最终可能由陷阱运动而不是相互作用设定，除非采用对振动激发完全不敏感的门方案。

“三种最优”不可混用

时间最优解只是参照系之一。Locher 等在 CRAB 框架下区分了三种”最优”：时间最优（最小化门时长 $T$）、里德堡时间最优（最小化里德堡总占据时间 $T_R={\int }_0^T{\sum }_j⟨n_{r,j}(t)⟩dt$）与参数最优（在给定 ansatz 内最小化参数数），三者互不等价。在理想阻塞 CZ 的具体算例中：时间最优解为 $({\Omega }_{0}T,{\Omega }_0{T}_R)=(7.611,2.958)$，而偏重衰减优化的六参数解为 $(7.725,2.936)$——门长了 $1.5\%$，里德堡占据反而略降。由于 $T_R$ 直接进入寿命误差 $\Gamma T_R$，“哪个最优”取决于误差预算里寿命与相干性谁是主导 Locher 2026。这与 NA077 的衰减权重优化（$\alpha \to 1300$）是同一权衡的两种参数化。

有限阻塞下的时间最优解结构

对阻塞型门（阻塞 CZ 门），Jandura 与 Pupillo 用 GRAPE 系统扫描无量纲门时长 $T{\Omega }_{\max }$，发现最小门时长 $T^{\ast }$ 附近门误差按二次律打开：

$$1-F=A(T^{\ast }-T{)}^2(T<T^{\ast }),1-F\approx 0(T\ge T^{\ast }),$$

拟合给出 $T^{\ast }{\Omega }_{\max }=7.612$、$A=0.0544$ Jandura 2022。这个二次打开规律是 QSL 逼近的典型景观形态：把 $T$ 压到极限附近的代价是缓增的，但越界后无解。他们进一步用庞特里亚金极大值原理（PMP，见Krotov 与 GOAT 方法的变分视角）证明：该时间最优脉冲不是 bang-bang 型，而是幅度恒定、相位连续光滑变化的解，可由 4 个（CZ）/6 个（C₂Z）参数的常微分方程完全重构 Jandura 2022。

误差–时间的显式权衡

里德堡寿命 $\Gamma$ 与有限阻塞强度 $B$ 给出可直接优化的门误差表达式 Jandura 2022：

$$1-F=\frac{\Gamma (T_R{\Omega }_{\max })}{{\Omega }_{\max }}+\frac{{\Omega }_{\max }^2\alpha }{B^2(T{\Omega }_{\max }{)}^2},$$

其中 $T_R$ 是里德堡占据时间积分、$\alpha$ 是脉冲的阻塞敏感系数。第一项随门变快而改善（占据时间短），第二项随门变快而恶化（需要更大的 ${\Omega }_{\max }$ 才能在更短时间内转完相位）——存在一个把两项平衡的最优 ${\Omega }_{\max }$。这说明”最快”与”最保真”并不重合：时间最优解是速度极限意义上的参照点，实际工作点由完整误差预算（见鲁棒最优控制）决定。

方法与实现

数值协议（NA055 路线）：以相位敏感的态到态目标函数做 Krotov 型优化 → 固定门时间收敛到 $\Delta F<10^{-4}$ → 步进压缩 $T$ 重复优化 → 收敛崩溃即到达该模型参数下的 QSL Goerz 2011。

半解析协议（NA077 路线）：先用 GRAPE（数百个分段常数参数，CZ 用 99 片、C₂Z 用 399 片）在 $B=\infty$ 下扫描 $T{\Omega }_{\max }$ 找到时间最优解；再用 PMP 把无限维控制景观约化成 4 维/6 维参数空间，用 Table 1 的初始协态即可不经 GRAPE 复现脉冲 Jandura 2022。关键结果：

• CZ：$T{\Omega }_{\max }=7.6114828$，$1-F=3.1\times 10^{-10}$（理想阻塞）；全局脉冲仅调相位、幅度恒定，比 Levine 等的方案（$T{\Omega }_{\max }=8.585$，见 Levine 2018）快约 $10\%$——而且即使放宽到允许单点寻址的脉冲，全局解仍是时间最优的：单点寻址对两比特门并无速度优势 Jandura 2022。
• C₂Z（三比特）：找到两个定性不同、时间相差 $<1\%$ 的时间最优全局脉冲，$T{\Omega }_{\max }=16.426439$ 与 $16.532211$，$1-F$ 分别为 $3.1\times 10^{-7}$ 与 $2.8\times 10^{-6}$ Jandura 2022；原生三比特门的时间下界研究见Toffoli/CCZ 门。
• 衰减鲁棒化：改以里德堡占据最小化为目标后，阻塞敏感系数 $\alpha$ 从时间最优值下降约 $30\%$ 收敛到 $\alpha \approx 1300$（$T{\Omega }_{\max }=60$ 处），以少量时间换寿命预算 Jandura 2022。

实验落地（Evered 等，87Rb）：把 Jandura 型时间最优思想参数化为相位轮廓 $\phi (t)=A\cos (\omega t-\varphi )+{\delta }_{0}t$，时间最优 CZ 的 $T=275\mathrm{ns}$（$A/2\pi =0.113$、$\omega /2\pi =1.038\mathrm{MHz}$）；成对并行 CZ 基准保真度 $99.54(2)\%$，60 比特并行执行 $99.48(2)\%$；原生 CCZ 门时间仅比时间最优 CZ 长 $44\%$，三比特 GHZ 保真度一致性给出 $F=97.9(2)\%$ Evered 2023。这是”QSL 参照点 → 少参数 ansatz → 全局脉冲并行执行”完整链条的实验演示。

时间最优 CZ 脉冲：恒定幅度包络与连续光滑相位轮廓 φ(t)，及其在 T 附近门误差 1-F=A(T*-T)^2 的打开曲线（图源：Jandura 2022 Fig. 1）*

参数与量级

量	典型值	来源
CZ 时间最优时长（理想阻塞）	$T{\Omega }_{\max }=7.611$（$1-F=3.1\times 10^{-10}$）	Jandura 2022
C₂Z 时间最优时长	$T{\Omega }_{\max }=16.43/16.53$（两条解）	Jandura 2022
最小门附近误差打开	$1-F=A(T^{\ast }-T{)}^2$，$T^{\ast }{\Omega }_{\max }=7.612$，$A=0.0544$	Jandura 2022
Levine 方案对照	$T{\Omega }_{\max }=8.585$（时间最优全局解比它快约 10%）	Levine 2018
Cs 里德堡寿命（300 K）	$1/\Gamma =540\mu \mathrm{s}$（$\Vert 107P_{3/2}⟩$）	Jandura 2022
实验时间最优 CZ	$T=275\mathrm{ns}$，保真度 $99.54(2)\%$	Evered 2023
晶格钙原子振动周期	$T_v\approx 2\mathrm{ns}$（$\nu \ge 250\mathrm{kHz}$，$d=200\mathrm{nm}$）	Goerz 2011
振动控制失效阈值	$C_3\gtrsim 4\times 10^{8}a.u.$ 时 $F_{00}\le 75\%$（$T=0.5-1\mathrm{ps}$）	Goerz 2011

与其他概念的关系

• QSL 的”操作型定义”（压缩 $T$ 直到优化不收敛）是量子最优控制框架的一部分；PMP 半解析约化与Krotov/GOAT的变分思想同源，GRAPE 实现见GRAPE。
• 时间最优解是阻塞 CZ 门与Toffoli/CCZ 门的参照下界；门相位的局域等价处理见门相位与局域等价。
• “最快”≠“里德堡占据最少”：时间最优/里德堡时间最优/参数最优的区分见CRAB/DCRAB的少参数 ansatz 现代实践。
• “最快”≠“最保真”：里德堡占据 $T_R$ 与阻塞泄漏的显式权衡把本条与鲁棒最优控制、误差预算与真实建模连接起来。
• 固定幅度、光滑相位这一时间最优解的函数形态，被解析整形脉冲与ORMD/ARMD的少参数 ansatz 继承；用随机基逼近极限附近解的另一路线见CRAB/DCRAB。
• 门时间进入寿命误差预算 $\Gamma T_R$ 的方式与里德堡原子的寿命物理直接相关。

参考文献

• M. H. Goerz, T. Calarco, C. P. Koch, “The quantum speed limit of optimal controlled phasegates for trapped neutral atoms”, J. Phys. B 44, 154011 (2011). NA055 · DOI 10.1088/0953-4075/44/15/154011
• 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
• H. Levine et al., “High-Fidelity Control and Entanglement of Rydberg-Atom Qubits”, Phys. Rev. Lett. 121, 123603 (2018). NA061 · DOI 10.1103/PhysRevLett.121.123603
• S. J. Evered et al., “High-fidelity parallel entangling gates on a neutral atom quantum computer”, Nature (2023). NA078 · DOI 10.1038/s41586-023-06481-y
• D. F. Locher et al., “Multiqubit Rydberg Gates for Quantum Error Correction”, PRX Quantum 7, 020354 (2026). NA086 · DOI 10.1103/j8fm-24cf