波动方程量子模拟

波动方程 ${\partial }_t^2\psi =c^2\Delta \psi$ 是量子计算机的”主场方程”：它的演化天然酉（能量守恒、模长不变），离散谱纯虚——不需要任何酉化技巧就能化成薛定谔方程。Liu–Lopes 在 2024 年给出了从对角化到真机的完整演示：QFT–对角相位–QFT 的三明治线路在 Quantinuum H1-1 离子阱上模拟了一维 Ricker 子波的传播Liu & Lopes 2024 (arXiv)。更早的 Costa–Love–Santos 量子行走方案（PRL 2019，本库未收录原文，以下数字转引自 QC090 的综述部分）把它推广到任意几何与混合边界。

从波动方程到薛定谔形式

一维波动方程 ${\partial }_t^2\psi =c^2{\partial }_x^2\psi$ 因式分解为左右行波（两个对流方程，见quantum-advection-equation）；等价地，引入辅助场

$$\varphi =i(-\Delta {)}^{-1/2}{\partial }_t\psi$$

则 $(\psi ,\varphi {)}^T$ 满足薛定谔方程

$$i{\partial }_t\left(\genfrac{}{}{0ex}{}{\psi }{\varphi }\right)=\left(\begin{array}{cc}0& \sqrt{-\Delta }\\ \sqrt{-\Delta }& 0\end{array}\right)\left(\genfrac{}{}{0ex}{}{\psi }{\varphi }\right)\equiv H_W\left(\genfrac{}{}{0ex}{}{\psi }{\varphi }\right).$$

波动哈密顿量 $H_W$ 厄米——这就是酉性的来源：第二时间导数方程化成两个耦合的一时间导数方程后，演化算符 $e^{-it{H}_W}$ 是酉的。物理直觉：$\varphi$ 就是”速度场”，$(\psi ,\varphi )$ 联合携带了守恒的能量内积。

QFT 对角化与三明治线路

周期边界 + 二阶中心差分下，离散 Laplacian 被量子傅里叶变换对角化：${\Delta }_a=\mathrm{QFT}D{\mathrm{QFT}}^{†}$（$D$ 为对角负特征值阵）。于是

$$e^{-it{H}_W}={\mathrm{QFT}}^{†}{e}^{-it{Z}_0\otimes \sqrt{-D}}\mathrm{QFT},$$

其中 $Z_0$ 是作用在 $\psi /\varphi$ 标签比特上的 Pauli-Z——中间的 $e^{-it{Z}_0\otimes \sqrt{-D}}$ 是对角算符，可拆成一串单/双比特旋转门（小角近似 $\sqrt{-D}$ 的本征值后，每个 Fourier 模式贡献相位 $t\sqrt{{\lambda }_k}$）。整个线路：变分砖墙拟设制备初态（深度 $\lfloor {\log }_2(n+1)\rfloor +1$）→ QFT → 对角演化 → 逆 QFT → 测量读出，$n+1$ 比特编码 $2^{n+1}$ 个格点Liu & Lopes 2024 (arXiv)。

图:近似模拟一维波动方程的量子线路(n+1=4 比特)——变分初态制备、QFT†-对角演化 exp(-itZ₀⊗√-D)-QFT、末态测量三段式(来源:QC090)

误差标度与真机结果

该方案的误差账目干净利落Liu & Lopes 2024 (arXiv)：

• 近似保真度损失随格点数 $ϵ\sim O(1/N^4)$（对角化+小角近似的联合标度）、随时间 $O(t^2)$；
• 双比特门总数 $O(n^2)$（编译到 H1-1 门集后与二次多项式拟合吻合）——多项式深度，NISQ 可承受。

真机演示：Quantinuum H1-1 离子阱上模拟 Ricker 子波（$n=6$ 比特、10000 shots），各时刻的波场采样与精确解吻合（误差棒为 Monte Carlo 采样误差 $\sigma /\sqrt{N_{\mathrm{shots}}}$）；无噪声仿真中保真度随比特数指数提升，含退极化噪声时出现”近似精度增益 vs 门噪声损失”的最优比特数。

图:Quantinuum H1-1 真机结果——不同时刻的 Ricker 子波波场,精确解(蓝实线)与 10000 shots 采样(橙虚线,含 ±ε_MC 误差棒)吻合(来源:QC090)

量子行走路线与联系

处理复杂几何（Dirichlet/Neumann 混合边界）的代表方案是 Costa–Love–Santos（PRL 2019）：把波动方程的时空几何编码成图，用量子行走（quantum walk）模拟图 Laplacian 的演化——据 QC090 综述部分转述，其空间复杂度 $O(d\log (\bar{l}/a))$（$d$ 维、域尺寸 $\bar{l}$、格距 $a$）、门复杂度 $\tilde{O}(t{d}^2/a)$，基于稀疏哈密顿量模拟；Suau 等随后的资源估计显示其渐近门复杂度下禁止性地大的前置常数。量子行走与 PDE 的一般联系（Childs 的图行走、离散行走与格点平移）见quantum-walk-pde——与 LBM 的流送步同构：都是”按离散速度平移振幅”的酉操作（见qlbm-streaming-circuit）。

定位与批判

波动方程的量子模拟是三重意义的标尺：理论上线性 + 酉 + 谱对角化三条件齐备，是复杂度最优的场景（快进定理对这类”量子动力学型”问题放行，见nonlinear-pde-quantum）；工程上它证明了 QFT–对角–QFT 结构在当前硬件可用；批判上它提醒读出墙仍在——波场全场重构需要 $O(N)$ 次测量，真机演示读的是采样分布本身（波场恰是概率幅的解释），这是波动方程独有的便利，不可外推到泊松或 NS。声学、电磁波传播、地震波模拟是其自然应用域（线性波动主导），非线性波动（激波、Burgers 型）则回到nonlinear-pde-quantum的困难清单。

延伸阅读

• quantum-advection-equation
• quantum-walk-pde
• hamiltonian-simulation
• quantum-spectral-method
• nisq-experiments