DRAG 与泄漏抑制

DRAG（Derivative Removal by Adiabatic Gate）的核心是一条解析规则：在主驱动包络之外，沿正交相位补一个与其时间导数成正比的分量，再把载频按包络平方微调——在绝热变换的框架下，这两个修正恰好逐阶消去流向泄漏能级的振幅，使短脉冲门不再被弱非谐性"出卖"。

问题背景

真实量子比特从来不是孤立二能级系统：超导 transmon 是弱非谐振子，里德堡原子在目标态 $|r⟩$ 附近还有一整族邻近 $n$ 态。当控制带宽与能级间隔（非谐性 $\Delta$）可比时，短脉冲必然把布居漏到计算子空间之外——这就是泄漏。压制泄漏的经典手段是拉长脉冲让频谱变窄，但门时间因此受”非谐性倒数”的下界约束，短门与低泄漏似乎不可兼得 Motzoi 2009。

Motzoi、Gambetta、Rebentrost 与 Wilhelm 于 2009 年给出的 DRAG 方案打破了这一权衡：不靠削窄频谱，而是主动补一个正交分量，让流向第三能级的振幅在每一时刻都相消。对超导比特的现实参数，它把门误差较当时最优的光滑包络再降一个数量级以上，且表明”弱非谐性原则上不是实现量子门的限制因素” Motzoi 2009。同一思想在里德堡门中的直接应用——用导数修正压制对邻近里德堡态的激发——见 Theis 2016（解析整形脉冲）。

理论模型

三能级模型与泄漏标度

考虑主跃迁 $|0⟩\leftrightarrow |1⟩$（拉比频率 $E_x(t)$）与泄漏跃迁 $|1⟩\leftrightarrow |2⟩$（耦合比 $\lambda ={\Omega }_{12}/{\Omega }_{01}$，非谐性 $\Delta$）。纯幅度调制的 Gaussian 脉冲门误差在短门极限下 $\propto E_{x,\max }^2/\Delta$；数值上 6 ns 的 $\pi$ 门已有约 $1\%$ 误差 Motzoi 2009。

绝热框架与逐阶相消

引入依赖脉冲的幺正变换（帧变换）

$$V(t)=\exp \left[-i\frac{E_x(t)}{2\Delta }\left({\hat{\sigma }}_y^{01}+\lambda {\hat{\sigma }}_y^{12}\right)\right],$$

其中 $E_x/\Delta$ 是绝热性参数；要求 $V(0)=V(t_g)=1$（即 $E_x(0)=E_x(t_g)=0$）使帧变换对门无净效应。在新帧中，有效哈密顿量除主驱动外多出”惯性项” ${\dot{E}}_x({\hat{\sigma }}_y^{01}+\lambda {\hat{\sigma }}_y^{12})/2\Delta$ 与交流斯塔克项——正是它们把布居甩进 $|2⟩$。选取正交分量与失谐

$$E_y(t)=-\frac{{\dot{E}}_x(t)}{\Delta },{\delta }_1(t)=\frac{({\lambda }^2-4)E_x^2(t)}{4\Delta },$$

即可把 $H_V$ 中的惯性项与相移误差逐项消去，泄漏被压到 $E_x^4/{\Delta }^3$ 阶；更高阶修正（含 $E_x$ 重标定、$E_y$ 的 $E_x^2{\dot{E}}_x$ 项等）可继续写出 Motzoi 2009。物理图像：DRAG 的正交分量让态矢量在布洛赫球上走一条”带倾斜进动轴”的路径，泄漏通道的激发振幅与主通道相位锁定而相消。

系统化的导数展开

Motzoi 与 Wilhelm 把 DRAG 推广为普适的导数跃迁压制框架：对三类离共振效应——(i) 同时被驱动的非目标跃迁（驱动串扰）、(ii) 未驱动的额外跃迁（能级梯子泄漏）、(iii) 波形边带（旋转波近似修正）——统一用”试函数（通常 Gaussian）的各阶导数展开”表示修正，并给出绝热误差与受控展开中各阶非绝热/受激误差的解析系数（一阶系数如 $b_1=-1/\Delta$），从而对密集谱系统系统地提升频率选择性 Motzoi 2013。适用边界也随之清晰：导数修正假设泄漏通道可按频点逐条处理；若门快到连”频点”都分辨不开，就应改走绝热路线（ARP 门、STIRAP 与暗态门）或正面数值优化（GRAPE）。

