3. 密度矩阵与量子信道
在理想量子电路模型中，我们习惯用态矢量 |ψ⟩ 描述量子态，用酉矩阵 U 描述演化。但真实量子处理器从来不是孤立的：量子比特会与衬底声子、控制线、环境电磁场、读出谐振腔等发生耦合，导致能量弛豫、相位退相干、门操作误差。此时态矢量不再够用，必须引入密度矩阵 ρ 与量子信道这一套语言。本节要回答三个问题：为什么开放系统必须用密度矩阵？什么样的映射才是物理上合法的量子信道？Kraus 算符表示到底在说什么？
一、从态矢量到密度矩阵：混合态的必然性
纯态 |ψ⟩ 对应的密度矩阵是 ρ = |ψ⟩⟨ψ|，它是一个秩为 1 的投影算符，满足 ρ² = ρ、Tr(ρ) = 1、ρ ≥ 0。当系统以概率 pi 处于某个纯态 |ψi⟩ 时，整体状态是混合态：
ρ = Σi pi |ψi⟩⟨ψi|，其中 pi ≥ 0，Σi pi = 1。
密度矩阵的三个基本性质是：厄米性 ρ† = ρ、迹为 1、半正定 ρ ≥ 0。半正定性保证任意测量结果的概率非负。判断纯态还是混合态，最方便的指标是纯度 P = Tr(ρ²)：纯态 P = 1，混合态 P < 1。例如单比特完全混合态 ρ = I/2，纯度 P = 1/2。
密度矩阵还有一个关键优势：偏迹。若复合系统 AB 处于 ρAB，只关心 A 时，对 B 求偏迹 ρA = TrB(ρAB)。即使 ρAB 是纯态（如贝尔态），ρA 也可能是完全混合态——这正是纠缠的体现，也是退相干把纯态变成混合态的数学根源。
二、量子信道：开放系统的演化映射
封闭系统的演化是酉变换 ρ → UρU†。开放系统的演化则是一个更一般的映射 ℰ：ρ → ℰ(ρ)。物理上合法的量子信道必须满足两个条件：
- 保迹（Trace-Preserving, TP）：Tr(ℰ(ρ)) = 1，概率守恒，不能凭空丢失或产生概率。
- 完全正定（Completely Positive, CP）：不仅 ℰ 本身把正定算符映为正定算符，而且对任意维度的辅助系统 R，ℰ ⊗ IR 也保持正定性。只要求“正定”是不够的，转置映射就是正定但不完全正定的经典反例。
CPTP 映射是量子信道的严格定义。它保证了：无论信道与多大的环境纠缠，作用后密度矩阵仍然合法。这一条件在实验上对应“不会出现负概率”，在数学上则直接导向 Kraus 表示。
三、Kraus 算符表示及其物理意义
Kraus 定理指出：映射 ℰ 是 CPTP 的，当且仅当存在一组算符 {Kk}，使得
ℰ(ρ) = Σk Kk ρ Kk†，且 Σk Kk†Kk = I。
这组 Kk 称为 Kraus 算符，个数 k 称为 Kraus 秩。物理意义非常清晰：每个 Kk 对应环境发生的一种“记录”或“分支”。系统与环境的联合酉演化 USE 作用在 |ψ⟩S|0⟩E 上，展开为 Σk Kk|ψ⟩S|k⟩E。若环境状态 |k⟩E 正交且我们不去观测环境，对 E 求偏迹就得到上面的 Kraus 形式。因此 Kraus 算符是“环境被追踪掉之后”留在系统上的有效作用。
Kraus 表示不唯一：{Kk} 与 {Σj ukjKj}（u 为酉矩阵）给出同一信道。这说明 Kraus 分解的“分支”不是物理上可直接区分的，只有信道整体 ℰ 才是可观测量。
四、典型信道的 Kraus 表示与参数
把抽象定义落到量子电路噪声上，最常见的三个单比特信道都有明确的 Kraus 算符。设 γ ∈ [0,1] 为噪声强度参数：
| 信道 | 物理过程 | Kraus 算符 | 关键参数 |
|---|---|---|---|
| 比特翻转 | X 错误随机发生 | K₀ = √(1−p) I，K₁ = √p X | p 为翻转概率 |
| 相位翻转 | Z 错误随机发生 | K₀ = √(1−p) I，K₁ = √p Z | p 为相位错误概率 |
| 振幅阻尼 | |1⟩ 衰减到 |0⟩（T₁ 过程） | K₀ = [[1,0],[0,√(1−γ)]]，K₁ = [[0,√γ],[0,0]] | γ = 1−e−t/T₁ |
| 相位阻尼 | 相位信息丢失（T₂ 过程） | K₀ = √(1−λ/2) I，K₁ = √(λ/2) Z | λ = 1−e−t/T_φ |
以振幅阻尼为例，把 K₀、K₁ 代入 ℰ(ρ) = K₀ρK₀† + K₁ρK₁†，对 ρ = [[ρ₀₀, ρ₀₁],[ρ₁₀, ρ₁₁]] 计算可得：
ℰ(ρ) = [[ρ₀₀ + γ·ρ₁₁, √(1−γ)·ρ₀₁],
[√(1−γ)·ρ₁₀, (1−γ)·ρ₁₁]]
对角元：ρ₁₁ 以速率 γ 衰减到 ρ₀₀，体现能量弛豫
非对角元：乘以 √(1−γ)，体现相位相干性损失这个结果直接对应实验可观测的 T₁ 衰减：制备 |1⟩ 后测量激发态布居 P₁(t) = e−t/T₁，与 γ = 1−e−t/T₁ 完全一致。而相位阻尼只压缩非对角元、不改变布居，对应 Ramsey 干涉条纹对比度按 e−t/T₂ 衰减。可见 Kraus 算符不是抽象符号，它的参数直接由 T₁、T₂、门错误率等实验量标定。
五、在噪声建模中的实际用法
在量子电路噪声建模流程中，密度矩阵与 Kraus 表示承担两个角色。第一，把器件物理参数翻译成信道参数：从 T₁、T₂、单门/双门错误率、读出保真度出发，为每个量子比特和每类门选择信道模型并标定 γ、p、λ。第二，把信道插入电路仿真：在理想门之后追加 ℰ，用密度矩阵模拟器（如 Qiskit Aer 的 density_matrix 方法）传播 ρ，得到含噪输出态，再计算保真度、纠缠度、期望值偏差。
需要注意，Kraus 表示假设噪声是马尔可夫的、无记忆的。真实超导量子比特存在 1/f 磁通噪声、准粒子中毒、串扰等非马尔可夫效应，此时需要更一般的量子过程层析或非马尔可夫主方程。但在当前 NISQ 阶段，CPTP + Kraus 仍是误差分析最实用、最可计算的起点。