鲁棒最优控制

鲁棒最优控制的要点是"不要等优化完再扫误差图"：在目标泛函里同时传播一组扰动哈密顿量、参数灵敏度态，或直接极小化保真度对误差的导数，让脉冲在名字被冠以"最优"之前就先对漂移免疫。

问题背景

量子最优控制框架给出的名义最优脉冲只对一个理想哈密顿量最优。真实实验里，双光子失谐随激光频率漂移、Rabi 频率随光强起伏（shot-to-shot 波动）、里德堡能级随杂散电场移动、相互作用强度随原子间距涨落——名义点上 $99.99\%$ 的门在 $\pm 1\%$ 的参数偏移下可能跌到 $99\%$。事后扫描一张”误差-参数”图并声称平坦，并不能改变脉冲本身的脆弱性；正确做法是把鲁棒性作为约束或惩罚项放进优化问题本身。

实现”免疫”有两条互补的路线。模型内鲁棒化：在仿真中扩展目标泛函或状态空间（ensemble、灵敏度、最坏情形），用GRAPE、Krotov/GOAT 或CRAB 等优化器求解；控制序列鲁棒化：不改泛函，而用复合脉冲、动力学解耦等滤波式设计把特定类型的误差在传播中抵消。两条路线最后都要在门标定闭环中接受独立验证。

理论模型

Ensemble 目标泛函

最直接的鲁棒化是”多哈密顿量并行传播”：对参数空间的扰动样本 $\{{\hat{H}}_n{\}}_{n=0}^{N-1}$（${\hat{H}}_0$ 为名义哈密顿量）同时演化全部计算基态，要求每个成员都到达目标门。Goerz 等对带单点寻址的里德堡 CPHASE 门给出的目标泛函为

$$J=1-\frac{1}{16N}\sum _{n=0}^{N-1}\sum _{k=1}^{16}\mathrm{Re}\mathrm{tr}[\hat{O}{\hat{\rho }}_k(0){\hat{O}}^{†}{\hat{\rho }}_{k,n}(T)]+\sum _j{\int }_0^T\frac{[{ϵ}_j(t)-{ϵ}_{j,\mathrm{ref}}(t){]}^2}{S(t)}\mathrm{d}t,$$

其中 $\hat{O}$ 是目标门，${\hat{\rho }}_{k,n}(T)$ 是第 $k$ 个基态在成员 $n$ 下的末态，第二项为 Krotov 型流量惩罚（$S(t)$ 为形状函数）Goerz 2014。传播在 Liouville 空间进行，因而可以显式包含中间态自发辐射等耗散——这与STIRAP 与暗态门中”衰减项与非绝热项竞争”的开放系统结构一致。ensemble 目标的代价是传播次数 $\times N$：这是鲁棒性换取数值成本的直接汇率。

灵敏度惩罚：把导数变成状态

若不确定参数为 $ϵ$，定义灵敏度态 $|s(t)⟩={\partial }_{ϵ}|\psi (t;ϵ)⟩{|}_{ϵ=0}$，对薛定谔方程求导得

$$i\frac{\mathrm{d}}{\mathrm{dt}}|s(t)⟩=H(t;0)|s(t)⟩+\left({\partial }_{ϵ}H\right){|}_0|\psi (t;0)⟩,$$

于是”对 $ϵ$ 的一阶不敏感”等价于压低 $\int ⟨s|s⟩\mathrm{d}t$：灵敏度像布居一样沿一条扩展的动力学方程传播，可作为约束或惩罚进入优化器 Mohan 2023。Mohan 等在 Sr-88 单光子里德堡门模型中用轨迹优化器 ALTRO 实现：把门相位、布居回归写成硬约束，把 $\dot{\Omega }$、$\dot{\Delta }$ 作为扩展控制并加以惩罚以偏好平滑波形。注意其边界：惩罚一阶导数有助于抑制快变，但不等于对真实电子学传递函数施加了严格频带限制；该工作主要处理 shot-to-shot 型准静态偏差，不能直接解读为对任意激光噪声谱都有效。

保真度易感性与普适鲁棒控制

不指定误差类型时，Poggi 等证明：一阶误差由保真度易感性（fidelity susceptibility）刻画，它可以写成超算符形式，且极小化时间平均算符的范数

$$J_V=\frac{1}{d}\Vert \bar{V}{\Vert }^2$$

