⚛️ quantum physics

Lie algebra-assisted quantum simulation and quantum optimal control via high-order Magnus expansions

本文介绍了一种可扩展的、由李代数辅助的方法，该方法能够高效地计算随时间变化的量子哈密顿量的高阶马格努斯展开，通过在中性原子平台上设计五比特相位门脉冲所展示的结果，显著加速了量子模拟和最优控制的应用。

原作者： R. F. dos Santos, S. J. J. M. F. Kokkelmans

发布于 2026-09-09

📖 1 分钟阅读🧠 深度阅读

原作者： R. F. dos Santos, S. J. J. M. F. Kokkelmans

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

在量子物理领域，物质的行为受控于与我们日常经验迥然不同的规则。这一领域的核心在于量子系统的概念，即由原子或电子等粒子组成的集合，这些粒子可以同时存在于多种状态之中。为了使这些系统能够应用于量子计算机等技术，科学家必须极其精确地引导其演化，通常使用激光或磁场等外部力量。这个过程被称为量子控制。然而，预测这些系统如何随时间变化的数学方法是出了名的困难。不同于在一个可预测弧度内摆动的简单单摆，一个受变化力驱动的量子系统涉及到一个复杂的相互作用网络，其中事件的顺序至关重要。如果顺序改变，结果也会改变。这种非对易性意味着，在尝试模拟即使是中等复杂度的系统时，标准的计算方法往往会失效或变得极其缓慢，导致研究人员缺乏关于其设备将如何表现的清晰蓝图。

埃因霍芬理工大学的一个研究小组开发了一种导航这一数学迷宫的新方法。他们创造了一种方法，能够快速且准确地模拟受时变控制驱动的量子系统，而这项任务此前需要巨大的计算能力。通过专注于一种被称为李代数（Lie algebra）的特定数学结构——它描述了量子系统中不同部分之间的关系——该团队找到了一种简化问题的方法。他们不再纠结于随着每一步而呈指数级变得愈发复杂的复杂高阶积分，而是将整个计算转化为一个多项式表达式。在这个新的框架下，量子系统的演化由一组仅取决于控制信号形状和实验持续时间的系数来描述。这种转变将计算量降低到了仅与控制函数的复杂度相关、而非与量子系统规模相关的水平，使得在微秒内计算出高阶近似成为可能。

研究人员通过将其应用于两个不同的挑战，展示了这一方法的威力。首先，他们测试了该方法在模拟从两个到九个量子比特（量子信息的基石）规模的系统动力学方面的效率。在这些测试中，新方法在不到一毫秒的时间内就评估了高达十二阶近似的系统演化。这一速度比以往的技术快了四个数量级，而以往的技术在超过四阶后便难以维持。结果表明，只要理解了底层的数学结构，该方法无论系统规模如何都能保持高度的准确性。这种能力对于量子模拟至关重要，因为在量子模拟中，科学家需要预测材料或复杂分子在特定条件下会如何表现，而无需进行实际的物理实验，因为后者可能无法实现或成本过高。

除了模拟之外，该团队还将他们的方法应用于中性原子量子计算机的控制脉冲设计，这是一种利用光束固定原子的平台。在这些系统中，原子根据彼此之间的距离进行相互作用，从而为生成量子门（量子计算机的逻辑运算）创造了丰富的环境。设计这些量子门通常涉及寻找操纵原子的完美激光脉冲序列。研究人员使用他们的新型多项式方法来优化一个五量子比特门的脉冲，这是一个需要同时精确协调多个原子的复杂操作。他们将控制脉冲视为一条平滑的连续曲线，而非一系列离散的步骤，从而允许他们同时调整脉冲的形状及其总持续时间。这种灵活性促成了具有高保真度运行的门的设计，证明了该方法能够处理现实世界量子硬件的复杂需求。

这项工作的意义在于它能够将一个计算上难以承受的问题转化为一个可控的问题。通过将量子动力学的困难积分转换为直观的多项式，研究人员提供了一个既快速又具有解析可微性的工具。这意味着该方法不仅可以预测结果，还可以引导优化过程本身，告诉研究人员如何通过微调控制信号来提高性能。虽然目前的结果是基于数值模拟和特定模型的，但该方法的扩展性表明它可能成为设计未来量子设备的标准工具。能够如此快速且精确地模拟和控制复杂的量子系统，消除了量子技术发展中的一个主要瓶颈，使强大的量子计算机和传感器的愿景离现实更近了一步。

技术摘要：基于李代数的量子模拟与高阶马格努斯展开的量子最优控制

