硅基波导非线性参数：n₂、有效非线性系数与损耗模型
做硅光计算芯片，非线性参数是绕不开的核心。说白了，光学激活函数能不能做出来，就看你对这几个参数摸得清不清。今天我就把非线性折射率 n₂、有效非线性系数 γ、还有损耗模型好好拆一遍。这里头有不少坑，我踩过几次，给你指条明路。
非线性折射率 n₂：材料的"脾气"
非线性折射率 n₂，是介质在强光下折射率变化的比例系数。公式很简单：Δn = n₂·I，I 是光强。硅的 n₂ 有多大？
- 晶体硅：典型值 4.7×10⁻¹⁸ m²/W（1550 nm 波长）
- 氮化硅：约 2.5×10⁻¹⁹ m²/W，比硅小一个数量级
- 非晶硅：受氢含量影响，范围 2～6×10⁻¹⁸ m²/W
嗯，要注意一点：这个数值在可见光和近红外波段差别挺大。我做过一个项目，用了不同批次的 SOI 片，测出来的 n₂ 差了一倍——后来发现是加工温度导致的晶格应力不同。你想想看，材料的"脾气"其实挺敏感的。
有效非线性系数 γ：波导的"综合战斗力"
有了 n₂，还得看波导结构。工程上更常用的是 γ（gamma），它把模场面积和材料参数打包了：
γ = 2π·n₂ / (λ·Aeff)
其中 Aeff 是有效模场面积。硅波导的 Aeff 可以做到 0.1～1 μm² 量级，γ 能达到 100～1000 /W/m。反观光纤，γ 才 1～10 /W/m。这就是为什么硅基非线性器件的长度可以做到毫米级。
有效模场面积怎么算？我一般不手算，直接用 FDTD 扫一遍。但有个经验：越小越好，但别过分压缩。光场被挤到波导边缘时，散射损耗会暴涨。我曾经为了追求高 γ，把波导做到 400×220 nm，结果损耗太大，非线性反而出不来。踩过一次就不再犯了。
| 波导类型 | n₂（×10⁻¹⁸ m²/W） | Aeff（μm²） | γ（/W/m） | 损耗（dB/cm） |
|---|---|---|---|---|
| 标准脊型（Si） | 4.7 | 0.2 | 960 | 2～3 |
| 狭缝波导（Si） | 4.7 | 0.08 | 2400 | 5～8 |
| 氮化硅脊型 | 0.25 | 0.5 | 20 | 0.5～1 |
| 悬浮膜波导 | 4.7 | 0.15 | 1280 | 10～15 |
看出来没？高 γ 往往伴随高损耗。这是个 trade-off，做芯片设计时得自己权衡。
损耗模型：光在波导里"折腾"了哪些能量
损耗会吃掉非线性效果，所以必须建模清楚。硅波导里损耗分为两大类：
- 线性损耗（传播损耗 αlin）：散射、材料吸收、弯曲损耗。通常 1～5 dB/cm。
- 非线性损耗：
- 双光子吸收（TPA）：两个光子同时被吸收，产生电子-空穴对。在 1550 nm，TPA 系数 βTPA ≈ 0.5～1.5 cm/GW。
- 自由载流子吸收（FCA）：TPA 产生的自由载流子吸收光子。尤其在高功率下，FCA 会严重限制非线性效果。
我一般用这样一个经验公式来估算总损耗：
α_total = α_lin + β_TPA · I + σ_FCA · N_carr其中 σ_FCA 是自由载流子吸收截面（约 1.45×10⁻²¹ m²），Ncarr 是载流子浓度。你怎么知道载流子浓度？通过载流子寿命模型算，或者直接测。
知识体系一览
下面这张 SVG 图，把 n₂、γ、损耗模型之间的关系串起来了。你顺着箭头看，就能知道参数怎么传导、设计时要抓哪些点。
为什么要把 γ 和损耗放在一起看？因为芯片级非线性效果最终取决于 γ / αlin 这个比值。比值越大，你才能在越短的长度上积累非线性相移。
参数提取与调试经验
我怎么测这些参数？分享一套标准流程：
- 测线性损耗：用 cut-back 法，测不同长度波导的透射率，斜率就是 αlin。
- 测 n₂：用 Z-scan 或者自相位调制（SPM）法。我个人偏爱 SPM，在波导一端打脉冲，看输出光谱展宽，反推 n₂。
- 测 γ：直接测四波混频（FWM）效率，配合已知的波导长度，用公式反算。
- 测 TPA 系数：测功率相关的非线性传输曲线，和理论模型拟合。
最后说个实用技巧：做光学激活函数时，非线性相移需要精确控制。你可以通过调节波导长度、输入功率、加的偏置电流来动态调 n₂。但记住——别让功率进入 TPA 饱和区，否则相位会大幅劣化。我一般在仿真阶段先用 self-focusing 公式算一遍，确认工作点在线性区。
好了，这部分参数你要反复实践才能手感。下次你调芯片时，如果发现非线性效果一直出不来，先查这几个数——八成是 γ 不够或者 FCA 在捣鬼。