1. 光的波动性与标量衍射理论

各位同学,咱们这门课叫《标量衍射理论指导下的DOE设计》。第一节课,得先把底子打好。说白了,就是搞清楚光到底是什么,以及我们怎么用数学去描述它传播时的那些“绕弯子”行为。

我刚开始接触衍射光学的时候,也犯过嘀咕:明明有麦克斯韦方程组那么精确的东西,为什么还要搞一套“标量近似”?后来在项目里吃过亏才明白——工程问题讲究的是“够用就好”。你想想看,如果每次设计DOE都要解一遍矢量波动方程,那项目周期得拖到猴年马月去。

1.1 光的电磁波本质

光是什么?这个问题物理学家吵了几百年。现在公认的结论是:光是一种电磁波。

麦克斯韦方程组告诉我们,变化的电场会产生磁场,变化的磁场又会产生电场,两者互相激发,就能在空间里传播出去。这就是电磁波。光波,就是频率在某个特定范围内的电磁波。

核心要点: 光波的传播,本质上就是电磁场在空间中的扰动传播。电场矢量 E 和磁场矢量 H 相互垂直,且都垂直于传播方向。

描述一个单色光波,我们通常用它的复振幅形式:

U(P, t) = A(P) * exp(j * ω * t) * exp(-j * k * r)

这里:

  • A(P) 是空间某点 P 的振幅
  • ω 是角频率
  • k = 2π/λ 是波数
  • r 是传播距离

在实际工程中,我们更关心的是光波在空间中的复振幅分布,也就是把时间因子 exp(jωt) 去掉后的部分。为什么?因为DOE只改变光的空间相位分布,时间项在干涉、衍射计算中会被消掉。

我的经验: 在Zemax或VirtualLab里做仿真时,记得确认光源设置里的“相干性”参数。我曾经因为默认设成了非相干光源,算出来的衍射效率怎么都对不上实验数据,折腾了两天才发现是这里的问题。

1.2 惠更斯-菲涅尔原理

好,现在我们知道光是波了。那波遇到障碍物会怎样?会绕过去。这就是衍射。

惠更斯在17世纪提出:波前上的每一点,都可以看作是一个新的子波源。这些子波源发出的球面波,在下一时刻的包络面,就是新的波前。

后来菲涅尔给这个想法加了点“料”:这些子波之间会发生干涉。也就是说,空间某一点的振动,是波前上所有子波源发出的球面波在该点叠加的结果。

惠更斯-菲涅尔原理的数学表达:

U(P) = (1/jλ) * ∬Σ [U(Q) * exp(jkr) / r] * cosθ * dS

其中:

  • U(P) 是观察点 P 的复振幅
  • U(Q) 是衍射屏上 Q 点的复振幅
  • r 是 Q 到 P 的距离
  • cosθ 是倾斜因子
  • Σ 是衍射屏上的开口区域

这个公式看着有点吓人,但说白了就是一句话:把衍射屏上的每个点都当成一个小喇叭,它们一起朝观察点喊话,你听到的是所有声音叠加后的结果。

1.3 基尔霍夫衍射公式

惠更斯-菲涅尔原理虽然物理图像清晰,但严格来说是个假设。基尔霍夫从波动方程出发,用数学严格推导出了衍射公式。

基尔霍夫做了两个关键假设:

  1. 在衍射屏开口处,场分布与没有屏时完全一样
  2. 在屏的遮挡部分,场分布为零

这两个假设听起来很合理,但严格来说并不完全正确。不过,在大多数工程问题中,它们足够精确。

基尔霍夫衍射公式长这样:

U(P) = (1/4π) * ∬Σ [U(Q) * ∂/∂n (exp(jkr)/r) - (exp(jkr)/r) * ∂U(Q)/∂n] * dS

嗯,比惠更斯-菲涅尔的公式更复杂了。但别怕,在实际应用中,我们很少直接用这个形式。

