技术摘要：量子斜交特征投影 (Quantum Oblique Eigenprojection)

问题陈述

本文解决了在量子计算机上执行斜交特征投影 (oblique eigenprojections) 的挑战。虽然针对厄米矩阵（其特征空间是正交的）的量子算法通过量子奇异值变换 (QSVT) 等框架已经非常成熟，但通用的非正规矩阵仍然是一个显著的障碍。对于非正规矩阵，希尔伯特空间分解为广义的、非正交的特征子空间。任务是构建一种量子算法，将输入态投影到与目标复谱区域（例如单位圆盘内的特征值）相关的特定不变子空间上，同时抑制与补区域相关的分量。

当谱为复数时，这个问题存在根本性的阻碍。与实谱情况不同，在复平面上，多项式无法近似指示函数（indicator function）而不产生显著误差，这使得解析函数（以及多项式）无法近似此类函数。因此，标准的 QSVT 和量子特征值变换 (QEVT) 技术——这些技术依赖于对解析函数的多项式逼近——无法直接解决此问题。

以往的方法（如使用柯西积分公式的留数积分法）可以处理非解析函数，即沿包围目标谱的轮廓线对矩阵留数 $(\zeta I - A)^{-1}$ 进行积分。然而，这些方法存在一个关键的低效性：它们会引入一个归一化因子 $\alpha$，其规模为 $\Theta(\|\Pi\|/\delta)$，其中 $\|\Pi\|$ 是投影的范数，$\delta$ 是谱间隙。这种膨胀的归一化因子显著增加了使用该投影作为子程序的任何算法的查询复杂度，特别是由于块编码的成功概率与 $\alpha^2$ 成反比。

方法论

作者提出了一种新颖的方法，通过利用由离散化留数积分生成的块编码的完整结构，改进了标准的留数积分法。

1. 离散化留数积分： 通过使用梯形法则对留数轮廓积分进行离散化，来近似斜交特征投影 $\Pi$。这产生了一个移位留数之和：$\Pi \approx \frac{1}{n} \sum_{j=0}^{n-1} (I - \omega^{-j}_n A)^{-1}$。
2. 块编码构建： 标准方法构造一个包含这些移位留数的块对角矩阵 $M$，并使用 QSVT 对其求逆。通过对辅助寄存器应用 Hadamard 变换，结果矩阵的左上角块近似于 $\Pi$。然而，剩余的块缺乏结构，迫使归一化因子由整个矩阵的最坏情况范数决定，即 $\Theta(1/\delta)$。
3. 离散傅里叶变换 (DFT) 分析： 本研究的核心洞察在于 $M^{-1}$ 的块并非任意的；它们具有循环结构 (circulant structure)。通过使用量子傅里叶变换 (QFT) 而不是 Hadamard 变换对 $M^{-1}$ 进行共轭，作者揭示了每一个块 $(k, l)$ 都对应于移位留数的特定离散傅里叶模 $S_m$。
   • 零阶模 $S_0$ 近似于特征投影 $\Pi$。
   • 高阶模 $S_m$ 在满足谱集条件（例如，如果矩阵是可对角化的或具有有界的数值范围）时呈几何级数衰减。
4. 双侧块预处理： 作者引入了一种在傅里叶基底下的双侧块预处理技术。他们应用一个在计算基底下的对角缩放算符 $T$ 来修改求逆前的矩阵。
   • 随着缩放参数 $t$ 的减小，预处理后的矩阵的第一行和第一列会增大。
   • 然而，这种预处理抑制了提取出的左上角块的逆间隙增长 (inverse-gap growth)，将其从 $\Theta(\|\Pi\|/\delta)$ 降低到 $O(c \|\Pi\|)$，其中 $c$ 是取决于谱集条件的常数（例如，对于数值范围，$c = 1+\sqrt{2}$；对于可对角化矩阵，则是条件数）。
   • 这种归一化是接近最优的，因为 $\|\Pi\|$ 是任何 $\Pi$ 的块编码的理论下界。

主要结果

论文确立了三个主要定理：

• 定理 1 (矩阵留数的 DFT)： 提供了在傅里叶基底下留数积分离散化的闭式表达式。它表明，这些块是由矩阵在内部和外部不变子空间上的限制决定的，并随移位的幂次几何级数衰减。
• 定理 2 (单位圆盘投影)： 证明了对于谱与单位圆不相交的矩阵 $A$，其指向内部子空间的斜交特征投影 $\Pi$ 可以以归一化因子 $\alpha_\Pi = O(c \|\Pi\|)$ 和查询复杂度 $O\left(\frac{c \alpha_A \|\Pi\|}{\delta} \log \frac{c \|\Pi\|}{\varepsilon} \log \frac{1}{\varepsilon}\right)$ 进行块编码。这与标准方法的查询复杂度相当，但极大地改善了归一化因子。
• 定理 3 (一般区域)： 将结果扩展到复平面中的任意目标区域。通过使用解析分离映射 $f$（例如，将目标区域映射到单位圆盘的有理函数），可以将问题转化为应用于 $f(A)$ 的单位圆盘情况。成本取决于 $f(A)$ 到单位圆的间隙。

应用

作者展示了这一改进的基元如何为以下三个特定应用带来显著加速：

1. 复谱特征态制备： 将 Low 和 Su 的 QEVT 算法扩展到具有复特征值的非正规矩阵。新方法制备特征态的查询复杂度缩放为 $O(\|\Pi\|^2 / (\eta \delta))$，相比于标准留数积分法的 $O(\|\Pi\|^2 / (\eta \delta^2))$ 实现了二次提升。这种提升是因为振幅放大轮数与间隙无关（缩放为 $O(\|\Pi\|/\eta)$），而标准方法需要 $O(\|\Pi\|/(\eta \delta))$ 轮。
2. 连续时间代数黎卡提方程 (Algebraic Riccati Equations)： 求解这些方程（源于最优控制和量子化学）等价于投影到非正规哈密顿量的不变子空间。与 Rodenas-Ruiz, Zhao, 和 Lee 之前的求解器相比，新方法在间隙依赖性上实现了三次提升，消除了总复杂度中的两个 $1/\delta$ 因子。
3. 普通 Sylvester 方程： 通过投影到块矩阵的不变子上来求解 $AX + XB = C$。该方法相对于 Wang 和 Liu 的直接增广方法，在间隙依赖性上实现了二次提升。

重要性与主张

本文声称，这项工作为涉及非正规矩阵的量子算法提供了一个基础基元 (fundamental primitive)。通过在保持与厄米矩阵相似的近线性间隙查询复杂度的同时，实现接近最优 $\|\Pi\|$ 的归一化因子，作者弥合了 QSVT 处理厄米矩阵的高效性与处理非厄米问题所需能力之间的鸿沟。

作者强调，其结果的意义在于：

• 它解决了长期困扰非解析函数留数相关方法的“归一化因子膨胀”问题。
• 它证明了非解析矩阵函数（特别是斜交特征投影）可以在量子计算机上实现，且其复杂度与处理厄米矩阵的 QSVT 相当，尽管无法直接使用多项式逼近。
• 通过 DFT 分析全块结构并应用双侧预处理的技术，为优化块编码提供了一种新的范式，这可能适用于其他非解析函数。

作者对研究范围保持了审慎态度，指出预处理依赖于留数特定的闭式 DFT，并不立即适用于所有构造（例如 Wang 和 Liu 的乘积留数法）。他们还指出，与 QSVT 中由目标函数的几何形状决定成本不同，本方法的成本受输入矩阵谱间隙的支配。