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⁻¹。

把这些指标组织成一张“噪声参数表”,是噪声感知编译的起点。下表给出一个示例(数值为合理假设,非实测):

参数符号示例值主要影响
能量弛豫时间T160 μs|1⟩ 衰减、振幅阻尼
相位相干时间T280 μs叠加态相位丢失
单比特门错误率ε₁3×10⁻⁴门序列累积错误
双比特门错误率ε₂8×10⁻³纠缠门主导错误
读出错误率ε_ro2×10⁻²测量结果偏差
ZZ 耦合强度ζ200 kHz条件相位、串扰

有了这些参数,编译器的成本函数才能从“门数量最少”升级为“预期错误率最低”。例如,一个双比特门错误率 8×10⁻³,而一个单比特门仅 3×10⁻⁴,那么把双比特门拆成多个单比特门在某些情况下反而更优——这正是噪声感知编译与纯逻辑综合的关键区别。

量子噪声分类与关键指标关系图 噪声来源 比特翻转 (X) 相位翻转 (Z) 振幅阻尼 热弛豫 退相干 串扰 (ZZ) 读出错误 物理机制 能量交换 相位随机化 环境耦合 控制脉冲泄漏 谐振腔耦合 测量链噪声 态制备不完美 关键指标与编译影响 T1:能量弛豫时间 → 缩短 |1⟩ 驻留,减少延迟 T2 / T2*:相位相干时间 → 插入回波、动态解耦 ε₁:单比特门错误率 → 控制单比特门数量 ε₂:双比特门错误率 → 优化纠缠门拓扑与数量 ε_ro:读出错误率 → 读出校准、选择测量比特 ζ:ZZ 耦合强度 → 避免同时驱动强耦合邻居

六、测量与判读的实操要点

要获得可用的噪声参数,建议按以下顺序执行。第一,先测 T1:准备 |1⟩ 态,扫描延迟 τ,读出 P1,拟合指数衰减。若拟合残差偏大,检查是否受到热激发或读出串扰影响。第二,测 Ramsey 与回波:比较 T2* 与 T2,若 T2* 远小于 T2,说明低频噪声显著,编译时应考虑插入回波脉冲。第三,用 RB 测单比特门错误率,用 XEB 或双比特 RB 测双比特门错误率,注意区分“每门错误率”与“每 Clifford 错误率”。第四,测读出混淆矩阵,用校准矩阵修正测量结果。第五,测 ZZ 耦合:让一个比特处于 |0⟩ 与 |1⟩,分别测邻居的 Ramsey 频率,频率差即 ZZ 强度。

判读原则:若 T2 ≈ 2T1,能量弛豫主导,优先缩短电路深度;若 T2 ≪ 2T1,相位噪声主导,优先插入动态解耦或回波;若双比特门错误率远高于单比特门,优先优化纠缠门拓扑与数量;若读出错误率超过 5%,先做读出校准再谈编译优化。