解析整形脉冲

解析整形脉冲用少数解析函数（Gaussian、Blackman 及其导数修正）替代方波与逐点自由波形：在频域里它们有指数衰减的旁瓣，在时域里首尾导数归零，于是有限阻塞与邻近里德堡态泄漏被压低数个数量级——而脉冲仍由几个解析参数完全刻画，标定与实验可行性与方波相当。

问题背景

早期里德堡阻塞门（见阻塞 CZ 门）使用恒定幅度的方波序列：控制原子 $\pi$ 脉冲–靶原子 $2\pi$ 脉冲–控制原子 $\pi$ 脉冲。方波在频域上是 sinc 型的，其旁瓣会共振驱动邻近里德堡态（如 Cs 中 $|n{p}_{3/2}⟩$ 附近的 $n^{\prime }p$、$(n\pm 1)p$ 家族），造成”主泄漏”与”阻塞泄漏”两类误差。对恒定幅度脉冲的详细分析给出，Rb/Cs 原子在 300 K 环境下的本征贝尔态保真度上限约 $F_B\sim 0.999$——这正是当年实验（$0.7-0.8$）与理论极限之间差距的结构性来源 Theis 2016。

出路有两条：一是数值最优控制（GRAPE、CRAB），保真度高但波形任意、难以标定；二是解析整形——Theis、Motzoi、Wilhelm 与 Saffman 证明，只需把方波换成平滑解析包络并做导数修正，同一门序列的贝尔保真度即可提升几个数量级：300 K 下 $F_B>0.9999$、门时间仅 $50\mathrm{ns}$（比此前方案快一个数量级），4 K 下 $F_B>0.99999$ Theis 2016。解析形状的优势在于：光滑（硬件带宽友好）、少参数（标定便宜）、可解释（每个系数对应一条被压制的跃迁），这些性质被后来的时间最优门与调制驱动方案继承（见量子速度极限与ORMD/ARMD）。

理论模型

频域压制判据

设初始脉冲 $f(t)$ 光滑且首端、末端前 $N=2m$ 阶导数为零，则跃迁到失谐为 $\delta$ 的泄漏态的激发幅度由有限傅里叶变换

$$S(f,\delta )={\int }_0^{T}f(t)e^{i\delta t}dt={\left(\frac{i}{\delta }\right)}^N{\int }_0^T{e}^{i\delta t}\frac{d^{N}f(t)}{d{t}^N}dt$$

给出（分部积分 $m$ 次）：泄漏幅度随失谐 ${\delta }^{-N}$ 快速衰减。若要在特定频率 ${\delta }_j$ 处实现完全相消（暗点），只需求解 $S({\epsilon }_x,{\delta }_j)=0$ Theis 2016。

DRAG 型导数展开

为此把控制波形写成初始光滑形状 ${\epsilon }_x^{(0)}(t)$（如 Gaussian）与其偶数阶导数的线性组合：

$${\epsilon }_x(t)={\epsilon }_x^{(0)}(t)+\sum _{k=1}^{N/2}{\alpha }_{2k}\frac{d^{2k}{\epsilon }_x^{(0)}(t)}{d{t}^{2k}},$$

由于 $S$ 对 $f$ 线性，压制 $m$ 条泄漏跃迁等价于线性方程组

$$1+\sum _{k=1}^{N/2}{\alpha }_{2k}(-i{\delta }_j{)}^{2k}=0,j=1,\dots ,m.$$

例如要同时压制两条失谐 ${\delta }_1,{\delta }_2$ 的跃迁，解出 ${\alpha }_2=({\delta }_1^2+{\delta }_2^2)/({\delta }_1^2{\delta }_2^2)$、${\alpha }_4=-1/({\delta }_1^2{\delta }_2^2)$。这些系数逐时刻消除泄漏振幅，可由严格的绝热变换迭代论证支撑 Theis 2016；该导数展开正是超导比特 DRAG 方法的里德堡推广（DRAG 原始机制见DRAG 与泄漏抑制）。

最优阻塞的解析设计

整形之外还有第二个解析自由度：阻塞移位 $B_0$ 本身可以选优。对 $n$ 相关的跃迁矩阵元 $\propto n^{-3/2}$ 与激发概率 $\propto {\Delta }^{-2}$，主泄漏通道的概率可写成 $B_0$ 的显式多项式（对各 $(n\pm 1)p_{3/2}$、$n^{\prime }p$、$n^{\prime \prime }p$ 家族求和），对 $B_0$ 求根即得最优阻塞。对 Cs 的设置 S2（$n=141$）解析最优值为 $B_0/2\pi \approx 0.68\mathrm{GHz}$，与数值扫描一致；设置 S1（$n=107$）的最优区间为 $0.7-2.7\mathrm{GHz}$ Theis 2016。这体现了”整形”的完整含义：波形与能级配置联合设计，而不是单边压波形。

方法与实现

脉冲构造流程（NA057）：选初始 Gaussian 包络 → 按 $\pi$/ $2\pi$ 脉冲分别取 $N=4$/$N=2$ 阶导数修正（$\pi$ 脉冲面对更多泄漏通道）→ 解线性方程组定 ${\alpha }_{2k}$ → 对靶原子再叠加”失谐 DRAG”修正相位（恒定失谐 $\Lambda$ 随 $1/{\tau }_t^2$ 标度，${\tau }_t=25\mathrm{ns}$ 时 $\Lambda /2\pi =124.07\mathrm{MHz}$）Theis 2016。

