热方程量子算法

热方程 ${\partial }_{t}u=\nu \Delta u$ 是扩散与耗散的原型——NS 的粘性项、湍流的能量级串终点都是它。它在量子算法里有双重身份：既是”非酉动力学如何模拟”的标准测试题（$e^{\nu t\Delta }$ 压缩范数，不是酉演化），又是 Schrödinger 化技术的第一演示对象——Jin–Liu–Yu 的原始论文Jin et al. 2024 (PRL)与 Hu–Jin–Liu–Zhang 的显式电路实现Hu et al. 2024 (Quantum)都以它为主算例；经 Cole–Hopf 变换它还是 Burgers 方程量子算法的隐藏内核Uchida et al. 2024 (arXiv)。

对角化：扩散在 Fourier 基下是乘法

对周期边界条件做 Fourier 展开 $u(x,t)={\sum }_k{\hat{u}}_k(t)e^{ikx}$，热方程逐模解耦：

$$\frac{d{\hat{u}}_k}{dt}=-\nu k^2{\hat{u}}_k⟹{\hat{u}}_k(t)=e^{-\nu k^{2}t}{\hat{u}}_k(0).$$

每个模式的演化只是幅度乘以 $e^{-\nu k^{2}t}$——差分离散后 Laplacian 被（量子）Fourier 变换对角化：${\Delta }_a=\mathrm{QFT}D{\mathrm{QFT}}^{†}$，$D$ 为对角特征值矩阵。问题恰恰出在 $e^{-\nu k^{2}t}<1$：演化算符 $e^{\nu t\Delta }$ 压缩范数、非酉，量子计算机不能直接执行。这是扩散类方程与波动方程（特征值实数符号翻转即得酉演化，见quantum-wave-equation）的本质分野，也是两条酉化技术的出发点。

酉化技术一：辅助比特 LCU

把非酉的 $e^{\nu t\Delta }$ 分解为酉算子的线性组合（linear combination of unitaries），在辅助比特的叠加分支上分别执行各酉分量，再以受控方式重组——成功概率随范数压缩而衰减，需后选择或幅度放大。Budinski 的 QLBM 碰撞步（$A=\frac{1}{2}(C_1+C_2)$ 的两酉分解）就是这一思路在格子玻尔兹曼语境下的实例（见budinski-qlbm）。LCU 路线通用，但每步后选择的成功率是反复时间步进的硬伤。

酉化技术二：Schrödinger 化

Jin–Liu–Yu 的 warped phase 变换给出系统解法：引入新变量 $w=e^{-p}u$（$p$ 为新增”扭曲”维），热方程 ${\partial }_{t}u=(\Delta -V)u$ 抬升为

$${\partial }_{t}w+(\Delta -V){\partial }_{p}w=0,$$

初值取偶延拓 $w(0,x,p)=\exp (-|p|)u_0(x)$；对 $p$ 做 Fourier 变换（变量 $\eta$）后得到解耦的实时间薛定谔方程组

$$i{\partial }_t\tilde{w}=\eta (-\Delta +V)\tilde{w},$$

每个 $\eta$ 模式是一个哈密顿量为 $\eta (-\Delta +V)$ 的酉演化——非酉的耗散被”搬进”了 $\eta$ 谱，量子侧只需标准 Hamiltonian 模拟。原解由 $u=\int \chi (p)wdp$ 恢复（细节见schrodingerization）Jin et al. 2024 (PRL)。

显式电路与高维优势条件

Hu–Jin–Liu–Zhang 把上述框架落成 Qiskit 门级电路：部件包括分步 LCU 的时间演化模块 $W_j(\gamma \tau ,\lambda )$、$V_0(\tau )$、热方程模块 $V_{\mathrm{heat}}(\tau )$，测量只需投影到 $p_k>0$ 的 Fourier 寄存器分支，含（逆）QFT。数值实验（状态向量 + 10000 shots）显示 $n_p=3,5,7$ 的 $p$ 寄存器随位数增大收敛于矩阵指数 $e^{At}$ 参考解，能量 $\Vert u{\Vert }^2$ 的测量与经典一致Hu et al. 2024 (Quantum)。

图:热方程数值解对比——矩阵指数算子 e^{At}(参考)、经典 Schrödinger 化实现与所提量子电路(n_p 位 p 寄存器,两种时间步)四者一致(来源:QC029)

复杂度对比（据该文分析）：经典有限差分到精度 $ϵ$ 需 $O(s\cdot 2^{n_{x}d}(T^2/ϵ+T/h^2))$（热方程 CFL 条件 $\tau =O(h^2)$ 是主导项）；量子电路对空间维度是对数的。结论是当 $h=O(ϵ)$ 且 $d>5$ 时量子取得对数级优势——扩散方程的量子价值在高维（多粒子相空间扩散、高维 Fokker–Planck），而不是一二维物理空间。这一定性结论与泊松/有限元分析一致（见quantum-poisson-solver、quantum-fem）。

接口：Cole–Hopf 与 Burgers

粘性 Burgers 方程 ${\partial }_{t}u+u{\partial }_{x}u=\nu {\partial }_x^{2}u$ 经 Cole–Hopf 变换 $u=-\nu {\partial }_x\ln \psi$ 精确地化为热方程——非线性被一次性变换消去，无需 Carleman 嵌入。Uchida 等的 Burgers 量子算法正是”量子热方程求解器 + 逆变换放在信息提取步”的组合：先演化 $\psi$ 的量子态，再在测量阶段直接提取多点函数 $P^{(n)}$ 与矩（两点关联、平坦度）等统计量（见quantum-burgers）Uchida et al. 2024 (arXiv)。反向地，变分路线中 Feynman–Kac 型抛物方程经 Wick 旋转与薛定谔方程相连，可用变分虚时演化求解（见feynman-kac）。

延伸阅读

• schrodingerization
• quantum-advection-equation
• quantum-burgers
• hamiltonian-simulation
• feynman-kac