4. 微环谐振器S参数建模

微环谐振器是硅光子系统中最常用的波长选择元件，它的S参数模型与直波导、MMI等器件有本质区别：微环的响应不是由几何长度直接决定，而是由环内一周的相位积累与耦合区的能量交换共同形成的谐振条件决定。因此，微环S参数建模的核心任务，是把耦合系数、环内损耗、相位这三个物理量映射为随波长快速变化的传输函数，再写成端口之间的S参数矩阵，供系统级仿真调用。

一、两种基本拓扑与端口定义

先约定参照系与端口编号，全文及图示保持一致。取输入光沿直波导从左向右传播，端口1为输入端口，端口2为直通输出端口。对于上下载型微环，另有一条与之平行的下载波导，端口3为下载端口（Drop），端口4为添加端口（Add）。所有场量采用归一化复振幅，耦合区用自耦合系数 t 与交叉耦合系数 κ 描述，理想无损耦合满足 |t|²+|κ|²=1。

全通型（All-Pass）只有一个耦合区，直通波导与环耦合，没有下载波导，因此只有端口1、2，表现为一个全通滤波器：直通端在谐振波长处出现陷波，其余波长几乎无损耗通过。上下载型（Add-Drop）有两个耦合区，光可以从端口1耦合进环、再耦合到下载波导从端口3输出，因此具有带通式的下载响应，同时直通端出现陷波。两者在S参数矩阵上的差别，正是端口数与矩阵维数的差别。

二、全通型传输函数推导

设环内一周的复传输因子为 a·e^jθ，其中 a 为环内单圈振幅透过率（含传播损耗与弯曲损耗），θ 为单圈相位积累。对全通型，输入场 E₁ 在耦合区分为直通分量与进入环的分量，环内场绕行一周后再次到达耦合区。设耦合区输入为 E₁ 与环内回返场 E_r，则直通输出与环内新注入场满足：

E2 = t·E1 + κ·Er
Ein_ring = -κ*·E1 + t*·Er
Er = a·e^{jθ}·Ein_ring

消去 E_r 与 E_{in_ring}，得到全通型直通端传输函数：

T_allpass(λ) = E2/E1 = (t - a·e^{jθ}) / (1 - t·a·e^{jθ})

谐振条件为 θ = 2mπ，此时 e^jθ=1，传输函数退化为实数 (t−a)/(1−ta)。当环内损耗很小即 a→1 时，该值趋近于零，形成深陷波；当 a 明显小于 1 时，陷波深度变浅，说明损耗直接限制了消光比。谐振波长由 θ=2πn_effL/λ=2mπ 决定，即 λ_m=n_effL/m，自由光谱范围 FSR≈λ²/(n_gL)。

三、上下载型传输函数与S参数矩阵

上下载型有两个耦合区，设上耦合区为 t₁、κ₁，下耦合区为 t₂、κ₂，环内单圈因子仍为 a·e^jθ。分别对两个耦合区列场方程并联立，可得直通端与下载端传输函数：

T_through = (t1 - t2*·a·e^{jθ}) / (1 - t1·t2*·a·e^{jθ})
T_drop    = (-κ1*·κ2·√a·e^{jθ/2}) / (1 - t1·t2*·a·e^{jθ})

谐振时 e^jθ=1，下载端幅度正比于 κ₁κ₂√a，直通端为 (t₁−t₂a)/(1−t₁t₂a)。当两个耦合区对称即 t₁=t₂=t 且 a→1 时，直通端趋近于零、下载端趋近于 1，实现完整的波长上下载。把上述四个端口的场关系写成矩阵形式，即得到上下载型微环的 4×4 S参数矩阵：

[E2]   [ S21  S23 ] [E1]
[E3] = [ S41  S43 ] [E4]

S21 = T_through,  S41 = T_drop
S23 = T_drop,     S43 = T_through

其中 S21 为输入到直通、S41 为输入到下载、S23 为添加端到直通、S43 为添加端到下载。由于结构互易，S23=S41、S43=S21，矩阵具有对称性。全通型可视为上下载型在 κ₂=0 时的退化情形，其S参数矩阵为 2×2，只有 S21 与 S12 两个非零项。

四、耦合系数、损耗与相位对谐振谱的影响

三个参数对谱形的影响各有侧重，建模时必须分别标定。耦合系数 κ 决定环与波导之间的能量交换强度：κ 增大使谐振线宽变宽、品质因数 Q 下降，同时下载端峰值升高；κ 过小则线宽极窄、对波长漂移过于敏感。环内损耗 a 决定谐振峰的深度与下载效率：a 越小，全通型陷波越浅、上下载型下载峰值越低，且谐振峰底部不再降到零。相位 θ 决定谐振峰的位置：θ 中的 n_eff 随波长与温度变化，直接导致谐振波长漂移，这也是热调谐与波长锁定的物理基础。

工程上常用临界耦合判断谱形：当 t=a 时全通型陷波达到最深，称为临界耦合；t>a 为过耦合，线宽较宽；t<a 为欠耦合，陷波变浅。上下载型则关注下载端是否达到最大峰值，通常要求两耦合区与损耗匹配。

五、用仿真数据拟合提取微环S参数

实际建模时，微环的 t、κ、a、n_eff 往往不能直接精确获知，需要从仿真或测试得到的谱线反推。以下给出一个可执行的拟合流程，示例参数为合理假设，用于说明方法而非实测结果。

1. 准备数据：用 FDTD 或传输矩阵法扫描波长，得到直通端与下载端的复传输谱，波长范围覆盖至少两个 FSR，采样点建议每 FSR 不少于 500 点，避免谐振峰欠采样。
2. 初值估计：从相邻谐振峰间距读出 FSR，由 λ_m=n_effL/m 估算 n_eff；由谐振峰 3 dB 带宽 Δλ 估算 Q=λ/Δλ，再由 Q 与 κ、a 的经验关系给出 t、a 初值。
3. 构造目标函数：把实测或仿真谱与传输函数计算谱的幅度差平方和作为残差，同时对直通端与下载端加权，避免只拟合单端口导致参数不唯一。
4. 非线性最小二乘迭代：以 t₁、t₂、a、n_eff、L 为待拟合量，用 Levenberg-Marquardt 迭代，约束 0<t<1、0<a<1。
5. 结果判读：检查拟合谱与数据在谐振峰顶、峰底、半高宽三处的偏差；若峰底偏差大，优先怀疑 a 或耦合区损耗；若峰位整体偏移，调整 n_eff 或 L；若 FSR 不符，检查环长 L 的取值。

下面给出一个可直接参考的拟合参数记录表，字段按微环S参数提取的实际需要设置。

完成上述提取后，把 t₁、t₂、a、n_eff、L 写回传输函数，即可生成该微环在目标波段内的 S 参数查找表或解析模型。验证时，用另一组未参与拟合的波长点计算直通与下载响应，与仿真数据对比，若峰位误差小于 0.05 nm、峰值误差小于 0.5 dB，可认为模型精度满足系统级仿真需求。