📚 课程目录第2章量子比特与密度矩阵基础

2. 量子比特与密度矩阵基础

在讨论噪声通道之前，必须先明确"被噪声作用的对象"用什么数学对象描述。量子比特（qubit）是量子信息的基本单元，而密度矩阵是描述量子态最通用的语言——它既能表示纯态，也能表示混态，还能自然处理噪声导致的经典不确定性。本节围绕态矢量、密度矩阵、纯态与混态、约化密度矩阵以及 Bloch 球表示展开，为后续噪声通道的输入输出建立统一的符号体系。

一、量子比特的态矢量表示

一个量子比特是二维复希尔伯特空间中的单位矢量。取计算基 \(|0\rangle=\begin{pmatrix}1\\0\end{pmatrix}\)、\(|1\rangle=\begin{pmatrix}0\\1\end{pmatrix}\)，任意纯态可写为

其中 \(\alpha,\beta\) 为复振幅，\(|\alpha|^2\) 与 \(|\beta|^2\) 分别是测量得到 0 和 1 的概率。由于全局相位不可观测，可令 \(\alpha=\cos(\theta/2)\)、\(\beta=e^{i\varphi}\sin(\theta/2)\)，于是纯态由两个实参数 \((\theta,\varphi)\) 唯一确定。这正是 Bloch 球坐标的来源：\(\theta\in[0,\pi]\) 为极角，\(\varphi\in[0,2\pi)\) 为方位角。

态矢量表示简洁，但只适用于纯态。一旦量子比特与外界发生纠缠或经历退相干，它就不再能用单个 \(|\psi\rangle\) 描述，必须升级到密度矩阵。

二、密度矩阵：纯态与混态的统一描述

密度矩阵 \(\rho\) 是作用在希尔伯特空间上的迹为 1 的半正定厄米算符：

\(\rho=\rho^\dagger,\quad \mathrm{Tr}(\rho)=1,\quad \rho\ge 0\)

纯态 \(|\psi\rangle\) 的密度矩阵为 \(\rho=|\psi\rangle\langle\psi|\)，它是秩 1 的投影算符，满足 \(\rho^2=\rho\)，等价判据是 \(\mathrm{Tr}(\rho^2)=1\)。混态则是若干纯态的经典概率混合：

\(\rho=\sum_i p_i|\psi_i\rangle\langle\psi_i|,\quad p_i\ge 0,\ \sum_i p_i=1\)

此时 \(\mathrm{Tr}(\rho^2)<1\)，该量称为纯度。纯度越低，说明态越"混合"，通常意味着噪声越强。判断纯态与混态最直接的方法就是计算 \(\mathrm{Tr}(\rho^2)\)：等于 1 为纯态，小于 1 为混态。

在计算基下，任意单量子比特密度矩阵可参数化为

ρ = 1/2 * ( I + x·σx + y·σy + z·σz )
  = 1/2 * [ [1+z,   x-iy],
            [x+iy,  1-z ] ]

其中 \(x,y,z\) 为实数，\((x,y,z)\) 就是 Bloch 矢量 \(\vec r\)，且 \(|\vec r|\le 1\)。这一参数化把抽象的算符与三维几何直接对应起来。

三、约化密度矩阵：子系统视角

当系统由多个量子比特组成时，我们往往只关心其中一个。设复合系统处于 \(\rho_{AB}\)，则子系统 A 的约化密度矩阵由对 B 求偏迹得到：

\(\rho_A=\mathrm{Tr}_B(\rho_{AB})\)

偏迹的物理含义是"忽略 B 的全部信息"。关键结论是：即使 \(\rho_{AB}\) 是纯态，只要 A 与 B 存在纠缠，\(\rho_A\) 就可能是混态。例如 Bell 态 \(|\Phi^+\rangle=(|00\rangle+|11\rangle)/\sqrt2\) 是纯态，但对 B 求偏迹后得到 \(\rho_A=I/2\)，即完全混态。这说明"混态"并不一定来自经典噪声，纠缠本身就会让局部观测呈现混合性。理解这一点对分析噪声通道至关重要：噪声既可能来自环境引入的经典不确定性，也可能来自与环境的纠缠。

