Pulser 把中性原子量子处理器抽象成"Device 约束 + Register 坐标 + Channel 通道 + Pulse 波形"四层对象:用户不再写抽象门电路,而是直接编排激光的振幅、失谐与相位随时间的轨迹——这是模拟与数字两种计算范式在中性原子平台上唯一的公共编程接口。

问题背景

可编程里德堡阵列具有两种截然不同的使用方式:模拟方式下,全局里德堡驱动直接实现 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:

  1. Device:硬件规格清单——里德堡能级选择、激光振幅/频率允许范围、原子间距上下限、可声明的通道集合;所有约束在构造序列时强制执行;2. Register:原子的坐标与 ID(量子比特映射),决定相互作用图;
  2. Channel:激光的抽象,按寻址(局域/全局)与跃迁类型(Rydberg/Raman)组织;
  3. Waveform:脉冲的振幅/失谐时间轨迹,可用预置形状(ramp、Blackman)或自定义函数;
  4. Pulse:振幅波形 + 失谐波形 + 固定相位,被加入序列时指明目标原子与执行通道,另有 delay(通道空转)与 target(把局域通道指向特定比特)等指令;
  5. Sequence:各通道上脉冲的调度表,绑定 Register 与 Device,是可发送给 QPU 的完整描述;
  6. Simulation:基于 QuTiP 的本地模拟入口,返回带后处理方法的结果对象。

读出同样在框架内显式建模:测量作用于 Register 中的全体原子,可在指定基(计算基或里德堡基)下进行,序列结束时声明的测量基决定返回的位串分布 Silvério 2022。

设计流程示例。 论文以 Bell 态 制备为例走通全流程 Silvério 2022:

  1. 电路改写为中性原子原生形式:CNOT 分解为 CZ + 单比特旋转(CNOT 非原生门),Hadamard 换成共振脉冲等效旋转,减少自由参数、降低出错面;
  2. 创建 Register:定义原子坐标与 ID;
  3. 加载目标 Device(示例用 Pasqal 研究设备 Chadoq2),以 Sequence(reg, Chadoq2) 绑定硬件约束;
  4. 声明通道:rydberg_global(全局里德堡)与 raman_local(局域 Raman);
  5. 用 Blackman 波形构造 CZ 脉冲,振幅上限由阻塞约束反推:max_val = Chadoq2.rabi_from_blockade(8);
  6. 末尾声明测量基并执行:发送真机,或调用 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)

与其他概念的关系

参考文献