⚛️ quantum physics

High-order splitting of non-unitary operators on quantum computers

该论文提出了一种利用复系数乘积公式将耗散动力学分解为实部和虚部时间演化的方法，成功在含噪量子处理器上实现了高精度的非幺正算子模拟，证明了高阶算子分裂方案在模拟耗散系统方面的准确性与实用性。

原作者： Peter Brearley, Philipp Pfeffer

发布于 2026-02-17

📖 1 分钟阅读🧠 深度阅读

原作者： Peter Brearley, Philipp Pfeffer

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

这篇论文讲述了一个关于如何在量子计算机上模拟“有损耗”的物理过程的突破性方法。

为了让你轻松理解，我们可以把这篇论文的核心内容想象成**“在量子世界里烹饪一道复杂的炖菜”**。

1. 背景：量子计算机的“洁癖”

想象一下，量子计算机是一个极其挑剔的厨师。它天生擅长做**“可逆”的菜肴（比如把面团揉成球，再完美地变回面团），这在物理上叫“幺正演化”**（Unitary Dynamics）。

但在现实世界中，绝大多数事情都是**“不可逆”的，也就是有“损耗”**的：

摩擦会让物体停下来（能量损耗）。
墨水在水中扩散（熵增）。
声波在空气中传播时会变弱（阻尼）。

这些过程在数学上叫**“耗散动力学”**。以前的量子算法很难处理这些“有损耗”的菜，因为一旦尝试模拟“反向”的损耗过程（比如让扩散的墨水自动聚拢），就像试图让打碎的鸡蛋自动复原一样，在数学上会导致数值爆炸，完全无法计算。

2. 问题：高阶方法的“副作用”

为了更精确地模拟物理过程，科学家们通常使用**“高阶分裂法”**。

低阶方法（如 1 阶、2 阶）：就像切菜，切得比较粗糙，但很稳。
高阶方法（如 4 阶、6 阶）：就像把菜切得极碎，能做出更美味的料理，精度极高。

但是，高阶方法有个致命缺点： 为了达到高精度，数学公式里必须包含**“负数时间步”**。

在模拟“可逆”过程时，负时间步只是“倒着走”，没问题。
在模拟“有损耗”过程时，负时间步意味着“让损耗变成增益”（比如让变冷的物体自动变热且更热），这会导致系统瞬间崩溃（数值不稳定）。

所以，以前大家不敢在量子计算机上用高阶方法模拟耗散过程，只能退而求其次，用精度很低的 2 阶方法。

3. 解决方案：给时间加点“魔法调料”

这篇论文的作者（Peter Brearley 和 Philipp Pfeffer）想出了一个绝妙的办法：使用“复数系数”的时间步。

这就好比我们在烹饪时，不再只用“正时间”（向前煮）或“负时间”（倒着煮），而是引入了**“虚时间”**（Imaginary Time）。

实部（Real Part）：代表正常的物理时间流逝，负责让能量耗散（比如让波慢慢变弱）。作者确保这部分系数是正的，保证过程是稳定的。
虚部（Imaginary Part）：代表一种“魔法旋转”。在量子力学里，虚时间演化其实等价于一种特殊的旋转（幺正演化）。

核心技巧：
作者设计了一套复杂的“配方”（乘积公式），把整个演化过程拆解成一系列简单的步骤：

实时间步骤：让波在真实世界里传播（幺正演化，量子计算机最擅长的）。
虚时间步骤：让波在“魔法维度”里旋转（也是幺正演化，量子计算机也能做）。
耗散步骤：利用复数系数的实部，让波在“收缩”（模拟能量损失）。

通过巧妙组合这些步骤，他们成功地把一个“会崩溃”的高阶耗散模拟，变成了一连串量子计算机能轻松处理的“简单旋转和收缩”操作。

4. 实验验证：在真实的量子芯片上“试菜”

为了证明这个方法真的有用，作者在IonQ Forte 1（一种真实的离子阱量子计算机）上进行了实验。

