3. Transmon 量子比特设计与参数提取

Transmon 是目前超导量子计算中最主流的量子比特类型。它的核心思想是把一个超导量子干涉仪（SQUID）与一个大电容并联，通过增大电容来抑制电荷噪声对能级的影响。要真正设计出一个可用的 Transmon，需要从版图几何出发，提取电容矩阵，计算约瑟夫森结能量与充电能量，最终确定比特频率、非谐性和可调性范围。本节围绕这条设计链路展开。

一、从版图到电容矩阵

Transmon 的电容网络由若干导体岛和它们之间的耦合电容构成。典型结构中包含：一个主岛（连接 SQUID 的一端）、一个接地平面、以及用于耦合的邻近结构（如读取谐振腔、耦合总线或其他比特）。设计的第一步是确定版图几何：叉指电容的指宽、指间距、指长、交叠长度，以及两个岛之间的间隙。

在电磁仿真中，我们把这些导体定义为端口，求解静电场，得到电容矩阵 C。对于一个包含 N 个导体的系统，电容矩阵是一个 N×N 的对称矩阵，对角元 C_ii 是导体 i 的自电容（即导体 i 与无穷远及其他导体之间的总电容），非对角元 C_ij 是导体 i 与导体 j 之间的互电容。实际设计中，我们更关心的是节点电容矩阵，即把接地平面作为参考节点后，各岛之间的等效电容。

以一个简化的两岛系统为例：主岛 A 与接地平面之间的电容为 C_A，主岛 A 与耦合岛 B 之间的电容为 C_AB，耦合岛 B 与地之间的电容为 C_B。在接地参考下，节点电容矩阵为：

C_node = [ C_A + C_AB    -C_AB      ]
         [   -C_AB     C_B + C_AB   ]

对于单比特设计，我们通常只关心主岛对地的总电容 C_Σ，它等于主岛与所有其他导体（包括地、耦合结构、控制线）之间电容的总和。这个 C_Σ 直接决定了充电能量 E_C。

二、约瑟夫森结能量与充电能量

约瑟夫森结是 Transmon 的非线性元件。一个结的约瑟夫森能量为 E_J = I_c Φ_0 / (2π)，其中 I_c 是结的临界电流，Φ_0 = h/(2e) 是磁通量子。对于 SQUID 结构，两个结并联，等效临界电流为 I_c(Φ_ext) = 2 I_c0 |cos(π Φ_ext / Φ_0)|，因此等效 E_J 可以通过外部磁通调节。

充电能量定义为 E_C = e² / (2 C_Σ)。Transmon 的工作区间要求 E_J / E_C ≫ 1，典型值在 50 到 100 之间。这个比值越大，电荷色散越小，对电荷噪声越不敏感；但比值过大时，非谐性会变小，导致比特与谐振腔之间的频率选择性变差。因此设计上需要在噪声免疫和非谐性之间取平衡。

以一个示例参数为例：假设通过仿真提取得到 C_Σ = 85 fF，则 E_C = e²/(2×85 fF) ≈ 0.226 GHz（以频率单位表示，E_C/h）。若设计目标 E_J/h = 18 GHz，则 E_J/E_C ≈ 80，处于合理区间。

三、频率与非谐性

Transmon 的能级在 E_J/E_C ≫ 1 时近似为：

E_0 ≈ -E_J + √(8 E_J E_C) / 2
E_1 ≈ -E_J + √(8 E_J E_C) × 3/2
E_2 ≈ -E_J + √(8 E_J E_C) × 5/2

由此得到比特频率 f_01 = (E_1 - E_0)/h ≈ (√(8 E_J E_C) - E_C)/h，非谐性 α = (E_2 - E_1)/h - (E_1 - E_0)/h ≈ -E_C/h。可以看到，非谐性主要由充电能量决定，而比特频率由 E_J 和 E_C 共同决定。

继续上面的示例：E_J/h = 18 GHz，E_C/h = 0.226 GHz，则 √(8 E_J E_C)/h ≈ √(8×18×0.226) ≈ 5.70 GHz，f_01 ≈ 5.70 - 0.226 ≈ 5.47 GHz，非谐性 α ≈ -0.226 GHz。这个非谐性足以把 |1⟩ → |2⟩ 跃迁与腔频率拉开，避免泄漏。

四、频率可调性与设计约束

通过 SQUID 的磁通调节，E_J 可以在 E_J_max 到接近 0 之间变化，因此比特频率可以在一个范围内连续调节。最大频率对应 Φ_ext = 0（或整数倍 Φ_0），最小频率对应 Φ_ext = Φ_0/2。实际设计中，我们通常把工作点设在 Φ_ext 略偏离 Φ_0/2 的位置，以保留足够的可调范围，同时避免磁通噪声敏感度过高。

可调性范围由 E_J_max 和 E_J_min 决定。若 E_J_max/h = 20 GHz，E_J_min/h = 2 GHz，则频率调节范围大约从 5.8 GHz 到 2.0 GHz。设计时需要确保：比特频率与读取谐振腔频率保持足够失谐（通常 1–2 GHz），两比特门工作点不与其他能级发生意外共振，且磁通噪声引起的频率抖动在可接受范围内。

五、参数提取的实操要点

在实际设计流程中，电容矩阵的提取通常借助有限元电磁仿真工具完成。操作上需要注意几点：第一，端口定义要覆盖所有可能影响 C_Σ 的导体，包括读取腔、耦合总线、磁通偏置线，甚至封装壁；第二，仿真收敛性要检查，网格加密后 C_Σ 的变化应小于 1%；第三，提取出的节点电容矩阵需要与电路模型对应，确认参考节点选择正确。

约瑟夫森结的参数则来自工艺线提供的临界电流密度 J_c。若 J_c = 1.5 μA/μm²，结面积 A = 0.04 μm²，则 I_c = 60 nA，E_J/h = I_c Φ_0 / (2π h) ≈ 30 GHz。SQUID 两个结并联后，最大等效 E_J/h 约为 60 GHz（若两结相同），但实际中会因结面积偏差而略有不对称，导致 E_J_min 不为零。这个不对称性需要在设计时预留，通常要求 E_J_max/E_J_min > 10 以保证足够的可调范围。

总结来说，Transmon 设计是一个从几何到电容、从电容到能量、从能量到频率的逐级映射过程。电容矩阵决定了 E_C，约瑟夫森结决定了 E_J，两者的比值决定了比特的基本工作区间，而 SQUID 的磁通调节能力则赋予频率可调性。理解这条链路中每个参数的来源和影响，是后续进行比特控制、读取和耦合设计的基础。