（$\bar{V}$ 为相互作用绘景中误差算符的时间平均）给出的解等价于让演化的一阶矩匹配 Haar 分布——构成一个 1-design，从而对任意未知的系统型误差同时鲁棒 Poggi 2024。与之呼应，Fromonteil 等按误差通道给里德堡门协议分级：纠缠保真度易感性 ${\chi }_C=0$ 称条件鲁棒，泄漏易感性 ${\chi }_P=0$ 称泄漏鲁棒，两者同时为零称完全鲁棒；把对强度误差鲁棒的 $C_{\pi /2}$ 序列执行两遍即可在保持条件鲁棒的同时获得对强度误差的泄漏鲁棒性 Fromonteil 2023。

最坏情形与区间目标

对指定的误差区间 $[-{ϵ}_0,{ϵ}_0]$，可以让优化器直接压平保真度曲线。Hou 等的”主动失谐鲁棒”（DER）目标是

$$J_{\mathrm{DER}}=[1-F(0){]}^2+[F_{\max }({ϵ}_{\delta })-F_{\min }({ϵ}_{\delta }){]}^2,{ϵ}_{\delta }\in [-{ϵ}_0,{ϵ}_0],$$

第一项守住名义点，第二项压制区间内起伏 Hou 2024。使用时要同时检查区间最低值：曲线平坦不等于保真度高。这类目标与本条其他目标一样，应当用与训练集不同的独立样本验证，避免”在采样点上完美”的过拟合。

滤波函数与复合脉冲：控制序列路线

时间相关噪声下的平均误差可用滤波函数框架估计：$1-F\simeq \int S(f)I(f)\mathrm{d}f$（噪声谱 $\times$ 协议响应），且 Magnus 展开可以推进到高阶——Green、Uys 与 Biercuk 给出了非对易控制下任意阶（至四阶）系综平均纠缠保真度的滤波函数表达式 Green 2012。复合脉冲本是 NMR 中纠正静态系统误差的设计，Kabytayev 等证明级联复合脉冲与动力学修正门对时间相关的非马尔可夫幅度/失谐噪声同样有一阶抑制能力，适用频率高达 Rabi 频率的约 $10\%$，并在 $1/f$（高频滚降 $1/f^2$）噪声谱下与数值模拟吻合 Kabytayev 2014。这一路线与误差预算与真实建模中 Tsai 等的保真度响应理论是同一数学的两种用法：前者主动设计 $I(f)$ 的零点，后者用它诊断现役波形。

方法与实现

• Krotov + ensemble（里德堡）：Goerz 等取 $N=24$ 个哈密顿量成员（里德堡能级移动 $\pm 300$ kHz、偶极耦合强度 $\pm 5\%$），以解析方案的双 $\pi$/ $2\pi$ 脉冲为初值，在 Liouville 空间（含中间态寿命 150 ns 的自发辐射）中优化，门时长 $T=800$ ns；对两类涨落同时鲁棒的解给出门误差 $1.04\times 10^{-4}$，并预期现实条件下 $99.90-99.99\%$ 的门保真度可行 Goerz 2014。优化器输出的脉冲出现第二频率分量，靠多路径干涉抑制中间态布居——鲁棒解往往不是”更平滑”，而是”更聪明”。
• ALTRO 轨迹优化（Sr）：Mohan 等比较两组约束：方案 A 取 $T=10/{\Omega }_{\max }$、$|\Delta (t)|\le 1.25{\Omega }_{\max }$；方案 B 取 $T=15/{\Omega }_{\max }$、$|\Delta (t)|\le 2.25{\Omega }_{\max }$——用更多时间与控制范围换取更宽的鲁棒窗口；该文的另一贡献是把里德堡主量子数 $n$（决定寿命、阻塞强度与极化率）与波形一起纳入选择 Mohan 2023。
• 谱约束 steering：Palao、Reich 与 Koch 表明，对脉冲频谱加滤波型约束（谱约束）可与单调收敛算法兼容，用于排除共振单光子通道、压制多余谱峰——这是在”解族”层面引导优化路径的手段，约束不同，收敛到的解族不同 Palao 2013。
• 鲁棒设计不必然赢：一个否定结果。Tsai 等在实测激光噪声谱下比较时间最优门、静态强度噪声鲁棒门（Fromonteil 型）与 CRAB 优化：Fromonteil 协议在零频处的强度响应确实为零，但在该实验的有限频率噪声与衰减共同作用下，时间最优门给出最佳预测性能 Tsai 2025。教训：鲁棒性只在被写进目标的那类误差上兑现；先量化误差预算（误差预算与真实建模），再决定为哪一项买单。

参数与量级

