4. 量子芯片电磁寄生参数提取的数值方法选型

量子芯片的互连结构横跨多个数量级：超导量子比特的约瑟夫森结与电容垫片在微米量级，共面波导传输线延伸到毫米量级，而封装引线与控制线又可能达到厘米量级。要在同一套模型中同时捕捉结区附近的准静态耦合与长线上的分布谐振，单一数值方法往往难以兼顾精度与效率。本节围绕有限元法（FEM）、矩量法（MoM）与部分元等效电路法（PEEC）三种主流方法，讨论它们在多尺度结构中的适用边界，并给出网格、边界与材料色散设置对提取精度的实际影响，最后形成按互连尺度选型的方法准则。

三种方法的物理建模视角差异

FEM 在体域内离散麦克斯韦方程，通过四面体或六面体网格求解电场或磁场的边值问题，天然适合处理非均匀介质、复杂三维几何与色散材料。它的未知量分布在体积内，网格量随结构尺寸的三次方增长，因此对毫米级以上的长线结构代价很高，但对微米级结区、空气桥、倒装焊球附近的场集中区域精度最好。

MoM 基于表面电流积分方程，只对导体表面与介质界面离散，未知量随表面积增长，特别适合开放边界下的辐射与长线耦合问题。由于格林函数已解析满足辐射条件，MoM 不需要人为截断外边界，在提取共面波导之间的互耦、封装腔体内的模式耦合时优势明显。但其系统矩阵稠密，直接求解复杂度约为 O(N³)，对含大量细小特征的量子芯片版图需要配合多层快速多极子或压缩算法。

PEEC 从部分电感、部分电容与部分电阻的电路等效出发，把导体离散为细丝或面片，生成可读入 SPICE 的集总或分布网络。它最贴近电路设计流程，能直接与约瑟夫森结模型、控制线阻抗、滤波器网络联合仿真，适合提取从直流到数 GHz 的互连寄生。PEEC 的准静态假设在结构尺寸接近波长十分之一时开始失效，因此对超导量子芯片中常见的 4–8 GHz 工作频段，毫米级长线需要谨慎评估延迟与辐射效应。

网格划分：精度与代价的平衡点

FEM 的网格质量直接决定寄生电容的收敛性。以超导量子比特的十字形电容垫片为例，垫片边缘的 fringe 场贡献了约 20%–30% 的对地电容，若网格在边缘处过粗，提取值会系统性偏低。实践中可采用自适应网格：先以 λ/20 的全局尺寸生成初始网格，再在电场能量密度梯度超过阈值的区域做 h 型加密，直到相邻两次迭代的电容变化小于 0.5%。对于 5 μm 宽的共面波导间隙，建议在间隙内至少布置 4–6 层网格，并在导体表面设置边界层单元以捕捉趋肤效应。

MoM 的网格关注的是表面电流的相位变化。导体表面剖分尺寸通常取 λ/10 以内，但对超导薄膜，还需考虑伦敦穿透深度 λ_L（铌约 40 nm，氮化铌约 200 nm）。当薄膜厚度与 λ_L 可比时，动生电感不可忽略，此时应在表面阻抗边界条件中引入复电导率，而不是单纯加密网格。PEEC 的细分则取决于最高关注频率对应的波长：若目标是 8 GHz，自由空间波长约 37.5 mm，在有效介电常数约 6 的衬底中约 15 mm，细丝段长度宜控制在 1 mm 以下，以保证每段延迟远小于周期。

边界条件：开放与封闭的取舍

量子芯片的电磁环境本质上是开放的，但 FEM 必须在有限计算域内截断。若把外边界设为理想电导体（PEC），会在封装腔体内引入虚假谐振，使 4–8 GHz 频段的传输参数出现不属于真实结构的尖峰。正确做法是使用理想匹配层（PML）或辐射边界，并把外边界置于距芯片表面至少 λ/4 处；对倒装封装结构，还需把封装盖板与焊球一并建模，否则盖板与芯片之间的腔体模式会被遗漏。MoM 由于积分方程本身满足辐射条件，只需在导体与介质界面离散，不存在截断误差，这也是它在开放结构互耦提取中更可靠的原因。PEEC 的边界处理体现在参考地的选取上：控制线若以封装地为回流路径，必须把地平面纳入部分电容网络，否则共模回流路径缺失，提取的串扰会偏小。

材料色散：超导芯片不可回避的修正

超导材料的表面阻抗具有频率依赖的复部，动生电感 L_k 与穿透深度相关。在 FEM 中可通过频变复电导率 σ(ω) = σ₁ − jσ₂ 实现，其中 σ₂ 对应库珀对的无耗响应；在 PEEC 中则把 L_k 作为附加的部分电感串联到对应细丝段。若忽略色散，在 6 GHz 以上提取的传输线特性阻抗会偏高约 3%–8%，相位常数偏差更明显。介质方面，硅衬底与 SiO₂ 界面层的损耗角正切也随频率变化，若只取单一常数，品质因数的提取会失真。建议在 1–10 GHz 内至少取 5 个频点做扫频，并对 L_k 与 tanδ 分别做线性或对数插值。

按互连尺度的方法选择准则

需要强调的是，上述判据中的数值（如 0.5%、λ/10、1 mm）属于工程经验参考值，实际项目中应结合具体版图、衬底参数与目标频段做收敛性测试后确定。方法选型不是一次性的，当工作频率上移或结构尺寸变化时，原先适用的准静态假设可能失效，应重新评估 PEEC 的适用性，或把长线部分交给 MoM、把结区交给 FEM 做分区提取。