包络选择与频域直觉

DRAG 并不排斥传统包络，而是与其互补。Gaussian 包络的频率响应带宽窄，在泄漏跃迁频率处的激发本就小；tangential（$\tanh$ 升降沿）包络则以明确的上升时间参数换取对边沿的精细控制。在足够长的门时里，光滑包络已经足够好；DRAG 的价值集中在短门区——那里包络频谱压不住 $\Delta$ 量级的泄漏频率，主动相消成为唯一出路 Motzoi 2009。数值优化的极限行为佐证了这一点：把 GRAPE（逐像素自由优化）约束到 1 ns 像素并以 Gaussian 为初值，收敛解的主分量仍是 Gaussian，只是长出了导数形状的正交分量——DRAG 是数值最优解的解析骨架 Motzoi 2009。

Magnus 展开视角

Ribeiro、Baksic 与 Clerk 提出基于 Magnus 展开的系统微扰方法：对”理想幺正 + 小误差项”的演化修正控制脉冲，使目标演化在展开的指定阶上得到复原。它覆盖绝热演化中的非绝热误差与任意杂散耦合泄漏，且不要求泄漏振幅逐时刻为零，只要求门末态正确——这放宽了 DRAG 的瞬时相消条件，在其失效的参数区仍然可用 Ribeiro 2017。

方法与实现

波形实现：DRAG 只需在 IQ 调制中加一条正交通道（$E_y\propto {\dot{E}}_x$）与一个随 $E_x^2$ 变化的微小的载频/相位修正，脉冲仍由 Gaussian 包络加几个解析系数描述。逐像素数值优化（GRAPE，见GRAPE）验证过：以 Gaussian 为初值、1 ns 像素约束下，优化收玫的 $E_y$ 形状正是 Gaussian 的导数——DRAG 是”把数值优化解解析化” Motzoi 2009。

量化收益：6 ns 高斯 $\pi$ 门的 $\sim 1\%$ 误差经 DRAG 修正降低数个数量级（同门时约 4 个数量级），且结果对包络形状不敏感——只要面积正确、首尾归零 Motzoi 2009。

泄漏的量化与基准：泄漏不是普通保真度能完全刻画的——它包含可逆的”漏出-回流”（seepage）与不可逆的离开。Wood 与 Gambetta 给出把泄漏率与回流率同平均门保真度并列的度量和三者的修正随机基准估计协议 Wood 2018；Chen 等在超导 Xmon transmon 上把随机基准与泄漏布居测量结合，系统表征单比特门：单纯 DRAG 存在”总保真度 vs 泄漏”的权衡，叠加额外的脉冲整形（含失谐微调）后门误差稳定低于 $10^{-3}$、泄漏率降至 $10^{-5}$ 量级，剩余泄漏主要来自非相干加热 Chen 2016。这套量化工具与门标定闭环衔接：泄漏率超标指向通道/漂移问题，回流率超标指向控制波形。

里德堡移植：里德堡门的泄漏通道是目标里德堡态附近的 $n\pm 1$、$n^{\prime }$ 家族。NA057 把”导数展开压频点”的 DRAG 逻辑推广到多条泄漏跃迁：$\pi$ 脉冲同时压制两条通道（$N=4$ 阶），$2\pi$ 脉冲压一条（$N=2$ 阶），对 Cs 将贝尔保真度推到 $F_B>0.9999$（50 ns 门、300 K）Theis 2016。与超导情形的差别在于泄漏通道不止一条、且”阻塞泄漏”经由里德堡相互作用的双激发通道进入，因此修正阶数与失谐配置需按能级图重新设计。

DRAG 波形示意：Gaussian 主分量 E_x(t) 与正交导数分量 E_y(t) ∝ dE_x/dt，及泄漏布居对比（图源：Theis 2016 Fig. 2）