量	典型值	来源
Ensemble 优化的里德堡 CPHASE	误差 $1.04\times 10^{-4}$（$T=800$ ns，$N=24$）	Goerz 2014
现实条件下的可行门保真度	$99.90-99.99\%$	Goerz 2014
Ensemble 采样窗口	能级移动 $\pm 300$ kHz，耦合 $\pm 5\%$	Goerz 2014
ALTRO 方案 A/B（Sr）	$T=10/{\Omega }_{\max }$、$|\Delta |\le 1.25{\Omega }_{\max }$；$T=15/{\Omega }_{\max }$、$|\Delta |\le 2.25{\Omega }_{\max }$	Mohan 2023
复合脉冲适用的噪声频率上限	$\sim 10\%$ Rabi 频率（时间相关非马尔可夫噪声）	Kabytayev 2014
滤波函数高阶展开	系综纠缠保真度至四阶（Magnus 展开）	Green 2012
DER 双 STIRAP CZ	$F\approx 0.9971$（含中间态衰减），精修后预测 $>0.9999$	Hou 2024
实测噪声下的对比	时间最优门优于鲁棒门（该实验噪声谱内）	Tsai 2025

与其他概念的关系

• 目标泛函、梯度与开放系统传播的公共底座见量子最优控制框架；GRAPE 与Krotov/GOAT 是承载 ensemble/灵敏度目标的两类主力优化器，CRAB 的少参数化则便于把鲁棒窗口当作外层参数扫描。
• 时间最优与鲁棒性是一对张力：量子速度极限约束了”多慢才够绝热”，见量子速度极限；解析方案的鲁棒性对照见解析整形脉冲、ARP 门与STIRAP 与暗态门。
• 滤波函数视角的另一半——用实测噪声谱写预算——在误差预算与真实建模；把鲁棒窗口交由实验闭环微调则进入门标定闭环。
• 复合脉冲与动力学解耦的误差抵消思想同样适用于单比特层：微波/Raman 驱动的标定与验证见中性原子单比特门，泄漏通道的处理见DRAG 与泄漏抑制。

参考文献

• M. H. Goerz, E. J. Halperin, J. M. Aytac, C. P. Koch, K. B. Whaley, “Robustness of high-fidelity Rydberg gates with single-site addressability”, Phys. Rev. A 90, 032329 (2014). Goerz 2014 · DOI 10.1103/PhysRevA.90.032329
• M. Mohan, R. de Keijzer, S. Kokkelmans, “Robust control and optimal Rydberg states for neutral atom two-qubit gates”, Phys. Rev. Research 5, 033052 (2023). Mohan 2023 · DOI 10.1103/PhysRevResearch.5.033052
• P. M. Poggi, G. De Chiara, S. Campbell, A. Kiely, “Universally Robust Quantum Control”, Phys. Rev. Lett. 132, 193801 (2024). Poggi 2024 · DOI 10.1103/PhysRevLett.132.193801
• T. J. Green, H. Uys, M. J. Biercuk, “High-Order Noise Filtering in Nontrivial Quantum Logic Gates”, Phys. Rev. Lett. 109, 020501 (2012). Green 2012 · DOI 10.1103/PhysRevLett.109.020501
• C. Kabytayev, T. J. Green, K. Khodjasteh, M. J. Biercuk, L. Viola, K. R. Brown, “Robustness of composite pulses to time-dependent control noise”, Phys. Rev. A 90, 012316 (2014). Kabytayev 2014 · DOI 10.1103/PhysRevA.90.012316
• J. P. Palao, D. M. Reich, C. P. Koch, “Steering the optimization pathway in the control landscape using constraints”, Phys. Rev. A 88, 053409 (2013). Palao 2013 · DOI 10.1103/PhysRevA.88.053409
• X. Fromonteil, D. Bluvstein, H. Pichler, “Protocols for Rydberg Entangling Gates Featuring Robustness against Quasistatic Errors”, PRX Quantum 4, 020335 (2023). Fromonteil 2023 · DOI 10.1103/PRXQuantum.4.020335
• Q.-L. Hou, H. Wang, J. Qian, “Active robustness against detuning error for Rydberg quantum gates”, Phys. Rev. Applied 22, 034054 (2024). Hou 2024 · DOI 10.1103/PhysRevApplied.22.034054
• R. B.-S. Tsai, X. Sun, A. Shaw, R. Finkelstein, M. Endres, “Benchmarking and Fidelity Response Theory of High-Fidelity Rydberg Entangling Gates”, PRX Quantum 6, 010331 (2025). Tsai 2025 · DOI 10.1103/PRXQuantum.6.010331