《光谱分析在AWG设计中的关键应用》

第3章 3 AWG的物理光学模型

3 · AWG的物理光学模型

做AWG设计这么多年,我越来越觉得——物理光学模型才是整个器件的灵魂。你光看光栅方程或者传输函数,总觉得隔了一层。今天咱们就把波动光学的底子彻底捋一遍,从输入波导到阵列波导的相位延迟,一个环节都不放过。

输入波导 自由传播区 (平板波导) 阵列波导 (ΔL) 自由传播区 (输出耦合) 输出波导 相位延迟 φ = β·ΔL (各臂等差递增) 波导/平板区 光传播方向

3.1 从波动光学说起

AWG本质上是一个多光束干涉仪。你想想看,光从输入波导耦合进平板波导,在自由传播区里像涟漪一样散开,然后被阵列波导的入口“接住”。每个阵列波导的长度都不一样,这就引入了相位差。最后在输出端的自由传播区再次干涉,不同波长的光就被分到了不同的输出波导。

我个人习惯把整个过程拆成三步:输入耦合 → 阵列相位积累 → 输出干涉。每一步都可以用波动光学的数学语言描述清楚。

核心思想: 阵列波导的等差长度 ΔL 产生等差相位延迟,等效于一个闪耀光栅。波长不同,干涉主极大的出射角就不同——这就是AWG的分波原理。

3.2 输入/输出波导的场分布

单模波导的模场一般用高斯近似。我记得有一次做硅基AWG,发现模场跟平板波导的耦合效率怎么都调不上去,后来发现是模场半径没匹配好。这里给出常用的表达式:

E_in(x, y) = A · exp[ -(x² + y²) / ω₀² ]

其中 ω₀ 是模场半径,A 是归一化常数。在平板波导入射端,这个场会衍射,形成一定的角谱分布。说白了,就是光在自由传播区里“散开”了,散开的角度由波长和波导尺寸决定。

3.3 平板波导(自由传播区)的传输模型

自由传播区可以看作一个二维的透镜系统。光从输入波导发出,经过一段距离 L_f 后,在阵列波导入口处形成一个空间分布。我常用菲涅尔衍射积分来描述:

U(x₁) = ∫ E_in(x₀) · exp[ j·k·(x₁ - x₀)² / (2L_f) ] dx₀

这里 k = 2π/λ,积分范围覆盖输入波导的模场。实际设计时,我们更关心的是阵列波导接收到的振幅和相位。嗯,这里要注意:平板波导的折射率通常跟阵列波导不同,边界处会有菲涅尔反射,不过一般通过抗反射设计压到很低了。

我在做一款25GHz通道间隔的AWG时,发现平板波导长度 L_f 稍微偏了0.5μm,中心波长就漂了0.3nm。所以自由传播区的长度必须严格按光栅方程校准,不能光靠仿真。

3.4 阵列波导的相位延迟分析

阵列波导是AWG的“心脏”。每根波导的长度依次增加 ΔL,对应的相位延迟为:

φ_m = β · (L₀ + m·ΔL)   (m = 0, 1, 2, ..., M-1)

其中 β = 2π·n_eff / λ 是传播常数,n_eff 是波导有效折射率。M 是阵列波导数,一般几十到几百根。为什么需要这么多?因为干涉的条纹锐度跟光束数量直接相关——阵列数越多,光谱通道的串扰越低。

我建议你记住一个经验值:阵列波导数 M 通常取 100~300。太少则串扰大,太多则版图面积过大且损耗增加。之前有个项目为了追求极致串扰,把M做到了400,结果插损大了1.5dB,得不偿失。

参数符号典型值对光谱的影响
波导长度差ΔL10~100 μm决定通道间隔
阵列波导数M100~300串扰 & 插损
有效折射率n_eff1.45~3.5中心波长 & 热漂移
平板波导长度L_f数百 μm~mm聚焦 & 通道均匀性

3.5 输出端的干涉与光谱形成

阵列波导的输出端再次进入自由传播区,各光束以不同的相位叠加。输出波导处的光场可以写成:

E_out(θ) = Σ A_m · exp[ j·(φ_m + k·d·sinθ) ]

这里 d 是阵列波导出口的间距,θ 是出射角。当相位差满足 2π 的整数倍时,出现干涉主极大。你想想看,波长一变,β 就变,主极大的方向也跟着变——这就是AWG的分波本质。

曾经我遇到一个案例:阵列波导的 ΔL 设计值完全正确,但流片回来光谱通道整体偏移了半个通道。排查了很久,发现是应力导致的有效折射率变化。所以物理模型里一定要留出折射率微调的余量,比如通过热光效应或紫外后处理来补偿。

3.6 模型中的关键近似与局限

上面这个模型做了几个近似:

  • 标量近似:忽略偏振耦合,实际中TE/TM模的有效折射率不同,会产生偏振相关频移。
  • 弱导近似:波导折射率差不大时,模场可用高斯近似;但对于高折射差平台(如Si₃N₄、SOI),需要用数值模场。
  • 等幅近似:假设所有阵列波导的振幅 A_m 相等,实际由于损耗和耦合不均匀,会有幅度锥削,影响串扰和通带形状。

不过话说回来,这个物理光学模型已经能准确预测AWG 90%以上的光谱行为。剩下的细节,靠仿真微调和工艺反馈就能搞定。

我个人习惯在初步设计阶段先用这个解析模型扫参数,确定 ΔL、M、L_f 的大致范围,再用FDTD或BPM做精细验证。这样效率最高,而且不容易陷入数值陷阱。