INTERCONNECT建模选型：紧凑模型与S参数模型的精度速度权衡

做光子集成的人几乎都绕不开这个问题INTERCONNECT里同一个模块你既能在元件库里拖一个理想模型跑通整条链路又愿意花几个小时去FDTD里把S参数磨出来再导进去用。前者快得让迭代停不下来后者准得让结果拿得出手但两条路之间的差距到底体现在哪里什么阶段该选哪边很多人其实没有系统梳理过。这篇文章就把 INTERCONNECT 的这两套建模方案放在一张台面上对比从各自原理、精度边界、耗时瓶颈到实际项目里怎么做取舍全都过一遍。无论你是刚开始接触 PIC 系统仿真的学生还是在做PDK开发、硅光芯片验证的工程师这篇内容应该都能帮你少走几步弯路。1. INTERCONNECT里的两种建模路径分别对应什么场景1.1 解析/紧凑模型链路架构阶段的主角INTERCONNECT 自带的标准元件库里绝大多数器件都是所谓的“解析模型”或者叫“紧凑模型”。这些模型不是从格子空间里求解麦克斯韦方程得到的而是基于器件在特定工作状态下的传输函数或者等效电路参数建立的。举个例子你在元件库里拖一个 Ring Resonator 元件出来它背后的物理本质是环腔耦合理论通常用自耦合系数、交叉耦合系数、环绕损耗、群折射率这几个参数就能描述。你设置半径、耦合系数、损耗它就能算出一套谐振谱线。同样一个 Directional Coupler 的模型也可以简化成两个传输系数 t 和 kMZM 调制器则用半波电压 Vπ 和消光比 ER 概括。这套模型的价值在于它把器件最核心的物理行为抽取成了几个互不干涉的参数数学形式简单频域、时域求解都非常快。你做一条包含几十个器件的复杂 WDM 链路用紧凑模型做一次时域仿真往往在毫秒到秒级改参数重跑几乎零成本。所以在项目最早期当你要验证链路架构是否合理、DSP算法能不能收敛、预算链路预算够不够的时候紧凑模型是唯一现实的选择。1.2 电磁提取的S参数模型性能收敛阶段的裁判第二套建模方案是走“物理级仿真 S参数导入”的路线。你先在 FDTD、MODE/Eigensolver、EME 这类器件级工具里精细建模得到器件在目标波长范围内的 S 参数矩阵然后以 Touchstone 文件.s2p、.s4p 之类或者 INTERCONNECT 自有的模型格式导入到系统级环境里作为一个黑盒器件参与后面的链路仿真。这套模型的精度来源是真实的电磁场求解。它不再需要你去猜耦合系数是多少因为 FDTD 已经把模式重叠积分、频散关系、模式失配损耗、弯曲辐射损耗这些东西全部算进去了。你关心的谐振波长偏移、串扰、反射、群延迟纹波都以幅度和相位信息的形式完整保存在 S 参数里。代价也很直观每一个器件都要先经过一轮甚至多轮电磁仿真才能拿到 S 参数而且在 INTERCONNECT 里一旦要做参数扫描最笨的办法就是所有参数组合全部重新跑一遍 FDTD。时间是分钟到小时级别地往上翻。1.3 两套模型在软件里的真实存在形式很多人第一次接触 INTERCONNECT 时容易混淆这两套模型是因为它们看起来“都是元件库里的一个图标”。这里需要说清楚紧凑模型是 INTERCONNECT 原生元件直接由 Internal 脚本或者内置 CMLCompact Model Library紧凑模型库支持参数面板上就是一堆数值输入框不需要外部文件。S参数模型通常以 GDS/物理结构为源头在 FDTD 后处理中导出为 Touchstone 文件然后在 INTERCONNECT 里用 S Parameter Model 元件导入。导入后它在功能上也是一个黑盒但内部数据是采样的频率响应表。还有一个容易搞混的选项是 co-simulation即把 FDTD 求解器动态挂到 INTERCONNECT 时域仿真里。这个是在两套方案之间的一种“实时物理验证”模式原理上不用提前准备 S 参数文件但每个时间步里都会去调用电磁仿真速度通常是灾难级别的工程上一般只在关键节点做少量验证用极少作为日常迭代工具。2. 紧凑模型到底牺牲了什么精度2.1 理想化假设一模式不匹配被完全忽略紧凑模型最大的问题不是“近似”而是它默认器件工作在某个理想工作点上。比如方向耦合器模型里它假设两根波导是平行且无限长的端口的模式完全对准耦合效率只由间距和长度决定。但实际版图里波导进耦合区之前可能有一个弯曲弯曲outgassing会把模式推向弯曲外侧导致进入耦合区时模式的横向分布已经不是标准的对称模式实际耦合强度会偏离解析值。这种模式不匹配在 TE0 到 TE0 的简单结构里可能只有几个百分点的影响但一旦涉及模式转换器、绝热锥形或者多模干涉结构误差会迅速加大。你如果用紧凑模型去预测 MMI 的额外损耗经常会出现“仿真说 0.2 dB测试出来 0.8 dB”这种落差不是因为你参数设置错了而是紧凑模型根本不具备表达模式重叠积分的维度。2.2 理想化假设二损耗被等效成固定常数紧凑模型里的波导损耗通常是一个单点数值比如 2 dB/cm。你设置好了之后它在整个波长范围内都用这个值去算。但实际上波导损耗是随波长变化的侧壁粗糙度导致的散射损耗在短波长更恶劣材料吸收边在长波方向也更恶劣弯曲波导的辐射损耗更是跟弯曲半径强相关。如果你在做的是一个宽带超低损耗传输链路或者多个弯曲半径混用的版图用固定损耗常数去预估整条链路的插损误差会一层层叠加。最典型的情况是紧凑模型算出来链路总损耗 4.3 dBS参数模型算出来 5.1 dB测试 5.4 dB。这个 0.8 dB 的差异对于系统链路预算来说可能刚好决定接收端灵敏度够不够。2.3 实际案例环形谐振器在两种模型下的谐振偏移我用一个实际的微环例子来说明偏差量级。假设直波导宽度 0.45 μm包层为氮化硅微环半径 5 μm耦合间隙 0.2 μm。紧凑模型里我设置群折射率 ng 2.1环周损耗 0.01 dB/round耦合系数按 0.3 左右给。跑出来的谐振波长在 1550.5 nmFSR 约 19 nm。随后我在 FDTD 里以 5 nm 的最小网格细化耦合区跑出 S 参数再导入 INTERCONNECT同样的结构得到的谐振波长变成了 1551.2 nm。一个 0.7 nm 的偏移来源于什么首先是弯曲区域模式频散比等效 ng 更复杂其次耦合区并不是“点耦合”而是有实际长度的分布式耦合它引入了一个额外的相位项。更头疼的是 FDTD 结果里还看到谐振谷的消光比只有 18 dB而紧凑模型默认能算出 35 dB。差距来自直通波导Taper与环之间模式失配带来的额外损耗以及耦合区高阶模的参与。这个案例很典型你用紧凑模型做出来的调谐曲线形状是对的趋势是对的但绝对频率位置、消光比、串联多个环之后的通带波纹这些细节全部不能直接用于最终设计。2.4 哪些场景可以继续放心用紧凑模型说了这么多紧凑模型的缺点并不是要否定它。至少有三类场景我依然会在项目里坚持用紧凑模型链路级架构探索。比如你先要在发射端、调制器、无源网络、接收端之间做功耗预算和带宽预算器件本身的非理想因素可以用“修正后的经验值”代替没必要上 S 参数。DSP 和算法开发。数字信号处理工程师需要的是一个能持续跑几十万比特高速信号的发射链路模型他们关心的是眼图质量趋势不是微环谐振波长 0.1 nm 的偏移。系统级蒙特卡洛分析。你给紧凑模型的每个参数加一个高斯分布扰动跑一千次系统仿真来找敏感性这种场景用 S 参数模型跑经济上完全不可行。关键是你要清楚什么时候紧凑模型的误差不会影响结论方向什么时候它会把你的结论带偏。3. S参数模型的精度从哪里来时间又花到哪里去3.1 S参数里藏着电磁场的完整因果链如果只看文件格式Touchstone 文件就是你扫描频率点上的一组复数矩阵没什么神秘的。它的优势在于S参数背后的数据是由电磁场数值解严格推导出来的里面自然嵌入了器件内部所有的模式演化过程。比如通过 FDTD 得到的方向耦合器 S 参数它已经包含了耦合区端面反射、模式转换、波长相关的耦合相位你在 INTERCONNECT 里用它做系统仿真时这些效应会自动参与计算不需要你额外建模。这里特别值得留意的是相位信息。紧凑模型计算群延迟往往靠解析色散公式而 S 参数里的群延迟来自真实的相位—频率关系。对滤波器设计、色散补偿、相干通信里的 IQ 时延匹配来说这个差异会直接决定链路性能评估是否可靠。不过“精度高”不等于“买了保险”S参数模型同样有重大的前提约束它只在导出 S 参数时的那个结构、那个波长范围里有效。你改了一个几何尺寸、换了一种包层材料、或者把单根波导靠近到另一个器件附近原来的 S 参数就不成立了。这就是为什么很多人说 S 参数模型的“泛化能力”差。3.2 时间开销一网格与仿真域要“甜点尺寸”S 参数模型的第一个时间瓶颈在器件电磁仿真。FDTD 的网格越细计算精度越高但时间步长也随之变小网格数量增大内存和 CPU 时间同步上升。问题是很多器件里真正需要细网格的区域并不大比如耦合间隙 200 nm 处可能需要 10 nm 甚至更细的网格而其他大部分波导区域 20–50 nm 的网格就足够。实操中的做法是开启 FDTD 的 mesh override region只在耦合区、尖锐拐角、模式剧烈变化的位置放置细网格其余区域用全局 mesh。这样做可以把 3D 仿真时间从几个小时压到几十分钟代价是你要对器件的物理场分布有一点预判。你如果整个仿真域都用 10 nm 网格一个半径 5 μm 的微环跑一次 1550 nm 波段仿真普通工作站可能连内存都吃紧。还需要权衡仿真区域大小。区域太小会把边界反射和模式截断误差引进来区域太大则白白增加网格数。一般建议在模式衰落到可以忽略的位置上设置 PML 边界留出至少半个波长的缓冲距离。这本身就是经验活只能靠一轮轮试出自己的“甜点尺寸”。3.3 时间开销二高Q结构需要长时域收敛FDTD 在 INTERCONNECT 建模链路里最“磨人”的一点是谐振结构的收敛时间。宽谱脉冲激励进入一个 Q 值很高的微环滤波器后腔内能量衰减非常慢你需要让仿真一直运行到瞬态能量衰减到足够低的水平否则提取的 S 参数会有明显抖动和纹波特别是谐振峰附近的谱线会“毛刺”得没法用。高 Q 微环Q 2万的时域收敛往往要跑到 100 万到几百万个时间步。虽然每一步的计算量不大但累计下来单次 FDTD 仿真轻松突破半小时。这在只算一个工作点的时候尚可接受一旦要做波长扫描或几何扫描次数会迅速膨胀。所以工程上遇到高 Q 结构我不会无脑迷信全波段 FDTD有时会改用 EME 或者频域本征模方法做分段的等效提取再在 INTERCONNECT 里手动拼接模型。3.4 时间开销三参数扫描会让成本成倍放大S 参数模型在 INTERCONNECT 里的性能瓶颈绝不只在建模阶段。很多人纠结“为什么我的微环参数扫描跑了一天还没出来”原因往往是每扫一个耦合间隙值就要重新跑一次 FDTD 再重新导入。FDTD 与 INTERCONNECT 之间是“离线”关系S 参数不会自动按几何差补。所以面对参数扫描必须预先规划采样矩阵而不是一次性把所有扫描目标和 S 参数生成绑在一起。比如你想扫耦合间隙 150/175/200/225/250 nm 五档不要在 INTERCONNECT 里做循环调用 FDTD先用 FDTD 并行把五个 S 参数文件导出来再一次性导入到 INTERCONNECT 里做参数化比较。一个文件一个文件地流水线跑浪费时间还容易在软件切换过程中出错。4. 精度速度权衡的实操策略我平时怎么选4.1 第一层按项目阶段分层我自己的项目流大致分成三个阶段每个阶段用的模型方案完全不同前端概念验证阶段所有器件都用紧凑模型目的只有一个把全链路跑通确认没有方案级错误确定系统参数的量级。这时候我不会花太多时间微调某个器件的精确响应因为架构一变大部分参数都要推倒重来。方案细化阶段开始对链路里最核心的几个器件做电磁提取。比如发射端的微环调制器、接收端的高阶滤波器这几个元件的 S 参数会直接决定系统能不能达标值得投入时间做高精度建模。其余外围器件继续保留紧凑模型。这个混合状态是 INTERCONNECT 项目里最常见的状态。最终验证阶段把所有关键部件的 S 参数更新到最终版版图对应的结果跑几轮端到端仿真并与测试数据对标。如果偏差大回头检查器件模型与版图一致性而不是直接在系统级仿真里加魔改系数。4.2 第二层混合建模只对关键器件开S参数混合建模听起来简单实际操作中最需要小心的是接口一致性。S 参数模型端口顺序、模式定义要和紧凑模型里的“等效端口”对齐否则你辛辛苦苦把微环模型换成 S 参数模型却发现链路通带变了——多半不是模型精度问题而是端口映射错了。我习惯在建 S 参数模型时先给它的端口命名和紧凑模型的端口一一对应比如 TA1、TA2、TB1、TB2。不要用一堆数字编号否则跨文件对接时根本分不清楚哪个是输入哪个是输出。再者混合模型里 S 参数器件两端连接的波导段最好额外加一小段“对接波导”来补偿端口参考面的位置。FDTD 提取 S 参数时参考面一般设置在波导末端但 INTERCONNECT 里元件之间的连线默认是零长度理想连接。如果你不把参考面位置偏差补回来测到的谐振波长会发生系统性偏移。这个细节最容易踩坑但又很少有人写进文档。4.3 第三层参数扫描前先做敏感性预筛选前面提到参数扫描会和 FDTD 次数绑定泛泛地说“我要对全部几何参数做扫描”是非常不理智的。一个微环相关几何参数可能就有环半径、波导宽度、波导厚度、耦合间隙、Taper长度五个参数各扫五个水平就是 3125 次 FDTD想都不要想。实际操作是先做一个敏感性预筛选用紧凑模型快速跑一遍 Sobol 或者拉丁超立方采样算清楚哪些参数对目标性能影响大。如果结果显示耦合间隙的影响占了 80%其他参数几乎无关紧要就让 FDTD 只扫耦合间隙其他参数用设计标称值。这样可以把 3125 次降到 5 次总计算量完全可控。这一步其实是用“数据驱动”的方式避免在哪个模型上都吃力不讨好我几乎每个项目都这么做。4.4 一个典型的PDK开发工作流如果是做硅光 PDK我个人推荐的固定工作流长这样用紧凑模型定义 PDK 标准元件的行为级描述定端口、定可调参数、定性能指标基线。对每个标准元件做几组电磁提取覆盖典型工艺角的边界生成 S 参数库。对提取的 S 参数做后处理比如平滑、插值、因果性校准再生成 INTERCONNECT 可加载的库文件。在 INTERCONNECT 里建一个“验证原理图”把库里的 S 参数模型和相应紧凑模型放在同一个测试台一键对比频响差异差异超过阈值就回头检查。最后一轮用 S 参数库跑系统验证甚至做小规模蒙特卡洛但只在敏感性筛选后保留的少数参数上做。这个流程不便宜前期投入很大但对后续每个定制芯片的验证都有复用价值。很多团队觉得 INTERCONNECT 精度和速度难以兼得其实是因为他们构建模型库的方式是“每次项目从零开始做S参数”从来没有积累和沉淀模型资产。5. 实测工程笔记采样、插值与模型管理5.1 频域采样点数怎么定才不被插值坑S 参数导入 INTERCONNECT 后并不会直接按原始文件逐点计算系统会做插值。插值算法的前提是原始数据足够密集能够还原真实的频率响应形态。对于透射峰很尖锐的滤波器如果采样点太稀疏插值会削掉峰顶甚至产生非物理的负群延迟。我的经验是保证每个谐振峰的 FWHM半高全宽内至少要有 4 到 8 个频点。比如 Q 值 1 万、波长 1550 nm 的谐振峰FWHM 大约 0.155 nm那么在 1540–1560 nm 扫描范围内至少要点 1000 个以上频点才安全。FDTD 输出 S 参数时可以自动按频点插值但别过度依赖默认设置一定要自己确认这些关键谐振峰的采样覆盖率。还有一个容易被忽略的问题频带边界效应。S 参数文件只覆盖你设定的波长范围超出这个范围后 INTERCONNECT 会外推外推通常是常数或线性延续极容易失真。如果你要在时域仿真里看短脉冲的传播频谱分量可能超出 S 参数覆盖范围结果会产生振铃和鬼影。所以时域仿真之前要有意识地确认信号带宽是否完全落在 S 参数的覆盖区间内。5.2 Touchstone文件的相位连续性与导出口径Touchstone 文件里的相位是以 [-180°, 180°] 或 [0°, 360°] 形式表示的在不同仿真器之间传输时相位卷绕可能导致插值错乱尤其在宽频带高色散器件中很明显。推荐导出时选择 unwrapped 相位或者至少保持和 INTERCONNECT 设置一致的相位口径。另外不同 FDTD 工具对端口阻抗、参考阻抗的定义也可能不同。标准 Touchstone 默认是 50Ω但波导端口的“阻抗”在光学语境下并不是电路里的 50Ω而是一个与模式功率和电压定义相关的等效值。如果参考阻抗定义不一致导入的 S 参数在级联计算中可能出现莫名其妙的反射。我一般在导出时把参考阻抗设为 1也就是归一化功率口避免不匹配问题。如果你拿到的是别人给的 S 参数文件第一件事就是用 INTERCONNECT 里自带的 RLC/Scope 观察一下 S 参数本身的回损和插损曲线确认没有异常跳跃和明显非因果性再继续往下做系统仿真。这一步 5 分钟能省下后面几小时的排查时间。5.3 模型版本管理和跨平台对接的建议S 参数模型本质上是文件这就不可避免地要面对版本管理问题。我给身边团队的建议是S 参数文件的命名必须包含器件名、几何尺寸、工艺角、提取工具、日期五个要素。比如“ring_r5_gap200_wpda_20250115.s4p”。缺乏这五个要素的模型三个月后没人能说清它是什么状态。在 INTERCONNECT 里用 S 参数模型做项目时建议和原始 FDTD 项目文件放在同一个工程目录下而不是复制散落在各个临时文件夹。因为跨平台对接时很多人会把 .s4p 文件拷走却忘了拷贝对应的 FDTD 源文件等回头想更新几何尺寸或网格精度时完全找不到源头。为保证模型可追溯我甚至会把 FDTD 项目的关键截图和网格设置记录一并归档防止以后返工。还有一个实用技巧S 参数库在 INTERCONNECT 里不要直接给几百个器件全部建 S Parameter Model 元件那样原理图会乱到没法看。建议把同类型、同工艺的 S 参数库封装成自定义元件在 CML 层面做好参数映射这样原理图上的表示仍然是拖一个“微环”图标但内部实现是 S 参数查找同时保留必要的可调参数。这些坑我基本都踩过一轮才总结出来。尤其是参考面和相位连续性那两点当年排查了整整两天才发现问题出在以为 S 参数“肯定是对的”上面。现在我做 INTERCONNECT 项目会刻意让两套模型并行存在草稿阶段靠紧凑模型快速出结果关键节点用 S 参数模型锁定细节两者互相印证才能在精度与速度之间拿到一个真正适合当前项目的平衡点。