📚 课程目录第4章泡利信道与量子噪声的算符和表示

4. 泡利信道与量子噪声的算符和表示

量子噪声在真实硬件上并不是“随机加一个错误”这么简单，它必须用完全正定且保迹的量子信道来描述。算符和表示（operator-sum representation）给出了最通用的描述方式：一个量子信道 \(\mathcal{E}\) 作用在密度矩阵 \(\rho\) 上，可以写成一组 Kraus 算符 \(\{E_k\}\) 的求和形式

\(\mathcal{E}(\rho)=\sum_k E_k \rho E_k^{\dagger}\)，并满足完备性条件 \(\sum_k E_k^{\dagger}E_k = I\)。

这个条件保证了信道是保迹的，即概率总和为 1。泡利信道是一类特别重要的信道：它的 Kraus 算符全部取自泡利矩阵 \(\{I, X, Y, Z\}\) 的倍数。由于泡利矩阵彼此正交且平方为 \(I\)，泡利信道的 Kraus 分解往往具有非常直观的“以某概率施加某个泡利错误”的物理含义。本节围绕三种最典型的泡利信道——比特翻转、相位翻转和去极化信道——展开它们的 Kraus 分解、矩阵形式以及参数含义。

比特翻转信道：以概率 p 施加 X

比特翻转信道（bit-flip channel）描述的是量子比特以概率 \(p\) 发生 \(|0\rangle \leftrightarrow |1\rangle\) 翻转，以概率 \(1-p\) 保持不变。它的 Kraus 算符为

\(E_0=\sqrt{1-p}\,I=\begin{pmatrix}\sqrt{1-p}&0\\0&\sqrt{1-p}\end{pmatrix},\quad E_1=\sqrt{p}\,X=\begin{pmatrix}0&\sqrt{p}\\\sqrt{p}&0\end{pmatrix}\)。

验证完备性：\(E_0^{\dagger}E_0+E_1^{\dagger}E_1=(1-p)I+pI=I\)，满足保迹条件。作用在任意单比特密度矩阵 \(\rho=\begin{pmatrix}\rho_{00}&\rho_{01}\\\rho_{10}&\rho_{11}\end{pmatrix}\) 上，得到

\(\mathcal{E}(\rho)=(1-p)\rho+p\,X\rho X=\begin{pmatrix}(1-p)\rho_{00}+p\rho_{11}&(1-2p)\rho_{01}\\(1-2p)\rho_{10}&(1-p)\rho_{11}+p\rho_{00}\end{pmatrix}\)。

从矩阵形式可以清楚看到：对角元（布居数）被以概率 \(p\) 交换，而非对角元（相干项）被因子 \(1-2p\) 衰减。当 \(p=0.5\) 时，非对角元完全消失，信道把任意态变成完全混合态 \(I/2\)，这正是最大比特翻转噪声的效果。在实验上，比特翻转常用来建模 \(T_1\) 弛豫中的一部分能量耗散过程，或经典比特翻转错误占主导的场景。

相位翻转信道：以概率 p 施加 Z

相位翻转信道（phase-flip channel）以概率 \(p\) 施加相位翻转 \(Z\)，以概率 \(1-p\) 保持不变。其 Kraus 算符为

\(E_0=\sqrt{1-p}\,I,\quad E_1=\sqrt{p}\,Z=\begin{pmatrix}\sqrt{p}&0\\0&-\sqrt{p}\end{pmatrix}\)。

作用在 \(\rho\) 上得到

\(\mathcal{E}(\rho)=(1-p)\rho+p\,Z\rho Z=\begin{pmatrix}\rho_{00}&(1-2p)\rho_{01}\\(1-2p)\rho_{10}&\rho_{11}\end{pmatrix}\)。

与比特翻转不同，相位翻转不改变布居数 \(\rho_{00},\rho_{11}\)，只把非对角元乘以 \(1-2p\)。也就是说，它破坏的是量子相干性而非能量分布。当 \(p=0.5\) 时，相干项被完全抹去，密度矩阵变成对角矩阵，量子叠加退化为经典混合。相位翻转信道是建模 \(T_2\) 退相干（纯退相位）的核心工具，在超导量子比特和自旋量子比特中都非常常见。

去极化信道：三种泡利错误等概率发生

去极化信道（depolarizing channel）是最具对称性的泡利信道。它以一个概率 \(p\) 随机施加 \(X\)、\(Y\)、\(Z\) 三种错误中的一种，每种概率为 \(p/3\)，以概率 \(1-p\) 保持不变。其 Kraus 算符为

\(E_0=\sqrt{1-p}\,I,\quad E_1=\sqrt{p/3}\,X,\quad E_2=\sqrt{p/3}\,Y,\quad E_3=\sqrt{p/3}\,Z\)。

完备性验证：\((1-p)I+\frac{p}{3}(X^{\dagger}X+Y^{\dagger}Y+Z^{\dagger}Z)=(1-p)I+\frac{p}{3}\cdot 3I=I\)。作用在 \(\rho\) 上，利用泡利矩阵的恒等式 \(\sum_{P\in\{I,X,Y,Z\}}P\rho P^{\dagger}=2\,\mathrm{tr}(\rho)I=2I\)（对单比特），可得

\(\mathcal{E}(\rho)=(1-p)\rho+\frac{p}{3}(X\rho X+Y\rho Y+Z\rho Z)=(1-\frac{4p}{3})\rho+\frac{4p}{3}\cdot\frac{I}{2}\)。

