3. Josephson结基本方程与静态特性

Josephson结之所以能成为超导量子计算、SQUID磁强计和超导电压基准的核心元件，根本原因在于它由两个简洁的方程支配。理解这两个方程，以及由它们导出的RCSJ模型和I-V曲线，是后续设计、表征和调试一切Josephson器件的基础。本节从直流效应出发，逐步引入交流效应，再通过RCSJ模型把结的静态与动态行为统一到一条可解析的I-V曲线上。

直流Josephson效应：零电压下的超导电流

当两块超导体通过弱连接（如薄绝缘势垒）耦合时，库珀对可以相干隧穿，形成无耗散的超导电流。这一现象由直流Josephson方程描述：

I_s = I_c sin φ

其中 I_s 是超导电流，I_c 是临界电流，φ 是两侧超导序参量的相位差。这个方程的关键含义是：在零电压状态下，结可以承载任意大小的直流电流，只要 |I| ≤ I_c；电流大小由相位差 φ 决定，而 φ 本身由外部电路条件约束。当电流超过 I_c 时，结无法继续维持零电压，必须切换到有阻态。

临界电流 I_c 不是任意常数，它由结的几何与材料参数决定。对于典型的隧道结，Ambegaokar-Baratoff 关系给出：

I_c = (π Δ(T) / (2 e R_N)) · tanh(Δ(T) / (2 k_B T))

其中 Δ 是超导能隙，R_N 是正常态电阻，e 是电子电荷，k_B 是玻尔兹曼常数。在低温极限 T ≪ T_c 下，tanh 项趋近于1，于是 I_cR_N ≈ πΔ/(2e)。这个乘积是表征结质量的重要指标：对于 Nb/Al-AlO_x/Nb 结，Δ ≈ 1.2 meV，理论 I_cR_N 约 1.9 mV。若实测值明显偏低，通常意味着势垒中存在针孔或界面退化。

交流Josephson效应：电压驱动相位演化

当结两端存在有限电压 V 时，相位差不再恒定，而是随时间线性演化。这就是交流Josephson方程：

dφ/dt = 2eV/ħ

其中 ħ 是约化普朗克常数。将两个方程结合，若 V 为常数，则 φ(t) = φ₀ + (2eV/ħ)t，代入直流方程得到：

I_s(t) = I_c sin(φ₀ + 2eVt/ħ)

这意味着在恒定电压下，超导电流以频率 f_J = 2eV/h 振荡。这个频率称为Josephson频率，换算关系为 f_J ≈ 483.6 GHz/mV。例如，当 V = 1 mV 时，振荡频率约为 483.6 GHz，落在亚毫米波段。这一效应是Josephson电压标准的基础：通过精确测量频率和电压，可以建立电压与基本常数的联系。

交流效应还解释了结在微波辐照下的Shapiro台阶现象：当外加微波频率 f_rf 与Josephson振荡发生锁相时，I-V曲线上会出现一系列等电压间隔的电流台阶，台阶电压满足 V_n = n h f_rf/(2e)。这是判断结是否正常工作的有力证据。

RCSJ模型：从理想结到真实器件的桥梁

实际Josephson结并非理想元件，它存在并联电容 C 和并联电阻 R。RCSJ（Resistively and Capacitively Shunted Junction）模型把结等效为理想Josephson元件与电阻、电容的并联组合。根据基尔霍夫电流定律，总电流为：

I = I_c sin φ + (V/R) + C dV/dt

利用交流Josephson方程将 V 替换为 (ħ/2e)dφ/dt，得到相位动力学方程：

I = I_c sin φ + (ħ/(2eR)) dφ/dt + (ħC/(2e)) d²φ/dt²

这个方程与机械摆的运动方程完全同构：I_c sin φ 对应重力恢复力，电阻项对应阻尼，电容项对应惯性。由此可以定义两个关键无量纲参数：

• McCumber参数 β_c = 2eI_cR²C/ħ，描述阻尼强弱。β_c ≪ 1 为过阻尼，β_c ≫ 1 为欠阻尼。
• 回滞判据：欠阻尼结在I-V曲线上表现出回滞——电流从零增大时在 I_c 处跳变到电压态，但电流减小时直到回滞电流 I_r ≈ 4I_c/(π√β_c) 才返回零电压态。

对于超导量子比特，通常希望 β_c 在 1 附近或略小于1，以兼顾可控性和避免过度回滞。对于SQUID磁强计，则常采用过阻尼设计以消除回滞，保证输出单值。

I-V曲线解析：从零电压态到有阻态

RCSJ模型给出的I-V曲线可以分为三个区域。在零电压态，|I| ≤ I_c，结电压为零，电流全部由超导分量承担。当电流超过 I_c 后，结进入有阻态，电压迅速上升，I-V曲线逐渐趋近于正常态电阻线 V = IR_N。在欠阻尼情况下，从有阻态回扫时会经过回滞路径，直到 I_r 处跳回零电压态。

实际测量中，我们通常用四线法采集I-V曲线：施加直流电流偏置，测量结两端电压。判断结是否正常，可以观察几个特征：零电压态电流是否平坦、跳变是否陡峭、回滞电流与临界电流的比值是否合理、以及高电压区是否接近 R_N。若零电压态出现斜率，说明存在并联电阻或测量引线电阻；若跳变圆滑，可能源于热涨落或微波干扰。

图中左侧为RCSJ等效电路，右侧为对应的I-V曲线。可以看到，零电压态在 I_c 处跳变到有阻态，回扫时在 I_r 处跳回，形成回滞环。回滞的大小直接由 β_c 决定：β_c 越大，回滞越明显。在微波辐照下，有阻态区域会出现等间距的Shapiro台阶，台阶电压由外加频率精确锁定。

参数提取与表征实践

从一条实测I-V曲线中，我们可以提取多个关键参数。首先读取临界电流 I_c，即零电压态跳变点的电流值。然后在高电压区拟合斜率，得到正常态电阻 R_N。回滞电流 I_r 与 I_c 的比值可用于估算 β_c：由 I_r/I_c ≈ 4/(π√β_c)，反解得到 β_c ≈ (4I_c/(πI_r))²。若已知结面积和势垒厚度，还可以从 R_N 和 I_c 反推势垒的比电阻和结电容。

需要注意的是，I_c 对温度敏感。在接近 T_c 时，I_c 按 (1-T/T_c)^1/2 或更复杂的规律下降。因此，表征时应在稳定的低温下进行，并记录温度。此外，环境电磁噪声会耦合到结上，导致I-V曲线出现台阶或跳变模糊。良好的屏蔽和滤波是获得干净曲线的必要条件。

总结来说，直流方程给出零电压态的电流-相位关系，交流方程给出电压驱动的相位演化，RCSJ模型把两者与电阻、电容结合，形成完整的相位动力学。I-V曲线则是这些方程在宏观测量中的直接体现。掌握从曲线中提取 I_c、R_N、β_c 的方法，是判断结质量、指导器件设计的第一步。