3. 光波导与无源器件的S参数提取
在硅光系统中，波导、弯曲波导和定向耦合器是最基础的无源构件。它们的S参数模型是否准确，直接决定了后续链路仿真的可信度。与有源器件不同，无源器件的S参数提取通常不依赖高速测试仪器，而是从几何结构与材料参数出发，通过模式求解得到传播常数，再构建传输矩阵，最终换算为频域S参数。本节以硅纳米线波导、弯曲波导和定向耦合器为例，完整走通这条提取路径。
一、从几何到模式：有效折射率提取
硅纳米线波导的典型截面为220 nm厚硅芯层，埋氧层2 μm，上包层为空气或二氧化硅。在1550 nm波段，单模条件通常对应宽度约450～500 nm。提取S参数的第一步，是求解该截面在给定波长下的本征模式。对于TE基模，可用有限差分或有限元方法求解矢量波动方程，得到模式的有效折射率 neff 与模场分布。
有效折射率并非固定值，它随波长变化，这一色散关系正是S参数频率依赖性的物理来源。以宽度500 nm的硅纳米线为例，在1500～1600 nm范围内，TE基模的 neff 大致从2.45降至2.30，群折射率 ng 约为4.2。若忽略色散而取常数，相位误差会随传播距离迅速累积。因此，实际提取时应至少在若干波长点上求解模式，再拟合为多项式或样条，供后续插值使用。
二、传输矩阵：把相位与损耗写进链路
得到 neff(λ) 后，一段长度为 L 的直波导可用传输矩阵描述。若定义参考面处的前向波为 a、反向波为 b，则直波导的传输关系为：
a_out = a_in · exp(-j·β·L) · exp(-α·L/2)
b_in = b_out · exp(-j·β·L) · exp(-α·L/2)
其中 β = 2π·n_eff/λ，α 为功率衰减系数（单位 1/m）这里采用“前向波沿 +z 传播、相位因子取 exp(-jβz)”的约定，全文与图示保持一致。对于低损耗硅波导，α 通常在1～3 dB/cm量级，短器件中可先忽略，但弯曲段和耦合区往往需要单独考虑辐射损耗与散射损耗。
弯曲波导不能简单用直波导相位替代。曲率会使模式向内侧偏移，有效折射率略高于直波导，同时引入弯曲损耗。工程上常用等效直波导近似：先求解弯曲截面模式得到 neff,bend，再用弧长 L=R·θ 代入传输矩阵。对于半径5 μm、转角90°的硅弯曲波导，弯曲损耗可控制在0.01 dB以下；半径缩小到2 μm时损耗显著上升，此时必须把弯曲损耗作为独立参数写入矩阵。
三、定向耦合器：从耦合模到S参数
定向耦合器由两条平行波导及中间间隙构成，是S参数提取中最具代表性的无源器件。设输入端口1、2，输出端口3、4，理想无损耦合器的S参数可写为：
S31 = cos(κ·L_c)
S41 = -j·sin(κ·L_c)
S32 = -j·sin(κ·L_c)
S42 = cos(κ·L_c)
其中 κ 为耦合系数，L_c 为耦合区长度κ 本身随波长和间隙宽度变化。间隙越小、波长越长，κ 越大。实际提取时，可在多个波长点求解耦合区超模（偶模与奇模），由两超模有效折射率之差得到 κ：
κ(λ) = π·(n_even - n_odd)/λ这一关系把模式求解与S参数直接连接起来。以间隙200 nm、耦合区长度10 μm的定向耦合器为例，在1550 nm附近若 neven−nodd≈0.08，则 κ≈0.162 μm-1，κ·L_c≈1.62 rad，对应交叉端口功率比约 sin²(1.62)≈0.998，接近完全耦合。若波长偏移到1600 nm，超模折射率差略增，κ·L_c 继续增大，耦合比会越过第一个极大值后回落，这正是定向耦合器波长敏感性的来源。
四、完整提取任务：定向耦合器S参数随波长变化
下面给出一个可复现的提取任务流程，参数均为教学示例，实际数值需以具体工艺和仿真设置为准。
- 确定几何与材料：硅芯层220 nm，波导宽500 nm，耦合间隙200 nm，耦合区长度10 μm，输入输出直波导各5 μm，弯曲半径5 μm。
- 波长扫描设置：1500～1600 nm，步长5 nm，共21个采样点。在每个波长点分别求解直波导TE基模、弯曲波导模式、耦合区偶模与奇模。
- 提取有效折射率：记录 neff(λ)、neff,bend(λ)、neven(λ)、nodd(λ)，并检查模式收敛性。
- 计算耦合系数：由 κ(λ)=π(neven−nodd)/λ 得到每个波长点的 κ，再计算 κ·L_c。
- 构建传输矩阵：将输入直波导、弯曲段、耦合区、输出直波导依次级联，得到总传输矩阵。
- 换算S参数：由总矩阵得到 S31、S41、S32、S42，绘制随波长变化曲线。
- 结果判读：检查1550 nm处交叉端口是否接近1、直通端口是否接近0；观察波长偏移时两端口功率的互补振荡；确认无源条件下 |S31|²+|S41|²≈1。
五、结果检查与参数敏感性
提取完成后，建议做三项检查。第一，能量守恒：在忽略损耗的示例中，|S31|²+|S41|² 应在1附近，偏差超过1%说明矩阵级联或端口定义有误。第二，互易性：S31 与 S13 应相等，S41 与 S14 应相等。第三，波长趋势：随着波长增大，κ 通常增大，κ·L_c 增大，交叉端口功率先升后降，直通端口互补变化。若曲线趋势相反，应检查 neven 与 nodd 是否标反。
参数敏感性方面，间隙宽度对 κ 的影响最强。间隙从200 nm减小到180 nm，κ 可能增大约30%，导致 κ·L_c 明显偏离原设计值，交叉端口功率下降。耦合区长度则直接线性影响 κ·L_c，是调节耦合比最直接的手段。实际设计中，常先固定间隙，再扫描 L_c 找到目标耦合比对应的长度。
| 参数 | 示例值 | 对S参数的主要影响 |
|---|---|---|
| 波导宽度 | 500 nm | 改变 n_eff 与模式限制，影响相位与色散 |
| 耦合间隙 | 200 nm | 强烈影响 κ，进而改变交叉/直通比 |
| 耦合区长度 | 10 μm | 线性改变 κ·L_c，决定工作点 |
| 弯曲半径 | 5 μm | 影响弯曲损耗与弯曲段相位 |
| 波长 | 1500～1600 nm | 通过色散改变 n_eff 与 κ，产生波长响应 |
把上述流程封装为脚本后，只需替换几何参数与材料折射率，即可批量提取不同器件的S参数。得到的S参数数据可直接写入Touchstone文件或供系统级仿真调用，完成从版图几何到频域模型的闭环。