1. 项目概述：光栅波导耦合器的经典仿真案例

这个Comsol案例虽然被作者称为"复古小案例"，但在光电子领域却是个历久弥新的基础课题。光栅波导耦合器作为集成光学中的关键元件，其耦合效率计算一直是研究生课程和工业研发中的必修内容。我十年前刚接触光子器件仿真时，第一个完整案例就是它，没想到现在依然是初学者最好的入门选择。

这类耦合器本质上是通过周期性光栅结构实现光波导与自由空间光束之间的能量转换。在硅光子学、光通信模块和各类传感器中，它的性能直接决定了整个系统的光学损耗。Comsol作为多物理场仿真标杆，其波动光学模块特别适合处理这类涉及模式耦合的复杂问题。

2. 核心物理场与边界条件设置

2.1 波动光学模块的选择依据

在Comsol启动界面选择"波动光学"模块而非"射频"或"光学射线"，是因为我们需要精确模拟近场衍射效应。光栅周期通常与波长相当（本例约1μm量级），此时几何光学近似完全失效。波动光学模块通过求解麦克斯韦方程组，能准确捕捉光栅区域的电磁场分布。

注意：Comsol 6.4版本后波动光学模块改进了周期性边界条件的处理方式，建议使用新版本来避免旧版中可能出现的收敛问题。

2.2 几何建模关键参数

创建一个二维模型即可满足基础仿真需求（三维计算量会指数级增长）。核心几何参数包括：

• 波导厚度：220nm（标准SOI晶圆参数）
• 光栅周期：Λ=630nm（由Bragg条件初步估算）
• 刻蚀深度：70nm（需权衡散射强度与工艺可行性）
• 占空比：0.5（齿宽/周期比）

% 参数化建模示例代码 lambda = 1550e-9; % 工作波长 n_eff = 2.8; % 波导有效折射率 theta = 10; % 目标耦合角度(度) Λ = lambda/(n_eff - sind(theta)); % 光栅周期计算

2.3 材料属性定义

硅波导的折射率需要设置为波长相关函数：

n_Si = 3.48 + 0.0001i @1550nm

二氧化硅衬底的折射率则简化为常数1.44。实际项目中建议使用Material Library中的实测数据，特别是当涉及温度变化时。

3. 网格划分的特殊技巧

3.1 手动网格控制要点

在光栅区域必须实施网格加密：

• 最大单元尺寸 ≤ λ/10/n_max = 1550/10/3.48 ≈ 45nm
• 使用边界层网格处理波导-光栅界面（3层，拉伸因子1.2）
• 远离光栅的区域可逐渐放宽到λ/5

常见报错解决：当出现"创建域的扫掠网格失败"时，检查是否有未被注意的微小几何交叉或重叠。

3.2 周期性边界条件实现

1. 在x方向设置Floquet周期边界
2. 顶部设置完美匹配层(PML)吸收辐射场
3. 波导端口使用模式分析边界，计算入射基模

边界条件设置要点： - 左/右侧：周期性条件，k_vector = (2*pi/Λ,0) - 顶部：PML厚度≥λ/2 - 底部：端口模式激励，选择TE0模

4. 耦合效率的计算方法

4.1 远场辐射计算

通过场计算器定义：

P_rad = integrate(emw.PoavZ, top_boundary) P_in = port_integrate(emw.Poavn, input_port) eta = P_rad/P_in

注意要排除PML区域的虚假能量吸收。

4.2 参数化扫描优化

典型优化变量：

1. 光栅周期（±50nm范围扫描）
2. 刻蚀深度（50-100nm）
3. 入射角（5-15度）

使用批处理模式运行参数组合，后处理时提取η-max对应的几何参数。

5. 常见问题排查手册

6. 高级技巧：耦合效率的影响因素分析

6.1 工艺误差的影响

通过参数扰动模拟实际制造偏差：

• 刻蚀深度±10nm → η变化约15%
• 侧壁角度偏离垂直 → 需启用"变形几何"接口
• 材料折射率波动 → 使用随机参数扫描

6.2 温度效应的耦合分析

添加固体传热和热膨胀模块：

1. 定义热源（如光吸收）
2. 计算温度分布
3. 通过"多物理场耦合"更新折射率
4. 观察η随温度漂移

7. 模型验证与实验对比

建议的验证步骤：

1. 收敛性测试：逐步加密网格直到η变化<1%
2. 解析解对比：简单情况下的Bragg条件验证
3. 文献比对：与已发表同类结构数据交叉验证

实测中常见差异来源：

• 表面粗糙度（未在模型体现）
• 材料吸收系数不准确
• 偏振纯度（模型假设理想偏振）

这个案例虽然基础，但掌握了它就能触类旁通处理各类周期性光学结构问题。我建议初学者可以尝试以下扩展练习：

1. 将二维模型扩展为三维倾斜光栅
2. 添加电光效应实现可调耦合
3. 结合Lumerical进行联合仿真