4. 电路建模基础与网表构建

超导量子比特耦合器的建模，最终要落到一组可以被仿真器读取的方程和一份可以被解析的网表上。本节从电路理论的两条基本定律出发，说明节点电压法与支路电流法在超导电路中的适用方式，再给出超导电路量子化的标准流程，最后落到 SPICE 网表的语法与元件定义，使读者能够独立写出一个可运行的耦合器电路网表。

节点电压法与支路电流法

节点电压法以节点电位为未知量，对每个独立节点列写基尔霍夫电流定律（KCL）：流入节点的电流代数和为零。对于含 N 个节点的连通电路，选定一个参考节点（地）后，可得到 N−1 个独立方程。写成矩阵形式即 Y·V = I，其中 Y 为节点导纳矩阵，V 为节点电压向量，I 为注入电流向量。节点电压法的优点是方程数与节点数对应，便于程序化生成，SPICE 类仿真器的改进节点分析（MNA）正是它的扩展形式。

支路电流法以支路电流为未知量，对每个独立回路列写基尔霍夫电压定律（KVL），再配合节点的 KCL 构成完备方程组。它的物理图像更直观：每条支路的电流直接对应一个物理量，适合支路数少、回路结构清晰的电路。对于超导量子电路，支路电流法还有一个额外好处——超导环中的环流本身就是量子化变量，用支路电流描述更贴近后续的哈密顿量构造。

两种方法在数学上等价，选择依据是电路拓扑与后续处理需求。节点数明显少于支路数时用节点电压法；需要显式追踪环流、构造磁通变量时用支路电流法。实际建模中常混合使用：先用节点电压法得到导纳矩阵，再把电感支路电流作为附加状态变量引入。

超导电路量子化流程

超导电路量子化的目标是把一个经典电路变成量子哈密顿量。标准流程可以概括为四步：

1. 确定电路拓扑与节点：标出所有节点、参考地、约瑟夫森结、电感与电容，明确每个元件连接在哪两个节点之间。
2. 选取广义坐标：通常取节点磁通 Φ_n（满足 V_n = dΦ_n/dt）作为广义坐标，或取支路电流作为坐标。约瑟夫森结的相位差与节点磁通直接相关：φ = 2πΦ/Φ₀，其中 Φ₀ = h/2e 为磁通量子。
3. 构造拉格朗日量：动能项来自电容储能 T = ½C(dΦ/dt)²，势能项来自电感储能 ½LI² 与约瑟夫森能量 −E_Jcos φ。写出 L = T − U，再由欧拉–拉格朗日方程得到运动方程。
4. 正则量子化：定义共轭动量 Q_n = ∂L/∂(dΦ_n/dt)，把 Φ 与 Q 提升为算符并施加对易关系 [Φ_n, Q_m] = iħδ_nm，得到量子哈密顿量。

以一个磁通耦合器为例：两个超导量子比特通过一个耦合电感 L_c 相连，每个比特用一个约瑟夫森结（临界电流 I_c）与结电容 C_J 并联近似。设节点磁通为 Φ₁、Φ₂，耦合支路电流为 I_c，则拉格朗日量可写为

L = ½C_J(Φ̇₁² + Φ̇₂²) + ½L_c I_c²
    + E_J1·cos(2πΦ₁/Φ₀) + E_J2·cos(2πΦ₂/Φ₀)
    − (Φ₁ − Φ₂ − L_c I_c)²/(2L_c)

这里 E_J = I_cΦ₀/2π。对 I_c 求偏导并消去该辅助变量，即可得到只含 Φ₁、Φ₂ 的有效哈密顿量，耦合项表现为 Φ₁Φ₂/(2L_c) 的形式。这一步是耦合器建模的核心：耦合强度由 L_c 决定，而比特非线性由 E_J 决定。

SPICE 网表语法与元件定义

SPICE 网表是纯文本，每行一条语句，以元件首字母标识类型，节点用整数或字符串编号，第一个节点通常约定为参考地（0 号节点）。超导电路建模常用的元件定义如下：

约瑟夫森结在标准 SPICE 中没有原生元件，通常用两种方式实现：一是用行为电流源 B 元件写出 I = I_c·sin(2πΦ/Φ₀)，配合一个积分环节得到相位；二是定义子电路，内部用非线性电阻或压控电流源近似。下面给出一个可直接运行的耦合器网表示例，两个比特通过电感耦合，参数为教学假设值：

* 超导量子比特耦合器示例网表
.param ic=3u cj=5f lc=2n
.subckt qubit n1 n0
CJ  n1 n0 {cj}
BJ  n1 n0 jjmod ic={ic}
.ends
XQ1 1 0 qubit
XQ2 2 0 qubit
LC  1 2 {lc}
IB1 0 1 dc 0.5u
IB2 0 2 dc 0.5u
.model jjmod jj(ic={ic})
.tran 0.01n 50n
.end

网表构建时有几个容易出错的点需要检查：节点编号必须与电路图一一对应，参考地统一用 0；单位后缀区分大小写，f 为飞、n 为纳、u 为微，M 与 m 含义不同；子电路调用时端口顺序必须与 .subckt 声明一致；行为源的表达式要保证量纲正确，相位以弧度为单位。

从网表到可验证的建模结果

网表写完后，判断建模是否正确的依据不是“仿真跑通了”，而是几个可检查的量。第一，检查直流工作点：在给定偏置电流下，约瑟夫森结两端的相位应满足 I_bias = I_c·sin φ，若偏置 0.5 μA、I_c = 3 μA，则 sin φ ≈ 0.167，φ ≈ 0.167 rad，可用 .op 语句输出节点电压核对。第二，检查瞬态波形：耦合电感两端电压应随两比特磁通差变化，若两比特初始条件相同，耦合支路电流应保持为零，这是对称性检验。第三，检查参数扫描：改变 L_c 从 1 nH 到 5 nH，观察耦合引起的频率劈裂是否随 L_c 单调变化，以此确认耦合项符号与量级正确。

需要强调的是，上述参数均为教学示例假设值，实际器件的临界电流、结电容与耦合电感需由工艺参数和版图提取得到。建模的价值在于：一旦网表结构正确，替换参数即可快速评估不同耦合器设计，而不必每次从方程重新推导。节点电压法给出方程骨架，支路电流法给出物理变量，量子化流程给出哈密顿量，SPICE 网表则把这一切固化为可复现、可版本管理的文本资产——这四者共同构成超导耦合器电路建模的完整链条。