2. 量子比特能级结构与泄漏通道分析
在超导量子计算中,Transmon 量子比特之所以能成为主流方案,关键在于它通过增大约瑟夫森结与分流电容的比值(即增大 EJ/EC)来抑制电荷噪声。但这一设计也带来了一个副作用:能级间距不再等距,且高能级(尤其是 |2⟩ 态)与计算子空间 {|0⟩, |1⟩} 之间的耦合始终存在。脉冲整形与泄漏抑制的所有物理依据,都建立在对这套能级结构和泄漏通道的准确理解之上。
一、Transmon 的能级结构:非谐性与计算子空间
把 Transmon 近似为弱非谐量子谐振子,其能级可写为
E_n ≈ -E_J + √(8 E_J E_C)·(n + 1/2) - (E_C/12)·(6n² + 6n + 3)
更常用的是用频率表示相邻能级跃迁:ω01 为基态到第一激发态的频率,ω12 为第一激发态到第二激发态的频率,二者之差即非谐性 α = ω01 − ω12。典型 Transmon 的 EC/h 在 200~350 MHz 量级,对应 α/h 约 200~300 MHz,而 ω01/2π 通常在 4~6 GHz。这意味着 |1⟩→|2⟩ 的跃迁频率比 |0⟩→|1⟩ 低约 200~300 MHz,正是这个“频率错位”让选择性驱动成为可能,也正是它让宽带脉冲容易同时激发两个跃迁。
计算子空间只包含 |0⟩ 和 |1⟩,|2⟩ 及以上统称泄漏态。需要强调的是,|2⟩ 并非“无关的高能级”:它通过驱动场与 |1⟩ 直接耦合,且其衰减通道 |2⟩→|1⟩ 会把泄漏能量重新注入计算子空间,形成相干与非相干混合的误差来源。因此泄漏分析必须同时关注“如何进入 |2⟩”和“进入后如何演化”两个方向。
二、泄漏进入 |2⟩ 的三类主要通道
通道一:驱动引起的多光子/非共振跃迁。 当施加频率为 ωd 的驱动时,若脉冲含有 ωd + ω12 附近的频率分量,或脉冲包络过快导致频谱展宽覆盖 ω12,就会在 |1⟩→|2⟩ 上产生非共振激发。更隐蔽的是,强驱动下 |0⟩→|2⟩ 的双光子过程也会被增强。以高斯脉冲为例,若总门时间 tg 过短,其频谱半宽约 0.44/tg;当 tg = 10 ns 时半宽约 44 MHz,尚不足以覆盖 200 MHz 的非谐性,但若叠加方波边沿或 DRAG 参数失配,高频分量会显著抬升。工程上常用 DRAG 校正:在正交分量上叠加与包络导数成正比的修正项,抵消 |1⟩→|2⟩ 路径的相位积累。
通道二:能级交叉与频率碰撞。 在多比特芯片中,每个 Transmon 的 ω01 会随磁通偏置或邻近比特状态变化。当两个比特的 |1⟩→|2⟩ 跃迁频率接近时(即 ω12(i) ≈ ω01(j)),会发生能级交叉或避免交叉,导致激发在比特间转移并伴随泄漏。这类通道对频率排布(frequency collision)极其敏感,通常在芯片设计阶段通过非谐性差异化来规避,但在可调耦合器架构中仍需在脉冲层面检查瞬时频率轨迹。
通道三:噪声诱导跃迁。 热噪声、磁通噪声和电荷噪声都会以一定速率把布居从计算子空间“踢”到 |2⟩。热光子数 nth 引起的 |1⟩→|2⟩ 激发速率近似正比于 nth·κ,其中 κ 为谐振腔或环境的耦合速率。磁通噪声则通过调制 EJ 改变瞬时频率,在偏置点附近形成等效的纵向驱动。这类通道的特点是速率较低但持续存在,往往在长时序或多脉冲序列中累积成可观测的泄漏。
三、从能级结构到脉冲设计约束
把上述通道翻译成脉冲设计语言,可以得到几条直接的工程约束。第一,脉冲带宽必须与 α 保持足够间隔:若门时间 tg 对应的频谱主瓣宽度接近 α,则 |1⟩→|2⟩ 的非共振激发不可忽略。经验上要求 1/tg ≪ α,例如 α/2π = 250 MHz 时,tg 取 20~30 ns 可使主瓣远小于非谐性。第二,任何幅度调制都会引入边带,方波边沿的 sinc 旁瓣衰减慢,因此实际脉冲应使用平滑包络(高斯、平方正弦等)并配合 DRAG。第三,多比特场景下必须检查驱动频率与邻近比特 ω12 的间隔,避免在频域上“撞车”。
以一个可复现的教学示例说明判断过程(以下数值为合理假设,非实测):设某 Transmon 的 ω01/2π = 5.00 GHz,α/2π = 250 MHz,则 ω12/2π = 4.75 GHz。若采用 16 ns 高斯脉冲实现 π 门,其频谱半宽约 0.44/16 ns ≈ 27.5 MHz,主瓣远小于 250 MHz,非共振泄漏在理想情况下很小;但若把脉冲截断为 8 ns 方波,频谱半宽升至约 55 MHz 且旁瓣拖尾严重,此时需在仿真中观察 |2⟩ 布居是否从 10⁻⁴ 量级抬升到 10⁻³ 量级。判断依据不是单看波形,而是看演化结束时 |2⟩ 的布居 P₂ 以及泄漏对门保真度的贡献。
四、泄漏通道的排查顺序
面对一个实测泄漏偏高的比特,合理的排查顺序是:先确认单比特门参数(tg、DRAG 系数 β)是否匹配当前 α,排除通道一;再用频谱扫描检查驱动频率附近是否存在邻近比特的 ω12 或耦合器模式,排除通道二;最后通过变温或改变偏置点观察泄漏速率是否随 nth、磁通噪声变化,定位通道三。这个顺序的依据是三类通道的可控性和出现频率不同:驱动参数最易调、频率碰撞次之、噪声最难改。