📚 课程目录第2章超导量子器件的物理基础

2. 超导量子器件的物理基础

要理解超导量子器件的频谱响应，必须先弄清三个层层递进的物理模型：约瑟夫森结提供非线性电感，LC 振荡电路提供谐振频率的“骨架”，而非线性谐振子模型则把二者统一起来，解释为什么量子比特的能级间隔会随激发态升高而变小。本节从器件物理出发，说明这些模型如何决定频谱响应与可测参数。

约瑟夫森结：非线性电感的来源

约瑟夫森结由两层超导体夹一层极薄绝缘势垒构成。当势垒足够薄时，库珀对可以隧穿通过，形成无耗散的超导电流。其核心关系是直流约瑟夫森方程：

I = I_c · sin(φ)
V = (Φ₀ / 2π) · dφ/dt
Φ₀ = h / 2e ≈ 2.0678 × 10⁻¹⁵ Wb

其中 I_c 是临界电流，φ 是结两端超导相位差，Φ₀ 是磁通量子。第一个式子说明电流与相位差呈正弦关系，而非线性；第二个式子说明电压正比于相位变化率。对相位求导可得结的小信号电感：

L_J(φ) = Φ₀ / (2π · I_c · cos φ)

当 φ 接近 0 时，L_J 最小；当 φ 接近 π/2 时，L_J 急剧增大。正是这种随相位（或电流）变化的电感，使约瑟夫森结成为可调的非线性元件。在典型 transmon 量子比特中，I_c 约为 20～100 nA，对应 L_J 约 3～15 nH，而结电容 C_J 约 1～10 fF，二者共同决定等离子体振荡频率。

LC 振荡电路：谐振频率的骨架

把约瑟夫森结与并联电容组合，就得到最基本的超导 LC 振荡电路。理想 LC 电路的谐振频率为：

ω₀ = 1 / √(L · C)
f₀ = ω₀ / (2π)

若取 L = 10 nH、C = 5 fF，则 f₀ ≈ 0.71 GHz。实际器件中，电容由版图几何决定，电感则来自约瑟夫森结。由于 L_J 随相位变化，谐振频率也会随驱动功率或磁通偏置而移动，这正是频谱测量中可观测的关键效应。

在读取腔设计中，常把约瑟夫森结阵列或超导传输线作为电感，与叉指电容或共面波导电容构成谐振器。读取腔的谐振频率通常设计在 5～8 GHz，与量子比特的工作频率保持一定失谐，以便通过色散读取获取比特状态信息。

非线性谐振子模型：从线性到量子

把 LC 电路量子化后，哈密顿量可写为：

H = 4E_C (n - n_g)² - E_J cos φ

其中 E_C = e²/(2C_Σ) 是充电能，E_J = Φ₀ I_c/(2π) 是约瑟夫森能。当 E_J ≫ E_C 时，相位涨落很小，可把 cos φ 展开到四阶：

H ≈ ω₀ a†a + (α/2) a†a†aa
ω₀ ≈ √(8 E_J E_C) - E_C
α ≈ -E_C

这里 α 是非谐性，即相邻能级间隔的差值。对 transmon，典型 E_C/h ≈ 200～350 MHz，E_J/h ≈ 10～20 GHz，因此 ω₀/2π 约 4～6 GHz，α/2π 约 -200～-350 MHz。负号意味着能级越高，间隔越小，这使我们可以用微波脉冲单独寻址 |0⟩↔|1⟩ 跃迁，而尽量不激发 |1⟩↔|2⟩。

模型如何决定频谱响应与可测参数

三个模型共同决定了器件在频域中的行为。约瑟夫森结的非线性电感 L_J(φ) 使谐振频率随相位变化，因此当外加磁通改变 φ 时，f₀ 会周期性移动，周期为 Φ₀。这一特性在频谱测量中表现为谐振峰位置的周期性漂移，是判断约瑟夫森结是否正常工作的直接证据。

LC 电路决定了谐振频率的量级。通过设计 C_J 和 I_c，可以把 f₀ 设定在目标频段。例如，若希望 f₀ = 5 GHz，取 C_J = 5 fF，则所需 L_J ≈ 2.0 nH，对应 I_c ≈ 160 nA。实际版图设计中还需考虑寄生电容和结阵列的串并联方式。

非线性谐振子模型则解释了能级结构。由于 α < 0，|1⟩↔|2⟩ 的跃迁频率低于 |0⟩↔|1⟩。在频谱测量中，若用固定频率微波驱动并扫描功率，可以观察到 |0⟩↔|1⟩ 跃迁先饱和，随后 |1⟩↔|2⟩ 跃迁出现，表现为双峰或峰分裂。通过拟合这些峰的位置和宽度，可以提取 ω₀₁、ω₁₂、α 以及退相干时间 T₁、T₂。

在频谱测量中，建议先做粗扫（宽频、低功率）找到谐振峰大致位置，再做细扫（窄频、多功率点）提取参数。粗扫时功率不宜过高，否则非线性效应会使峰形畸变，难以判断真实 f₀。

从模型到参数提取的实践要点

实际测量中，矢量网络分析仪输出的 S₂₁ 参数包含幅度和相位信息。对谐振峰进行拟合时，常用输入-输出模型：

S₂₁(f) = 1 - (Q_l / Q_c) / (1 + 2j Q_l (f - f₀)/f₀)

其中 Q_l 是负载品质因数，Q_c 是耦合品质因数。拟合得到 f₀ 和 Q_l 后，结合版图提取的 C_J，可反推 L_J 和 I_c。若测得的 f₀ 与设计值偏差超过 5%，通常需要检查结面积、氧化条件或寄生电容。

对于非线性参数 α，可通过双音谱或功率扫描提取。在低功率下，谐振峰近似洛伦兹形；随着功率升高，峰向低频移动（对应 α < 0），移动量 Δf ≈ -α · n̄，其中 n̄ 是腔内平均光子数。通过测量不同功率下的峰位，拟合 Δf 与 n̄ 的关系，即可得到 α。

磁通偏置下的 f₀ 调谐曲线也是重要诊断工具。若 f₀ 随磁通呈标准余弦形，说明约瑟夫森结行为正常；若曲线不对称或周期异常，可能提示结阵列中存在不对称或额外电感。这些判断依据都直接来自本节三个模型的组合分析。

注意：约瑟夫森结的临界电流 I_c 对温度敏感。在稀释制冷机中，若实际工作温度高于设计值，I_c 会下降，导致 L_J 增大、f₀ 降低。频谱测量前应确认器件处于基础温度（通常 < 20 mK），否则提取的参数会偏离真实值。

总结来看，约瑟夫森结提供非线性电感，LC 电路设定谐振频率，非线性谐振子模型给出能级间隔和非谐性。三者共同决定了频谱响应中的峰位、峰形和功率依赖关系，而 f₀、Q、g、α 等可测参数正是从这些响应中拟合提取的。理解模型与参数之间的映射关系，是后续进行频谱测量和参数提取的基础。