4. 电学仿真方法与工具链

硅光调制器的光学响应，最终由硅波导中的自由载流子浓度变化驱动。无论是载流子注入型（PIN）还是载流子耗尽型（PN 结反偏）调制器，其核心电学问题都可以归结为：在给定电极电压下，波导有源区内的电子与空穴浓度如何重新分布，以及这种分布如何随频率变化。要回答这个问题，必须回到半导体器件的漂移扩散（Drift-Diffusion, DD）模型，并借助 TCAD 工具链完成从结构到电学特性的求解。

漂移扩散模型：电学仿真的物理内核

漂移扩散模型由三组方程耦合而成。第一组是泊松方程，描述电势 ψ 与空间电荷的关系：

∇·(ε∇ψ) = -q(p - n + N_D⁺ - N_A⁻)

其中 ε 为介电常数，q 为元电荷，n、p 分别为电子和空穴浓度，N_D⁺、N_A⁻ 为电离杂质浓度。第二、三组是电子与空穴的连续性方程，配合电流密度表达式：

J_n = qμ_n n E + qD_n ∇n
J_p = qμ_p p E - qD_p ∇p
∂n/∂t = (1/q)∇·J_n - R
∂p/∂t = -(1/q)∇·J_p - R

这里 μ 为迁移率，D 为扩散系数（由爱因斯坦关系 D = μkT/q 联系），R 为复合率。在硅光调制器的典型工作条件下，掺杂浓度在 10¹⁷～10²⁰ cm⁻³ 量级，电场强度可达 10⁵ V/cm，因此迁移率必须考虑掺杂浓度与电场强度的依赖关系，常用 Caughey-Thomas 模型或 Lombardi 表面迁移率模型修正。复合项则需同时考虑 SRH 复合、俄歇复合和辐射复合，其中俄歇复合在重掺杂区不可忽略。

载流子输运的数值求解策略

DD 方程是强非线性耦合偏微分方程组，无法解析求解，必须离散化后迭代。主流 TCAD 工具（Sentaurus Device、Atlas、Minimos-NT 等）采用有限体积法或有限元法，在网格节点上离散化。求解流程通常为：先做平衡态求解得到初始电势与载流子分布，再逐步施加偏压，每一步用 Newton-Raphson 迭代求解耦合方程组。对于硅光调制器，网格划分需要特别关注 PN 结耗尽区、波导侧壁和氧化层界面，这些位置浓度梯度大，网格需加密到 1～5 nm 量级。

收敛性是实际仿真中最常遇到的问题。常见原因包括：掺杂分布突变导致初始猜测差、偏压步长过大、迁移率模型在强场下不连续。工程上的做法是先用小步长（如 0.05 V）扫描到目标偏压，再在目标点做精细求解；必要时启用阻尼 Newton 或 Gummel 迭代作为预条件。

PN 结与 MOS 电容：两类核心电学结构

硅光调制器中最常见的电学结构有两类。第一类是横向或纵向 PN 结，工作在反偏或零偏附近，通过耗尽区宽度随电压变化来改变波导中的自由载流子。以横向 PN 结为例，当反偏电压从 0 V 增加到 -3 V 时，耗尽区展宽，波导核心区的载流子被抽走，有效折射率实部增大、虚部减小。仿真中需要提取的关键量是：不同偏压下的耗尽区边界位置、波导截面内的载流子浓度积分（即 ΔN_e 和 ΔN_h），以及结电容 C_j(V)。

第二类是 MOS 电容结构，常见于积累型调制器。它由硅波导、超薄氧化层（如 5～10 nm HfO₂ 或 SiO₂）和多晶硅/金属栅构成。当栅压变化时，硅表面积累或耗尽，载流子集中在波导表面几个纳米范围内。MOS 结构的仿真难点在于量子效应和界面态：当氧化层厚度低于 5 nm 时，反型层载流子分布偏离经典 DD 模型，需要引入密度梯度或 Schrödinger-Poisson 修正。此外，界面陷阱电荷会改变平带电压，仿真中通常用固定电荷面密度 Q_f 来等效。

从电学结果到光学调制：数据传递链

TCAD 电学仿真输出的不是最终的光学指标，而是波导截面上的载流子浓度分布。要得到调制效率，需要把 ΔN_e(x,y) 和 ΔN_h(x,y) 映射到光学模式上，计算有效折射率变化：

Δn_eff = Γ_e · Δn_e + Γ_h · Δn_h
Δα = Γ_e · Δα_e + Γ_h · Δα_h

其中 Γ 为模式与有源区的重叠因子，Δn 和 Δα 由自由载流子色散效应（Soref-Bennett 关系）给出。这一步通常由光学模式求解器（如 Lumerical MODE、COMSOL Wave Optics）完成，电学与光学工具之间通过网格插值传递数据。工程上要特别注意网格匹配：电学网格在结区加密，光学网格在模式中心加密，两者插值误差会直接影响 Δn_eff 的精度。

工具链与实操要点

完整的电学仿真工具链通常包括：结构建模（Sentaurus Structure Editor / Atlas 命令语言）、掺杂定义（高斯或 erfc 分布）、物理模型选择（迁移率、复合、隧穿）、求解设置（偏压扫描、频率扫描）和结果提取（IV、CV、载流子分布）。对于硅光调制器，建议在仿真中至少完成以下检查：

• 平衡态下耗尽区宽度与理论值（√(2εV_bi/qN)）偏差是否在 10% 以内；
• 反偏漏电流是否由 SRH 或 BTB 隧穿主导，量级是否合理；
• C-V 曲线在反偏区是否呈现 1/C² 线性关系，以验证掺杂分布假设；
• 小信号 AC 仿真中，3dB 频率是否与 R_s·C_j 估算一致；
• 载流子浓度分布是否出现非物理负值或振荡，若有则需调整网格或阻尼因子。

实际工程中，电学仿真往往不是一次完成的。典型迭代路径是：先用简化结构（无量子修正、无界面态）快速扫描，确认掺杂和几何参数的大致范围；再逐步加入迁移率退化、俄歇复合、量子修正等模型，观察关键指标的变化趋势；最后用实验数据（如 C-V、IV、眼图）做标定。每一步都要记录模型开关和参数值，否则后续无法复现。对于硅光调制器，最关键的标定量通常是结电容和串联电阻，因为它们直接决定调制带宽和驱动电压。