对流方程量子模拟

一维对流方程 ${\partial }_{t}u+{\partial }_{x}u=0$ 的解是精确平移 $u(x,t)=u_0(x-t)$——在离散格点上，一个时间步就是循环移位矩阵作用在振幅编码的态上，而移位是置换、置换是酉的。对流方程因此是全部 PDE 中与量子硬件”血缘最近”的一个：不需要任何酉化技巧就能走完一个时间步。它同时是波动方程的因子（见quantum-wave-equation）、LBM 流送步的宏观极限（见qlbm-streaming-circuit）和对流扩散方程的极限情形（见advection-diffusion-equation）。

平移算子的酉结构

把 $N=2^n$ 个格点上的 $u$ 振幅编码为 $|u⟩={\sum }_j{u}_j|j⟩$。周期边界下向前平移一步是循环移位

$$P|j⟩=|(j+1)\mathrm{mod}N⟩,P^{†}P=I,$$

显然酉。更有用的是它在 Fourier 基下对角：$P=\mathrm{QFT}D{\mathrm{QFT}}^{†}$，$D=\mathrm{diag}(e^{2\pi ik/N})$ 的相位全在单位圆上——与热方程的实衰减谱（见quantum-heat-equation）对照，对流算子的谱是纯虚的，演化 $e^{-ikt}$ 天然酉。这就是量子计算对”输运主导”问题的结构性契合：多步对流 = 幂 $P^t$，仍是酉，可以无测量地连续推进；Todorova–Steijl 的无碰撞流送线路正是把这一性质按 LBM 的每个离散速度方向各实现一遍（$\hat{S}={⨁}_i{\hat{S}}_i$ 的受控 increment 门，见qlbm-streaming-circuit）Sanavio & Succi 2024 (IntechOpen)。

迎风格式的量子化：非厄米的处理

纯平移是理想化；带变系数、迎风离散后的对流算子矩阵非厄米（虽仍可近似酉）。Hu–Jin–Liu–Zhang 的 Schrödinger 化电路首次处理了迎风格式的非厄米矩阵：对迎风差分产生的算符做 warped phase 抬升，分步 LCU 实现演化模块，数值实验（Qiskit、10000 shots）显示解收敛于 $e^{At}$ 参考且在间断处非振荡——优于此前中心差分文献的振荡行为Hu et al. 2024 (Quantum)。

图:迎风格式对流方程的数值解对比——矩阵指数参考、经典 Schrödinger 化与量子电路实现一致,间断处无振荡(来源:QC029)

复杂度上对流方程比热方程更宽松：CFL 条件为 $\tau =O(h)$（热方程是 $O(h^2)$），据该文分析经典有限差分成本 $O(s\cdot 2^{n_{x}d}(T^2/ϵ+T/h))$，量子电路取得对数级优势的条件放宽到 $d>4$。注意这个优势声明同样受 CFL–时间步数耦合的制约（显式推进的步数随网格多项式增长，见nonlinear-pde-quantum的端到端批判）。

NISQ 演示：一维波动方程

Liu–Lopes 把”对流的酉结构”推成了一个完整的 NISQ 真机演示：一维波动方程 ${\partial }_t^2\psi =c^2{\partial }_x^2\psi$ 因式分解为左右行波（两个对流方程），引入 $\varphi =i(-\Delta {)}^{-1/2}{\partial }_t\psi$ 后化为薛定谔形式 $i{\partial }_t(\psi ,\varphi {)}^T=H_W(\psi ,\varphi {)}^T$，波动哈密顿量

$$H_W=\left(\begin{array}{cc}0& \sqrt{-\Delta }\\ \sqrt{-\Delta }& 0\end{array}\right)$$

厄米；周期边界下 QFT 对角化 ${\Delta }_a$，于是 $e^{-it{H}_W}$ 实现为 ${\mathrm{QFT}}^{†}\to e^{-it{Z}_0\otimes \sqrt{-D}}\to \mathrm{QFT}$——对角算子夹在两个 QFT 之间，$n+1$ 比特即可模拟，双比特门数 $O(n^2)$。他们在 Quantinuum H1-1 离子阱真机上模拟了 Ricker 子波的传播（6 比特、10000 shots），保真度随比特数指数提升、随时间 $O(t^2)$ 退化Liu & Lopes 2024 (arXiv)（完整推导见quantum-wave-equation）。这个演示对对流类问题的启示是：QFT–对角相位–QFT 的三明治结构在当前硬件上已经可行。

对流扩散与 LBM 量子化

加入扩散项后算符谱离开单位圆，纯酉性破坏——这正是 Budinski 对流扩散量子 LBM 的处理对象：D1Q2/D2Q5 格子上，流送（对流）用受控移位保持酉，碰撞（含松弛/扩散）非酉经 LCU 分解 + 后选择重缩放实现；且若做适当重缩放，可连续多个时间步而无须测量重初始化Budinski 2021 (QIP)（细节见budinski-qlbm）。速度场固定时对流算子可在步间复用；一旦速度场随时间更新（非线性流），复用即失效——非线性的代价再次从侧面显现Sanavio & Succi 2024 (IntechOpen)。

延伸阅读

• quantum-wave-equation
• qlbm-streaming-circuit
• budinski-qlbm
• quantum-heat-equation
• advection-diffusion-equation