光与物质的周期性耦合：布拉格光栅（FBG）的模态分析与工程实践

1. 引言：波导中的微扰引入

对于光学工程师而言，光纤布拉格光栅（FBG）不仅仅是一个滤光片，它本质上是圆柱波导中的**一维光子带隙（1D Photonic Bandgap）**结构。通过在光纤纤芯引入周期性的介电常数微扰（$\Delta \varepsilon$），我们在波矢空间（k-space）人为制造了一个“禁带”。

本文将跳过基础概念，直接从波动方程出发，利用耦合模理论解析 FBG 的光谱响应及其在色散管理和极限传感中的物理机制。

2. 理论框架：耦合模方程（Coupled Mode Equations）

在弱导波近似下，光栅的作用是将前向传输的基模（Forward Fundamental Mode, amplitude $A(z)$）能量耦合到后向传输模（Backward Mode, amplitude $B(z)$）。

2.1 介电常数的周期性微扰

光纤纤芯的折射率调制可表示为：
$$ n(z) = n_{eff} + \overline{\Delta n}(z) + \Delta n(z) \cos\left(\frac{2\pi z}{\Lambda} + \phi(z)\right) $$
其中 $\overline{\Delta n}$ 是直流平均折射率变化，$\Delta n(z)$ 是交流调制幅度（决定耦合强度），$\phi(z)$ 是相位啁啾项。

2.2 耦合模方程推导

根据麦克斯韦方程组，引入微扰极化项，通过正交归一化条件，可导出描述 $A(z)$ 和 $B(z)$ 演化的微分方程组：

$$ \frac{dA}{dz} = -i\kappa B e^{i2\delta z} - \alpha A $$
$$ \frac{dB}{dz} = i\kappa^* A e^{-i2\delta z} + \alpha B $$

其中关键参数定义如下：

• $\kappa$（交流耦合系数）： $\kappa = \frac{\pi \Delta n}{\lambda} \eta$，表征单位长度的反射能力，$\eta$ 为模场重叠积分因子。
• $\delta$（失谐量 Detuning）： $\delta = \beta - \frac{\pi}{\Lambda} = 2\pi n_{eff} \left( \frac{1}{\lambda} - \frac{1}{\lambda_B} \right)$。它描述了当前波长与布拉格中心波长 $\lambda_B$ 的偏离程度。

2.3 反射率的解析解

对于均匀光栅（$\kappa$ 为常数，无啁啾），解上述方程组可得反射率 $R$：

$$ R(\lambda, L) = \left| \frac{B(0)}{A(0)} \right|^2 = \frac{\sinh^2(\sqrt{\kappa^2 - \delta^2}L)}{\cosh^2(\sqrt{\kappa^2 - \delta^2}L) - \frac{\delta^2}{\kappa^2}} $$

在中心波长处（$\delta = 0$），公式简化为工程中最常用的近似：
$$ R_{max} = \tanh^2(\kappa L) $$

工程推论： 要获得高反射率（$R>99.9%$），必须提高 $\kappa L$ 乘积。对于短光栅（如激光器腔镜），需要极高的 $\Delta n$（强光敏光纤或载氢处理）；对于长光栅，则需保证刻写的均匀性。

3. 旁瓣抑制与切趾技术（Apodization）

从傅里叶变换的角度看，有限长度 $L$ 的均匀光栅，其空间域是矩形窗函数（Rectangular Window），因此其频域响应必然伴随着Sinc函数的振荡，即旁瓣（Side-lobes）。

在密集波分复用（DWDM）系统中，旁瓣会导致信道串扰。为了消除旁瓣，必须对耦合系数 $\kappa(z)$ 进行切趾（Apodization），即让折射率调制深度在光栅两端平滑衰减至零（如高斯分布或余弦平方分布）。

$$ \Delta n(z) = \Delta n_{max} \exp\left(-\frac{4 \ln 2 \cdot z^2}{FWHM^2}\right) $$