模拟对象：模拟**“阻尼波”**（就像在充满空气的房间里敲击音叉，声音会逐渐变小）。
对比测试：他们分别用了 1 阶、2 阶、4 阶和 6 阶的方法。
结果惊人：
- 1 阶和 2 阶：虽然电路简单，但模拟出来的波形很不准，就像用粗勺子舀汤，漏了很多细节。
- 4 阶：虽然电路变复杂了（步骤更多），但在嘈杂的量子硬件上，它模拟出的波形最接近真实情况，甚至比 2 阶更准！
- 6 阶：虽然理论上精度最高，但因为电路太深（步骤太多），量子计算机的噪声（就像厨房里的杂音）把结果给淹没了，反而不如 4 阶。

结论：在目前的“含噪声”量子计算机上，4 阶方法是性价比最高的选择。它既利用了高阶方法的精度优势，又没有被噪声彻底击垮。

5. 总结与意义

这篇论文就像是为量子计算机颁发了一张**“高阶烹饪许可证”**。

以前：因为怕“负时间”导致崩溃，我们只能做简单的低精度模拟。
现在：通过引入“复数时间”和巧妙的拆分，我们可以用高阶方法稳定地模拟摩擦、扩散、阻尼等现实世界中最常见的物理过程。

这意味着，未来量子计算机不仅能模拟完美的、理想化的物理世界，还能更真实、更准确地模拟充满摩擦和损耗的现实世界。这对于流体力学、材料科学、甚至金融建模等领域，都是一次巨大的飞跃。

一句话总结：
作者发明了一种“魔法配方”，让量子计算机能够用高精度的方法，稳定地模拟那些会“漏气”、“变慢”或“消失”的物理现象，并且已经在真实的量子芯片上证明了它的威力。

这是一份关于论文《High-order splitting of non-unitary operators on quantum computers》（量子计算机上非幺正算子的高阶分裂）的详细技术总结。

1. 研究背景与问题 (Problem)

核心挑战：量子计算机天然适合模拟幺正（Unitary）动力学，因为量子门模型本质上是可逆的。然而，大多数物理过程（如摩擦、粘滞、扩散、边界通量）以及开放量子系统的演化都是耗散（Dissipative）和不可逆的，表现为非幺正动力学。
现有局限：
- 传统的算子分裂（Operator Splitting）方法（如 Lie-Trotter 和 Strang 分裂）在模拟幺正动力学时非常有效。
- 为了获得高阶精度（高于 2 阶），传统的乘积公式（Product Formulas）必须引入负系数。
- 对于幺正动力学，负系数意味着模拟“时间反向演化”，这在物理上是可行的（因为 $e^{-iHt}$ 也是幺正的）。
- 对于耗散动力学，负系数意味着模拟“时间反向的耗散”，这在数值上是不稳定的，且无法在不重新缩放的情况下直接编码为幺正算子。
现状：目前量子模拟耗散动力学的方法主要局限于 2 阶精度，缺乏高阶、稳定且实用的量子算法。

2. 方法论 (Methodology)

作者提出了一种基于复系数乘积公式的高阶算子分裂方法，用于在量子计算机上稳定地模拟非幺正动力学。

算子分解：
- 将演化生成元 $M$ 分解为厄米部分（ $H_1$ ）和反厄米部分（ $iH_2$ ），即 $M = H_1 + iH_2$ 。
- $H_1$ 对应耗散部分（实部）， $iH_2$ 对应幺正部分（虚部）。
复系数乘积公式：
- 利用满足特定约束的复系数 $a_j$ 和实系数 $b_j$ 构建乘积公式。
- 约束条件：耗散子步的系数 $a_j$ 必须满足 $\Re(a_j) > 0$ （实部为正），以确保耗散演化的收缩性（数值稳定）；幺正子步的系数 $b_j$ 必须满足 $b_j > 0$ 。
- 演化机制：
  - 实时间演化：利用 $H_1$ 的实部进行耗散演化，利用 $H_2$ 进行幺正演化。
  - 虚时间演化：利用 $H_1$ 的虚部（ $\Im(a_j)$ ）进行演化。由于系数是复数，虚部对应的演化 $e^{i H_1 \Im(a_j) \Delta t}$ 本质上是幺正的（表现为相位门）。
