4. 耦合光学基础与光束传播分析

光纤阵列（FA）与微透镜的耦合设计，本质上是一个高斯光束被透镜变换后与目标波导（或另一段光纤）模场相互匹配的问题。要判断一个设计方案能走多远，必须先把三件事说清楚：光在自由空间里怎么传播、透镜对波前做了什么变换、以及两个模场重叠多少才算“耦合上了”。本节围绕这三条主线展开，并给出用光束传播法（BPM）或有限差分方法观察光场分布的思路，最后讨论耦合效率的理论上限从何而来。

一、高斯光束：用两个参数描述整条光路

单模光纤出射的光近似为基模高斯光束。在傍轴近似下，它沿传播方向 z 的场分布由束腰半径 w₀ 和波长 λ 唯一确定。定义瑞利长度 z_R = πw₀²/λ，则任意位置的光斑半径为

w(z) = w₀ · sqrt(1 + (z/z_R)²)
R(z) = z · (1 + (z_R/z)²)      // 波前曲率半径
θ  ≈ λ / (π w₀)                 // 远场半发散角

这三个量构成一个整体：束腰越小，发散越快；束腰越大，准直距离越长。以 1310 nm 波长、模场半径 w₀ = 4.6 μm（对应约 9.2 μm 模场直径）为例，z_R ≈ π×4.6²/1.31 ≈ 50.7 μm。也就是说，光离开端面仅约 50 μm 后光斑就扩大到 √2 倍。这解释了一个常见工程现象：FA 与微透镜之间的空气间隙若达到数百微米，光束已明显发散，必须靠透镜重新收束，而不能指望“对准就行”。

需要强调符号约定：本文统一取光传播方向为 +z，波前曲率半径 R(z) 在发散段为正、会聚段为负，q 参数写作 1/q = 1/R − jλ/(πw²)。后续矩阵运算与图示均遵循此约定。

二、模场直径匹配：耦合效率的第一道门槛

两段光路对接时，耦合效率由两模场的归一化重叠积分决定。若两侧都是高斯模、且波前曲率已通过透镜校正到近似平面，则功率耦合效率可写成

η = [ 2·w₁·w₂ / (w₁² + w₂²) ]² · exp( −2d² / (w₁² + w₂²) )

其中 d 为横向偏移。第一项是模场半径失配的惩罚，第二项是横向对准误差的惩罚。代入一组示例参数：若 w₁ = 4.6 μm、w₂ = 5.2 μm，仅半径失配一项就有 η ≈ 0.994；但若横向偏移 d = 2 μm，指数项降到约 0.83，总效率约 0.82。这说明在小模场系统里，对准精度往往比模场匹配更早成为瓶颈。工程上常把“模场直径匹配”和“容差分配”放在一起考虑：先让 w₁/w₂ 落在 0.9～1.1 区间，再把剩余预算留给横向与轴向误差。

三、ABCD 矩阵：把整条光路压缩成一次乘法

傍轴光学中，光线状态用 (r, r') 表示，任何一段光学元件都对应一个 2×2 矩阵。自由空间距离 L、薄透镜焦距 f 的矩阵分别为

M_free(L) = [1  L]      M_lens(f) = [ 1   0]
            [0  1]                  [-1/f 1]

对 q 参数则有 q₂ = (A q₁ + B)/(C q₁ + D)。FA—微透镜—目标波导的典型链路可写成 M = M_free(L₂)·M_lens(f)·M_free(L₁)，其中 L₁ 是光纤端面到透镜的距离，L₂ 是透镜到目标面的距离。设计时通常固定 L₁、L₂ 的机械约束，反解 f 使目标面处的光斑半径等于目标模场半径、且波前曲率接近零。以 w₀ = 4.6 μm、L₁ = 200 μm 为例，若希望目标面光斑仍为 4.6 μm，需要透镜把发散光束重新聚焦，f 大致落在数百微米量级——这正是微透镜而非普通透镜的原因。

四、用 BPM 与有限差分看“光场到底长什么样”

ABCD 矩阵只给出光斑半径和曲率，无法反映透镜边缘、端面菲涅尔反射、轴向像差等细节。这时需要数值求解傍轴亥姆霍兹方程。BPM 的思路是把传播方向切成 Δz 小步，每一步在频域乘一个相位因子、在空域乘透镜相位，交替推进；有限差分（FD-BPM）则用 Crank–Nicolson 格式对横向二阶导做隐式离散，稳定性更好，适合折射率对比度较大的微透镜结构。

一个可操作的示例流程如下（参数为教学假设，非实测）：计算窗口 60 μm × 400 μm，横向网格 0.1 μm，纵向步长 0.5 μm，波长 1310 nm；入射场取 w₀ = 4.6 μm 的高斯分布；在 z = 200 μm 处放置焦距 f = 250 μm 的薄透镜相位 exp(−jk r²/(2f))；记录每个 z 截面的 |E|²。观察三个判据：目标面处光斑半径是否接近 5.2 μm、波前曲率是否接近零、以及场分布是否出现明显旁瓣。若旁瓣超过主峰 10%，通常意味着透镜孔径截断过强，需要增大透镜口径或缩短 L₁。

五、耦合效率的理论极限从哪里来

把前面的关系串起来，理论极限由三个因子相乘决定：模场半径失配、波前曲率失配、以及横向/轴向对准误差。即便对准完美、曲率校正完美，只要 w₁ ≠ w₂，效率就封顶在 [2w₁w₂/(w₁²+w₂²)]²。例如 4.6 μm 对 5.2 μm 封顶约 99.4%，而 4.6 μm 对 8 μm 则降到约 84%。这就是为什么在 FA 与微透镜设计中，先做模场直径的“粗匹配”，再用透镜做“细校正”，而不是反过来指望透镜把任意模场都变完美。

实际设计时建议按以下顺序推进：先用 ABCD 矩阵扫参，得到 (f, L₁, L₂) 的可行域；再用 BPM 在可行域内验证光场细节与孔径截断；最后把横向容差、轴向容差代入效率公式，确认在 ±1 μm 偏移下效率仍满足指标。这样得到的不是单一“最优点”，而是一个有工程余量的设计窗口，这才是耦合光学分析真正的交付物。