2. 量子比特与量子门基础回顾

噪声建模的第一步不是写代码,而是把“量子态如何被描述、门如何作用、测量如何读出”这三件事统一到同一套数学语言里。只有先固定表示法,后面讨论退相干、门误差、读出误差时,才能明确每个误差究竟作用在态矢量、密度矩阵还是测量算符上。本节按“态—门—测量”的顺序梳理这套语言,并给出可直接用于后续噪声建模的约定。

一、态矢量:纯态的表示与相位约定

单个量子比特的纯态写成二维复向量:

|ψ⟩ = α|0⟩ + β|1⟩,  α, β ∈ ℂ,  |α|² + |β|² = 1

计算基取 |0⟩ = (1, 0)ᵀ,|1⟩ = (0, 1)ᵀ。α、β 是概率幅,|α|² 与 |β|² 分别是测量得到 0 和 1 的概率。整体相位 e^{iγ} 不可观测,因此 |ψ⟩ 与 e^{iγ}|ψ⟩ 描述同一物理态;但相对相位(α 与 β 之间的相位差)可观测,它决定 X、Y 基下的干涉结果。这一点在噪声建模中很关键:退相位噪声破坏的正是相对相位,而不是整体相位。

把态矢量对应到布洛赫球更直观。定义极化矢量 r = (⟨X⟩, ⟨Y⟩, ⟨Z⟩),纯态满足 |r| = 1,位于球面。例如 |0⟩ 对应 r = (0, 0, 1),|+⟩ = (|0⟩+|1⟩)/√2 对应 r = (1, 0, 0)。后续把噪声看作对 r 的收缩与旋转,比直接操作复数幅更易理解。

二、密度矩阵:从纯态到混合态

真实器件中的量子比特几乎从不是纯态,而是纯态的统计混合。密度矩阵把“量子叠加”和“经典不确定性”统一起来:

ρ = Σ_i p_i |ψ_i⟩⟨ψ_i|,  p_i ≥ 0,  Σ_i p_i = 1

纯态是特例,此时 ρ = |ψ⟩⟨ψ|,且 ρ² = ρ、Tr(ρ²) = 1。混合态满足 Tr(ρ²) < 1。常用纯度 P = Tr(ρ²) 和冯·诺依曼熵 S(ρ) = −Tr(ρ log₂ρ) 来量化“混”的程度。单比特密度矩阵可参数化为:

ρ = (I + r_x X + r_y Y + r_z Z) / 2

其中 r = (r_x, r_y, r_z) 是布洛赫矢量,|r| ≤ 1。纯态 |r| = 1,完全混合态 r = 0 对应 ρ = I/2。这个参数化是噪声建模的核心工具:退相干使 |r| 缩短,相干旋转使 r 绕某轴转动,振幅阻尼则把 r 拉向 −z 或 +z 方向。测量算符的期望值也可直接写成 ⟨O⟩ = Tr(ρO)。

三、量子门:幺正演化与矩阵表示

理想量子门是幺正矩阵 U,满足 U†U = I。作用在态矢量上是 |ψ'⟩ = U|ψ⟩,作用在密度矩阵上是 ρ' = UρU†。单比特常用门及其矩阵为:

X = [[0,1],[1,0]]      Y = [[0,-i],[i,0]]     Z = [[1,0],[0,-1]]
H = (1/√2)[[1,1],[1,-1]]
S = [[1,0],[0,i]]      T = [[1,0],[0,e^{iπ/4}]]
Rx(θ) = cos(θ/2)I − i sin(θ/2)X

双比特门以 CNOT 为例,控制位为高时翻转目标位:

CNOT = [[1,0,0,0],
        [0,1,0,0],
        [0,0,0,1],
        [0,0,1,0]]

矩阵基序按 |00⟩, |01⟩, |10⟩, |11⟩ 排列。任何单比特门都可写成绕布洛赫球某轴的旋转,例如 Rx(θ) 使 r 绕 x 轴转 θ。噪声建模中,门误差常建模为“理想 U 之后叠加一个小的错误旋转”,即 U_actual = U_error · U_ideal,其中 U_error 由过旋转角、串扰或控制脉冲失真引起。

四、测量算符:从可观测量到读出

可观测量的谱分解为 O = Σ_m m P_m,P_m 是投影算符。对态 ρ 测量得到结果 m 的概率为 p(m) = Tr(P_m ρ),测量后态坍缩为 P_m ρ P_m / p(m)。计算基测量对应 P_0 = |0⟩⟨0|、P_1 = |1⟩⟨1|,于是:

p(0) = Tr(|0⟩⟨0| ρ) = ρ₀₀
p(1) = Tr(|1⟩⟨1| ρ) = ρ₁₁

也就是说,计算基测量直接读出密度矩阵对角元。若在 X 基测量,则先施加 H 门再计算基测量,等价于测量算符 X。期望值统一为 ⟨O⟩ = Tr(ρO)。真实读出存在误差:本应读 0 却读成 1 的概率记为 p(1|0),反之记为 p(0|1),这可用读出混淆矩阵描述,是后续读出误差建模的入口。

量子比特数学语言:态 → 门 → 测量 态矢量 |ψ⟩(纯态) |ψ⟩ = α|0⟩ + β|1⟩ |α|² + |β|² = 1 整体相位不可观测 相对相位可观测 布洛赫球:|r| = 1 r = (⟨X⟩, ⟨Y⟩, ⟨Z⟩) 密度矩阵 ρ(混合态) ρ = Σ p_i |ψ_i⟩⟨ψ_i| Tr(ρ) = 1, ρ = ρ† ≥ 0 ρ = (I + r·σ)/2 纯度 P = Tr(ρ²) ≤ 1 |r| < 1 表示混合 r = 0 → ρ = I/2 测量算符 O O = Σ m P_m p(m) = Tr(P_m ρ) 坍缩:P_m ρ P_m / p(m) 期望:⟨O⟩ = Tr(ρO) 计算基:p(0)=ρ₀₀ 读出误差 p(1|0), p(0|1) 统计混合 读出 量子门:幺正演化 U(U†U = I) 态矢量:|ψ'⟩ = U|ψ⟩ 密度矩阵:ρ' = UρU† 布洛赫球:r 绕轴旋转 X, Y, Z, H, S, T, Rx(θ), CNOT … 噪声建模的注入点(本节语言对应) 退相位:|r| 横向收缩 ρ 非对角元衰减 门误差:U_actual = U_error · U_ideal 读出误差:混淆矩阵 p(1|0), p(0|1)

五、把四件工具串成噪声建模的入口

现在把上述语言落到一个可操作的小例子上。设某超导量子比特经 Ramsey 实验后测得布洛赫矢量 r = (0.62, 0.18, 0.75),则密度矩阵为:

ρ = (I + 0.62X + 0.18Y + 0.75Z)/2
  = [[0.875, 0.31 − 0.09i],
     [0.31 + 0.09i, 0.125]]

对角元 0.875 与 0.125 给出计算基测量概率;非对角元模长 |ρ₀₁| ≈ 0.323 反映相干性。若理想纯态应有 |r| = 1,而此处 |r| = √(0.62² + 0.18² + 0.75²) ≈ 0.99,说明相干性略有损失。把这一测量重复在不同延迟时间 t 上,拟合 |ρ₀₁(t)| 的指数衰减即可提取 T₂;拟合 ρ₁₁(t) 的弛豫即可提取 T₁。这就是从“态的描述”过渡到“噪声参数提取”的最小闭环。