2. 超导量子比特的电磁建模基础

要对超导量子比特做电磁仿真与参数提取,第一步不是打开仿真器,而是把“量子比特”翻译成一套电路与电磁量都能描述的对象。超导量子比特本质上是一个非线性谐振电路:线性部分由电容、电感(含传输线、谐振腔)构成,非线性部分由约瑟夫森结提供。仿真的目标,是把版图几何、材料与结参数,映射成哈密顿量里的 E_C、E_J、ω_q、α 等量;反过来,设计时又要把目标频率与非线性,翻译回电容矩阵和结面积。理解这条双向映射链,是后续所有参数提取工作的基础。

一、从 LC 谐振子到约瑟夫森结非线性

一个理想 LC 谐振子的能量是二次型,能级等间距,无法单独寻址。约瑟夫森结由两段超导体夹一层薄绝缘势垒构成,其电流—相位关系为 I = I_c sin φ,对应的势能为 U(φ) = −E_J cos φ,其中 E_J = Φ₀ I_c / 2π,Φ₀ = h/2e ≈ 2.068×10⁻¹⁵ Wb 为磁通量子。把 cos φ 展开,二次项给出等效电感,四次项给出非线性。定义约瑟夫森电感

L_J = Φ₀ / (2π I_c cos φ)   (小信号,φ≈0 时 L_J0 = Φ₀ / (2π I_c))
E_J = Φ₀ I_c / (2π) = (Φ₀/2π)² / L_J0

这就是“非线性电感”的来源:它不是常数,而是随相位(即随电流)变化,因此能级间距不再相等,出现负的非谐 α = ω₁₂ − ω₀₁ < 0,从而可以用微波脉冲单独驱动 0→1 跃迁。

二、Transmon 的等效电路与参数映射

Transmon 是在电荷比特基础上并联一个大电容 C_B,把 E_C 压到远小于 E_J(典型 E_J/E_C ≈ 30–100),从而对电荷噪声不敏感。其等效电路为:一个约瑟夫森结(L_J 与 C_J)与大电容 C_B 并联,再通过耦合电容 C_g 接到读腔或驱动线。哈密顿量写作

H = 4E_C (n − n_g)² − E_J cos φ
E_C = e² / (2C_Σ),  C_Σ = C_B + C_J + C_g + ...
ω₀₁ ≈ √(8 E_J E_C) − E_C
α ≈ −E_C

这里 C_Σ 是节点总电容,由版图几何决定;E_J 由结面积与临界电流密度 J_c 决定。电磁仿真的任务,就是用静电场或射频场求解器提取 C_Σ、C_g,再结合工艺给出的 J_c 算出 E_J,最终得到 ω₀₁ 与 α。例如目标 ω₀₁/2π = 5.0 GHz、α/2π = −330 MHz,则 E_C/h ≈ 330 MHz,C_Σ ≈ e²/(2·E_C) ≈ 58 fF;再由 ω₀₁ ≈ √(8E_JE_C) − E_C 反解 E_J/h ≈ 9.5 GHz,对应 I_c ≈ 19.3 nA。这些数值属于教学示例,实际设计需按具体工艺 J_c 与版图重新提取。

Transmon 等效电路 JJ (E_J, C_J) C_B C_g C_Σ = C_B + C_J + C_g Fluxonium 等效电路 L_a (超电感) C_q (小) E_C 小、E_L 小、E_J 大 → 重通量子 几何/工艺 → 电磁量 → 哈密顿量参数 映射链 版图几何 C_B, C_g, L_a 工艺参数 J_c, 结面积 电磁提取 C_Σ, L_J, E_J 哈密顿量参数 E_C, E_J, ω₀₁, α 注:图中数值为教学示例,实际需按具体工艺与版图重新提取。

三、Fluxonium 的等效电路与参数特点

Fluxonium 在约瑟夫森结上并联一个超电感 L_a 和一个小电容 C_q,形成“重通量子”模式。其哈密顿量为

H = 4E_C n² − E_J cos φ + (1/2) E_L (φ − φ_ext)²
E_L = (Φ₀/2π)² / L_a

与 Transmon 不同,Fluxonium 的 E_J 可以很大(几十 GHz),E_C 和 E_L 都很小(约 0.5–1 GHz),因此波函数在相位上高度局域,对电荷噪声和磁通噪声都有较好的保护。其频率对 φ_ext 敏感,通常工作在半个磁通量子附近。电磁建模时,L_a 往往由超导纳米线或结阵列实现,需要用射频仿真提取其等效电感与寄生电容;C_q 则决定 E_C,进而影响非谐与驱动强度。

四、参数提取的实际流程与判断依据

以 Transmon 为例,一次完整的参数提取通常包含以下环节。首先在版图中确定电容盘尺寸、结面积与耦合间隙;然后用静电场求解器(如有限元或边界元)提取电容矩阵,得到 C_Σ 与 C_g。接着根据工艺给出的 J_c 与结面积 A 计算 I_c = J_c·A,再得到 E_J。最后代入公式计算 ω₀₁ 与 α,并与目标值比较。

判断结果是否合理,可以看三个量:E_J/E_C 是否落在 30–100 之间(太小则电荷噪声敏感,太大则非谐过小、驱动变慢);α/2π 是否在 −200 至 −400 MHz 之间(决定门速度与泄漏);C_g 是否与读腔耦合强度匹配(影响读取保真度与 Purcell 限制)。若 ω₀₁ 偏离目标,优先调结面积或 J_c;若 α 不合适,则调电容盘尺寸改变 C_Σ。这些调整在仿真中迭代,直到参数落入可接受窗口。

建模时注意:约瑟夫森结的 C_J 通常不可忽略,尤其在小结面积下可能占总电容的 10%–20%;提取 C_Σ 时应把结的寄生电容一并计入,否则 E_C 会被高估。

五、从电磁量到哈密顿量的双向映射

正向映射是“设计→仿真→参数”:给定几何与工艺,算出 E_C、E_J、ω₀₁、α。反向映射是“目标→参数→几何”:给定目标频率与非线性,反解 C_Σ 与 I_c,再回到版图尺寸。实际工作中,反向映射往往先给出初值,再用正向仿真迭代修正。理解这条双向链,才能在仿真结果偏离时快速定位是电容提取、结参数还是耦合设计出了问题。

对于 Fluxonium,映射关系更复杂:ω₀₁ 同时依赖 E_J、E_C、E_L 和外部磁通 φ_ext,且对 φ_ext 的敏感度很高。仿真时通常先固定 φ_ext 在工作点,再扫描 E_L 与 E_C,观察频率与非谐的变化,找到对噪声不敏感的区域。这一过程需要把电磁提取得到的 L_a、C_q 与结参数一起代入,才能得到可信的能谱。

总结来说,超导量子比特的电磁建模,核心是把“几何—电容—电感—结—哈密顿量”这条链打通。Transmon 的关键是 C_Σ 与 E_J 的比值,Fluxonium 的关键是 E_L 与 E_J 的配合。掌握这些映射关系,后续的谐振腔耦合、读取链设计和参数提取才有可靠的物理依据。