3. 约瑟夫森结室温电阻与临界电流换算

室温电阻到低温临界电流的换算链

约瑟夫森结的临界电流 Ic 无法在室温下直接测量，但结在正常态下的室温电阻 Rn 可以稳定测得。Ambegaokar-Baratoff（AB）关系把两者联系起来：在零温、BCS 弱耦合近似下，

Ic(0) = (π·Δ(0) / (2·e·Rn)) · tanh(Δ(0) / (2·k_B·T))

当 T ≪ Tc 时 tanh 项趋近 1，工程上常用简化式 Ic ≈ π·Δ/(2·e·Rn)。对铝（Al）结，Δ(0) ≈ 180 μeV，代入常数可得 Ic·Rn ≈ 0.28 V（更精确地，πΔ/(2e) ≈ 0.283 V）。因此只要测得 Rn，就能快速估算 Ic：Ic[A] ≈ 0.283 / Rn[Ω]。

换算的实操要点

测量 Rn 时应在远高于 Tc 的温度（如 4.2 K 以上或室温）用四线法，避免引线电阻和接触电阻混入。若在室温测量，需注意 Rn 随温度有轻微变化，通常以接近 Tc 以上的正常态平台值为准。换算时先确认结的材料体系：铝结取 Ic·Rn ≈ 0.28 V，铌结 Δ 更大，系数相应提高，不能混用。

下面给出一个可直接运行的换算与离散性评估脚本，输入为各结的实测 Rn 与设计面积，输出 Ic 估算值、均值与相对标准差。

import numpy as np

# 材料常数：铝结 Ic*Rn 乘积（V），示例取 0.283
IC_RN_PRODUCT = 0.283

# 示例数据：结面积 A (um^2) 与实测室温电阻 Rn (kOhm)
areas = np.array([0.04, 0.04, 0.04, 0.09, 0.09, 0.09, 0.16, 0.16, 0.16])
rn_kohm = np.array([12.1, 12.5, 13.4, 5.4, 5.6, 5.9, 3.0, 3.1, 3.2])

rn_ohm = rn_kohm * 1e3
ic = IC_RN_PRODUCT / rn_ohm          # 单位 A
ic_na = ic * 1e9                     # 单位 nA

for a in np.unique(areas):
    mask = areas == a
    ic_g = ic_na[mask]
    rn_g = rn_kohm[mask]
    print(f"A={a:.2f} um^2 | Rn={rn_g.mean():.2f} kOhm "
          f"| Ic={ic_g.mean():.1f} nA "
          f"| Ic相对离散={ic_g.std(ddof=1)/ic_g.mean()*100:.1f}%")

示例输出（由上述输入直接计算得到）：

A=0.04 um^2 | Rn=12.67 kOhm | Ic=22.4 nA | Ic相对离散=5.4%
A=0.09 um^2 | Rn=5.63 kOhm | Ic=50.3 nA | Ic相对离散=4.5%
A=0.16 um^2 | Rn=3.10 kOhm | Ic=91.3 nA | Ic相对离散=3.2%

离散性评估与设计判断

从上表可见，Rn 的绝对波动（约 0.1–0.9 kΩ）在小面积结上占比更大，导致 Ic 相对离散随面积减小而上升。若 Transmon 频率对 Ic 的容差要求为 ±3%，则 0.04 μm² 的结在示例数据下已接近或超出边界，需要收紧工艺均匀性或改用更大面积结。评估时应以同一批次、同一晶圆上的多点测量为样本，计算组内标准差，而不是用单点值代表整片。

实际设计中，可把目标 Ic 与允许的相对离散作为输入，反推所需的 Rn 均值与工艺控制窗口：Rn_target = 0.283 / Ic_target，再要求 σ_Rn / Rn_target 不超过 Ic 的相对容差。这样就把版图面积选择与工艺统计直接关联起来。