3. 量子芯片版图与材料参数准备

在量子芯片电磁仿真的完整工作流中,版图与材料参数准备是决定后续仿真可信度的第一道关口。版图决定了“几何在哪里”,材料参数决定了“电磁波如何传播与损耗”,两者缺一不可。本节围绕四个关键动作展开:版图导入、层定义、超导薄膜与介质材料电磁参数设定、界面损耗初步设定,并给出可直接参考的参数配置示例。

一、版图导入:从设计文件到仿真几何

量子芯片版图通常以 GDSII(.gds/.gdsii)格式交付,部分团队也会使用 OASIS(.oas)或 DXF。导入时需重点关注三件事:单位、坐标系与布尔运算。GDSII 默认数据库单位为 1 nm,而仿真软件常以 μm 或 m 为内部单位,若单位换算错误,一个 5 μm 宽的传输线可能被识别为 5 nm 或 5 mm,直接导致仿真失效。建议在导入后立即用“测量工具”核对已知结构尺寸,例如读取一个设计值为 10 μm 的共面波导中心导体宽度,确认读回值一致。

坐标系方面,需确认版图原点与仿真坐标原点是否对齐。若版图原点位于芯片左下角,而仿真中把原点设在芯片中心,则所有端口位置都会偏移。常见做法是在导入时指定一个参考点(如芯片几何中心或某个标记点),并记录该点在版图中的坐标,后续端口设置均以此为基准。

布尔运算用于处理版图中重叠或相邻的图形。例如,超导薄膜层与刻蚀窗口层可能部分重叠,导入后需执行“差集”或“合并”操作,得到实际金属区域。若忽略这一步,仿真中可能出现金属与真空重叠的非法几何,导致网格划分失败或结果异常。

二、层定义:把版图层映射为物理对象

GDSII 中的每一层(layer)只是一个编号,本身没有物理含义。层定义的任务就是建立“版图层号 → 物理对象”的映射关系。以典型的超导量子芯片为例,常见层包括:

  • Layer 1:铌(Nb)或铝(Al)超导薄膜,用于传输线、谐振腔、量子比特电容/电感结构;
  • Layer 2:约瑟夫森结区域,通常由 Al/AlOx/Al 叠层构成;
  • Layer 3:空气桥或跨线层,用于抑制寄生槽线模式;
  • Layer 4:刻蚀窗口或介质开孔层,用于定义介质去除区域。

在仿真软件中,每一层需要指定材料、厚度和优先级。厚度方向通常沿 z 轴,需与后续端口和边界条件的方向约定一致。例如,若规定 z 轴向上为正,则超导薄膜厚度可设为 200 nm,位于 z=0 到 z=0.2 μm 之间。层优先级决定重叠区域的归属,一般金属层优先级高于介质层,刻蚀窗口优先级高于金属层。

三、超导薄膜与介质材料电磁参数

超导薄膜在微波频段的电磁行为不能简单用理想导体近似。在远低于临界温度(如 10 mK)时,超导体的表面阻抗由两部分组成:

Z_s = R_s + jωL_k

其中 R_s 为表面电阻,反映残余损耗;L_k 为动态电感,来源于库珀对的质量惯性。对于铌薄膜,在 5 GHz 下典型表面电阻约为 1–10 μΩ,动态电感约为 0.1–0.5 pH/□。这些参数直接影响谐振频率和品质因数,必须在材料设置中显式输入。

介质材料方面,量子芯片常用的有硅(Si)、蓝宝石(Al₂O₃)、二氧化硅(SiO₂)和氮化硅(SiN)。关键电磁参数包括相对介电常数 ε_r 和损耗角正切 tanδ。示例参数如下:

材料ε_r(示例)tanδ(示例,10 mK)典型用途
高阻硅11.71×10⁻⁶衬底
蓝宝石9.3(各向异性)5×10⁻⁷衬底
SiO₂3.91×10⁻³介质层/界面
SiN7.02×10⁻⁴钝化层

注意:上述数值为教学示例,实际值随薄膜生长工艺、温度和后处理变化,应以测量或文献数据为准。蓝宝石的介电常数具有各向异性,ε_r 在 c 轴和 a 轴方向不同,仿真中若使用各向同性近似,需在报告中注明。

