📚 课程目录第3章 RCSJ等效电路模型

3. RCSJ等效电路模型

在真实器件中，Josephson结并不是一个孤立的理想超导弱连接。两块超导电极之间隔着极薄的势垒层，天然构成一个平板电容器；势垒层存在的微小漏电通道、引线与电极的欧姆损耗，又构成一条并联的电阻支路。把这三条物理通道并联起来，就得到本节要讨论的RCSJ模型（Resistively and Capacitively Shunted Junction，电阻-电容并联分流的结模型）。它是分析结的直流I-V特性、开关行为与动态响应的最基础、也最常用的等效电路。

一、电路拓扑与三条支路

RCSJ模型把结抽象为三条并联支路，共用同一个结电压 V（即两端电极间的电压）：

超导支路：承载库珀对隧道电流，满足第一Josephson关系 I_s = I_c sin φ，其中 I_c 为临界电流，φ 为结两端超导波函数的相位差。
电阻支路：承载正常电子隧道电流与漏电流，I_n = V/R_n，R_n 为正常态电阻。
电容支路：承载位移电流，I_d = C dV/dt，C 为结电容。

根据第二Josephson关系，结电压与相位差满足 dφ/dt = 2eV/ħ。把三条支路电流相加并代入该关系，就得到RCSJ模型的核心电路方程：

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

这个方程在形式上与“受迫阻尼摆”完全同构：sinφ 项相当于重力回复力矩，一阶导数项相当于阻尼，二阶导数项相当于惯性。理解这一点，后面所有关于I-V曲线和开关行为的讨论都会变得直观。

二、归一化参数：把方程变成无量纲形式

直接使用有量纲的方程不便于比较不同结的性能。工程上通常引入三个归一化量：归一化电流 i = I/I_c、归一化电压 v = V/(I_cR_n)、归一化时间 τ = t/τ_J，其中 τ_J = ħ/(2eI_cR_n) 称为Josephson时间常数。代入后方程化为：

i = sinφ + dφ/dτ + β_c·d²φ/dτ²

这里出现了一个决定结行为“性格”的关键无量纲参数——McCumber阻尼参数：

β_c = (2e/ħ)·I_c·R_n²·C = ω_c·R_n·C

其中 ω_c = 2eI_cR_n/ħ 为特征频率。β_c 的物理含义是电容充放电时间常数 R_nC 与Josephson时间常数之比，它直接衡量结的“惯性”相对于“阻尼”有多强：

β_c ≪ 1（过阻尼）：电容效应可忽略，结近似为纯电阻分流，I-V曲线回滞极小，开关行为可逆。
β_c ≫ 1（欠阻尼）：电容储能显著，结在撤去偏置后仍能持续振荡，I-V曲线出现明显回滞，是数字逻辑与RSFQ电路常用的工作区。
β_c ≈ 1：临界阻尼附近，回滞刚好出现或消失的过渡区。

工程上常把 β_c 作为设计自由度：给定 I_c 与 R_n 后，通过选择结面积（决定 C）把 β_c 调到目标区间。例如某Nb/AlO_x/Nb结若 I_c = 100 μA、R_n = 50 Ω、C = 0.5 pF，则 β_c ≈ 2e/ħ × 100μA × 2500Ω² × 0.5pF ≈ 0.38，属于接近临界阻尼的轻度欠阻尼结。

三、直流I-V特性的两种形态

RCSJ模型的直流I-V特性由方程中 dφ/dτ 与 d²φ/dτ² 项是否可忽略决定，呈现两种典型形态。

过阻尼极限（β_c → 0）：方程退化为 i = sinφ + dφ/dτ。当偏置电流 i < 1 时，存在稳定的静态解 sinφ = i，结电压为零，对应I-V曲线上的零电压态；当 i > 1 时，sinφ 无法平衡偏置，φ 持续增长，结进入电压态，平均电压满足 v = √(i² − 1)。因此过阻尼结的I-V曲线是：零电压支路一直延伸到 i = 1，然后沿 v = √(i² − 1) 平滑上升，无回滞。

欠阻尼情形（β_c 较大）：当电流从零增大到 I_c 时，结仍停留在零电压态；一旦超过 I_c，结跳到电压态，且由于电容储能，即使把电流降回 I_c 以下，结仍维持在电压态，直到电流降到回滞电流 I_r 才跳回零电压态。回滞电流可由近似式估算：

I_r / I_c ≈ 4 / (π·√β_c)   （β_c ≫ 1 时适用）

例如 β_c = 4 时，I_r/I_c ≈ 0.64，意味着电流要降到临界值的六成左右才返回零电压态。这一回滞窗口正是SFQ/RSFQ逻辑中“存储”与“触发”的基础，也是参数提取时必须准确测得的量。

四、参数提取的实操路径

从实测I-V曲线反推 I_c、R_n、C 是RCSJ模型最直接的工程用途。典型流程如下：

读临界电流：在零电压态向电压态跳变的拐点读取 I_c，注意区分正负扫描方向，取平均以消除扫描速率引起的偏差。
读正常态电阻：在远高于能隙电压的电压态区段（如 V > 2Δ/e）对I-V曲线做线性拟合，斜率即 1/R_n。
读回滞电流：从电压态回扫，记录跳回零电压态时的电流 I_r，由 I_r/I_c 反解 β_c。
算电容：由 β_c = (2e/ħ)I_cR_n²C 解出 C = β_cħ/(2eI_cR_n²)。

上述提取假设结满足理想RCSJ行为。实际器件中若存在并联的额外分流电阻、势垒不均匀或测量引线电阻，会使拟合出的 R_n 偏小、β_c 偏大。建议先用多组不同结面积的器件交叉验证 C 是否与面积成正比，以判断提取结果是否可信。

把RCSJ方程、β_c 与I-V回滞三者联系起来看，就能理解为什么同一个“结”在不同偏置条件下会表现出截然不同的开关特性：本质上是阻尼与惯性的竞争。掌握这一竞争关系，后续讨论结的高频响应、噪声引起的误翻转以及参数提取的误差来源时，都能回到这个统一的物理图像上。