4、基尔霍夫衍射理论:基尔霍夫积分定理、边界条件假设、菲涅尔-基尔霍夫衍射公式的推导思路

基尔霍夫衍射理论是标量衍射理论的基石,它将惠更斯-菲涅尔原理从一种物理直觉提升为严格的数学表述。该理论的核心在于利用格林定理,将空间中某一点的复振幅场与其边界上的场值及其法向导数联系起来。本章将系统梳理基尔霍夫积分定理的数学基础、边界条件假设的物理意义,以及如何从该定理出发推导出菲涅尔-基尔霍夫衍射公式,并讨论其适用边界与常见误区。

基尔霍夫衍射理论推导路径与边界条件 基尔霍夫积分定理 格林定理 + 亥姆霍兹方程 U(P) = (1/4π)∮[U∂G/∂n - G∂U/∂n]dS 自由空间格林函数: G = exp(ikr)/r 适用于任意封闭曲面 选择 曲面 基尔霍夫边界条件 ① 开孔Σ上:U, ∂U/∂n 等于无屏时的入射场 ② 挡光屏背面Σ₁: U = 0, ∂U/∂n = 0 ⚠ 数学上不相容,但物理近似 (基尔霍夫近似) 代入 化简 菲涅尔-基尔霍夫 衍射公式 U(P) = (A/iλ)∬Σ [exp(ikr)/r]· [cos(n,r)+cos(n,l)]/2 dS 倾斜因子 K(θ) = (cosθ₀+cosθ)/2 自动包含,无需人为假设 理论成立的关键假设与限制 标量近似 忽略矢量偏振,孔径≫λ 观察距离远 r ≫ λ,r ≫ 孔径尺寸 边界条件近似 屏上U和∂U/∂n同时为零 单色波假设 时间因子exp(-iωt) 推导核心:将积分定理应用于衍射屏,利用边界条件消去屏后贡献, 得到仅对开孔区域的面积分表达式 适用场景:夫琅禾费衍射、菲涅尔衍射(近场/远场) 不适用:孔径尺寸接近λ、极短距离、强偏振效应

4.1 基尔霍夫积分定理:从格林定理到衍射积分

基尔霍夫积分定理是连接波动方程与衍射积分的桥梁。其数学起点是亥姆霍兹方程,即单色光波在均匀介质中传播时,复振幅 \( U(P) \) 满足:

∇²U + k²U = 0

其中 \( k = 2\pi/\lambda \) 为波数。为了求解空间任意一点 \( P \) 的场值,基尔霍夫引入了辅助函数——自由空间格林函数 \( G(P, Q) = \frac{\exp(ikr)}{r} \),其中 \( r = |\mathbf{r}_P - \mathbf{r}_Q| \) 表示源点 \( Q \) 到观察点 \( P \) 的距离。该格林函数本身也满足亥姆霍兹方程(除 \( r=0 \) 的奇点外)。

利用格林第二恒等式,将 \( U \) 和 \( G \) 代入,并对包围 \( P \) 点的任意封闭曲面 \( S \) 积分,可以得到:

U(P) = (1/4π) ∮_S [ U ∂G/∂n - G ∂U/∂n ] dS

这就是基尔霍夫积分定理。它表明:封闭曲面 \( S \) 上任意一点的场值 \( U \) 及其法向导数 \( ∂U/∂n \),共同决定了曲面内部任意一点 \( P \) 的场值。这一结论在数学上是严格的,不依赖于任何物理近似。然而,实际衍射问题中,我们通常只知道有限孔径上的场分布,而无法获知整个封闭曲面上的完整信息,因此必须引入边界条件假设。

核心理解:基尔霍夫积分定理本质上是一个数学恒等式,它将体积分转化为面积分。其物理意义在于:曲面上的每一点都可以看作一个次级波源,而观察点的场是所有这些次级波源贡献的相干叠加。这与惠更斯-菲涅尔原理的物理图像完全吻合,但基尔霍夫给出了严格的数学形式。

4.2 边界条件假设:基尔霍夫近似的合理性与局限

在实际衍射问题中,我们面对的是带有开孔的不透明屏。设屏位于 \( z=0 \) 平面,开孔区域为 \( \Sigma \),屏的挡光部分为 \( \Sigma_1 \)。为了应用基尔霍夫积分定理,需要选择一个封闭曲面,通常由三部分组成:开孔 \( \Sigma \)、屏的背面 \( \Sigma_1 \)、以及一个半径趋于无穷大的半球面 \( \Sigma_2 \)。