四、界面损耗的初步设定

量子芯片的主要损耗来源往往不是体材料,而是界面。典型界面包括:金属-衬底界面、金属-空气界面、衬底-空气界面。这些界面处存在几纳米厚的非晶层或氧化层,其损耗角正切可比体材料高 2–3 个数量级。

初步设定时,可采用“薄层等效”方法:在金属与衬底之间插入一层厚度 2–5 nm、ε_r≈4、tanδ≈1×10⁻³ 的界面层。该层不参与网格细化,而是通过表面阻抗边界条件或等效介质参数引入损耗。另一种做法是直接指定界面损耗因子,例如在仿真软件中设置“surface loss tangent”,典型初始值取 1×10⁻³ 至 5×10⁻³。

界面损耗的初步设定不必追求一次精确。建议先以保守值(如 tanδ=2×10⁻³)运行一次本征模仿真,观察品质因数 Q 的量级。若 Q 远高于实验值,再逐步增大界面损耗;若 Q 远低于实验值,则检查是否误将体材料损耗设得过高。

五、参数配置示例与检查清单

以下给出一个可直接参考的配置片段(以某仿真工具的参数文件风格为例):

# 层定义
layer 1 name=Nb  material=niobium  thickness=0.2um  priority=10
layer 2 name=JJ  material=alox     thickness=0.002um priority=20
layer 3 name=AirBridge material=vacuum thickness=1.0um priority=5
layer 4 name=EtchWindow material=vacuum thickness=0.2um priority=30

# 超导薄膜参数
niobium: Rs=5uOhm  Lk=0.3pH_per_square  Tc=9.2K

# 介质参数
silicon: epsilon=11.7  tanD=1e-6
sapphire: epsilon=9.3  tanD=5e-7  anisotropy=uniaxial
sio2: epsilon=3.9  tanD=1e-3

# 界面损耗
interface metal_substrate: thickness=3nm epsilon=4.0 tanD=2e-3
interface substrate_air: thickness=2nm epsilon=1.0 tanD=1e-3

完成上述设置后,建议按以下清单逐项核对:

  1. 版图单位与仿真单位是否一致?用已知尺寸结构验证。
  2. 所有 GDS 层是否都已映射到物理对象?有无遗漏层?
  3. 超导薄膜的 R_s 和 L_k 是否已输入?是否与工作频率和温度匹配?
  4. 介质材料的 ε_r 和 tanδ 是否已设置?各向异性材料是否指定了方向?
  5. 界面损耗层是否已插入?其厚度和 tanδ 是否在合理范围?
  6. 层优先级是否会导致重叠区域归属错误?
量子芯片版图与材料参数准备流程 ① 版图导入 GDSII / OASIS 单位·坐标系·布尔运算 ② 层定义 Layer→物理对象映射 材料·厚度·优先级 ③ 材料电磁参数 Rs·Lk·εr·tanδ 超导薄膜+介质 ④ 界面损耗 薄层等效 tanδ 初设 超导薄膜参数(示例) 铌(Nb)薄膜: 表面电阻 Rs ≈ 5 μΩ @ 5 GHz 动态电感 Lk ≈ 0.3 pH/□ 临界温度 Tc ≈ 9.2 K 介质材料参数(示例) 高阻硅:εr=11.7,tanδ=1×10⁻⁶ 蓝宝石:εr=9.3(各向异性),tanδ=5×10⁻⁷ SiO₂:εr=3.9,tanδ=1×10⁻³ SiN:εr=7.0,tanδ=2×10⁻⁴ 界面损耗初步设定(薄层等效) 超导金属(Nb) 界面层 2–5 nm,tanδ≈2×10⁻³ 衬底(Si / Sapphire) 界面层不参与网格细化, 通过表面阻抗边界或等效介质引入损耗。 初始值建议:tanδ = 1×10⁻³ ~ 5×10⁻³。

完成以上四步后,版图与材料参数准备即告一段落。此时应保存一份参数快照,记录所有层映射、材料参数和界面损耗设定,以便后续仿真结果出现异常时回溯。下一阶段将进入端口设置与网格划分,但前提是本节的所有参数已通过检查清单验证。