- 通过这种方式，整个非幺正动力学被分解为一系列简单的哈密顿量演化（实时间和虚时间交替），所有子步均可在量子计算机上高效实现。
高阶公式选择：
- 采用了 Castella 等人提出的4 阶公式（5 个阶段）和 Bernier 等人提出的6 阶公式（16 个阶段）。
- 这些公式均满足 $\Re(a_j) > 0$ 和 $b_j > 0$ 的约束，确保了算法的稳定性。

3. 关键贡献 (Key Contributions)

理论突破：首次展示了利用复系数乘积公式在量子计算机上实现耗散动力学的高阶（4 阶和 6 阶）算子分裂。解决了高阶分裂中负系数导致的数值不稳定问题。
算法设计：提出了一种将耗散动力学分解为实时间（耗散）和虚时间（幺正/相位）演化的具体方案。该方法直接生成了耗散演化算子的未缩放块编码（Block Encoding）。
电路实现：推导了针对经典阻尼波方程（Damped Wave Equation）的量子电路，包括：
- 基于傅里叶谱方法的离散化。
- 利用受控旋转门（Controlled $R_Y$ ）实现波传播。
- 利用辅助量子比特和后选择（Post-selection）或振幅放大实现阻尼项。
硬件验证：在 IonQ Forte 1 离子阱量子处理器上进行了端到端的实验验证，涵盖了从 1 阶到 6 阶的不同精度方案。

4. 实验结果 (Results)

数值模拟（Statevector Simulations）：
- 在 9 量子比特系统上模拟了阻尼波方程。
- 结果显示，所有高阶公式（4 阶和 6 阶）均达到了理论预期的误差缩放比例（ $O(\Delta t^p)$ ）。
- 在相同的误差容限下，高阶公式所需的 CNOT 门数量显著少于低阶公式，特别是在高精度需求下效率提升明显。
硬件实验（IonQ Forte 1）：
- 配置：使用 16 个网格点（4 个数据量子比特），模拟一个分裂步。
- 精度对比：
  - 1 阶：结果明显不准确。
  - 2 阶和 4 阶：表现出高度准确性，在精度和电路深度之间取得了最佳平衡。4 阶方法在噪声环境下仍比 2 阶略优（误差范数更小）。
  - 6 阶：由于电路深度过大（232 个 CNOT 门），噪声效应占主导地位，导致精度反而低于 1 阶和 2 阶方案。
- 后选择概率：实验测量的后选择成功概率与理论值高度吻合，证明了算法在硬件上的可行性。
结论：尽管电路深度增加，但在含噪中等规模量子（NISQ）设备上，4 阶分裂方案能够提供比低阶方案更高的实际精度，证明了该方法在近中期量子硬件上的实用性。

5. 意义与展望 (Significance)

扩展量子模拟能力：将高阶算子分裂的灵活性和高精度从幺正动力学扩展到了耗散动力学，填补了量子计算在模拟物理耗散过程（如流体力学、热传导、化学反应）方面的空白。
资源效率：研究表明，对于给定的误差容限，高阶方法（如 4 阶）比低阶方法（如 2 阶）所需的资源（CNOT 门数）更少，或者在相同资源下精度更高。
NISQ 时代的适用性：实验证明了在当前的噪声硬件上，适度增加阶数（如从 2 阶提升到 4 阶）是有益的，但过高的阶数（如 6 阶）会因噪声累积而失效。这为未来算法设计提供了重要的工程指导。
未来工作：
- 寻找更高阶（>6 阶）且满足稳定性约束的系数。
- 探索更高效的辅助量子比特压缩技术（如压缩门 gadget）以减少量子比特需求。
- 将该方法推广到更广泛的工业应用，如湍流模拟和复杂材料科学。

总结：该论文提出并验证了一种利用复系数乘积公式在量子计算机上稳定模拟耗散动力学的高阶方法。通过在离子阱处理器上的成功演示，证明了 4 阶分裂方案在噪声环境下仍能超越低阶方案，为未来在量子计算机上解决广泛的非幺正物理问题奠定了坚实基础。