基尔霍夫提出了以下边界条件假设:

  • 开孔区域 \( \Sigma \) 上:场值 \( U \) 和法向导数 \( ∂U/∂n \) 与不存在屏时完全相同,即等于入射波的场分布。
  • 挡光屏背面 \( \Sigma_1 \) 上:场值 \( U = 0 \) 且法向导数 \( ∂U/∂n = 0 \)。
  • 无穷远半球面 \( \Sigma_2 \) 上:由于索末菲辐射条件,当 \( R \to \infty \) 时,该面上的积分贡献趋于零,可以忽略。

这些假设合称为基尔霍夫边界条件。其中,屏背面场及其导数同时为零的假设在数学上是不自洽的——如果 \( U \) 和 \( ∂U/∂n \) 在某个有限区域内同时为零,那么由亥姆霍兹方程的唯一性定理,整个空间场都应为零。然而,基尔霍夫近似在物理上却取得了巨大成功,原因在于:对于孔径远大于波长的情形,屏边缘的衍射效应相对较弱,屏背面确实几乎没有光场穿透,因此该近似在积分意义上误差很小。

边界条件项 数学表述 物理含义 适用条件
开孔 \( \Sigma \) \( U = U_{inc} \), \( ∂U/∂n = ∂U_{inc}/∂n \) 入射波不受屏扰动 孔径尺寸 ≫ λ
屏背面 \( \Sigma_1 \) \( U = 0 \), \( ∂U/∂n = 0 \) 屏完全不透光 理想导体屏近似
无穷远面 \( \Sigma_2 \) 积分贡献 → 0 辐射条件 观察点距屏有限远
重要提醒:基尔霍夫边界条件在数学上是不自洽的,它同时规定了 \( U \) 和 \( ∂U/∂n \) 为零,这违反了亥姆霍兹方程解的唯一性。但该近似在孔径远大于波长时误差很小。若孔径尺寸接近波长或小于波长,必须采用更严格的矢量衍射理论(如瑞利-索末菲理论)。

4.3 菲涅尔-基尔霍夫衍射公式的推导思路

将基尔霍夫边界条件代入积分定理,经过一系列代数运算,可以得到菲涅尔-基尔霍夫衍射公式。推导的关键步骤如下:

  1. 选择格林函数:采用自由空间格林函数 \( G = \exp(ikr)/r \),并计算其法向导数 \( ∂G/∂n = \cos(\mathbf{n}, \mathbf{r}) \cdot (ik - 1/r) \cdot \exp(ikr)/r \)。在远场近似下(\( r \gg \lambda \)),\( 1/r \) 项可忽略,得到 \( ∂G/∂n \approx ik \cos(\mathbf{n}, \mathbf{r}) \cdot \exp(ikr)/r \)。
  2. 代入入射波:设入射波为从点源 \( S \) 发出的球面波 \( U_{inc} = A \exp(ikl)/l \),其中 \( l \) 是源点到孔径上某点的距离。同样计算其法向导数 \( ∂U_{inc}/∂n \approx ik \cos(\mathbf{n}, \mathbf{l}) \cdot A \exp(ikl)/l \)。
  3. 化简积分表达式:将上述表达式代入基尔霍夫积分定理,并利用边界条件消去屏背面和无穷远面的贡献,仅保留对开孔 \( \Sigma \) 的面积分。
  4. 引入倾斜因子:整理后得到:
U(P) = (A / iλ) ∬_Σ [ exp(ik(l+r)) / (l·r) ] · [ cos(n,l) + cos(n,r) ] / 2 dS

其中 \( \cos(\mathbf{n}, \mathbf{l}) \) 是入射光线与屏法线的夹角余弦,\( \cos(\mathbf{n}, \mathbf{r}) \) 是衍射光线与屏法线的夹角余弦。两者之和的一半即为倾斜因子 \( K(\theta) = (\cos\theta_0 + \cos\theta)/2 \)。

该公式的物理意义非常清晰:观察点 \( P \) 的场是孔径上所有次级波源贡献的叠加,每个次级波源的振幅正比于入射波振幅,并乘以倾斜因子;相位由传播距离 \( l+r \) 决定。倾斜因子的引入是基尔霍夫理论的重要贡献——它自动给出了惠更斯-菲涅尔原理中需要人为假设的倾斜因子,且当 \( \theta = 0 \)(正前方)时 \( K = 1 \),当 \( \theta = \pi \)(反向)时 \( K = 0 \),解释了为什么不存在反向传播的次级波。

推导技巧:在推导过程中,注意保留 \( 1/r \) 和 \( 1/l \) 项时,要区分哪些项在远场条件下可以忽略。通常保留 \( 1/(lr) \) 作为振幅项,而将 \( ik \) 项保留在指数前。另外,法向导数的计算要特别注意方向——\( \mathbf{n} \) 指向观察点一侧为正方向。

4.4 案例拆解:点源照射圆孔的菲涅尔衍射

