Qubitization 与二维信号子空间

形式陈述

令 $A=H/\alpha$ 为Hermitian收缩， $\|A\|\le1$ 。先假设给出的信号酉 $S$ 还满足

S=S^\dagger,\qquad S^2=I, \qquad(\langle0^a|\otimes I)S(|0^a\rangle\otimes I)=A.

取辅助零子空间的投影 $\Pi=|0^a\rangle\langle0^a|\otimes I$ ，以及反射 $R=2\Pi-I$ 。Qubitization行走定义为

W=RS.

对 $A$ 的任一特征值 $\lambda\in(-1,1)$ ， $W$ 在一个明确的二维不变子空间内有矩阵

W_\lambda=\begin{pmatrix} \lambda&\sqrt{1-\lambda^2}\\ -\sqrt{1-\lambda^2}&\lambda \end{pmatrix}, \tag{1}

因而特征相位为 $\pm\arccos\lambda$ 。 $\lambda=\pm1$ 的情形退化为一维，须单独处理。[1, Corollary 9]

这里固定 $W=RS$ 的次序和下文第二个基向量的正号约定；换成 $SR$ 或改变基向量符号会改写旋转方向，但不能在同一推导中随意混用 $\arccos$ 与 $\arcsin$ 。

直觉

信号块告诉我们：从辅助零子空间出发，经 $S$ 后返回的振幅就是 $\lambda$ 。酉性迫使剩余幅度大小为 $\sqrt{1-\lambda^2}$ 。

如果 $S$ 本身是一个反射，第二次作用会沿同一个平面回来；再乘辅助反射，就成为该平面中的旋转。每个特征值各有自己的平面，但一套电路同时作用于它们，不必先知道特征向量。

两个信号平面分别对应lambda等于正负二分之一，旋转相位为正负pi/3和正负2pi/3；旁边标出端点正负1退化为一维。

例子与边界

二维基到底是哪两个向量

由Hermitian谱分解，取 $A|u_\lambda\rangle=\lambda|u_\lambda\rangle$ ，令

|g_\lambda\rangle=|0^a\rangle|u_\lambda\rangle, \qquad |b_\lambda\rangle= \frac{S|g_\lambda\rangle-\lambda|g_\lambda\rangle} {\sqrt{1-\lambda^2}}.

信号块条件给 $\Pi S|g_\lambda\rangle=\lambda|g_\lambda\rangle$ ，所以 $\Pi|b_\lambda\rangle=0$ ；它与 $|g_\lambda\rangle$ 正交。再用 $S$ 酉，可算得分子范数平方为 $1-\lambda^2$ ，故 $b_\lambda$ 归一化。

定义本身给

S|g_\lambda\rangle=\lambda|g_\lambda\rangle+s|b_\lambda\rangle, \qquad s=\sqrt{1-\lambda^2}.

使用 $S^2=I$ ，对第二个向量作用得

S|b_\lambda\rangle=s|g_\lambda\rangle-\lambda|b_\lambda\rangle.

因此在基 $(g_\lambda,b_\lambda)$ 上， $S$ 是 $\begin{pmatrix}\lambda&s\\s&-\lambda\end{pmatrix}$ ， $R$ 是 $\operatorname{diag}(1,-1)$ ，二者相乘正好得到式 (1)。这同时证明该平面不变，不能只从一个正确的左上元素猜出整个二维矩阵。

正负特征值与退化端点

若 $H=\operatorname{diag}(\alpha/2,-\alpha/2)$ ，第一平面的矩阵为

\begin{pmatrix}1/2&\sqrt3/2\\-\sqrt3/2&1/2\end{pmatrix},

相位为 $\pm\pi/3$ ；第二平面则把两个对角项换成 $-1/2$ ，相位为 $\pm2\pi/3$ 。负特征值没有被取绝对值，也没有丢失符号。

当 $\lambda=1$ 或 $-1$ ，返回振幅已经具有单位绝对值，所以没有正交泄漏， $Sg_\lambda=\lambda g_\lambda$ 。此时 $b_\lambda$ 的公式不能使用，因为分母为0；直接得到 $Wg_\lambda=\lambda g_\lambda$ ，相位分别为0或 $\pi$ 。

一般块编码为什么要先处理成反射

只知道某个酉 $U$ 左上块为Hermitian $A$ ，并不保证 $U^2=I$ 。缺少这个条件时，上面对 $Sb_\lambda$ 的推导失效， $RU$ 未必停在那一个二维平面。

一种显式修正是增加一位 $c$ ，定义

S'=|0\rangle\langle1|_c\otimes U +|1\rangle\langle0|_c\otimes U^\dagger.

它满足 $S'^\dagger=S'$ 、 $S'^2=I$ 。用新的信号投影 $\Pi'=|+\rangle\langle+|_c\otimes\Pi$ ，其压缩块为

\frac{A+A^\dagger}{2}=A.

于是可对 $R'=2\Pi'-I$ 和 $S'$ 应用前面的证明。[1, Lemma 10]

实现 $S'$ 需要在标志控制下选择 $U$ 或 $U^\dagger$ ，可由一次受控 $U$ 、一次受控 $U^\dagger$ 加一个标志翻转构造。已知电路时可编译这些门；只有未知不带控制的酉黑盒时，不能把新增控制当作自动免费权限。

推论与应用

旋转幂为何产生Chebyshev响应

取整数 $k\ge0$ ，令 $\theta=\arccos\lambda$ 。式 (1) 的 $k$ 次方有左上元素 $\cos(k\theta)$ ，所以返回信号子空间的块为

(\langle0^a|\otimes I)W^k(|0^a\rangle\otimes I) =T_k(A),

其中Chebyshev 多项式 $T_k(x)=\cos(k\arccos x)$ 。这是一条对所有特征值同时成立的算子关系；实现它仍需 $k$ 次行走调用，不能把幂符号当作一步电路。

量子信号处理在这些二维旋转之间插入可选相位，把固定的 $T_k$ 响应推广为满足实现约束的多项式。该用途与标记集合搜索不同：这里要处理的是信号块的谱，而不是仅由某个图的连通性或搜索谱隙决定运行。

归一化和成本仍须保留

行走相位编码的是 $H/\alpha$ 的特征值，不是未缩放 $H$ 。若后续想模拟 $e^{-iHt}$ ，目标函数应是 $\lambda\mapsto e^{-i\alpha t\lambda}$ ；漏掉 $\alpha$ 就模拟了错误的时间尺度。

在已经给出反射型 $S$ 的模型中，一步 $W$ 用一次 $S$ 和一次辅助反射；一般块编码扩张还要按前节计算受控/逆调用。若辅助信号态不是计算基零态而是由PREP制备的 $|G\rangle$ ，反射还需PREP、零反射与PREP逆，不能只数核心信号门。

本页给出的是把矩阵谱变成可处理旋转的接口。态制备、最终读出以及所需多项式的次数仍在后续算法里计费；二维分析简化了证明，不代表实际数据空间只有两维。

参考资料

[1] Guang Hao Low and Isaac L. Chuang, Hamiltonian Simulation by Qubitization, Quantum 3, 163, 2019，§4 Lemma 8、Corollary 9、Lemma 10及Theorem 2：不变平面、反射条件与一般信号酉扩张。本页固定自己的平面基与 $RS$ 次序，直接计算式 (1)。

形式陈述 ​

二维基到底是哪两个向量 ​

正负特征值与退化端点 ​

一般块编码为什么要先处理成反射 ​

旋转幂为何产生Chebyshev响应 ​