注意: 基尔霍夫公式在边界上的自洽性存在问题。当观察点靠近衍射屏边缘时,计算结果会与实际情况有偏差。我在设计小尺寸DOE(比如几十微米量级)时,就遇到过这个问题——仿真和实测差了10%以上。后来改用瑞利-索末菲公式才把误差降下来。

1.4 菲涅尔-基尔霍夫衍射公式的标量近似

终于到重点了。为什么我们要做标量近似?

因为光波是矢量波,电场和磁场各有三个分量,严格计算需要解六个耦合的偏微分方程。这在工程上几乎不可行。

标量近似的核心思想: 把光波当成一个标量场来处理,忽略其矢量特性。也就是说,我们只关心电场的一个分量(比如 x 方向的分量),并假设其他分量可以类似处理。

这样做需要满足两个条件:

  • 衍射孔径远大于波长(通常要求孔径 > 10λ)
  • 观察距离远大于波长

在标量近似下,菲涅尔-基尔霍夫公式简化为:

U(x, y) = (exp(jkz) / jλz) * ∬ U(ξ, η) * exp{ jk/2z * [(x-ξ)² + (y-η)²] } * dξ dη

这个公式,就是菲涅尔衍射积分。它是我们设计DOE最常用的工具之一。

为什么这个公式重要?

因为它把衍射计算变成了一个卷积或傅里叶变换问题。而傅里叶变换,我们有FFT算法,计算速度飞快。这就是现代DOE设计软件能在几秒钟内算出复杂衍射图案的原因。

如果观察距离更远,满足 z >> k(ξ²+η²)/2,菲涅尔积分还可以进一步简化为夫琅禾费衍射:

U(x, y) = (exp(jkz) / jλz) * exp{ jk/2z * (x²+y²) } * ∬ U(ξ, η) * exp{ -j2π/λz * (xξ + yη) } * dξ dη

这个公式就更漂亮了——它直接就是衍射屏复振幅的傅里叶变换。

我的建议: 刚开始学DOE设计时,先别急着用复杂的迭代算法。先用夫琅禾费近似手算几个简单例子(比如矩形孔、圆孔的衍射图案),把物理图像建立起来。我当年就是靠这个“笨办法”,才真正理解了为什么DOE能把一束光变成任意形状。

1.5 本章知识结构

下面这张图,把本章的核心逻辑串起来了。建议你保存下来,后面学完整个课程再回来看,会有更深的理解。

标量衍射理论知识结构图 标量衍射理论核心逻辑 光的电磁波本质 惠更斯-菲涅尔原理(子波干涉) 基尔霍夫衍射公式(严格数学推导) 标量近似(孔径 > 10λ,观察距离远大于波长) 菲涅尔衍射(近场) 夫琅禾费衍射(远场) 工程 常用 FFT

1.6 标量近似的适用范围

标量近似不是万能的。我整理了一个表格,帮你快速判断什么时候能用,什么时候不能用。

条件 标量近似适用 需用矢量理论
孔径尺寸 > 10λ < 10λ
观察距离 远大于波长 近场(亚波长距离)
衍射角 < 30° > 30°
偏振效应 可忽略 显著(如偏振分束DOE)
典型应用 常规DOE、光束整形 超表面、亚波长光栅

避坑指南: 我曾经设计过一个用于激光雷达的DOE,要求衍射角达到40°。我用标量近似算出来的效率是85%,结果实际做出来只有62%。后来一查,问题就出在大角度衍射时标量近似失效了。所以,大角度设计一定要用矢量方法验证。

1.7 小结

这一章我们走完了从“光是电磁波”到“标量衍射公式”的完整逻辑链。你不需要记住所有公式的推导细节,但一定要理解:

  • 光为什么能绕过障碍物(衍射的本质)
  • 惠更斯-菲涅尔原理的物理图像(子波干涉)
  • 标量近似在什么条件下成立(孔径 > 10λ,观察距离远大于波长)
  • 菲涅尔和夫琅禾费衍射的区别(近场 vs 远场)