2. 超导量子比特的核心物理基础

要理解超导量子比特为什么能“同时处于0和1”,必须回到两个最基本的电路元件:LC谐振电路和约瑟夫森结。前者提供了量子比特的“频率骨架”,后者提供了让能级变得非线性的关键机制。本节从这两个对象出发,逐步推导出电荷量子比特、磁通量子比特和Transmon三种主流结构,并说明它们的参数如何决定量子特性。

一、LC谐振电路:线性但不可用作量子比特

一个理想LC谐振电路由电感L和电容C并联构成,其谐振角频率为 ω = 1/√(LC),对应频率 f = 1/(2π√(LC))。在量子化之后,LC电路的能量本征态是等间距的谐振子能级:En = ħω(n + 1/2)。这意味着相邻能级间隔恒为 ħω,从|0⟩激发到|1⟩所需的能量与从|1⟩激发到|2⟩完全相同。

问题正出在这里。量子计算要求我们把量子比特限制在|0⟩和|1⟩两个能级上,但如果|1⟩→|2⟩的跃迁频率与|0⟩→|1⟩相同,任何用于驱动比特的微波脉冲都会同时激发更高能级,导致信息泄漏。因此,纯线性LC电路无法作为量子比特,必须引入非线性元件使能级间隔不再相等。

二、约瑟夫森结:非线性的来源

约瑟夫森结由两层超导体夹一层极薄绝缘势垒构成(S-I-S结构)。当绝缘层足够薄时,库珀对可以隧穿通过,形成无耗散的超导电流。其两个核心关系为:

  • 电流-相位关系:I = Ic sin φ,其中Ic为临界电流,φ为结两端超导波函数的相位差。
  • 电压-相位关系:V = (Φ0/2π)·dφ/dt,其中Φ0 = h/2e ≈ 2.07×10-15 Wb为磁通量子。

约瑟夫森结的势能为 EJ(1 − cos φ),其中EJ = IcΦ0/2π称为约瑟夫森能量。这个余弦势阱不是抛物线,因此量子化后的能级间距随能级升高而减小,天然打破了等间距结构。这正是量子比特所需的非线性。

工程上常用“约瑟夫森结电感”LJ = Φ0/(2πIccos φ)来近似描述结在小相位下的等效电感。结的物理尺寸(面积)直接决定Ic,进而决定EJ,这是版图设计中最重要的可调参数之一。

三、三种主流量子比特结构

电荷量子比特(Cooper Pair Box)由一个超导岛通过约瑟夫森结与 reservoir 耦合,岛上的库珀对数目n受门电压控制。其哈密顿量由 charging energy EC = (2e)2/(2CΣ) 和约瑟夫森能量EJ竞争决定。当EC ≫ EJ时,电荷态|n⟩是良好本征态,比特对电荷噪声极其敏感,相干时间通常只有纳秒到微秒量级。这是早期电荷比特的致命弱点。

磁通量子比特(Flux Qubit)由一个超导环路中嵌入1~3个约瑟夫森结构成,信息编码在环路持续电流的方向上(顺时针/逆时针)。其势能是余弦项的叠加,形成双势阱结构。两个势阱之间的隧穿由EJ控制,势阱深度由环路电感决定。磁通比特对磁通噪声敏感,但通过调节外部磁通偏置可以工作在“最优点”(sweet spot),此时一阶导数dω/dΦ = 0,对磁通噪声不敏感。

Transmon是电荷比特的改进版:通过增大并联电容CΣ,使EJ/EC从电荷比特的约1提高到20~100。此时电荷色散关系变得平坦,比特对电荷噪声的敏感度呈指数下降。代价是能级非谐性(anharmonicity)α = (E12 − E01)/h 减小到约 −200~−300 MHz,但仍足以区分|1⟩→|2⟩跃迁。Transmon的典型参数为:f01 ≈ 4~6 GHz,EC/h ≈ 200~300 MHz,EJ/h ≈ 10~20 GHz,T1可达100 μs以上。它是目前超导量子处理器中最主流的比特类型。

从LC谐振到三种超导量子比特的结构演化 LC谐振电路(线性) Eₙ = ħω(n+½) 能级等间距 → 无法选态 ω = 1/√(LC) 引入 约瑟夫森结(S-I-S) 超导 绝缘 超导 I I = I꜀ sin φ E = Eⱼ(1 − cos φ) 余弦势阱 → 能级非等间距 组合 三种量子比特结构 电荷比特 E꜀ ≫ Eⱼ,电荷噪声敏感 磁通比特 双势阱,持续电流方向编码 Transmon Eⱼ/E꜀ ≈ 20~100,电荷噪声不敏感 关键参数如何决定量子特性 E꜀ = (2e)²/2CΣ 充电能(电容决定) Eⱼ = I꜀Φ₀/2π 约瑟夫森能(结面积决定) f₀₁ ≈ √(8EⱼE꜀)/h 比特跃迁频率 α ≈ −E꜀ 非谐性(选态能力) T₁, T₂ 相干时间 设计权衡: 增大CΣ → E꜀减小 → 电荷噪声敏感度指数下降(Transmon路线),但非谐性α同步减小,|1⟩→|2⟩泄漏风险上升; 增大结面积 → Eⱼ增大 → f₀₁升高,但结寄生电容与工艺均匀性约束需同步评估。

四、参数如何决定量子特性:一个设计示例

假设我们要设计一个目标频率f01 ≈ 5.0 GHz的Transmon。取EC/h = 250 MHz(对应总电容CΣ ≈ 78 fF),则非谐性α ≈ −250 MHz。由近似关系f01 ≈ (√(8EJEC) − EC)/h,反解可得EJ/h ≈ 13.0 GHz,对应EJ/EC ≈ 52,处于Transmon的典型工作区间。若结的临界电流密度为Jc = 2 μA/μm²,则所需结面积约为Ic/(Jc) ≈ 13.0 GHz × h/(2e·Φ0/2π) / Jc,计算得Ic ≈ 26 nA,结面积约0.013 μm²,对应约200 nm × 65 nm的隧道结。这些数值属于教学示例假设,实际设计需结合具体工艺的Jc均匀性和结电阻测量结果迭代修正。

判断设计是否合理的三个观察点:① 频率是否落在微波馈线带宽内且避开两能级缺陷态密度较高的频段;② 非谐性|α|是否大于驱动脉冲带宽的3~5倍,以保证选态选择性;③ EJ/EC是否足够大,使电荷色散曲线在工作点附近足够平坦。

五、从结构到特性的因果链

三种比特的差异本质上来自同一个参数比值EJ/EC的不同取值。电荷比特取EJ/EC ≪ 1,电荷态是良好量子数,但噪声敏感;磁通比特通过环路几何引入双势阱,用磁通偏置而非电荷偏置来调谐;Transmon取EJ/EC ≫ 1,牺牲部分非谐性换取对电荷噪声的指数级抑制。理解这条因果链,就能在版图设计阶段预判:改变结面积会同时移动f01和α,改变并联电容主要移动EC从而同时影响α和噪声敏感度。这正是“结构决定量子特性”在工程上的具体含义。