量化收益：在相同门时间下，Gaussian 包络把泄漏布居误差较方波压低最多 $2.5$ 个数量级；加上 DRAG 导数修正再降约 $1.5$ 个数量级。计入全部泄漏通道与 300 K 里德堡衰变后，贝尔态保真度 $F_B>0.9999$（门时间 $50\mathrm{ns}$，含两个 $\pi$ 脉冲与一个 $2\pi$ 脉冲）；4 K 环境 $F_B>0.99999$ Theis 2016。

解析包络家族的谱系：

• Gaussian / DRAG 系列：频域指数衰减，首尾归零，是 NA057 与 DRAG 与泄漏抑制的主体；文献中把高斯两翼专门削平的变体（winglet 形脉冲）也属此家族。
• Blackman 包络：窗函数型 tapering，频谱旁瓣衰减极快，被中性原子单比特微波门 Wang 2016、Pulser 平台内置波形 Pulser 2022 与时间最优门实验的 AOM 建模 Evered 2023 广泛采用。
• 对称绝热解析脉冲：Saffman 给出 $\tanh$ 型升降沿加解析失谐的对称 CZ 方案，绝热跟踪双原子暗态，属”解析可写、绝热鲁棒”路线 Saffman 2020，见ARP 门与STIRAP 与暗态门。

标定接口：解析形状的少参数性天然适配闭环标定——包络由几个系数唯一决定，实验上只需扫描/拟合这几个量（Gaussian 宽度、导数修正系数、失谐修正量），而不必面对 GRAPE 式的成百上千自由度；这一”少参数波形 + 闭环微调”的工作方式与门标定闭环（ORBIT/RB/Ad-HOC）直接衔接。需要注意的边界是：光滑包络压的是频域泄漏，它不自动消除失谐漂移、光移与多普勒等通道误差，后者仍需鲁棒最优控制或硬件手段处理。

与现代少参数方案的衔接：NA057 的”初始形状 + 导数修正 + 少量系数”结构是各类低维 ansatz 的原型——时间最优门的相位 ansatz（量子速度极限）、ORMD 的 Bernstein 多项式（ORMD/ARMD）与 CRAB 的谐波展开（CRAB/DCRAB）都遵循”解析、光滑、少参数”的同一条工程路线。

方波/Gaussian/DRAG 三种波形的频谱对比与贝尔失真度随门时间曲线（图源：Theis 2016 Fig. 1）

参数与量级

量	典型值	来源
贝尔态保真度（整形脉冲，300 K）	$F_B>0.9999$，门时间 $\sim 50\mathrm{ns}$	Theis 2016
贝尔态保真度（4 K）	$F_B>0.99999$	Theis 2016
方波本征上限（300 K）	$F_B\sim 0.999$	Theis 2016
Gaussian/DRAG 泄漏压制	方波→Gaussian 约 2.5 个数量级；再加 DRAG 约 1.5 个数量级	Theis 2016
最优阻塞	$B_0/2\pi \approx 0.68\mathrm{GHz}$（Cs $n=141$）；$0.7-2.7\mathrm{GHz}$（$n=107$）	Theis 2016
里德堡态寿命	${\tau }_n=538/969\mu \mathrm{s}$（Cs $n=107/141$，300 K）	Theis 2016
失谐 DRAG 修正	$\Lambda /2\pi =124.07\mathrm{MHz}$（${\tau }_t=25\mathrm{ns}$，$\propto 1/{\tau }_t^2$）	Theis 2016
压制通道数	$\pi$ 脉冲 $N=4$（两条通道），$2\pi$ 脉冲 $N=2$（一条通道）	Theis 2016

与其他概念的关系

• 导数展开的底层机制（正交分量补偿、绝热框架）在DRAG 与泄漏抑制中专门展开；两者的差别是超导三能级梯子 vs 里德堡多泄漏通道。
• 解析形状换取的频域压制，与方波门的泄漏上限共同进入阻塞 CZ 门的误差预算；全口径分解见误差预算与真实建模。
• “少参数解析 ansatz”在时间最优门（量子速度极限）、调制驱动（ORMD/ARMD）与随机基优化（CRAB/DCRAB）中反复出现；带宽受限优化的系统化处理见鲁棒最优控制。
• 绝热极限下的解析脉冲（暗态跟踪）见ARP 门与STIRAP 与暗态门；微波域的 Blackman 整形单比特门见中性原子单比特门。
• 压制的对象是里德堡相互作用能级图中的邻近 $n$ 态：能级选择（$n$、$B_0$）与波形必须联合优化。

参考文献

• 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
• M. Saffman, “Symmetric Rydberg controlled-Z gates with adiabatic pulses”, Phys. Rev. A 101, 062309 (2020). Saffman 2020 · DOI 10.1103/PhysRevA.101.062309
• Y. Wang, A. Kumar, T.-Y. Wu, D. S. Weiss, “Single-qubit gates based on targeted phase shifts in a 3D neutral atom array”, Science (2016). NA058 · DOI 10.1126/science.aaf2581
• 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
• H. Silvério, S. Grijalva, C. Dalyac, L. Leclerc, et al., “Pulser: An open-source package for the design of pulse sequences in programmable neutral-atom arrays”, Quantum 6, 629 (2022). NA177 · DOI 10.22331/q-2022-01-24-629