2. 量子噪声的数学表征
在理想模型中，量子算法用纯态 \(|\psi\rangle\) 和酉变换 \(U\) 描述。但真实量子处理器上，量子比特会与环境耦合，发生退相干、能量弛豫和门操作误差。要分析噪声如何影响算法，必须先把“噪声”从物理直觉提升为可计算的数学对象。本节围绕密度矩阵、Kraus 算符、量子信道，以及三类最常用的噪声模型展开，并说明它们如何改变量子态的演化轨迹。
一、密度矩阵：从纯态到混合态的统一描述
纯态 \(|\psi\rangle\) 的密度矩阵定义为 \(\rho = |\psi\rangle\langle\psi|\)。它满足 \(\rho^\dagger=\rho\)、\(\mathrm{Tr}(\rho)=1\)、\(\rho^2=\rho\)。当系统以概率 \(p_i\) 处于 \(|\psi_i\rangle\) 时，混合态写作：
ρ = Σ_i p_i |ψ_i⟩⟨ψ_i|, Σ_i p_i = 1, p_i ≥ 0此时 \(\mathrm{Tr}(\rho^2)<1\)，纯度下降。密度矩阵的关键价值在于：测量算符 \(M\) 的期望值为 \(\langle M\rangle=\mathrm{Tr}(\rho M)\)，因此无论纯态还是混合态，都能用同一套线性代数工具处理。对单量子比特，\(\rho\) 可展开为布洛赫球形式：
ρ = (I + r_x X + r_y Y + r_z Z) / 2其中 \((r_x,r_y,r_z)\) 是布洛赫向量。纯态对应 \(|r|=1\)，混合态对应 \(|r|<1\)。噪声对量子态的影响，直观上就是让布洛赫向量缩短、旋转或向某方向偏移。
二、Kraus 算符与量子信道
封闭系统的演化是酉变换 \(\rho' = U\rho U^\dagger\)。开放系统的演化更一般，由量子信道 \(\mathcal{E}\) 描述。任何完全正定且保迹的映射都可写成 Kraus 表示：
E(ρ) = Σ_k K_k ρ K_k†, 且 Σ_k K_k† K_k = I其中 \(\{K_k\}\) 称为 Kraus 算符，个数由环境维度决定。保迹条件 \(\sum_k K_k^\dagger K_k=I\) 保证概率守恒。若只有一个 Kraus 算符且满足 \(K^\dagger K=I\)，信道退化为酉演化。Kraus 表示不唯一，但对应的信道作用唯一。实际计算中，我们通常先写出 Kraus 算符，再作用于密度矩阵，得到噪声后的状态。
三、泡利信道：比特翻转、相位翻转与去极化
泡利信道用泡利算符 \(I,X,Y,Z\) 构造，是最常用的离散噪声模型。单比特泡利信道写作：
E(ρ) = (1-p_x-p_y-p_z) ρ + p_x XρX + p_y YρY + p_z ZρZ当 \(p_x=p_y=p_z=p/3\) 时，称为去极化信道：
E_depol(ρ) = (1-p) ρ + (p/3)(XρX + YρY + ZρZ)去极化信道把布洛赫向量按因子 \((1-4p/3)\) 均匀缩短，不改变方向。若只保留 \(X\) 项，就是比特翻转信道；只保留 \(Z\) 项，就是相位翻转信道。泡利信道适合描述门误差、测量误差和部分退相干，因为其 Kraus 算符就是泡利矩阵的平方根形式，便于在稳定子形式、纠错码和蒙特卡洛轨迹模拟中使用。
四、振幅阻尼与相位阻尼
振幅阻尼描述能量弛豫，例如超导量子比特从 \(|1\rangle\) 衰减到 \(|0\rangle\)。其 Kraus 算符为：
K0 = [[1, 0], [0, √(1-γ)]]
K1 = [[0, √γ], [0, 0]]其中 \(\gamma=1-e^{-t/T_1}\)，\(T_1\) 是能量弛豫时间。作用在 \(|1\rangle\) 上：
E(|1⟩⟨1|) = (1-γ)|1⟩⟨1| + γ|0⟩⟨0|可见激发态概率按 \(e^{-t/T_1}\) 衰减，同时相干项 \(\rho_{01}\) 乘以 \(\sqrt{1-\gamma}\)。振幅阻尼既改变布洛赫向量的 \(z\) 分量，也压缩横向分量。
相位阻尼只破坏相干性，不交换能量。其 Kraus 算符为：
K0 = [[1, 0], [0, √(1-λ)]]
K1 = [[0, 0], [0, √λ]]作用后，\(\rho_{01}\) 和 \(\rho_{10}\) 乘以 \(\sqrt{1-\lambda}\)，而 \(\rho_{00}\)、\(\rho_{11}\) 不变。相位阻尼对应 \(T_2\) 过程，布洛赫向量的 \(x,y\) 分量按 \(e^{-t/T_2}\) 收缩，\(z\) 分量不变。若同时存在振幅阻尼和相位阻尼，横向相干通常按 \(e^{-t/T_2}\) 衰减，而 \(T_2\le 2T_1\)。
五、噪声对量子态演化的影响
把上述信道作用到具体量子态上，可以清楚看到噪声如何改变算法结果。以单比特叠加态 \(|\psi\rangle=(|0\rangle+|1\rangle)/\sqrt{2}\) 为例，其密度矩阵为：
ρ = 1/2 [[1, 1],
[1, 1]]经过相位阻尼后：
ρ' = 1/2 [[1, √(1-λ)],
[√(1-λ), 1]]对角项不变，非对角项衰减。这意味着测量得到 \(|0\rangle\) 和 \(|1\rangle\) 的概率仍各为 1/2，但干涉效应减弱。若算法依赖相位干涉（如量子傅里叶变换、相位估计），相位阻尼会直接降低干涉条纹可见度，导致结果分布变宽。
经过振幅阻尼后：
ρ' = 1/2 [[1+γ, √(1-γ)],
[√(1-γ), 1-γ]]此时 \(|1\rangle\) 概率从 1/2 降为 \((1-\gamma)/2\)，\(|0\rangle\) 概率升为 \((1+\gamma)/2\)，布洛赫向量向 \(+z\) 方向偏移。若算法依赖 \(|1\rangle\) 态或激发态布居，振幅阻尼会引入系统性偏差。
对于多比特系统，噪声还会通过纠缠传播。例如两比特 Bell 态 \(|\Phi^+\rangle=(|00\rangle+|11\rangle)/\sqrt{2}\) 在其中一个比特经历振幅阻尼后，密度矩阵的非对角项 \(\rho_{00,11}\) 乘以 \(\sqrt{1-\gamma}\)，纠缠度下降。若两个比特都经历去极化，保真度按 \((1-4p/3)^2\) 量级衰减。这解释了为什么深度电路和长演化时间会显著降低算法成功率。
六、从数学表征到噪声适应
实际建模时，建议按以下顺序推进：第一，根据硬件标定数据确定 \(T_1\)、\(T_2\)、门错误率和读出错误率；第二，把单比特门误差近似为去极化信道，把空闲等待近似为振幅阻尼加相位阻尼；第三，写出每个时刻的 Kraus 算符或等效的 Pauli 转移矩阵；第四，用密度矩阵模拟或蒙特卡洛轨迹方法传播噪声；第五，比较理想分布与噪声分布的保真度、总变差距离或期望值偏差，判断噪声是否已超出算法容忍阈值。
理解密度矩阵、Kraus 算符和三类典型信道，是后续噪声适应与鲁棒优化的前提。只有把噪声写成可计算的映射，才能定量分析它对量子态演化的影响，并据此选择纠错、缓解或鲁棒控制策略。