2. 硅光基础与仿真对象建模

硅光芯片的仿真不是从“画一个波导”开始，而是从“确定要仿真的是什么物理对象”开始。同一段硅波导，在光域仿真中是一条折射率对比约 2.0 的矩形介质通道，在热域仿真中是一块热导率约 148 W/(m·K) 的硅梁，在电域仿真中则可能是一条需要注入载流子的 PIN 结。若建模时把这三者混为一谈，后续的光电热联合求解必然出现边界条件错配。因此本节的核心任务是：把四类核心器件拆解为几何、材料、物理场三个可独立校核的建模单元。

一、硅波导：几何与模式约束

以 220 nm 顶层硅的 SOI 平台为例，单模条形波导的典型截面为宽 500 nm、高 220 nm，埋氧层 2 μm，上包层为空气或 SiO₂。这个尺寸不是随意选取的：在 1550 nm 波长下，500×220 nm 的硅条对 TE₀ 模的有效折射率约为 2.44，而对 TE₁ 模的截止条件恰好被压制。仿真时几何建模必须保留埋氧层和上包层，因为模式泄漏与衬底耦合损耗直接取决于包层厚度。

材料参数方面，硅在 1550 nm 的折射率取 3.476，热光系数 dn/dT 约 1.86×10⁻⁴ K⁻¹，热导率 148 W/(m·K)；SiO₂ 折射率 1.444，热导率 1.4 W/(m·K)。这些数值必须作为温度的函数或至少作为分段常数写入材料库，否则热光调谐的仿真结果会系统性偏小。波导的网格划分建议在硅芯区采用 10 nm 量级的四面体或矩形网格，包层区可放宽至 50 nm，以在模式求解精度与内存之间取得平衡。

二、调制器：电光耦合的建模重点

硅光调制器主流为载流子耗尽型马赫-曾德尔结构，其核心是嵌入波导两侧的 PIN 结。几何上需要明确：掺杂区距波导中心的横向偏移、结深、掺杂浓度梯度。例如一个典型设计为 P 区距中心 300 nm、N 区距中心 300 nm、本征区宽 600 nm，掺杂浓度从 10¹⁸ cm⁻³ 渐变至 10²⁰ cm⁻³ 接触区。

电学仿真需要求解泊松方程与漂移-扩散方程，输出载流子浓度分布；光学仿真则需将载流子浓度映射为折射率与吸收系数的扰动，即等离子体色散效应。这一步是调制器建模的关键接口：Δn 与 Δα 必须由载流子浓度插值得到，而非直接指定。热学仿真则关注结区焦耳热引起的温升，它会反向改变折射率，形成热-电-光三场耦合。仿真时应先做直流偏置下的静态求解，再扫描反向偏压从 0 V 到 −4 V，观察有效折射率变化量是否达到约 10⁻³ 量级，对应 π 相移所需臂长。

三、探测器：吸收与收集的协同

硅本身在 1550 nm 不吸收，因此硅光探测器通常采用 Ge 吸收层或表面等离激元结构。以 Ge PIN 探测器为例，几何上需定义 Ge 区的长度、宽度、厚度以及与硅波导的对接方式。典型参数为 Ge 长 10 μm、宽 500 nm、厚 500 nm，本征区掺杂 10¹⁵ cm⁻³，P 区和 N 区分别 10¹⁹ cm⁻³。

建模时，光学仿真给出光场在 Ge 区的重叠积分，电学仿真给出耗尽区宽度与电场分布，二者共同决定响应度。仿真中应检查：在 −1 V 偏压下，耗尽区是否覆盖大部分 Ge 吸收区；若未覆盖，则响应度会明显低于 0.8 A/W 的预期值。热学方面，暗电流引起的温升通常较小，但在高功率入射下仍需评估 Ge 与 Si 界面处的热应力。

四、光栅耦合器：周期结构与相位匹配

光栅耦合器用于将光纤模场耦合进硅波导，其几何核心是周期、刻蚀深度、占空比和入射角。以 1550 nm 为例，若采用 220 nm 硅层、70 nm 刻蚀深度、周期 630 nm、占空比 0.5，则垂直耦合效率在仿真中通常可达到约 40%～50%。仿真时需使用严格耦合波分析或时域有限差分，扫描周期与刻蚀深度，观察耦合效率峰值对应的相位匹配条件。

材料参数上，光栅区的硅与 SiO₂ 折射率与波导一致，但刻蚀引入的侧壁粗糙度会带来散射损耗。仿真中可通过在侧壁添加随机扰动来评估损耗量级。热学方面，光栅耦合器对温度不敏感，但若与调制器集成，仍需考虑整体热串扰。

五、建模顺序与校核要点

四类器件的建模应遵循统一顺序：先建立几何模型并确认关键尺寸，再赋予材料参数并检查温度依赖性，然后分别进行单物理场求解以验证边界条件，最后通过耦合接口进行联合仿真。以调制器为例，若电学求解得到的载流子浓度在波导中心处低于 10¹⁶ cm⁻³，则说明耗尽区未充分覆盖光场，此时应调整掺杂偏移或偏压，而非直接修改光学结果。

校核时建议关注三个量：光域的有效折射率与损耗、电域的耗尽区宽度与电容、热域的温升与热时间常数。这三个量分别对应器件的相位效率、调制带宽和热调谐速度。若其中任一量偏离预期，应回到几何或材料参数中查找原因，而不是在耦合求解中强行修正。只有单物理场结果可信，多物理场联合仿真才有意义。