这个结果非常漂亮：去极化信道把任意输入态以概率 \(4p/3\) 替换为完全混合态 \(I/2\)，以概率 \(1-4p/3\) 保留原态。注意这里的“有效去极化概率”是 \(4p/3\)，而不是 \(p\)，这是初学者最容易混淆的地方。当 \(p=3/4\) 时，有效概率达到 1，输出恒为 \(I/2\)。去极化信道在量子纠错阈值分析、纠缠提纯和量子容量计算中都是标准模型。

三种信道的对比与参数关系

把三种信道放在一起看，可以发现它们共享同一个结构：以概率 \(1-p\) 保留原态，以概率 \(p\) 施加某个泡利错误。区别只在于施加的是单个 \(X\)、单个 \(Z\)，还是 \(X,Y,Z\) 三种等概率混合。比特翻转只交换布居数，相位翻转只衰减相干项，而去极化信道同时影响布居数和相干项，并且对任意输入态都给出相同的“收缩”效果。

在参数关系上，去极化信道的有效去极化概率是 \(4p/3\)，这意味着当 \(p\) 较小时，去极化信道对态的破坏速度是比特翻转或相位翻转的 \(4/3\) 倍。如果希望去极化信道与比特翻转信道在相干项衰减上等效，需要满足 \(1-2p_{\text{bit}}=1-4p_{\text{dep}}/3\)，即 \(p_{\text{dep}}=3p_{\text{bit}}/2\)。这个换算关系在把不同文献中的噪声参数统一到同一模型时非常实用。

Kraus 分解的矩阵形式与 Bloch 球图像

对单比特态，密度矩阵也可以用 Bloch 矢量 \(\vec{r}=(r_x,r_y,r_z)\) 表示：\(\rho=(I+\vec{r}\cdot\vec{\sigma})/2\)。泡利信道对 Bloch 矢量的作用可以写成线性收缩：

比特翻转：\(r_x \to (1-2p)r_x,\ r_y \to (1-2p)r_y,\ r_z \to r_z\)；
相位翻转：\(r_x \to (1-2p)r_x,\ r_y \to (1-2p)r_y,\ r_z \to r_z\)（注意相位翻转在 Bloch 球上表现为绕 Z 轴的收缩，但 \(r_z\) 不变）；
去极化：\(\vec{r} \to (1-4p/3)\vec{r}\)，即整个 Bloch 球均匀收缩。

这里需要特别说明：相位翻转信道在 Bloch 球上把 \(r_x,r_y\) 同时乘以 \(1-2p\)，而 \(r_z\) 保持不变，因此 Bloch 球被压扁成一个椭球。去极化信道则把整个球均匀缩小，保持球对称性。比特翻转信道把 \(r_x,r_y\) 乘以 \(1-2p\)，\(r_z\) 不变，但注意 \(X\) 错误实际上交换 \(|0\rangle\) 和 \(|1\rangle\)，在 Bloch 球上对应绕 X 轴旋转 180°，因此其完整作用并不是简单的分量收缩，而是先旋转再收缩。上面的分量公式是在“错误以概率 p 随机发生”的平均意义上成立的，读者在推导时应以密度矩阵的 Kraus 求和为准，Bloch 矢量公式作为辅助直观。

实际计算示例

假设输入态为 \(|+\rangle=(|0\rangle+|1\rangle)/\sqrt{2}\)，其密度矩阵为 \(\rho=\frac{1}{2}\begin{pmatrix}1&1\\1&1\end{pmatrix}\)。取比特翻转信道参数 \(p=0.1\)，则

\(\mathcal{E}(\rho)=0.9\rho+0.1 X\rho X=0.9\cdot\frac{1}{2}\begin{pmatrix}1&1\\1&1\end{pmatrix}+0.1\cdot\frac{1}{2}\begin{pmatrix}1&1\\1&1\end{pmatrix}=\frac{1}{2}\begin{pmatrix}1&0.8\\0.8&1\end{pmatrix}\)。

可以看到非对角元从 1 衰减到 0.8，因子正好是 \(1-2p=0.8\)。若改用去极化信道且希望达到相同的相干项衰减，需要 \(1-4p/3=0.8\)，解得 \(p=0.15\)。这说明在相同“名义错误概率”下，去极化信道对相干性的破坏更强；反过来，若两种信道要产生相同的相干衰减，去极化信道的参数应取比特翻转的 1.5 倍。

在量子纠错码的性能评估中，通常把物理错误统一折算为去极化信道参数 \(p_{\text{dep}}\)，因为去极化模型对 X、Y、Z 错误对称，便于与纠错码的阈值比较。若实验测得的是比特翻转率 \(p_X\) 和相位翻转率 \(p_Z\)，可以近似取 \(p_{\text{dep}}\approx p_X+p_Z\) 作为保守估计，再根据具体码型调整。

注意区分“Kraus 算符个数”与“错误类型数”。比特翻转和相位翻转各有两个 Kraus 算符，去极化有四个。Kraus 算符的个数并不等于信道“复杂度”，而是取决于所选分解方式；同一个信道可以有不同个数的 Kraus 表示，但最小个数由信道 Choi 矩阵的秩决定。去极化信道的最小 Kraus 个数为 4，不能压缩到 2 个。

总结来说，泡利信道的算符和表示为量子噪声建模提供了统一而简洁的框架：比特翻转用 \(\{I,X\}\)，相位翻转用 \(\{I,Z\}\)，去极化用 \(\{I,X,Y,Z\}\)，每个 Kraus 算符前的系数平方对应错误发生的概率。掌握这三种信道的 Kraus 分解与矩阵形式，是后续分析量子纠错、退相干和噪声阈值的基础。