Fast Hamiltonian engineering from cut polytope geometry

量子计算机有望解决当今机器无法解决的问题，但它们是出了名的脆弱。为了使它们发挥作用，科学家必须通过编程让量子设备模拟特定系统的行为，例如分子或新材料，从而模拟复杂量子系统的行为。这一过程被称为哈密顿量工程（Hamiltonian engineering），它涉及将一个天然执行某些相互作用的原生机器，通过一系列控制脉冲来塑造其行为，从而模仿另一种所需的相互作用。挑战在于，这些控制脉冲的执行需要时间，而量子系统活跃的时间越长，就越容易将其脆弱的量子态丢失给噪声和误差。因此，最关键的目标是找到尽可能快的脉冲序列来实现所需的模拟，从而最大限度地减少机器暴露在环境中的时间。

一个研究团队开发了一种全新的统一方法，用于为包括量子比特（qubits）、高维量子比特（qudits）和费米子（fermions）在内的各种量子系统寻找这些最优脉冲序列。通过将问题视为高维空间中形状的几何谜题，他们创建了专门针对特定系统和目标模拟定制控制脉冲的算法。他们的方法始终能找到接近理论最优解的方案，性能显著优于以往依赖随机猜测的方法。在复杂模型的测试中，他们的技术将所需时间降低到了一个不随系统规模增长的水平，而旧方法则会随着系统的增大而变得更慢。

问题的核心在于量子系统是如何相互作用的。想象一下，一个量子设备天然允许粒子以某种特定方式进行相互作用，但科学家想要模拟的是另一种类型的相互作用。为了弥补这一差距，科学家施加了一层层的控制操作，即脉冲，这些脉冲会扭转系统的状态。目标是找到这些“扭转”的正确组合，使得系统有效地表现得就像是在遵循新的规则一样。研究人员意识到，对于许多重要的量子系统，天然相互作用与控制脉冲之间的关系遵循一个简单的规则：脉冲仅通过特定的相位偏移来改变相互作用，就像将旋钮转到一个特定的角度。这一观察使他们能够将寻找最佳脉冲的复杂任务转化为一个几何问题。

他们将这个问题可视化为一束光线从广袤的多维空间中的一个起点射出。目标相互作用定义了这条光线的方向。通过可用脉冲可以创造的所有可能相互作用构成了一个特定的几何形状，研究人员称之为“多胞体”（polytope）。最快的模拟对应于这条光线第一次接触到该形状表面的点。如果光线很快撞击到形状，则模拟速度快；如果必须行进很远，则模拟速度慢。研究人员证明，对于大型系统，要在合理的时间内完美解决找到这个精确点的数学问题是不可能的。这是计算机科学中一个已知的难题，这意味着对于任何实际应用，人们必须接受一个非常好的近似解而非完美答案。

为了克服这一点，该团队设计了一个巧妙的变通方法。他们没有试图直接撞击精确的形状，而是将问题放宽到一个包围原始形状的更平滑、更简单的形状上。然后，他们使用一种数学技术来弯曲搜索光线的路径，以便当他们最终从这个放宽后的形状中选取特定脉冲时，结果能恰好落在他们需要的位置。这个过程生成了他们所说的“信息驱动型”（informed）脉冲——这些控制序列不是随机选择的，而是根据系统的具体细节计算出来的。这些信息驱动型脉冲随后被输入到一个标准的优化工具中，以确定最终的时序和顺序。

研究人员在三种不同类型的量子系统上测试了该方法。首先，他们研究了作为目前大多数量子计算机基础的标准量子比特系统。他们将这种新方法与使用脉冲随机采样的现有技术进行了对比。在这些测试中，他们这种信息驱动的方法始终能找到比旧方法快得多的解决方案，通常能达到理论最佳时间的百分之几以内。相比之下，旧的随机方法往往需要更多的时间，特别是当目标模拟的复杂度增加时。

接下来，他们将该方法应用于量子比特（qudits），这是一种具有两个以上状态的量子单元，提供了更丰富的可能性。在这里，挑战更大，因为涉及的相互作用包含复数。他们的算法成功应对了这些复杂性，找到了近乎最优的脉冲序列。结果显示，他们的方法可以适应目标的特定结构，而随机方法即使在硬件允许更精细控制的情况下也无法改进。

最后，他们在费米子系统中测试了该方法，这些系统被用于模拟材料中的电子。这是一个特别困难的情况，因为其中的相互作用涉及不能占据相同状态的粒子。他们使用了一个被称为霍夫施塔特模型（Hofstadter model）的模型，该模型描述了在磁场中运动的电子。在这种场景下，他们的方法与旧的随机方法之间的差异是巨大的。随机方法所需的模拟时间随网格规模呈线性增长；随着网格变大，模拟时间也成比例增加。然而，他们的信息驱动方法找到的解决方案中，时间保持不变，无论网格规模如何。这意味着对于大规模材料模拟，他们的方法可以比以往快几个数量级，使以前难以实现的模拟突然变得可行。

这些发现的意义不仅在于速度。在量子模拟中，系统的运行时间直接关系到它积累了多少噪声。更快的模拟意味着更少的噪声和更准确的结果。此外，在模拟相互作用粒子的情况下，模拟的速度决定了可以构建的相互作用强度。更快的算法可以模拟比硬件天然支持的更强的相互作用，从而为研究新的物相打开了大门。研究人员还表明，他们的方法具有鲁棒性；即使控制脉冲不完美或执行需要有限的时间，算法也可以通过调整来抑制这些误差，而不损失其速度优势。

通过在单一几何框架下统一处理量子比特、量子比特（qudits）和费米子，这项工作为量子模拟器的自动编程提供了一个强大的新工具。它使该领域从试错法或随机猜测转向了系统化、有数学依据的方法。虽然由于完美解在计算上难以实现，该方法依赖于近似处理，但结果表明这些近似是非常紧密的。算法始终能提供近乎最优的性能，这表明可以模拟的理论极限比我们之前认为的要接近现实。这一进展让我们离使用量子设备解决化学和材料科学等现实世界问题的梦想又近了一步。