2. 量子噪声基础与表征
量子比特不是理想二能级系统。它与环境耦合、被控制脉冲扰动、被相邻比特影响,最终表现为计算结果的概率性偏差。噪声感知编译的第一步,不是急着优化电路,而是把“噪声是什么、从哪里来、用什么指标刻画”说清楚。本节围绕比特翻转、相位翻转、退相干、振幅阻尼、热弛豫、串扰、读出错误七类典型噪声,以及T1、T2、门错误率等关键指标展开,并给出可操作的测量与判读方法。
一、从“错误作用在哪个自由度”分类
单比特纯态可写成 |ψ⟩ = α|0⟩ + β|1⟩。比特翻转(bit flip)交换 |0⟩ 与 |1⟩ 的振幅,对应 Pauli-X 型错误;相位翻转(phase flip)保持 |0⟩ 不变、给 |1⟩ 加一个负号,对应 Pauli-Z 型错误;两者同时发生即 Y 型错误。若把布洛赫球作为参照系,X 错误使态矢量绕 x 轴旋转,Z 错误绕 z 轴旋转,Y 错误绕 y 轴旋转。这一参照约定在全文与图示中保持一致:z 轴对应 |0⟩/|1⟩ 的能量轴,x、y 轴对应叠加态的相位方向。
需要强调的是,真实噪声很少是单一 Pauli 错误,而是连续旋转加随机相位。但量子纠错与噪声感知编译常把噪声“Pauli 化”近似,因为任意单比特错误都能分解为 I、X、Y、Z 的线性组合,这样便于用错误率参数驱动编译决策。
二、退相干、振幅阻尼与热弛豫
退相干(decoherence)描述量子叠加相位信息的丢失。它分两类:能量弛豫使 |1⟩ 衰减到 |0⟩,相位弛豫使叠加态的相对相位随机化。振幅阻尼(amplitude damping)是能量弛豫的典型信道模型,其 Kraus 算符为
E0 = [[1, 0], [0, sqrt(1-γ)]]
E1 = [[0, sqrt(γ)], [0, 0]]
其中 γ = 1 − exp(−t/T1) 表示在时间 t 内发生衰减的概率。当 γ 很小时,E1 近似把 |1⟩ 映射到 |0⟩,这正是比特翻转错误的一个来源。热弛豫(thermal relaxation)是振幅阻尼的推广:环境并非零温,|0⟩ 也会以一定概率被激发到 |1⟩。设激发速率 Γ↑、衰减速率 Γ↓,则
1/T1 = Γ↓ + Γ↑
p_excited = Γ↑ / (Γ↓ + Γ↑)
在超导量子比特中,典型工作温度约 10 mK,热激发概率通常很小,但并非为零。若 T1 = 60 μs、门时长 40 ns,则单次门时间内衰减概率约为 1 − exp(−0.04/60) ≈ 6.7×10⁻⁴,这是评估门错误率下限的重要量级。
三、T1、T2 与 T2* 的区别
T1 是能量弛豫时间,通过反转恢复实验测量:先把比特激发到 |1⟩,等待可变延迟 τ,再读出 |1⟩ 的剩余概率 P1(τ),拟合 P1(τ) = A·exp(−τ/T1) + B。T2 是相位相干时间,通过 Ramsey 干涉测量:施加两个 π/2 脉冲,中间插入延迟 τ,扫描相位后拟合振荡包络 exp(−τ/T2)。若存在低频噪声导致频率抖动,则观测到的是 T2*(去相位时间),它通常短于 T2;用 Hahn 回波(插入一个 π 脉冲)可部分抵消低频噪声,测得更接近本征的 T2。
三者关系可近似写为 1/T2 = 1/(2T1) + 1/Tφ,其中 Tφ 是纯相位弛豫时间。当 T2 明显小于 2T1 时,说明相位噪声占主导,编译时应优先考虑动态解耦或回波插入;当 T2 接近 2T1 时,能量弛豫是主要限制,此时缩短电路深度、减少 |1⟩ 态驻留时间更有效。
四、串扰与读出错误
串扰(crosstalk)指对一个比特施加操作时,相邻比特受到非预期影响。常见形式包括:控制线间的微波泄漏导致 ZZ 耦合,使相邻比特频率随对方状态移动;同时驱动多个比特时的残余耦合;以及读出谐振腔之间的相互干扰。ZZ 耦合会引入条件相位错误,表现为相邻比特的 Ramsey 频率随对方状态变化,典型量级在几十 kHz 到几 MHz。编译时可通过避免同时激发强耦合邻居、插入补偿脉冲或选择低串扰的比特对来缓解。
读出错误(readout error)发生在测量阶段:本应是 |0⟩ 却判为 |1⟩,或反之。它用 2×2 混淆矩阵描述:
P(0|0)=0.985 P(1|0)=0.015
P(0|1)=0.030 P(1|1)=0.970
读出错误不改变量子态本身,但会污染统计结果。可通过读出校准矩阵求逆做修正,或在编译时把测量操作安排在错误率较低的比特上。
五、门错误率与噪声参数汇总
门错误率(gate error rate)是编译最直接的成本输入。单比特门错误率通常用随机基准测试(RB)测量:施加随机 Clifford 序列,拟合存活概率随序列长度的指数衰减,得到每门平均错误率。双比特门错误率更高,常用交叉熵基准测试(XEB)或双比特 RB 测量。典型超导器件中,单比特门错误率约 10⁻⁴~10⁻³,双比特门约 10⁻³~10⁻²,读出错误约 10⁻²~10⁻¹。
把这些指标组织成一张“噪声参数表”,是噪声感知编译的起点。下表给出一个示例(数值为合理假设,非实测):
| 参数 | 符号 | 示例值 | 主要影响 |
|---|---|---|---|
| 能量弛豫时间 | T1 | 60 μs | |1⟩ 衰减、振幅阻尼 |
| 相位相干时间 | T2 | 80 μs | 叠加态相位丢失 |
| 单比特门错误率 | ε₁ | 3×10⁻⁴ | 门序列累积错误 |
| 双比特门错误率 | ε₂ | 8×10⁻³ | 纠缠门主导错误 |
| 读出错误率 | ε_ro | 2×10⁻² | 测量结果偏差 |
| ZZ 耦合强度 | ζ | 200 kHz | 条件相位、串扰 |
有了这些参数,编译器的成本函数才能从“门数量最少”升级为“预期错误率最低”。例如,一个双比特门错误率 8×10⁻³,而一个单比特门仅 3×10⁻⁴,那么把双比特门拆成多个单比特门在某些情况下反而更优——这正是噪声感知编译与纯逻辑综合的关键区别。
六、测量与判读的实操要点
要获得可用的噪声参数,建议按以下顺序执行。第一,先测 T1:准备 |1⟩ 态,扫描延迟 τ,读出 P1,拟合指数衰减。若拟合残差偏大,检查是否受到热激发或读出串扰影响。第二,测 Ramsey 与回波:比较 T2* 与 T2,若 T2* 远小于 T2,说明低频噪声显著,编译时应考虑插入回波脉冲。第三,用 RB 测单比特门错误率,用 XEB 或双比特 RB 测双比特门错误率,注意区分“每门错误率”与“每 Clifford 错误率”。第四,测读出混淆矩阵,用校准矩阵修正测量结果。第五,测 ZZ 耦合:让一个比特处于 |0⟩ 与 |1⟩,分别测邻居的 Ramsey 频率,频率差即 ZZ 强度。