问题陈述
在随时间变化的驱动下，量子系统的演化对于量子计算、量子传感和量子通信等新兴技术至关重要。然而，量子动力学的高维性和非交换特性使得模拟和优化这些系统变得极具挑战性。虽然马格努斯展开（Magnus expansion, ME）提供了一个通过常数“有效哈密顿量”来近似随时间变化系统有效动力学的解析框架，但计算高阶项在计算上是非常昂贵的。现有方法依赖于对对易子进行高阶时间积分，这限制了其实际应用仅限于低阶（通常为四阶），或者需要复杂的解析表达式，且扩展性较差。这一瓶颈阻碍了高效量子操作的设计，特别是在噪声中规模量子（NISQ）设备上的多比特门设计。

方法论
作者提出了一种可扩展的方法来评估马格努斯展开，该方法将计算量降低到仅取决于随时间变化的控制函数的自由度，而非整个系统的规模。该方法针对形式为 $H(t) = A + d(t)B$ 的哈密顿量，其中 $A$ 是常数漂移项， $B$ 是由标量函数 $d(t)$ 調制的受控项。

该方法的核心包含三个关键步骤：

控制函数的多项式表示：将控制函数 $d(t)$ 表示为最高次数为 $m$ 的多项式。这使得马格努斯展开中涉及的高阶时间积分可以通过解析法求解。
利用动力学李代数：有效哈密顿量被证明属于由 $A$ 和 $B$ 生成的动力学李代数。通过在这一李代数的基底中表示哈密顿量，作者预先计算了展开所需的对易关系（结构常数）。
多项式系数推导：高阶积分和对易子被系统地简化为多项式表达式。马格努斯展开表示为 $M^{(k_M)}(t, \vec{d}) = \sum_\mu a_\mu^{(k_M)}(t, \vec{d}) L_\mu$ ，其中 $L_\mu$ 是李代数的基元素，而 $a_\mu$ 是取决于时间 $t$ 和控制参数 $\vec{d}$ 的多项式系数。这些系数封装了复杂的动力学过程，并针对给定的控制模型进行一次性计算。

核心贡献

可扩展的高阶评估：该方法将计算马格努斯展开的复杂度降低到对控制参数的多项式求值。这使得评估高达十二阶的展开可以在不到一毫秒内完成，比以往的技术（例如依赖于链式积分的技术）快了几个数量级。
解析可微性：生成的多项式表达式对于控制参数是解析可微的。这一特性使得可以直接应用基于梯度的优化算法进行量子最优控制。
GrAPE 的扩展：作者利用截断的马格努斯展开扩展了梯度上升脉冲工程（GrAPE）算法以处理连续脉冲。与通常使用常数阶跃函数的标准 GrAPE 不同，该方法优化的是连续脉冲形状（以 Hermite 插值样条表示），并同时优化脉冲持续时间。

结果

数值效率：在稀疏（最近邻伊辛模型）和稠密（长程伊辛模型）模型上的基准测试表明，只要截断阶数 $k_M$ 固定，评估马格努斯展开多项式的时间随系统规模增加几乎保持不变（范围在微秒到 $\sim 100$ $\mu$ s 之间）。动力学系数的生成随量子比特数量呈弱扩展。
误差缩放：该方法保留了马格努斯展开的时间对称性质。数值测试确认，只要时间截断阶 $\Gamma$ 足够大，误差将按 $O(t^{k_M+3})$ 进行缩放（对于偶数阶截断 $k_M$ ），这与理论预期一致。
最优控制应用：该方法被用于为中性原子平台上的里德堡原子设计五比特相位门族 ( $C_4P(\phi)$ $C_{4} P (ϕ)$ )。该系统使用 $n=k+1$ $n = k + 1$ 个原子，其中 $k$ $k$ 是控制量子比特的数量；因此， $C_4P$ $C_{4} P$ 门涉及 5 个原子，这与摘要中描述的“五比特相位门”一致。研究还包括一个三比特受控相位门 ( $C_2Z$ $C_{2} Z$ )。
- 对于三比特受控相位门 ( $C_2Z$ )，该方法成功优化了一个持续时间为 $\Omega T = 11.5$ 的脉冲，这比理论上的时间最优脉冲慢 6%，归因于对激光失谐而非相位的优化。
- 对于五比特门族 ( $C_4P(\phi)$ )，作者采用了“热启动”策略，即使用针对较低相位角的优化脉冲作为较高相位角的初始猜测。该方法成功设计了整个门族直至 $C_4Z$ 的脉冲，其持续时间与最先进的结果相当，同时直接将脉冲持续时间作为一个参数进行优化。

意义与主张
论文声称，该方法提供了一种可以比现有技术显著更快进行评估的马格努斯展开多项式表达式，从而能够广泛应用于量子模拟和量子最优控制。通过将计算负担转移到控制自由度，该方法实现了此前难以实现的的高阶模拟。作者强调，该表达式的解析可微性促进了脉冲设计的新方法，特别是展示了在高精度下设计中性原子平台多比特控制门的能力。这项工作将该方法定位为一种既能快速模拟随时间变化的量子系统，又能作为研究其动力学（特别是在 Floquet 系统和量子材料背景下）的解析工具。

类似论文