技术摘要：用于微分方程的有效哈密顿量量子求解器

问题陈述

微分方程（DEs）是科学建模的基础，但对于高维、非线性或多尺度问题，其计算需求仍然非常高。虽然量子计算在处理线性代数和模拟方面具有潜在优势，但现有的量子微分方程求解器在状态准备、读取成本和收敛保证方面面临权衡。一个特定的挑战在于“有效哈密顿量”（effective-Hamiltonian）方法，即通过将微分方程编码，使得其解对应于一个正半定哈密顿量的基态。该方法的先前迭代版本（例如 Wu 等人的研究）面临三个主要局限性：

1. 非零条件： 难以直接将通用的非零初始、边界或数据条件纳入齐次形式（$A|f\rangle = 0$）中，而不丢失解的物理尺度。
2. 非线性简并： 对于非线性微分方程，将状态空间提升到张量积（例如 $|f\rangle^{\otimes D}$）会引入大规模的基态简并。零能量子空间中的大多数状态都不具备所需的重复乘积结构，从而无法表示物理解，导致求解器失效。
3. 表示限制： 依赖全局谱基（如切比雪夫多项式）对于具有局部特征或不连续性的解可能效率低下，而替代的网格化表示则需要不同的算符构建方式。

方法论

本文通过四个主要的工程技术发展扩展了有效哈密顿量框架：

1. 通过齐次化处理非零条件

作者提出将多个非零初始、边界和数据条件转换为齐次的、尺度不变的约束。通过利用一个已知的非零参考条件 $f^{(m)}(x_s) = f_s \neq 0$，任何其他条件 $f^{(n_l}(x_l) = f_l$ 都可以重写为相对约束：
$$\left( \phi(x_l)^T D^{n_l} - \frac{f_l}{f_s} \phi(x_s)^T D^m \right) |f\rangle = B_l |f\rangle = 0$$
这使得非零数据可以直接纳入标准的 $A|f\rangle = 0$ 形式中，而无需使用增广系统。

2. 另一种右端项形式 ($A|f\rangle = |b\rangle$)

引入了一种基于线性系统 $A|f\rangle = |b\rangle$ 的替代构建方式，其中源项和非零数据被编码在 $|b\rangle$ 中。有效哈密顿量定义为：
$$H_b = A^\dagger (I - |b\rangle\langle b|) A$$
该算符是正半定的，其基态（零能量）对应于满足 $A|f\rangle$ 落在 $|b\rangle$ 张成空间内的解。这种方法自然地处理了非齐次方程，而不依赖于特定的非零参考条件进行缩放，尽管它要求 $|b\rangle \neq 0$。

3. 解决非线性简并问题

对于需要提升空间 $|f\rangle^{\otimes D}$ 的非线性微分方程，作者分析了由于零能量态不具备所需乘积结构而导致的非物理简并。他们提出了两种缓解策略：

• Ansatz 层级的限制： 在变分算法（如 VQE）中，将试探态限制为形式 $|\Psi(\theta)\rangle = |u(\theta)\rangle^{\otimes D}$，可以确保搜索空间仅包含物理有效的乘积态，从而在优化过程中有效地绕过简并。
• 哈密顿量约束： 本文研究了通过添加对称性约束（例如强制执行 $f_j f_k = f_k f_j$）来减少简并，并指出可能需要针对特定问题的约束来完全打破简并。

4. 网格上的振幅编码

该框架被扩展到了振幅编码，其中量子比特的振幅代表离散网格上的函数值，而非谱系数。这需要构建特定的算符：

• 导数算符 ($D$)： 通过有限差分矩阵实现。
• 乘法算符 ($M_a$)： 用于已知函数的对角矩阵。
• 非线性乘积算符 ($N$)： 用于重新排列张量积振幅以表示逐点乘积的算符。
  这种方法避免了基函数截断误差，但引入了网格离散化误差，并且需要 $O(2^n)$ 的经典输出来进行完整的状态重建。

关键结果

所提方法通过 Python（OpenQARP）在无噪声态矢量模拟中，在四个微分方程示例上进行了验证：

1. RC 电路充电（线性）：

   • 证明了对于零初始条件（$\beta=0$），标准的 $A_0$ ($A|f\rangle=0$) 构建由于缺乏非零参考而失效，而 $A_b$ ($A|f\rangle=|b\rangle$) 构建仍然有效。
   • 对于非零条件，两种构建方式都能得到准确的解，尽管它们根据编码方式（切比雪夫 vs 振幅）表现出不同的谱间隙。

2. Airy 方程（线性，振荡型）：

   • 表明在固定量子比特数（$N=4$）下，切比雪夫编码在处理平滑、振荡解方面优于有限差分振幅编码，这归功于谱收敛性。
   • 确认了 $A_0$ 和 $A_b$ 构建均能成功处理多个非零初始条件。

3. 复合壁热传导（线性，导数不连续）：

   • 强调了全局切比雪夫展开无法捕捉材料界面处的导数不连续性，从而导致错误的通量约束。
   • 域分解切比雪夫编码和振幅编码成功重现了解析解，证明了根据解的正则性调整表示形式的必要性。

4. Fisher-KPP 行进波前（非线性）：

   • 应用了针对二次非线性的张量积提升。
   • 验证了通过将 VQE ansatz 限制为重复乘积态（$|u\rangle^{\otimes 2}$），可以成功恢复行进波前解，从而避免了非物理基态。
   • 观察到完整的提升哈密顿量仍然具有高度简并性（例如，切比雪夫编码的简并度为 34–47，振幅编码为 211–240），因此必须使用 ansatz 限制或数值延拓技术来引导优化。

意义与主张

本文声称通过以下方式扩大了基于基态的量子微分方程求解器的适用范围：

• 泛化约束： 使将多个非零条件直接纳入齐次形式成为可能，并提供了一种不依赖于特定参考条件的替代非齐次形式 ($A_b$)。
• 澄清权衡： 提供了一项对比分析，表明选择谱编码（切比雪夫）还是网格（振幅）编码，很大程度上取决于解的平滑度及其边界条件。
• 解决非线性： 为处理提升空间中的非线性提供了一条切实可行的路径（ansatz 限制），同时承认哈密顿量的基本简并问题仍是一个开放性挑战。

作者保持了审慎的态度，指出虽然无噪声模拟验证了哈密顿量构建和解的恢复，但尚未建立高效的量子基态准备或明确的量子优势。文中指出，未来的工作需要在变分训练保证、贫瘠高原（barren plateaus）以及处理非多项式非线性的可扩展性方面进行深入研究。