切趾的代价是有效光栅长度缩短，导致反射谱的主峰展宽（FWHM 增加）。这是工程设计中必须权衡的带宽-串扰博弈。

4. 色散管理：啁啾光栅（Chirped FBG）

当光栅周期 $\Lambda(z)$ 沿轴向线性变化时，不同波长的光满足布拉格条件的物理位置 $z_{Bragg}$ 不同。

$$ \lambda_B(z) = 2 n_{eff} \Lambda(z) $$

短波长在光栅近端反射，长波长在光栅远端反射。这引入了与频率相关的群时延（Group Delay, $\tau_g$）：

$$ \tau_g(\lambda) = \frac{2 n_{eff} z(\lambda)}{c} $$

其斜率 $D = \frac{d\tau_g}{d\lambda}$ 即为色散量（ps/nm）。

• 应用场景： 在高速光通信（>10Gb/s）中，通过设计特定的 $\Lambda(z)$ 分布，使 CFBG 的色散正好与传输光纤（SMF-28）的色散符号相反、大小相等，实现全光链路的色散补偿。

5. 高阶传感模型：张量与交叉敏感

在传感应用中，简单的线性公式往往不足以描述复杂工况。我们需要引入光弹效应张量（Photo-elastic Tensor）。

5.1 应变光效应的严格表述

轴向应变 $\epsilon_z$ 不仅改变光栅周期，还通过泊松效应改变光纤截面尺寸和折射率分布。波长漂移的精确表达为：

$$ \frac{\Delta \lambda_B}{\lambda_B} = \epsilon_z - \frac{n_{eff}^2}{2} [p_{12} - \nu(p_{11} + p_{12})] \epsilon_z $$

其中 $p_{ij}$ 为石英玻璃的弹光系数分量，$\nu$ 为泊松比。通常，有效弹光系数 $P_e \approx 0.22$，综合灵敏度系数 $1 - P_e \approx 0.78$。

5.2 温度与应变的解耦

由于温度和应变都会导致 $\Delta \lambda_B$，在实际工程（如航空复合材料监测）中必须解耦。常用的方法是构建传输矩阵：

$$ \begin{bmatrix} \Delta \lambda_1 \ \Delta \lambda_2 \end{bmatrix} = \begin{bmatrix} K_{\epsilon 1} & K_{T 1} \ K_{\epsilon 2} & K_{T 2} \end{bmatrix} \begin{bmatrix} \epsilon \ \Delta T \end{bmatrix} $$

通过引入第二个参考光栅（仅感温不感应变）或利用双折射效应（测量 PANDA 光纤的两个偏振模的差分漂移），求解上述矩阵的逆矩阵即可实现解耦。

6. 光纤激光器中的谐振腔设计

在高功率光纤激光器（Ytterbium-doped Fiber Laser）中，FBG 充当了全反镜（HR）和输出耦合镜（OC）。

• 线宽控制： HR 光栅通常设计为高反射率（>99%）、宽带宽（1-2nm），以覆盖增益谱；OC 光栅设计为低反射率（5-10%）、窄带宽，以选模并压窄输出线宽。
• 受激拉曼散射（SRS）抑制： 在 kW 级激光器中，为了抑制 SRS，OC 光栅的设计波长和带宽必须严格匹配，确保能量集中在基频光，并通过大模场面积（LMA）光纤光栅降低功率密度。

7. 结语

从耦合模方程的解析解到非均匀光栅的数值模拟（Transfer Matrix Method），布拉格光栅的设计已经进入了高度数学化和定制化的阶段。

未来的挑战在于极端环境下的材料改性（如飞秒激光刻写的 Type-II 光栅耐受 1000°C 高温）以及在集成光路（PIC）波导上实现复杂的片上滤波器。对于光学工程师而言，理解 $\kappa$ 与 $\Lambda(z)$ 的相空间演化，是掌握这把“光之尺”的关键。

---