2. 超导与约瑟夫森结基础

SQUID 之所以能测到远低于单个磁通量子的磁场变化,根源在于两件事:电子在超导态下结成库珀对并形成宏观量子相位,以及两个超导体被弱连接隔开时发生的约瑟夫森隧穿。理解这两点,才能理解后面 SQUID 的 I-V 曲线、调制特性和噪声来源。

超导能隙与库珀对

在低于临界温度 Tc 时,某些材料的电子之间通过晶格振动产生净吸引作用,两两配对成库珀对。配对使单电子激发谱打开一个能隙 Δ,即拆散一个库珀对所需的最小能量约为 2Δ。以铌(Nb)为例,Tc ≈ 9.2 K,Δ ≈ 1.5 meV,对应 2Δ/kBTc ≈ 3.5,这是 BCS 理论给出的普适比值。能隙的存在意味着:只要温度 kBT ≪ Δ,热激发几乎无法拆散库珀对,超导电流可以无阻流动。

更关键的是,所有库珀对凝聚到同一个量子态,可以用一个宏观波函数 Ψ = |Ψ|eiθ 描述。振幅 |Ψ| 的平方给出库珀对密度,相位 θ 在整个超导体内部保持相干。正是这个相位,成为约瑟夫森效应的“控制变量”。

约瑟夫森方程

当两块超导体通过绝缘层、微桥或点接触形成弱连接时,两侧相位 θ1、θ2 的差值 φ = θ2 − θ1 决定了两条基本关系:

直流约瑟夫森关系:  I = I_c · sin(φ)
交流约瑟夫森关系:  dφ/dt = 2e·V / ħ = 2πV / Φ₀

其中 Φ₀ = h/2e ≈ 2.07×10⁻¹⁵ Wb 为磁通量子

第一式说明:即使两端电压为零,也能流过不超过临界电流 Ic 的超导电流,且电流大小由相位差决定。第二式说明:一旦结上有直流电压 V,相位差就会以角频率 ω = 2eV/ħ 持续旋转,产生频率 f = 2eV/h 的高频振荡。例如 V = 1 μV 时,f ≈ 483.6 MHz。这两式是 SQUID 工作的物理内核:SQUID 通过外部磁通改变相位差,从而改变临界电流,最终把磁通变化转成可测的电压变化。

结电容与电阻参数

真实的约瑟夫森结不是理想元件。两块超导电极之间隔着绝缘层,构成一个平板电容,结电容 C 通常为 fF~pF 量级。同时,在有限温度或电压下存在准粒子隧穿通道,可用电阻 R 表示。于是结的等效电路是“理想约瑟夫森元件 + 电容 C 并联 + 电阻 R 并联”,这就是 RCSJ 模型(Resistively and Capacitively Shunted Junction)。

RCSJ 模型把结的总电流写成三支路之和:

I = I_c·sin(φ) + (V/R) + C·dV/dt
且 V = (ħ/2e)·dφ/dt

归一化后: i = sin(φ) + dφ/dτ + β_c·d²φ/dτ²
其中 β_c = 2e·I_c·R²·C / ħ 为 Stewart-McCumber 参数

βc 是判断结工作模式的核心参数。βc ≪ 1 时,电容被电阻充分阻尼,I-V 曲线无迟滞,称为过阻尼结;βc ≫ 1 时,结在回扫时出现迟滞回线,称为欠阻尼结。SQUID 通常希望工作在过阻尼或近临界阻尼状态,以避免磁滞导致的多值输出和额外噪声。

RCSJ 模型对 SQUID 行为的影响

把两个约瑟夫森结并联成一个超导环路,就得到直流 SQUID。环路中总相位差受外部磁通 Φ 约束:φ1 − φ2 = 2πΦ/Φ0。因此临界电流随磁通周期性调制:

I_c(Φ) = 2I_c0 · |cos(πΦ/Φ₀)|

当 Φ = nΦ₀ 时,I_c 最大;
当 Φ = (n+1/2)Φ₀ 时,I_c 最小。

在恒流偏置下,SQUID 两端电压随磁通呈近似正弦的周期变化,周期正好是一个 Φ₀。这就是 SQUID 能作为磁通—电压转换器的原因。RCSJ 模型进一步告诉我们:IcR 乘积决定电压调制幅度,R 和 C 决定响应速度和阻尼程度,而热噪声会使 I-V 曲线圆滑、降低调制深度。因此设计 SQUID 时,需要根据目标灵敏度和带宽,在 Ic、R、C 之间折中。

约瑟夫森结 RCSJ 模型与 SQUID 磁通—电压转换链 约瑟夫森结结构 超导体 1(相位 θ₁) 绝缘势垒层(C) 超导体 2(相位 θ₂) φ = θ₂ − θ₁ RCSJ 等效电路 JJ I_c·sinφ R(准粒子) C(结电容) 直流 SQUID 环路 结1 结2 Φ 临界电流与电压的磁通调制(周期 Φ₀) Φ I_c Φ₀ 2Φ₀ 3Φ₀ 4Φ₀ 5Φ₀ I_c 最小 I_c 最大 β_c ≪ 1:过阻尼,无迟滞 β_c ≫ 1:欠阻尼,有迟滞 I_c·R 决定电压调制幅度 热噪声使曲线圆滑

参数选择与工程判断

在实际设计或选用 SQUID 时,可以从 RCSJ 模型出发做几个快速判断。第一,看 βc:若 βc 大于 1,I-V 曲线会出现迟滞,恒流偏置下可能跳到电压态,需要改用恒压偏置或增加并联电阻阻尼。第二,看 IcR 乘积:它近似决定电压调制深度 ΔV,典型 Nb/AlOx/Nb 结的 IcR 在 0.1~1 mV 量级,太小则信号弱,太大则容易受热涨落影响。第三,看 C:结面积越大,C 越大,βc 越高,响应越慢,因此高频 SQUID 常采用小面积结。第四,看温度:kBT 与 IcΦ0/2π 的比值决定热噪声对调制曲线的展宽程度,液氦温度(4.2 K)下通常远小于 1,而液氮温度(77 K)下高温超导 SQUID 的噪声明显更大。

理解顺序建议:先抓住“相位差 φ 是控制量”,再记住两个约瑟夫森方程,然后把 C 和 R 并联进去得到 RCSJ,最后看 βc 和 IcR 如何决定 SQUID 的调制深度、迟滞和噪声。这样后续分析 SQUID 的偏置点、调制曲线和噪声谱时,就不会只停留在公式记忆层面。