四、Bloch 球表示

Bloch 球把单量子比特态映射到单位球内：球面上的点（\(|\vec r|=1\)）对应纯态，球内的点（\(|\vec r|<1\)）对应混态，球心 \(\vec r=0\) 对应最大混态 \(I/2\)。常用参考点如下：

\(|0\rangle\)：北极，\(\vec r=(0,0,1)\)
\(|1\rangle\)：南极，\(\vec r=(0,0,-1)\)
\(|+\rangle=(|0\rangle+|1\rangle)/\sqrt2\)：x 轴正向，\(\vec r=(1,0,0)\)
\(|+i\rangle=(|0\rangle+i|1\rangle)/\sqrt2\)：y 轴正向，\(\vec r=(0,1,0)\)

Bloch 球的最大价值在于把噪声通道的作用可视化。以退极化通道为例，它把 Bloch 矢量按比例 \(1-p\) 收缩：\(\vec r\mapsto(1-p)\vec r\)。当 \(p=0\) 时态不变，\(p=1\) 时所有态都被压到球心，成为完全混态。比特翻转通道则把 Bloch 矢量绕 x 轴旋转，相位翻转通道绕 z 轴旋转。因此，观察 Bloch 矢量在通道作用前后的长度与方向变化，就能直观判断噪声类型与强度。

五、从密度矩阵到噪声通道的接口

噪声通道在数学上是一个完全正定保迹映射（CPTP map）\(\mathcal{E}\)，它把输入密度矩阵 \(\rho_{\text{in}}\) 变为输出 \(\rho_{\text{out}}=\mathcal{E}(\rho_{\text{in}})\)。用算符和表示（Kraus 表示）可写为

\(\mathcal{E}(\rho)=\sum_k E_k\rho E_k^\dagger,\quad \sum_k E_k^\dagger E_k=I\)

Kraus 算符 \(E_k\) 刻画了噪声的具体作用方式。例如退极化通道的 Kraus 算符为 \(E_0=\sqrt{1-3p/4}\,I\)、\(E_1=\sqrt{p/4}\,\sigma_x\)、\(E_2=\sqrt{p/4}\,\sigma_y\)、\(E_3=\sqrt{p/4}\,\sigma_z\)。把它作用到 \(\rho=\frac12(I+\vec r\cdot\vec\sigma)\) 上，可直接算出 \(\vec r\mapsto(1-p)\vec r\)，与 Bloch 球的几何图像完全一致。

因此，分析噪声通道的典型流程是：写出输入态的密度矩阵（或 Bloch 矢量）→ 确定通道的 Kraus 算符或 Bloch 映射规则 → 计算输出密度矩阵 → 用纯度、保真度或 Bloch 矢量长度判断噪声影响。保真度 \(F(\rho_{\text{in}},\rho_{\text{out}})=\left(\mathrm{Tr}\sqrt{\sqrt{\rho_{\text{in}}}\,\rho_{\text{out}}\sqrt{\rho_{\text{in}}}}\right)^2\) 是常用的量化指标，纯态输入时可简化为 \(\langle\psi|\rho_{\text{out}}|\psi\rangle\)。

实用判据：给定单量子比特输入 \(|0\rangle\)（Bloch 矢量 (0,0,1)），若通道输出 Bloch 矢量为 (0,0,1−p)，则说明该通道等效于以概率 p 发生比特翻转或退极化；若输出为 (0,0,1) 但纯度下降，则需检查是否存在相位噪声。先看方向、再看长度，是快速定位噪声类型的有效顺序。

本节建立的符号体系——态矢量、密度矩阵、纯态/混态判据、约化密度矩阵、Bloch 球与 Kraus 算符——将贯穿后续所有噪声通道的建模与分析。掌握这些工具后，无论面对退极化、比特翻转、相位翻转还是振幅阻尼通道，都可以用同一套语言写出输入、作用与输出，并用纯度与保真度定量评估噪声强度。