📚 课程目录第3章量子系统噪声的数学表征

3. 量子系统噪声的数学表征

在量子测量误差分析中，仅用态矢量 |ψ⟩ 描述系统是不够的。真实量子系统与环境耦合后，会从纯态退化为混合态，测量结果表现为概率分布而非确定值。要定量刻画这种“噪声如何改变量子态”，需要一套完整的数学语言：密度矩阵描述状态本身，量子信道描述状态演化，Kraus算符给出信道的可计算形式，主方程描述连续时间演化，而完全正映射则规定了什么样的演化在物理上是合法的。

密度矩阵：混合态的通用描述

密度矩阵 ρ 是描述量子态最通用的工具。纯态 ρ=|ψ⟩⟨ψ| 满足 ρ²=ρ 且 Tr(ρ²)=1；混合态则 Tr(ρ²)<1。以单量子比特为例，任意密度矩阵可写成布洛赫球形式：

ρ = (I + r·σ)/2,  r = (rx, ry, rz),  |r| ≤ 1
纯态: |r| = 1    混合态: |r| < 1    完全混态: r = 0

布洛赫矢量 r 的长度直接反映态纯度。噪声的作用往往表现为 r 的收缩：退相干使横向分量 rx、ry 衰减，退相位使 rz 保持而横向分量衰减，振幅阻尼则同时改变 rz 并收缩横向分量。测量误差分析中，我们常通过量子态层析重建 ρ，再观察 r 的变化来判断噪声类型与强度。

量子信道：状态演化的映射

量子信道是密度矩阵到密度矩阵的映射 ℰ: ρ → ℰ(ρ)，要求满足两个条件：迹保持 Tr(ℰ(ρ))=1，保证概率归一；完全正性，即对任意辅助系统维度都保持正定。完全正性排除了转置这类“正但非完全正”的映射，因为转置会与纠缠态结合后产生负本征值，物理上不可实现。

常见的单比特信道包括：

比特翻转信道：以概率 p 施加 X 门，ℰ(ρ)=(1−p)ρ+pXρX；
相位翻转信道：以概率 p 施加 Z 门，ℰ(ρ)=(1−p)ρ+pZρZ；
退极化信道：ℰ(ρ)=(1−p)ρ+(p/3)(XρX+YρY+ZρZ)，布洛赫矢量整体收缩为 (1−4p/3)r；
振幅阻尼信道：描述能量耗散，Kraus 算符为 E₀=diag(1,√(1−γ))，E₁=√γ|0⟩⟨1|。

Kraus算符：信道的可计算表示

任何完全正、迹保持映射都可写成 Kraus 表示：

ℰ(ρ) = Σ_k E_k ρ E_k†,   约束条件: Σ_k E_k† E_k = I

Kraus 算符 {E_k} 的个数称为信道的 Kraus 秩。以振幅阻尼为例，取 γ=0.1：

E0 = [[1, 0], [0, sqrt(0.9)]] ≈ [[1,0],[0,0.9487]]
E1 = [[0, sqrt(0.1)], [0, 0]] ≈ [[0,0.3162],[0,0]]

验证: E0†E0 + E1†E1 = [[1,0],[0,0.9]] + [[0,0],[0,0.1]] = I ✓

若输入态为 |1⟩⟨1|，则 ℰ(|1⟩⟨1|)=E₀|1⟩⟨1|E₀†+E₁|1⟩⟨1|E₁†=0.9|1⟩⟨1|+0.1|0⟩⟨0|，即系统以 10% 概率衰减到 |0⟩。这正是量子测量中“读出错误”或“能量弛豫”的数学来源。Kraus 算符的实用价值在于：只要写出 E_k，就能直接计算任意输入态的演化结果，无需解微分方程。

主方程：连续时间噪声演化

当噪声连续作用时，用 Lindblad 主方程描述：

dρ/dt = −(i/ħ)[H, ρ] + Σ_k (L_k ρ L_k† − ½{L_k†L_k, ρ})

第一项是幺正演化，第二项是耗散项。L_k 称为跳跃算符。对振幅阻尼，取 L=√γ σ₋，可得 dρ₁₁/dt=−γρ₁₁，即激发态布居按 e^{−γt} 衰减。对退相位，取 L=√(γ_φ/2) σ_z，则非对角元按 e^{−γ_φ t} 衰减，而布居不变。主方程把“噪声强度”与“时间常数”直接联系起来，是分析测量误差随等待时间增长的核心工具。

完全正映射的物理意义与判据

完全正性（CP）是量子信道合法性的底线。一个映射 ℰ 完全正，当且仅当其 Choi 矩阵 J(ℰ)=Σ_{ij}|i⟩⟨j|⊗ℰ(|i⟩⟨j|) 半正定。Choi 矩阵的秩等于 Kraus 秩，其本征值谱反映信道的“噪声通道结构”。在实验上，通过量子过程层析得到 ℰ 的完整作用，再构造 Choi 矩阵并检查最小本征值：若为负，说明数据拟合出的映射不满足完全正性，通常源于测量误差、态制备错误或模型选择不当。

实用判断流程：先由过程层析得到 ℰ 的矩阵表示 → 构造 Choi 矩阵 J(ℰ) → 检查 J(ℰ)≥0 且 Tr(ℰ(ρ))=1 → 若通过，再做 Kraus 分解得到 {E_k} → 用 {E_k} 预测任意输入态的测量分布，并与实验对比。

从数学表征到测量误差分析

在实际量子测量实验中，这套数学工具的使用路径通常是：先通过量子过程层析获得信道的 Choi 矩阵，检查其半正定性；若通过，则提取 Kraus 算符，建立可计算的噪声模型；再用主方程拟合噪声随时间的演化曲线，提取 γ、γ_φ 等参数；最后用该模型预测不同输入态下的测量结果分布，与实验直方图对比。若预测与实验偏差较大，说明噪声模型不完整，需要引入更多跳跃算符或考虑非马尔可夫效应。

需要特别注意的是，Kraus 表示不唯一：{E_k} 与 {U_{kj}E_j}（U 为任意幺正矩阵）给出同一信道。因此从实验数据提取 Kraus 算符时，应固定规范或直接使用 Choi 矩阵进行比较，避免因表示不同而误判噪声类型。此外，完全正性只是必要条件，实际信道还需满足迹保持；若过程层析得到的映射迹不保持，通常意味着测量存在系统偏差或态制备不完整。

常见误区：将“正映射”等同于“完全正映射”。转置映射是正的但不是完全正的，若在分析中误用，会得到负概率或负密度矩阵。判断时务必构造 Choi 矩阵并检查最小本征值，而不是只验证单态正定性。

总结来说，密度矩阵给出状态描述，量子信道给出演化框架，Kraus 算符给出可计算形式，主方程给出连续时间动力学，完全正映射给出物理合法性判据。五者构成一个从静态到动态、从抽象到可计算的完整体系，是后续噪声参数估计、误差缓解和测量不确定度分析的数学基础。