为了加深理解,我们考察一个经典场景:单色点源 \( S \) 照射一个半径为 \( a \) 的圆孔,观察点 \( P \) 位于圆孔正后方距离 \( z \) 处。设源距孔径为 \( l_0 \),观察点距孔径为 \( r_0 \)。

初始条件:波长 \( \lambda = 632.8 \) nm(He-Ne激光),孔径半径 \( a = 1 \) mm,源距 \( l_0 = 1 \) m,观察距离 \( r_0 = 1 \) m。此时菲涅尔数 \( N_F = a^2 / (\lambda z) \approx 1.58 \),属于菲涅尔衍射区。

分析过程:由于圆孔具有轴对称性,采用极坐标 \( (\rho, \phi) \) 表示孔径上的点。入射距离 \( l = \sqrt{l_0^2 + \rho^2} \approx l_0 + \rho^2/(2l_0) \),衍射距离 \( r = \sqrt{r_0^2 + \rho^2} \approx r_0 + \rho^2/(2r_0) \)。倾斜因子在傍轴条件下近似为 1。代入菲涅尔-基尔霍夫公式,得到轴上点的场为:

U(P) = (A / iλ) · exp(ik(l₀+r₀)) / (l₀·r₀) · ∫₀^a exp[ikρ²(1/l₀+1/r₀)/2] · 2πρ dρ

令 \( u = \rho^2 \),则 \( d\rho = du/(2\rho) \),积分变为:

U(P) = (A / iλ) · exp(ik(l₀+r₀)) / (l₀·r₀) · π ∫₀^{a²} exp[iku(1/l₀+1/r₀)/2] du

该积分可直接计算,结果为:

U(P) = (A / iλ) · exp(ik(l₀+r₀)) / (l₀·r₀) · π · [exp(ik a²(1/l₀+1/r₀)/2) - 1] / [ik(1/l₀+1/r₀)/2]

结果分析:当 \( a^2(1/l_0+1/r_0)/2 = m\lambda \)(\( m \) 为整数)时,轴上点出现亮斑;当 \( a^2(1/l_0+1/r_0)/2 = (m+1/2)\lambda \) 时,出现暗斑。这解释了菲涅尔衍射中轴上点的亮暗交替现象。若 \( a \to \infty \)(无屏),则积分趋于 \( 2\pi \delta \) 函数,恢复自由传播的球面波。

复盘经验:该案例展示了菲涅尔-基尔霍夫公式的实际计算流程。关键在于:① 正确写出 \( l \) 和 \( r \) 的傍轴近似;② 利用对称性简化积分;③ 注意积分上下限与孔径边界的关系。实际工程中,若孔径形状不规则,通常需要数值积分(如快速傅里叶变换法)。

4.5 常见误区与总结复盘

在学习基尔霍夫衍射理论时,初学者常陷入以下误区:

  • 误区一:认为基尔霍夫边界条件是严格成立的。事实上,屏背面 \( U=0 \) 和 \( ∂U/∂n=0 \) 同时成立在数学上不可能,它只是物理近似。在孔径边缘附近,该近似误差较大,会产生边缘衍射波。
  • 误区二:混淆菲涅尔-基尔霍夫公式与瑞利-索末菲公式。两者对边界条件的处理不同:瑞利-索末菲理论只规定 \( U \) 或 \( ∂U/∂n \) 之一为零,避免了数学不自洽,但两者在远场和傍轴条件下结果几乎一致。
  • 误区三:忽略倾斜因子的作用。在非傍轴条件下,倾斜因子 \( K(\theta) \) 对衍射图样的影响不可忽略,尤其是大角度衍射时。
  • 误区四:将标量理论用于亚波长孔径。当孔径尺寸与波长可比时,矢量效应(偏振)变得显著,标量理论失效。
对比维度 基尔霍夫理论 瑞利-索末菲理论
边界条件 U 和 ∂U/∂n 同时为零 仅 U=0 或仅 ∂U/∂n=0
数学自洽性 不自洽 自洽
倾斜因子 (cosθ₀+cosθ)/2 cosθ 或 cosθ₀
适用场景 傍轴、远场、孔径≫λ 更广,包括非傍轴
工程常用度 教学常用,物理直观 数值计算更稳定
总结:基尔霍夫衍射理论的核心贡献在于:① 用严格的数学推导将惠更斯-菲涅尔原理转化为可计算的积分公式;② 自动给出倾斜因子,解释了次级波传播的方向性;③ 建立了从边界条件到衍射积分的完整方法论。尽管存在数学上的不自洽,但在绝大多数光学工程场景(孔径远大于波长、傍轴近似成立)中,该理论给出了与实验高度吻合的结果。后续章节中,我们将在此基础上进一步推导菲涅尔近似和夫琅禾费近似,它们将大大简化计算,并适用于更广泛的工程场景。