4. 光波导与无源器件电路级建模

在硅光电路级仿真中，光波导与无源器件不能像晶体管那样直接调用工艺库模型，必须先把电磁场问题“压缩”成电路可用的等效参数。核心思路是：用有效折射率法把三维波导降为二维等效平板，再用本征模展开得到传播常数与模场分布，最后把弯曲、交叉、色散、温度等效应折算成 S 参数或等效折射率随波长、温度的变化关系，形成紧凑模型。

一、有效折射率法与模式求解

对截面为 500 nm × 220 nm 的硅条形波导，包层为 SiO₂，先按 TE 基模求解。有效折射率法把矩形波导分解为两个方向的一维限制问题：水平方向由宽度 W 决定，垂直方向由厚度 H 决定。以 1550 nm 为例，硅材料折射率约 3.476，SiO₂ 约 1.444。用一维平板波导色散方程求解，可得到水平等效折射率 neff,x ≈ 2.85，垂直方向 neff,y ≈ 2.60，再组合得到整体 neff ≈ 2.44。这个值直接决定波导的相位常数 β = 2πneff/λ，是后续所有电路级参数的基础。

实际提取时，建议用本征模展开（EME）或有限差分模式求解器做一次“标定”：扫描宽度 400～600 nm，记录 neff 与群折射率 ng。例如宽度从 500 nm 增至 520 nm，neff 约上升 0.012，ng 约上升 0.03。这些灵敏度系数是紧凑模型里“宽度偏差”端口的系数来源。

硅波导截面模式求解与电路级参数映射 Si 500×220 nm SiO₂ 包层 SiO₂ 埋氧 截面：W=500 nm, H=220 nm 模式求解 本征模展开结果 n_eff ≈ 2.44 @1550 nm n_g ≈ 4.20 β = 2πn_eff/λ 模场直径 ≈ 450 nm （示例值，随工艺浮动） 映射 电路级紧凑参数 S21 相位 = -βL 损耗 α (dB/cm) n_g 决定群时延 dn_eff/dT 温度系数 → 波长/温度端口 色散效应（示例） λ=1530 nm → n_eff≈2.452 λ=1550 nm → n_eff≈2.440 λ=1570 nm → n_eff≈2.428 斜率 ≈ -0.0006 /nm 温度效应（示例） T=25°C → n_eff≈2.440 T=65°C → n_eff≈2.452 dn_eff/dT ≈ 1.8×10⁻⁴ /K → 热光调谐基础

二、弯曲与交叉损耗的紧凑建模

直波导模型只解决传播问题，实际版图中大量存在弯曲和交叉。弯曲损耗来自辐射泄漏，半径越小、折射率对比越低，损耗越大。对 500×220 nm 硅波导，弯曲半径 R=5 μm 时 TE 模弯曲损耗约 0.02 dB/90°，R=3 μm 时升至约 0.15 dB/90°，R=2 μm 时可能超过 1 dB。紧凑模型通常写成：

α_bend(R) = a · exp(-b·R) + c
# 示例拟合：a=2.5, b=0.9, c=0.005 (dB/90°, R单位μm)
# R=5 → 2.5·exp(-4.5)+0.005 ≈ 0.033 dB/90°

交叉损耗则与交叉角度和模式重叠有关。小角度交叉（<30°）会激发高阶模和辐射模，损耗显著；大角度交叉（>60°）损耗可降到 0.01～0.05 dB。建模时把交叉当作一个二端口网络，插入损耗用角度查表或拟合式给出，串扰则用相邻端口耦合系数表示。若交叉处有锥形过渡或模式转换器，需单独提取 S 参数。

三、色散与温度效应的参数提取

色散决定不同波长下的相位和群时延，温度效应决定器件工作点漂移。提取时建议做两组扫描：

  • 波长扫描：1530～1570 nm，步长 1 nm，记录 neff(λ)。用二次多项式拟合 neff(λ)=n₀+n₁(λ-λ₀)+n₂(λ-λ₀)²，其中 n₁ 对应色散斜率，n₂ 对应色散曲率。
  • 温度扫描：25～85°C，步长 10°C，记录 neff(T)。硅的热光系数约 1.8×10⁻⁴ /K，SiO₂ 约 1×10⁻⁵ /K，因此波导整体 dneff/dT 主要由硅区限制因子决定。

把两组数据合并，得到紧凑模型的核心表达式：

n_eff(λ,T) = n₀ + n₁(λ-λ₀) + n₂(λ-λ₀)² + k_T(T-T₀)
# 示例：n₀=2.440, n₁=-6e-4 /nm, n₂=1e-6 /nm², k_T=1.8e-4 /K
# λ₀=1550 nm, T₀=25°C

该式可直接用于计算相位变化 Δφ = 2π·Δneff·L/λ，进而得到马赫-曾德尔干涉仪、微环等无源器件的温度调谐和波长漂移。例如 L=100 μm 的波导，温度升高 40 K，Δneff≈0.0072，相位变化约 2π×0.0072×100/1.55 ≈ 2.9 rad，对应微环谐振峰漂移约 0.8 nm。

四、紧凑模型构建与验证要点

一个可用的波导紧凑模型至少包含以下端口和参数：光输入/输出端口、波长端口、温度端口、宽度偏差端口。内部用 neff(λ,T,W) 计算相位和损耗，再转成 S 参数。验证时，先用直波导长度扫描确认相位斜率与 ng 一致，再用弯曲半径扫描确认损耗拟合误差小于 20%，最后用温度扫描确认谐振峰漂移方向与 dneff/dT 符号一致。

参数提取时，建议把“材料色散”和“结构色散”分开记录：材料色散来自硅和 SiO₂ 的折射率随波长变化，结构色散来自波导截面限制随波长变化。两者在紧凑模型中可合并为 n₁、n₂，但分开记录有助于定位工艺偏差来源。
弯曲损耗拟合式只在特定截面和偏振下成立。若波导宽度从 500 nm 改为 450 nm，或从 TE 切换到 TM，必须重新提取 a、b、c 系数，不能直接套用。

实际工程中，常把上述参数写成 JSON 或 CSV 参数表，由仿真器读取后生成 S 参数。例如：

{
  "waveguide": {
    "width_nm": 500,
    "height_nm": 220,
    "n0": 2.440,
    "n1_per_nm": -0.0006,
    "n2_per_nm2": 1e-6,
    "k_T_per_K": 0.00018,
    "loss_dB_cm": 2.0,
    "bend": {"a": 2.5, "b": 0.9, "c": 0.005}
  }
}

这样，电路级仿真就能在保持精度的同时，把电磁场计算量降到最低，为后续系统验证提供可复用的无源器件模型。

← 上一章 📖 返回目录 下一章 →