4、衍射理论基础：惠更斯-菲涅尔原理、基尔霍夫衍射公式、菲涅尔近似与夫琅禾费近似
衍射是波动光学中标志性的物理现象，它揭示了光在传播过程中偏离几何光学直线传播的本质。从工程应用的角度看，无论是光学镜头的分辨率极限、光栅光谱仪的色散行为，还是激光束的整形与传输，其底层物理模型都建立在衍射理论之上。本章将沿着“物理思想→数学严格化→工程近似”的脉络，系统梳理惠更斯-菲涅尔原理、基尔霍夫衍射公式以及菲涅尔与夫琅禾费两种近似条件，帮助学习者建立从现象到公式、从公式到工程判据的完整认知链。
上图展示了衍射理论从物理原理到工程近似的完整框架。惠更斯-菲涅尔原理提供了“次波相干叠加”的物理图像，基尔霍夫公式给出了严格的数学表达，而菲涅尔与夫琅禾费近似则是在不同传播距离下的工程化简化。下面我们逐一深入剖析每个环节的数学内涵与适用边界。
4.1 惠更斯-菲涅尔原理：衍射的物理图像
4.1.1 核心原理
惠更斯-菲涅尔原理是理解衍射现象的基石。它包含两个关键要素：次波假设与相干叠加。惠更斯提出，波前上的每一点都可以看作一个新的球面波源（次波源）；菲涅尔则补充了最关键的一点——这些次波在空间某点相遇时会发生相干叠加，其复振幅的矢量叠加决定了该点的总光场。这一补充将几何光学中“光线”的概念升级为波动光学的“波前传播”。
数学上，设衍射屏上的孔径透射函数为 \( U_0(\xi, \eta) \)，观察平面上任一点 \( P(x, y) \) 的光场可以表示为：
U(x,y) = (1/jλ) ∬ U₀(ξ,η) · [exp(jkr)/r] · K(θ) dξ dη其中 \( r = \sqrt{(x-\xi)^2 + (y-\eta)^2 + z^2} \) 是次波源到观察点的距离，\( K(\theta) \) 为倾斜因子，表示次波在不同方向上的辐射强度差异。菲涅尔假设 \( K(\theta) = (1+\cos\theta)/2 \)，当观察距离远大于孔径尺寸时，\( \theta \approx 0 \)，倾斜因子近似为1。
4.1.2 判断方法：何时使用该原理
惠更斯-菲涅尔原理是衍射理论的通用起点，适用于任何衍射场景。但在实际工程计算中，直接使用该积分往往面临数学困难。判断是否需要从该原理出发，主要看以下两点：
- 孔径与波长的相对尺度：当孔径尺寸与波长可比拟时（如亚波长光栅），必须严格使用矢量衍射理论，标量惠更斯-菲涅尔原理失效。
- 观察距离的远近：距离越远，积分中的二次相位因子变化越平缓，越容易做近似处理。
4.1.3 案例拆解：单缝衍射的物理图像
以单缝衍射为例，缝宽 \( a = 0.1 \text{ mm} \)，波长 \( \lambda = 632.8 \text{ nm} \)（He-Ne激光），观察距离 \( z = 1 \text{ m} \)。将缝内每一点视为次波源，所有次波在观察屏上叠加。在缝中心附近，各次波光程差较小，干涉加强；在偏离中心的位置，光程差增大，干涉减弱甚至相消。这一过程直观展示了“次波相干叠加”的物理内涵。实际计算时，我们并不需要逐点叠加，而是利用傅里叶变换工具快速求解。
4.1.4 常见误区
- 误区一：将惠更斯-菲涅尔原理等同于几何光学中的“光线传播”。实际上，该原理强调的是波的叠加，而非直线传播。
- 误区二：忽略倾斜因子的作用。在近场或大角度衍射场景中，倾斜因子对衍射图样的边缘强度有明显影响。
- 误区三：认为次波源是真实存在的物理源。次波是数学构造，用于描述波前的传播行为，并非实际光源。
4.2 基尔霍夫衍射公式：数学严格化
4.2.1 核心原理
基尔霍夫利用格林定理，将惠更斯-菲涅尔原理从“物理假设”提升为“数学定理”。他通过设定衍射屏的边界条件（屏上开孔处场不为零，屏后阴影处场为零），推导出严格的衍射积分公式：
U(P) = (1/4π) ∬ [U₀ · ∂G/∂n - G · ∂U₀/∂n] dS其中 \( G \) 为格林函数，\( \partial/\partial n \) 表示沿法线方向的偏导数。在傍轴近似下，该公式可简化为：
U(x,y) = (exp(jkz)/jλz) · ∬ U₀(ξ,η) · exp{jk[(x-ξ)²+(y-η)²]/2z} dξ dη这一形式被称为菲涅尔衍射积分，是连接严格理论与工程近似的桥梁。
4.2.2 判断方法：基尔霍夫公式的适用边界
基尔霍夫公式的推导基于两个关键假设：标量场和边界条件。因此，其适用条件为：
- 衍射孔径的尺寸远大于波长（通常要求 \( a > 10\lambda \)）；
- 观察点距离衍射屏足够远，满足傍轴条件；
- 衍射屏为理想黑屏，即屏上开孔处场不受扰动。
在实际光学系统中，这些条件通常可以满足。但当孔径接近波长量级（如近场光学显微镜的探针）或涉及偏振效应时，必须使用矢量衍射理论（如严格耦合波分析RCWA）。
4.2.3 案例拆解：矩形孔径的菲涅尔衍射
考虑一个 \( 2\text{ mm} \times 2\text{ mm} \) 的方形孔径，波长 \( \lambda = 500 \text{ nm} \)，观察距离 \( z = 100 \text{ mm} \)。利用菲涅尔衍射积分，可以数值计算观察面上的光强分布。由于菲涅尔数 \( N_F = a^2/(\lambda z) = (1\text{ mm})^2 / (500\text{ nm} \times 100\text{ mm}) = 20 \)，属于典型的近场衍射。此时衍射图样呈现清晰的边缘振荡条纹，不能简单地用夫琅禾费近似处理。工程中，这类计算常用于微透镜阵列的焦深分析或MEMS光开关的衍射损耗评估。
4.2.4 常见误区
- 误区一：认为基尔霍夫公式是精确解。实际上，它仍是标量近似，忽略了电磁场的矢量特性。
- 误区二：混淆格林函数与传播子的概念。格林函数是数学工具，传播子（如 \( \exp(jkr)/r \)）是物理量。
- 误区三：在近场区域强行使用基尔霍夫公式而不检查边界条件。当观察距离小于波长时，该公式不再适用。
4.3 菲涅尔近似与夫琅禾费近似
4.3.1 核心原理：两种近似的数学本质
菲涅尔衍射积分中的二次相位因子 \( \exp\{jk[(x-\xi)^2+(y-\eta)^2]/2z\} \) 是计算的难点。根据观察距离的不同，可以对该相位因子做不同程度的简化：
| 近似类型 | 相位因子处理 | 数学形式 | 物理意义 |
|---|---|---|---|
| 菲涅尔近似 | 保留二次项 \( (x^2+y^2)/2z \) 和交叉项 \( (x\xi+y\eta)/z \) | 卷积积分 | 球面波前在观察面上的投影 |
| 夫琅禾费近似 | 忽略二次项，仅保留线性项 \( (x\xi+y\eta)/z \) | 傅里叶变换 | 平面波前在远场的投影 |
夫琅禾费近似的核心条件是 \( z \gg k(\xi^2+\eta^2)_{\max}/2 \)，即观察距离远大于孔径尺寸的平方除以波长。此时，衍射图样与孔径的傅里叶变换直接相关，这为光学信息处理（如空间滤波、图像识别）提供了理论基础。
4.3.2 判断方法：如何选择近似类型
工程中，我们通常通过以下步骤快速判断：
- 计算菲涅尔数 \( N_F = a^2/(\lambda z) \)；
- 若 \( N_F \ll 1 \)（例如小于0.1），采用夫琅禾费近似；
- 若 \( N_F \geq 1 \)，采用菲涅尔近似；
- 若 \( N_F \) 在0.1~1之间，需根据精度要求决定，必要时直接数值求解菲涅尔积分。
4.3.3 案例拆解：光栅光谱仪的衍射级次
以反射式光栅光谱仪为例，光栅常数 \( d = 1/1200 \text{ mm} \)，入射角 \( \alpha = 30^\circ \)，波长 \( \lambda = 500 \text{ nm} \)。根据光栅方程 \( d(\sin\alpha + \sin\beta) = m\lambda \)，可以计算各衍射级次的衍射角。在夫琅禾费近似下，衍射图样是孔径函数的傅里叶变换，光栅的周期结构导致离散的衍射级次。实际设计中，我们关心的是光谱分辨率 \( R = mN \)（\( N \) 为光栅刻线总数），这直接由夫琅禾费衍射的角宽度决定。若误用菲涅尔近似，会得到错误的衍射角分布，导致光谱仪设计失败。
4.3.4 常见误区与总结复盘
- 误区一：将菲涅尔数作为唯一判据。实际上，孔径形状、照明方式（平面波或球面波）也会影响近似的有效性。
- 误区二：忽略透镜的傅里叶变换作用。在4F系统中，即使物理距离很短，焦平面上的分布仍是夫琅禾费型。
- 误区三：混淆“远场”与“夫琅禾费区”。远场是物理概念，夫琅禾费近似是数学条件，两者并不完全等价。