问题背景
可编程里德堡阵列具有两种截然不同的使用方式:模拟方式下,全局里德堡驱动直接实现 Ising 型哈密顿量的时间演化,用于量子多体模拟与优化问题;数字方式下,激光脉冲序列组合出量子门以执行线路。云平台通常只暴露”门电路”这一高层 API,把脉冲细节封死在设备专属编译之后;但对量子控制研究(波形整形、最优控制、鲁棒性设计)与混合模拟-数字算法的开发者而言,脉冲层的可编程性才是刚需 Silvério 2022。Pulser 填补的正是这一层:一个描述”物理层面发生了什么”的开源 Python 框架,同一份脉冲序列既可以发给真实 QPU,也可以在本地小系统上模拟。
理论模型:驱动、相互作用与两种范式
双能级驱动。 Pulser 对每个通道建模为对目标跃迁 的经典驱动,含三个时间依赖的控制量:Rabi 频率 (振幅)、失谐 与相位 ,在旋转波近似下单原子哈密顿量为
在 Bloch 球图像中,该哈密顿量每时每刻描述绕轴 的旋转,角速度为有效 Rabi 频率 Silvério 2022。共振脉冲(、固定 )实现绕 x-y 平面某轴的旋转 ;由于绕 z 轴旋转等效于相位参考平移(虚拟 Z),任意单比特旋转可由 型构造完成,无需真实移动激光频率 Silvério 2022。
模拟范式的相互作用。 在基态–里德堡”Ising”构型下(Pulser 的 ground-rydberg 通道),范德瓦尔斯相互作用(见里德堡相互作用)把多体哈密顿量补全为含横场的量子 Ising 形式
其中 为原子间距、耦合 ——阵列坐标本身成为可编程参数,这正是”场可编程”的含义(见FPQA 计算模型)。XY 构型( 的交换型相互作用)在论文发表时尚未实现。
数字范式的门。 双比特门由里德堡阻塞(见里德堡阻塞)介导:对控制-靶对的脉冲序列使 之外的基矢获得条件相位,即 CZ 门;配合单比特旋转即可通用。数字基用基态超精细态 (Raman 通道寻址),模拟基直接用 ;当通道同时涉及两个跃迁时, 原子系统占据 维希尔伯特空间 Silvério 2022。
方法与实现
代码架构:七类对象。 Pulser 的核心组件及关系为 Silvério 2022:
- Device:硬件规格清单——里德堡能级选择、激光振幅/频率允许范围、原子间距上下限、可声明的通道集合;所有约束在构造序列时强制执行;2. Register:原子的坐标与 ID(量子比特映射),决定相互作用图;
- Channel:激光的抽象,按寻址(局域/全局)与跃迁类型(Rydberg/Raman)组织;
- Waveform:脉冲的振幅/失谐时间轨迹,可用预置形状(ramp、Blackman)或自定义函数;
- Pulse:振幅波形 + 失谐波形 + 固定相位,被加入序列时指明目标原子与执行通道,另有
delay(通道空转)与target(把局域通道指向特定比特)等指令; - Sequence:各通道上脉冲的调度表,绑定 Register 与 Device,是可发送给 QPU 的完整描述;
- Simulation:基于 QuTiP 的本地模拟入口,返回带后处理方法的结果对象。
读出同样在框架内显式建模:测量作用于 Register 中的全体原子,可在指定基(计算基或里德堡基)下进行,序列结束时声明的测量基决定返回的位串分布 Silvério 2022。
设计流程示例。 论文以 Bell 态 制备为例走通全流程 Silvério 2022:
- 电路改写为中性原子原生形式:CNOT 分解为 CZ + 单比特旋转(CNOT 非原生门),Hadamard 换成共振脉冲等效旋转,减少自由参数、降低出错面;
- 创建 Register:定义原子坐标与 ID;
- 加载目标 Device(示例用 Pasqal 研究设备 Chadoq2),以
Sequence(reg, Chadoq2)绑定硬件约束; - 声明通道:
rydberg_global(全局里德堡)与raman_local(局域 Raman); - 用 Blackman 波形构造 CZ 脉冲,振幅上限由阻塞约束反推:
max_val = Chadoq2.rabi_from_blockade(8); - 末尾声明测量基并执行:发送真机,或调用 Simulation 在本地采样。
CZ 脉冲序列对各计算基矢累积的相位在论文 Table 1 中逐行列出,配合相位与局域幺正的重定义即可识别出受控 Z(这一”相位–局域等价”的分析思路见门相位与局域等价)。其中里德堡阻塞条件为 ,这也是 Device 对象用”最小原子间距”约束脉冲合法性的物理依据。
数字与模拟的公共接口。 同一 Sequence 对象里,全局 Rydberg 通道跑模拟段(Ising 演化、绝热斜坡),局域 Raman 通道跑数字段(单/双比特门),这种”混合数字-模拟算法”的表述能力是 Pulser 相对门电路语言(如 OpenQASM)的根本差异 Silvério 2022。
模拟能力边界。 本地模拟用 QuTiP 精确演化, 维空间下个人计算机可高效模拟的系统约为 个原子;更大规模需真机执行或张量网络等专门方法 Silvério 2022。
两个官方用例。 论文第 4 节用内置模拟器复现了实验结果,展示模拟/数字两种范式的典型用法 Silvério 2022:
- 绝热制备反铁磁态:原子排成间距 的方阵,全局 Rydberg 通道上施加”上升–扫描–下降”三段脉冲,对应量子 Ising 相图中从顺磁相进入反铁磁穹顶的一条路径(示例参数 、 从 扫到 ,);结果用自旋-自旋关联函数 定量评估。这类 Ising 态在中性原子实验中已做到数百比特,超出经典模拟的精度范围;
- QAOA 变分算法:把图问题(单位圆盘图)映射为 Register 几何,构建参数化脉冲序列并迭代优化角度——脉冲层接口让变分参数直接作用于波形,这是门标定闭环与变分优化共享的基础设施。
采样接口形如 res.sample_final_state(N_samples=1e4),返回有限采样噪声下的位串计数分布;seq.draw() 可视化序列(通道目标、测量基等标注其上)。
生态位置。 Pulser 面向的用户群横跨量子控制(波形设计,如解析整形脉冲与ORMD/ARMD 调制驱动的原型载体)、模拟物理(Ising 模型、MIS 优化)与数字线路(脉冲级编译目标);其”Device 规格约束脉冲合法性”的设计让编译器(见编译器公平评测与可视化)与控制研究者共享同一硬件模型。
参数与量级
| 量 | 典型值 | 来源 |
|---|---|---|
| 编程层级 | 脉冲级(振幅/失谐/相位波形) | Silvério 2022 |
| 通道组织 | 局域/全局寻址 × Rydberg/Raman 跃迁 | Silvério 2022 |
| 波形库 | 自定义、ramp、Blackman 等 | Silvério 2022 |
| 本地模拟规模 | 原子(QuTiP 精确演化) | Silvério 2022 |
| 希尔伯特空间维度 | 数字+模拟混合时 | Silvério 2022 |
| 示例目标设备 | Chadoq2(Pasqal 研究 QPU 规格) | Silvério 2022 |
| 分发方式 | 开源 Python 包(pip install pulser) | Silvério 2022 |

Pulser 主要类关系图:Sequence 居中,连接 Device/Channel/Waveform/Pulse/Register/Simulation(图源:Silvério 2022 Fig. 5)
与其他概念的关系
- 全局 Rydberg 通道的模拟范式即里德堡模拟优化与 MIS的编程入口;相互作用的物理来源见里德堡相互作用与 Förster 共振。
- 数字范式的核心资源是阻塞 CZ 门与多比特阻塞门;脉冲形状如何抑制阻塞误差见解析整形脉冲与ORMD/ARMD,最优控制的框架化描述见量子最优控制框架。
- Device 对象编码的硬件约束(间距、振幅范围、寻址能力)是FPQA 计算模型的形式化来源之一;标定这些参数的闭环流程见门标定闭环。
- 脉冲合法性受激光系统与光镊物理约束,背景见光镊阵列与魔强度陷阱与中性原子量子计算硬件全景。
参考文献
- 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 · arXiv:2104.15044