2. 硅基光波导的热光效应与温漂物理机制
硅光器件的波长漂移，根源不在“光”本身，而在“热”。硅材料的热光系数（thermo-optic coefficient，TOC）约为 1.8×10⁻⁴ /K，比二氧化硅（约 1×10⁻⁵ /K）高出一个数量级。这意味着温度每变化 1 K，硅波导的有效折射率 n_eff 就会变化约 1.8×10⁻⁴。对于工作在 1550 nm 附近的微环谐振器，这个变化足以让谐振峰移动约 0.1 nm 量级——而微环的 3 dB 带宽往往只有 0.05～0.2 nm，温漂因此成为波长锁定必须首先量化的物理量。
一、热光效应的物理来源
热光效应描述折射率随温度的变化，定义为：
dn/dT = (∂n/∂T)_ρ + (∂n/∂ρ)_T · (∂ρ/∂T)第一项是温度直接改变电子极化率带来的贡献，第二项是热膨胀导致密度 ρ 下降、进而改变折射率的间接贡献。在硅中，第一项占主导，且符号为正：温度升高，折射率增大。这与许多聚合物（负热光系数）相反，也是硅光温漂方向一致、便于统一补偿的原因。
需要区分两个层次：材料热光系数 dn/dT 描述体材料；波导有效折射率温度系数 dn_eff/dT 描述具体波导结构。后者不仅取决于材料，还取决于模式在芯层与包层中的功率分布。对于硅芯/二氧化硅包层波导，模式限制在硅中的比例越高，dn_eff/dT 越接近硅的体值；若包层改用氮化硅或聚合物，则会被“稀释”。
定量估算的起点：对 1550 nm 波长，硅的 dn/dT ≈ 1.86×10⁻⁴ /K，二氧化硅 dn/dT ≈ 1.0×10⁻⁵ /K。若波导中约 80% 模式能量位于硅芯，则 dn_eff/dT ≈ 0.8×1.86×10⁻⁴ + 0.2×1.0×10⁻⁵ ≈ 1.51×10⁻⁴ /K。
二、从折射率变化到波长偏移
谐振波长满足 λ_res = n_eff · L / m，其中 L 为环周长，m 为整数阶数。对温度求导：
dλ_res/dT = (λ_res / n_eff) · (dn_eff/dT) + (λ_res / L) · (dL/dT)第二项来自热膨胀，硅的热膨胀系数约 2.6×10⁻⁶ /K，比热光项小约两个数量级，工程估算中常先忽略，再作为修正项。于是简化为：
dλ_res/dT ≈ (λ_res / n_g) · (dn_eff/dT)这里用群折射率 n_g 替代 n_eff 更准确，因为谐振峰移动由群速度决定。对典型硅波导，n_g ≈ 4.2，dn_eff/dT ≈ 1.5×10⁻⁴ /K，λ=1550 nm 时：
dλ_res/dT ≈ (1550 / 4.2) × 1.5×10⁻⁴ ≈ 0.055 nm/K也就是说，温度升高 1 K，微环谐振峰红移约 0.055 nm。若芯片温度波动 ±2 K，波长漂移可达 ±0.11 nm，远超许多 DWDM 通道的容差。
三、微环与马赫-曾德尔结构的温漂差异
微环和 MZI 都受热光效应影响，但表现形式不同。微环是谐振型器件，温漂表现为谐振峰整体平移，漂移量由 λ/n_g · dn_eff/dT 决定；MZI 是干涉型器件，温漂表现为输出相位变化：
Δφ = (2π/λ) · Δn_eff · ΔL其中 ΔL 为两臂长度差。MZI 的温漂灵敏度与 ΔL 成正比，而微环与周长 L 相关。关键差异在于：微环的 FSR 通常只有几纳米，温漂 0.055 nm/K 在几 K 内就会跨过一个 FSR，造成“跳模”式失锁；MZI 的相位是连续累积的，温漂表现为输出功率余弦式起伏，周期由 FSR 决定，但不会出现阶跃跳变。因此，微环对温度更敏感、锁定更难，MZI 则更适合做宽范围粗调。
四、温漂-波长偏移的定量估算方法
把上面的物理关系整理成可复用的估算流程，只需四步：
- 确定材料与模式限制因子。 查硅的
dn/dT，用模式仿真得到硅芯中的功率占比 Γ。若没有仿真，可按 Γ≈0.75～0.85 估算。 - 计算
dn_eff/dT。dn_eff/dT ≈ Γ·(dn_Si/dT) + (1−Γ)·(dn_SiO2/dT)。 - 获取群折射率
n_g。 由波导色散曲线或 FSR 反推：n_g = λ²/(FSR·L)。 - 代入
dλ/dT = λ/n_g · dn_eff/dT。 得到每开尔文的波长漂移，再乘以允许温升得到总漂移。
以一个教学示例说明：某微环半径 R=10 μm，周长 L≈62.8 μm，实测 FSR≈3.2 nm，λ=1550 nm。由 n_g = λ²/(FSR·L) = 1550²/(3.2×62800) ≈ 4.28。取 Γ=0.8，则 dn_eff/dT ≈ 0.8×1.86×10⁻⁴ + 0.2×1.0×10⁻⁵ ≈ 1.51×10⁻⁴ /K。于是 dλ/dT ≈ 1550/4.28 × 1.51×10⁻⁴ ≈ 0.0547 nm/K。若芯片温控精度为 ±0.1 K，则波长不确定度约 ±0.0055 nm，对 50 GHz（0.4 nm）通道间隔而言仅占 1.4%，可以接受；若温控只有 ±1 K，则漂移 0.055 nm，已接近通道间隔的 14%，必须补偿。
估算时注意符号约定：本文以温度升高为正方向，硅的 dn/dT 为正，因此谐振峰向长波方向（红移）移动。若使用负热光包层或聚合物，符号相反，需在公式中保留符号。
五、热膨胀项的修正
前面忽略了热膨胀。严格来说，dλ/dT 还应加上 λ·α_sub/ n_g 量级的贡献，其中 α_sub 为衬底热膨胀系数。对硅衬底 α≈2.6×10⁻⁶ /K，修正量约 1550×2.6×10⁻⁶/4.28 ≈ 9.4×10⁻⁴ nm/K，仅为热光项的 1.7%。在 ±5 K 范围内，修正量约 0.005 nm，通常可忽略；但在高精度锁定（<0.01 nm）或大温变场景下，应纳入模型。
更完整的做法是把 dn_eff/dT 和 dL/dT 一起写入温度-波长传递函数，用实测数据拟合出综合系数 K_T = dλ/dT，再用于前馈补偿。这样既保留了物理可解释性，也吸收了工艺偏差。