2. 超导量子比特的电路量子电动力学基础

要理解超导量子比特为什么能像“人造原子”一样工作，需要把三个物理对象串起来看：LC振荡器提供线性谐振子的经典图像，约瑟夫森结提供无耗散的非线性电感，二者组合成非线性振子，最终在低温量子极限下表现为分立能级。本节按这条线索展开，并给出可用于数值分析的参数化描述。

一、LC振荡器：线性谐振子的经典与量子图像

理想LC回路由电感L和电容C并联构成，谐振角频率为 ω₀ = 1/√(LC)。在经典图像中，能量在电感磁能与电容电能之间往复交换，相空间轨迹是椭圆；在量子图像中，哈密顿量写成 H = Q²/(2C) + Φ²/(2L)，其中Q为节点电荷、Φ为节点磁通，二者满足对易关系 [Φ, Q] = iℏ。引入升降算符后 H = ℏω₀(a†a + 1/2)，能级等间距分布：Eₙ = ℏω₀(n + 1/2)。

等间距带来一个关键后果：相邻能级跃迁频率完全相同，都是 ω₀。若用单一频率的微波驱动，会同时激发 0→1、1→2、2→3 等所有跃迁，无法把系统当作只使用最低两个能级的二能级系统。这正是LC谐振器不能直接做量子比特的根本原因——它缺少非谐性（anharmonicity）。

工程上常用传输线谐振器或集总元件谐振器实现，典型参数例如 L ≈ 1 nH、C ≈ 1 pF，则 ω₀/2π ≈ 5 GHz，对应 ℏω₀ ≈ 20 μeV，远小于室温热能 k_BT ≈ 26 meV，因此必须在稀释制冷机中把温度降到 10 mK 量级，使 k_BT ≪ ℏω₀，否则热激发会淹没量子能级结构。

二、约瑟夫森结：无耗散的非线性电感

约瑟夫森结由两片超导体夹一层极薄绝缘势垒构成，其核心关系为：超导电流 I = I_c sin φ，结两端电压 V = (Φ₀/2π) dφ/dt，其中 φ 是两侧超导序参量的相位差，Φ₀ = h/2e ≈ 2.07×10⁻¹⁵ Wb 为磁通量子，I_c 为临界电流。对电流关系求导可得结的微分电感：

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

当 φ 接近 0 时 L_J ≈ Φ₀/(2πI_c)，是一个确定值；但随 φ 增大，L_J 随 cos φ 变化，呈现非线性。更重要的是，结的势能 U(φ) = −E_J cos φ，其中 E_J = Φ₀I_c/2π 称为约瑟夫森能量。这个余弦势不是抛物线，正是它打破了能级的等间距。

需要注意，约瑟夫森结本身不引入欧姆耗散（在低于能隙的电压下准粒子电流被抑制），这使它能作为量子电路中的低损耗非线性元件。实际器件中，结的并联电容 C_J 与外部并联电容共同决定充电能 E_C = e²/(2C_Σ)，这是另一个关键能量尺度。

三、非线性振子：从余弦势到可分辨能级

把约瑟夫森结并联一个电容，就得到最简单的超导量子比特雏形——库珀对盒或传输子（transmon）。其哈密顿量为：

H = 4E_C n² − E_J cos φ

其中 n 为库珀对数目算符，与 φ 满足 [φ, n] = i。当 E_J ≫ E_C 时（transmon 工作区，典型 E_J/E_C ≈ 20～100），相位涨落很小，可把余弦势在 φ = 0 附近展开：

−E_J cos φ ≈ −E_J + (E_J/2)φ² − (E_J/24)φ⁴ + …

二次项给出谐振频率 ℏω₀₁ ≈ √(8E_JE_C)，四次项则引入负的非谐性修正。按微扰论可得近似结果：

E₀₁ ≈ √(8E_JE_C) − E_C
α = E₁₂ − E₀₁ ≈ −E_C

也就是说，transmon 的非谐性 α 约等于 −E_C，典型值在 −200 MHz 到 −350 MHz 之间。这个量级远大于线宽和谐振器耦合强度，因此可以用微波脉冲单独驱动 0↔1 跃迁，而几乎不激发 1↔2 跃迁，从而把器件当作二能级量子比特使用。

四、能级量子化的物理图像与数值分析要点

能级量子化的直接证据来自势阱中波函数的离散本征值。对 transmon，φ 的涨落均方根约为 (8E_C/E_J)^(1/4)，E_J/E_C 越大，波函数越局域在余弦势底部，能级越接近谐振子但保留固定非谐性。数值分析时通常采用以下两条路径：

• 相位基展开：在 φ ∈ [−π, π] 上取 N 个格点（如 N = 101），构造动能矩阵 4E_C n² 的离散傅里叶表示与势能对角矩阵 −E_J cos φ，对角化得到本征值 E₀、E₁、E₂，再计算 ω₀₁ 与 α。
• 电荷基展开：以库珀对数目 n 为基，势能项 −E_J cos φ 耦合相邻 n，形成三对角块矩阵，截断到 |n| ≤ 10 即可获得千赫兹量级收敛精度。

判断数值结果是否可信，可检查三点：一是截断维数加倍后 ω₀₁ 变化小于 1 MHz；二是 α 的符号为负且量级接近 −E_C；三是波函数在 φ = ±π 边界处幅值接近零。若边界幅值不可忽略，说明 E_J/E_C 过小或截断不足，需要扩大相位范围或增加基矢。

五、从物理图像到数值分析的工作流

把上述物理图像落到数值分析，通常按以下顺序推进：先由版图或设计值确定 E_J 与 E_C（E_J 由结面积和临界电流密度决定，E_C 由总电容决定）；再选定基矢与截断；然后对角化哈密顿量，提取 E₀、E₁、E₂；最后计算 ω₀₁、α 并与设计目标比对。若 ω₀₁ 偏离目标，优先调整 E_J（即结面积）；若 α 偏离目标，优先调整 E_C（即并联电容），因为 α ≈ −E_C 对电容更敏感。

需要强调的是，以上数值均为教学示例中的合理假设，真实器件还需考虑结的寄生电容、谐振器耦合、磁通偏置点以及材料损耗。但无论细节如何，LC谐振器—约瑟夫森结—非线性振子—分立能级这条主线不变：线性谐振子给出等间距能级的起点，余弦势打破等间距，非谐性使最低两个能级可以被单独寻址，从而构成超导量子比特的物理基础。