参数与量级

量	典型值	来源
未修 Gaussian 短门误差	$\sim 1\%$（$\pi$ 门，$t_g=6\mathrm{ns}$），$\propto E_x^2/\Delta$	Motzoi 2009
DRAG 后同门时误差	降低约 4 个数量级（$t_g=6\mathrm{ns}$）	Motzoi 2009
正交分量	$E_y=-{\dot{E}}_x/\Delta$；失谐 ${\delta }_1=({\lambda }^2-4)E_x^2/4\Delta$	Motzoi 2009
泄漏残余阶数	一阶修正后 $\sim E_x^4/{\Delta }^3$	Motzoi 2009
一阶导数系数（频选框架）	$b_1=-1/\Delta$	Motzoi 2013
泄漏率/回流率 + 平均门保真度	修正 RB 协议可同时稳健估计三者	Wood 2018
里德堡 π/2π 脉冲压制阶数	$\pi$：$N=4$（双通道）；$2\pi$：$N=2$（单通道）	Theis 2016
里德堡门收益	$F_B>0.9999$（50 ns，300 K，Cs）	Theis 2016

与其他概念的关系

• DRAG 的”初始形状 + 导数修正”模板在里德堡多通道版本即解析整形脉冲；后者还叠加了最优阻塞能级设计。
• 逐时刻相消并非唯一目标：只约束门末态的 Magnus 系统化方法与放宽约束的优化目标属于鲁棒最优控制的谱系；绝热极限的对照方案见ARP 门与STIRAP 与暗态门。
• 泄漏率/回流率的度量与修正随机基准属于门标定闭环（ORBIT/RB/Ad-HOC）的工具箱；泄漏进入误差预算的全口径分解见误差预算与真实建模。
• 当泄漏通道是相互作用导致的阻塞泄漏时，问题与阻塞 CZ 门、多比特阻塞门的门设计耦合。
• 波形解析系数的普适构造思想（对参数求导、逐阶展开）与GOAT 生成轨道法同源。
• 泄漏谱的物理来源（弱非谐梯子 vs 里德堡 $n$ 族）与里德堡原子的能级结构、硬件全景中的控制链路带宽共同决定该用几阶修正。
• 与”哪类误差可被闭环诊断”的问题相关：可累积的相干误差指向校准，计算子空间外布居指向泄漏通道与回收动力学——诊断逻辑见门标定闭环。

参考文献

• F. Motzoi, J. M. Gambetta, P. Rebentrost, F. K. Wilhelm, “Simple Pulses for Elimination of Leakage in Weakly Nonlinear Qubits”, Phys. Rev. Lett. 103, 110501 (2009). Motzoi 2009 · DOI 10.1103/PhysRevLett.103.110501
• F. Motzoi, F. K. Wilhelm, “Improving frequency selection of driven pulses using derivative-based transition suppression”, Phys. Rev. A 88, 062318 (2013). Motzoi 2013 · DOI 10.1103/PhysRevA.88.062318
• H. Ribeiro, A. Baksic, A. A. Clerk, “Systematic Magnus-Based Approach for Suppressing Leakage and Nonadiabatic Errors in Quantum Dynamics”, Phys. Rev. X 7, 011021 (2017). Ribeiro 2017 · DOI 10.1103/PhysRevX.7.011021
• C. J. Wood, J. M. Gambetta, “Quantification and characterization of leakage errors”, Phys. Rev. A 97, 032306 (2018). Wood 2018 · DOI 10.1103/PhysRevA.97.032306
• Z. Chen et al., “Measuring and Suppressing Quantum State Leakage in a Superconducting Qubit”, Phys. Rev. Lett. 116, 020501 (2016). Chen 2016 · DOI 10.1103/PhysRevLett.116.020501
• L. S. Theis, F. Motzoi, F. K. Wilhelm, M. Saffman, “High-fidelity Rydberg-blockade entangling gate using shaped, analytic pulses”, Phys. Rev. A 94, 032306 (2016). NA057 · DOI 10.1103/PhysRevA.94.032306