2. Transmon量子比特频率与色散参数计算

从结参数到比特频率：一条可复算的参数链

Transmon 比特的核心参数只有两个：约瑟夫森结的临界电流 I_c 与并联电容 C_Σ。它们先合成约瑟夫森能量 E_J 与充电能量 E_C，再决定能级间距、非谐性和色散频移。本节用一组明确标注为“教学示例”的参数走完整条链路，并给出容差分析的可执行做法。

计算约定：取 Φ₀ = h/(2e) ≈ 2.0678×10⁻¹⁵ Wb，h ≈ 6.626×10⁻³⁴ J·s，e ≈ 1.602×10⁻¹⁹ C。能量统一用频率单位 GHz 表示（即 E/h），这样所有结果可直接与频谱仪读数对照。

第一步：由 I_c 与 C_Σ 得到 E_J、E_C

约瑟夫森能量与充电能量分别为：

E_J = Φ₀ · I_c / (2π)          [J]
E_C = e² / (2 · C_Σ)           [J]
换算成频率：E_J/h、E_C/h        [Hz]

示例参数取 I_c = 25 nA、C_Σ = 95 fF。代入得 E_J/h ≈ 8.22 GHz、E_C/h ≈ 0.203 GHz，比值 E_J/E_C ≈ 40.5，处于典型 Transmon 区间（20–60），说明电荷噪声敏感度已被有效压制。

第二步：能级频率、非谐性与色散频移

Transmon 的能级近似为 E_n ≈ √(8 E_J E_C)·(n+1/2) − E_C/12·(6n²+6n+3)。由此得到：

  • 跃迁频率 f₀₁ ≈ (√(8 E_J E_C) − E_C)/h ≈ 5.62 GHz；
  • 非谐性 α = f₁₂ − f₀₁ ≈ −E_C/h ≈ −0.203 GHz，负号表示能级间距随 n 增大而收缩；
  • 与谐振腔耦合时，色散频移 χ ≈ −(g²/Δ)·α/(Δ+α)，取 g/2π = 80 MHz、Δ/2π = 1.2 GHz，得 χ/2π ≈ +1.0 MHz（符号取决于失谐方向，此处 Δ 取腔频高于比特频）。

注意 χ 与 α 同号且量级远小于 Δ，这正是可分辨比特态、又不破坏比特相干的前提。

约瑟夫森结 I_c = 25 nA 并联电容 C_Σ = 95 fF E_J/h ≈ 8.22 GHz E_J/E_C ≈ 40.5 E_C/h ≈ 0.203 GHz 电荷噪声已压制 f₀₁ ≈ 5.62 GHz 比特工作频率 α ≈ −0.203 GHz 非谐性（负） χ/2π ≈ +1.0 MHz g=80 MHz, Δ=1.2 GHz 容差敏感度：±5% 参数扰动对 f₀₁ 的影响（示例估算） −5% 标称值 +5% f₀₁ I_c ±5% → f₀₁ ∓2.5% C_Σ ±5% → f₀₁ ±2.5% 目标窗口 ±1%

第三步：容差如何吃掉频率目标

由 f₀₁ ∝ √(I_c/C_Σ) 可知，I_c 与 C_Σ 对频率的敏感度都是半阶：任一参数 ±5% 只带来约 ∓2.5% 的频率偏移。但芯片上真正难控的是结面积与氧化层厚度，I_c 的片内离散常达 ±3%~±8%，而 C_Σ 由版图决定、精度可做到 ±1% 以内。因此频率误差主要来自 I_c。

若目标是把 f₀₁ 控制在 ±1%（约 ±56 MHz）以内，则要求 I_c 偏差不超过 ±2%。单靠工艺做不到时，工程上通常两条路并行：其一，在结旁并联一组可激光修调的电容（trim cap），把 C_Σ 向上微调以拉回频率；其二，设计时把 f₀₁ 目标定在腔频失谐 Δ 的容差带中央，给 χ 留出余量。下表给出可直接套用的敏感度系数。

参数标称值扰动f₀₁ 变化α 变化χ 变化
I_c25 nA+5%+2.5%（+140 MHz）≈0≈0
I_c25 nA−5%−2.5%（−140 MHz）≈0≈0
C_Σ95 fF+5%−2.5%（−140 MHz）−5%（−10 MHz）−5%
C_Σ95 fF−5%+2.5%（+140 MHz）+5%（+10 MHz）+5%
做容差分析时，先算 ∂f₀₁/∂ln I_c = +f₀₁/2、∂f₀₁/∂ln C_Σ = −f₀₁/2，再乘参数相对偏差，比反复重算完整能级公式快得多，且误差在 ±10% 扰动内可忽略。
色散频移 χ 的符号由失谐方向决定：腔频高于比特频（Δ>0）时 χ 为正，反之为负。读腔频谱时若把符号搞反，会把比特态 |0⟩、|1⟩ 的峰位判反。计算前务必先固定 Δ 的符号约定，并在正文、表格与版图标注中保持一致。

可执行的计算脚本

下面这段 Python 不依赖第三方库，直接复现上述全部数值，并输出容差扫描表，可作为设计初期的快速核算工具。

import math

PHI0 = 2.067833848e-15   # Wb
H    = 6.62607015e-34    # J·s
E    = 1.602176634e-19   # C

def transmon(ic_nA, c_fF, g_MHz=80.0, delta_GHz=1.2):
    ic = ic_nA * 1e-9
    c  = c_fF * 1e-15
    EJ = PHI0 * ic / (2 * math.pi)      # J
    EC = E**2 / (2 * c)                 # J
    f01 = (math.sqrt(8 * EJ * EC) - EC) / H / 1e9   # GHz
    alpha = -EC / H / 1e9                            # GHz
    g = g_MHz / 1e3
    chi = -(g**2 / delta_GHz) * alpha / (delta_GHz + alpha)  # GHz
    return f01, alpha, chi, EJ / EC

f01, alpha, chi, ratio = transmon(25.0, 95.0)
print(f"f01={f01:.3f} GHz  alpha={alpha:.3f} GHz  chi={chi*1e3:.2f} MHz  EJ/EC={ratio:.1f}")

for d in (-0.05, 0.0, 0.05):
    f_ic, _, _, _ = transmon(25.0 * (1 + d), 95.0)
    f_c,  _, _, _ = transmon(25.0, 95.0 * (1 + d))
    print(f"扰动 {d:+.0%}: dI_c -> {f_ic:.3f} GHz, dC -> {f_c:.3f} GHz")

预期输出（示例参数下的计算结果，非实测）：f01≈5.62 GHz、alpha≈−0.203 GHz、chi≈+1.0 MHz、EJ/EC≈40.5；扰动行显示 I_c 与 C_Σ 各 ±5% 时频率分别向相反方向偏移约 0.14 GHz。若你的输出与表中数值偏差超过 1%，优先检查单位换算（nA→A、fF→F）与 Φ₀ 是否误用了 h/e 而非 h/2e。

结果判读与设计闭环

  • E_J/E_C 落在 20–60：电荷噪声压制充分，非谐性仍在 150–300 MHz 量级，可用纳秒级门操作区分 |1⟩→|2⟩ 泄漏。
  • f₀₁ 与腔频失谐 Δ 应保持在 0.8–1.5 GHz：太小则 χ 增大但比特受腔噪声反作用，太大则读取信号变弱。
  • 若容差扫描显示 f₀₁ 超出 ±1% 窗口，优先在版图中加入 trim cap 阵列，而不是放宽 I_c 的工艺规格——后者会同时恶化 χ 与读取保真度。
← 上一章 📖 返回目录 下一章 →