通过薛定谔化进行非线性偏微分方程的哈密顿模拟

这项研究为何重要

我们在自然中看到的许多图案——从培养皿中化学物质的反应和扩散到细胞在组织中的扩散——都由既复杂又非线性的方程描述。随着系统规模增长，在现有经典计算机上精确模拟这些系统的代价迅速变得难以承受。本研究提出了一种利用量子计算机处理这类重要方程的方法，为分析远超经典能力的巨型真实系统打开了可能的途径。

从混乱的现实到可求解的方程

科学家和工程师常用偏微分方程来描述诸如温度、浓度或形变等量在空间和时间上的变化。当这些方程是线性时，已有量子算法可以非常高效地模拟它们的行为。然而，许多真实现象——如湍流、大尺度材料变形和反应-扩散过程——由非线性方程控制，其中变化的规律依赖于当前状态本身。这些非线性恰恰使得系统表现丰富，但也让经典与量子方法都难以求解。

将非线性问题转化为线性问题

作者集中研究一种常用于材料与生物模式建模的特定非线性方程——反应-扩散方程。他们的第一步是应用一种称为 Carleman 线性化的数学技术。概念上，该方法通过跟踪基本变量以及它们到某一阶的所有乘积，将原始非线性系统替换为一个更大的线性系统。实际上，这个无限层次在可控水平被截断，产生一个规模很大但纯粹线性的系统，仍然近似原始的非线性动力学。这一步使问题更契合量子硬件的特性，后者天然处理线性演化。

让耗散动力学看起来像量子动力学

即便在线性化之后，得到的方程通常描述的是耗散系统，其中量可能不可逆地衰减或扩散。相比之下，量子演化是守恒的：它保全整体概率，并在数学上由酉算子表示。为弥合这一差距，作者使用了一种基于扭曲相位变换的薛定谔化方法。他们引入一个额外的人工变量，并重新表述线性系统，使其演化可以写成与支配量子力学的薛定谔方程相同的数学形式。在这个扩展的空间中，时间演化是酉的，原则上可以直接作为量子哈密顿模拟来实现。

在模型问题上测试该方法

为了评估他们的方法，研究人员将组合程序——Carleman 线性化加上薛定谔化（简称 CLS）——应用到一个经典的非线性反应-扩散模型，即 KPP–Fisher 方程。他们对空间和辅助变量进行离散化，然后使用标准数值技术模拟时间演化，并将基于 CLS 的结果与更传统的差分方法进行比较。模拟波形的形状和运动在各方法间高度一致，详细的误差分析表明精度如何依赖于诸如线性化截断阶数和空间网格精细程度等选择。研究发现误差按可控且可预测的方式表现，主要由标准数值近似支配，而非 CLS 变换本身的根本缺陷。

这对未来量子模拟意味着什么

简言之，这项工作表明，一大类非线性方程可以系统地重新表述为量子计算机设计来处理的形式。虽然本研究使用经典模拟来验证方案，但相同步骤可以被转化为实现相应哈密顿演化的量子电路。如果未来的量子设备能在规模上实现这些电路，CLS 有望实现对极其复杂系统的高效模拟——例如大规模化学反应网络或复杂的相分离过程——在这些情况下经典方法变得不可行。主要结论是，物理世界中的非线性行为并不必然排斥基于强大哈密顿方法的量子算法；相反，通过合适的数学桥梁，它可以被引入量子领域。

引用: Sasaki, S., Endo, K. & Muramatsu, M. Hamiltonian simulation for nonlinear partial differential equation by Schrödingerization. Sci Rep 16, 11743 (2026). https://doi.org/10.1038/s41598-026-44920-8

关键词: 量子计算, 哈密顿模拟, 非线性偏微分方程, 反应-扩散, Carleman 线性化