4. MZI工作原理:光程差与相位差、干涉输出公式推导

好,咱们进入正题。马赫曾德尔干涉仪(MZI)这东西,说白了就是一个光的“分分合合”游戏。光进去,先被分成两路,走不同的路径,然后再合到一起。合的时候,两路光如果“步调一致”,就亮;如果“步调相反”,就暗。就这么简单。

但工程上,我们得用数学把它说清楚。我当年刚接触MZI时,总觉得公式推导很绕。后来发现,抓住光程差这个核心,一切都顺了。

4.1 从光程差到相位差

光在介质里走,真正影响它“相位”的,不是几何距离,而是光程(Optical Path Length, OPL)。

核心公式:

光程 = 几何路径长度 × 介质折射率

光程差 ΔL = n₁L₁ - n₂L₂

相位差 Δφ 和光程差的关系是:

Δφ = (2π / λ₀) × ΔL

其中 λ₀ 是真空中的波长。你看,只要光程差是波长的整数倍,相位差就是 2π 的整数倍,两路光完全同相,干涉增强。如果光程差是半波长的奇数倍,相位差就是 π 的奇数倍,完全反相,干涉相消。

我在项目里调试一个光纤MZI传感器时,就遇到过这个问题。当时输出总是不稳定,查了半天,发现是光纤的折射率随温度变了,导致光程差漂移。嗯,从那以后,我设计任何MZI系统,第一件事就是评估环境对光程的影响。

4.2 两束光的干涉输出公式

现在我们来推导输出光强。假设输入光振幅为 E₀,经过第一个分束器(50:50),两路振幅都变成 E₀/√2。

两路光分别经历相位 φ₁ 和 φ₂ 后,在第二个分束器处合束。合束后的输出振幅为:

E_out = (E₀/2) × [exp(iφ₁) + exp(iφ₂)]

输出光强 I_out 正比于振幅的模平方:

I_out = |E_out|² = (E₀²/4) × |exp(iφ₁) + exp(iφ₂)|²
      = (I₀/2) × [1 + cos(Δφ)]

其中 I₀ = E₀² 是输入光强,Δφ = φ₁ - φ₂ 是两臂的相位差。

我个人习惯:把这个公式记成“1 + cos(Δφ)”形式。它太常用了,几乎所有MZI的分析都从这里出发。你想想看,当 Δφ = 0 时,I_out = I₀,全部输出;当 Δφ = π 时,I_out = 0,完全相消。

4.3 干涉输出的关键特性

从公式里,我们能读出几个重要信息:

  • 余弦响应:输出光强随相位差呈余弦变化,不是线性的。
  • 周期为 2π:相位差每变化 2π,输出就循环一次。
  • 消光比:最大光强与最小光强的比值,理想情况下是无穷大,实际受分束比和损耗影响。
相位差 Δφ 输出光强 I_out 物理状态
0 I₀ 完全相长干涉
π/2 I₀/2 中间状态
π 0 完全相消干涉
3π/2 I₀/2 中间状态
2π I₀ 回到相长干涉

我曾经踩过的坑:在设计一个高精度MZI传感器时,我以为只要保证两臂几何长度相等就行。结果发现,即使长度相等,如果两臂的波导材料折射率有微小差异,光程差照样不为零。后来我学乖了,设计时直接用光程差来算,而不是几何长度。

4.4 MZI干涉过程可视化

下面这张图,把整个干涉过程画出来了。你可以看到光怎么分、怎么走、怎么合,以及相位差如何决定输出。

MZI干涉过程与输出关系 输入光 E₀ BS1 上臂 L₁, n₁ 相位 φ₁ 下臂 L₂, n₂ 相位 φ₂ BS2 输出 I_out Δφ = φ₁ - φ₂ 输出光强公式 I_out = (I₀/2)[1 + cos(Δφ)] 关键说明 1. 光经过BS1后,分成等振幅的两路:上臂和下臂。 2. 两臂的光程差 ΔL = n₁L₁ - n₂L₂ 决定相位差 Δφ。 3. 在BS2处合束,输出光强随 Δφ 呈余弦变化。 4. 当 Δφ = 0, 2π, 4π... 时输出最大;Δφ = π, 3π, 5π... 时输出为零。 5. 实际系统中,分束比、损耗、偏振都会影响消光比。

4.5 实际应用中的注意事项

公式推导完了,但工程上还有几个细节要注意:

  1. 分束器不是理想的50:50——实际分束比会有偏差,导致消光比下降。
  2. 偏振态要匹配——两路光的偏振方向如果不一致,干涉效率会降低。
  3. 温度稳定性——折射率随温度变化,光程差会漂移,输出不稳定。

我的建议:如果你在做MZI实验,先别急着调精度。先把两臂的光程差调到接近零,确认能看到清晰的干涉条纹。然后再一步步引入你要测量的变化。这样排查问题会快很多。