基爾霍夫衍射公式

在光學裏，菲涅耳-基爾霍夫衍射公式（Fresnel-Kirchoff's diffraction formula）可以應用於光波傳播的理論分析模型或數值分析模型。從菲涅耳-基爾霍夫衍射公式，可以推導出惠更斯－菲涅耳原理，並且解釋一些惠更斯－菲涅耳原理無法解釋的物理現象與結果。菲涅耳-基爾霍夫衍射公式常被稱為「基爾霍夫衍射公式」（Kirchoff's diffraction formula）。

從基爾霍夫積分定理，在假定一些近似之後，可以推導出菲涅耳-基爾霍夫衍射公式。

## 惠更斯原理

惠更斯原理是克里斯蒂安·惠更斯於1678年提出的關于波傳播的理論。惠更斯原理表明，假設在時間$t=t_0$由主波源Q0發射出的球面波，在時間$t=t_1$傳播到波前$\mathbb{S}$，那麼位於波前$\mathbb{S}$的每一個面元素向量$\mathrm{d}\mathbf{S}$都可以被視為一個次波源，所有從這些次波源發射出的次波，在之後時間$t=t_2$波前的包絡面就是主波源Q0所發射出的球面波在時間$t=t_2$的波前。

## 惠更斯-菲涅耳原理

波動有兩個基本屬性：

- 它是波擾的傳播
- 它具有時空周期性，能夠相干疊加。

惠更斯原理只闡述了前一條屬性，奧古斯丁·菲涅耳將惠更斯提出的次波的概念加以延伸，提出用「次波相干疊加」的點子來解釋衍射現象，這就是惠更斯-菲涅耳原理。這原理表明，波前$\mathbb{S}$的每個面元素向量$\mathrm{d}\mathbf{S}'$都可以視為次波源，它們會發射出次波，在空間任意一點P的波擾是所有這些次波在該點P的相干疊加。設定位於波前$\mathbb{S}$的任意一點Q，它在點P貢獻的複振幅為$\mathrm{d}\psi(\mathbf{r},\mathbf{r}')$；其中，$\mathbf{r}$、$\mathbf{r}'$分別為點P、點Q的位置。在點P的總波擾為

$\psi(\mathbf{r})=\oint_{\mathbb{S}}\,\mathrm{d}\psi(\mathbf{r},\mathbf{r}')$。

為了將這公式具體化，菲涅耳憑藉直覺對$\mathrm{d}\psi(\mathbf{r},\mathbf{r}')$作出了如下假設：

$\mathrm{d}\psi(\mathbf{r},\mathbf{r}')\propto\mathrm{d}S'$。

$\mathrm{d}\psi(\mathbf{r},\mathbf{r}')\propto\psi(\mathbf{r}')$。

- 次波源發射出的次波應是球面波，其中$k$是波數：

$\mathrm{d}\psi(\mathbf{r},\mathbf{r}')\propto\frac{e^{ikR}}{R}$；其中，$\mathbf{R}=\mathbf{r}-\mathbf{r}'$是從點Q到點P的位移向量。

- 次波源發射出的次波是各向異性的。假設$\hat{\mathbf{n}}$是與面元素向量$\mathrm{d}\mathbf{S}'$同方向的單位向量，$\chi$是$\hat{\mathbf{n}}$與$\hat{\mathbf{R}}$之間的夾角，則傾斜因子$K(\chi)$與$\mathrm{d}\psi(\mathbf{r},\mathbf{r}')$的關係為

$\mathrm{d}\psi(\mathbf{r},\mathbf{r}')\propto K(\chi)$。

根據以上假設可以得到如下菲涅耳衍射積分公式

$\psi(\mathbf{r})=c\oint_{\mathbb{S}}\,\psi(\mathbf{r}')K(\chi)\frac{e^{ikR}}{R}\,\mathrm{d}S'$。

其中，$c$是比例常數。

## 菲涅耳-基爾霍夫衍射公式

在菲涅耳衍射積分公式提出六十餘年後，古斯塔夫·基爾霍夫用嚴格的數學理論推導出菲涅耳-基爾霍夫衍射公式：

$\psi(\mathbf{r})=-\ \frac{i\psi_0}{2\lambda}\oint_{\mathbb{S}} \left(\frac{e^{ik(r'+R)}}{r'R} \right) [\cos\alpha+\cos\chi] \,\mathrm{d}S'$。

其中，$\alpha$、$\chi$分別是$\hat{\mathbf{r}'}$、$\hat{\mathbf{R}}$與$\hat{\mathbf{n}}$之間的夾角。

推論從點光源Q0發射的單色光波，其波擾的數值大小與傳播距離成反比，在位置$\mathbf{r}'$以方程式表達為$\psi(\mathbf{r}') = \psi_0 e^{ikr'}/r'$。又在其發射出的球面波的波前任意位置，$\hat{\mathbf{r}'}$與$\hat{\mathbf{n}}$同向，夾角$\alpha=0$。設定比例常數$c=-i/\lambda$。

$K(\chi)=(1+\cos\chi)/2$，則可得到菲涅耳衍射積分公式。

## 嚴格導引

基爾霍夫積分定理應用格林第二恆等式來推導出齊次波動方程式的解答，這解答是以波動方程式在任意閉合曲面$\mathbb{S}$的每一個點的解答和其一階導數來表達。

對於單頻率波，解答為

$\psi(\mathbf{r})=\frac{1}{4\pi}\oint_\mathbb{S} \left[\psi(\mathbf{r}')\nabla'\left( \frac{e^{ikR}}{R} \right) -\left( \frac{e^{ikR}}{R} \right)\nabla'\psi(\mathbf{r}')\right] \cdot\,\mathrm{d}\mathbf{S}'$。

或者

$\psi(\mathbf{r})=\frac{1}{4\pi}\oint_\mathbb{S} \left[\psi(\mathbf{r}')\frac{\partial}{\partial n'}\left( \frac{e^{ikR}}{R} \right) -\left( \frac{e^{ikR}}{R} \right)\frac{\partial \psi(\mathbf{r}')}{\partial n'}\right] \,\mathrm{d}S'$。

其中，$\mathbf{r}$、$\mathbf{r}'$分別是從點Q0到點P、點Q的位移向量，$\psi(\mathbf{r})$是在點P的波擾，$\mathbf{R}=\mathbf{r}- \mathbf{r}'$是從點Q到點P的位移向量，$R$是其數值大小，$k$是波數，$\nabla'$是對於源位置$\mathbf{r}'$的梯度，$\mathrm{d}\mathbf{S}'$是從閉合曲面$\mathbb{S}$向外指出的微小面元素向量，$\frac{\partial}{\partial n'}$是閉合曲面$\mathbb{S}$的法向導數。

在推導基爾霍夫衍射公式的過程中，基爾霍夫做了以下假定：

- 點波源與孔隙之間的距離$r'$超大於波長$\lambda=2\pi/k$。
- $R$超大於波長$\lambda$。

### 點波源

從點波源Q0發射的單頻率波，其能量與傳播距離平方成反比，波擾的數值大小與傳播距離成反比，在點Q的波擾以方程式表達為

$\psi(\mathbf{r}') = \psi_0\frac{e^{ikr'}}{r'}$。

其中，$\psi_0$是複值波幅。

假設點P在閉合曲面$\mathbb{S}$之外，應用基爾霍夫積分定理的方程式，可以得到在點P的波擾：

$\psi(\mathbf{r}) =\frac{\psi_0}{4\pi}\oint_{\mathbb{S}} \left[\left(\frac{e^{ikr'}}{r'} \right)\nabla'\left( \frac{e^{ikR}}{R} \right) -\left( \frac{e^{ikR}}{R} \right)\nabla'\left(\frac{e^{ikr'}}{r'} \right)\right] \cdot\,\hat{\mathbf{n}}\mathrm{d}S'$。

其中，$\hat{\mathbf{n}}$是與$\mathrm{d}\mathbf{S}'$同方向的單位向量。

注意到球面出射波的梯度為

$\nabla'\left( \frac{e^{ikR}}{R} \right) =-\left( \frac{e^{ikR}}{R} \right)\left(ik-\ \frac{1}{R}\right)\hat{\mathbf{R}}$。

$\nabla'\left(\frac{e^{ikr'}}{r'} \right) =\left( \frac{e^{ikr'}}{r'} \right)\left(ik-\ \frac{1}{r'}\right)\hat{\mathbf{r}'}$。

從基爾霍夫所做的假定，$k\gg 1/R$、$k \gg 1/r'$（例如，假設距離大約為1mm，則對於波長在0.4μm至0.7μm之間的可見光，可以做這假定；但對於波長在1mm至1m之間的微波，這假定不適用），則上述兩個公式近似為

$\nabla'\left( \frac{e^{ikR}}{R} \right) =-ik\left( \frac{e^{ikR}}{R} \right)\hat{\mathbf{R}}$。

$\nabla'\left(\frac{e^{ikr'}}{r'} \right) =ik\left( \frac{e^{ikr'}}{r'} \right)\hat{\mathbf{r}'}$。

所以，在點P的波擾

$\begin{align}\psi(\mathbf{r}) & =-\ \frac{i\psi_0}{2\lambda}\oint_{\mathbb{S}} \left(\frac{e^{ik(r'+R)}}{r'R} \right)(\hat{\mathbf{r}'}\cdot\hat{\mathbf{n}}+\hat{\mathbf{R}}\cdot\hat{\mathbf{n}}) \,\mathrm{d}S' \\ & =-\ \frac{i\psi_0}{2\lambda}\oint_{\mathbb{S}} \left(\frac{e^{ik(r'+R)}}{r'R} \right) (\cos\alpha+\cos\chi) \,\mathrm{d}S' \\ \end{align}$。

其中，$\alpha$、$\chi$分別是$\hat{\mathbf{r}'}$、$\hat{\mathbf{R}}$與$\hat{\mathbf{n}}$之間的夾角。

這就是菲涅耳-基爾霍夫衍射公式，或基爾霍夫衍射公式。

### 傾斜因子

如右圖所示，假設閉合曲面$\mathbb{S}$是圓球面，點波源Q0與圓球面$\mathbb{S}$的圓心同點。在圓球面$\mathbb{S}$的任意位置，$\hat{\mathbf{r}'}$與$\hat{\mathbf{n}}$同向，所以。

$\cos\alpha=1$。

注意到$r'$是圓球面$\mathbb{S}$的半徑，對於這積分，$r'$值不變，可以從積分裏提出。在點P的波擾為

$\psi(\mathbf{r}) =-\ \frac{i\psi(\mathbf{r}')}{\lambda}\oint_{\mathbb{S}} \left(\frac{e^{ikR}}{R}\right)K(\chi)\,\mathrm{d}S'$。

其中，$K(\chi)=\frac{1+\cos\chi}{2}$為傾斜因子。

應用惠更斯－菲涅耳原理，所得到在點P的波擾的方程式，就是這方程式。但是，惠更斯－菲涅耳原理無法解釋相位差與傾斜因子的物理原因。傾斜因子使得次波的波幅會因為傳播方向而不同；朝著主波方向，波幅較大；逆著主波方向，波幅較小。這解釋了為甚麼波動只會朝著前方傳播的物理現象。

### 惠更斯－菲涅耳原理

仔細詮釋惠更斯－菲涅耳原理的方程式：從點波源Q0發射的波幅為$\psi_0$的球面波，在點Q的波擾為$\psi(\mathbf{r}')=\psi_0 e^{ikr'}/r'$；而從點Q發射的次波，將傾斜因子與相位差納入考量，所貢獻出的波擾，在點P為

$-\ \frac{i\psi(\mathbf{r}')}{\lambda}\left(\frac{e^{ikR}}{R}\right)K(\chi)$。

總合所有與點Q同波前的點次波源在點P所貢獻出的波擾，就可以得到$\psi(\mathbf{r})$。

換另一種直接方法來詮釋，從點波源Q0發射的球面波，在點P的波擾為

$\psi(\mathbf{r}) =\psi_0 \frac{e^{ikr}}{r}$。

假若這兩種詮釋都正確，則從這兩種$\psi(\mathbf{r})$的表達式分別計算出的結果，應該可以被核對為相等：

$\psi_0 \frac{e^{ikr}}{r}=-\ \frac{i}{\lambda}\oint_{\mathbb{S}_1} \left(\frac{\psi_0 e^{ikr'}}{r'}\right) \left(\frac{e^{ikR}}{R}\right)K(\chi)\,\mathrm{d}S'$。

為了簡易計算，假設$r\gg r'$，則以下近似成立：

$R\approx r-r' \cos(\theta)$。

$R^2\approx r^2-2rr' \cos(\theta)$。

$K(\chi)=\frac{1+\cos(\chi)}{2}\approx \frac{1+\cos(\theta)}{2}$。

其中，$\theta$為$\mathbf{r}'$與$\mathbf{r}$之間的夾角。

所以，在點P的波擾可以近似為

$\begin{align}\psi(\mathbf{r}) & \approx -\ \frac{i\psi_0}{2\lambda} \frac{e^{ik(r'+r)}}{r'r}\int_0^{\pi} e^{-ikr'\cos(\theta)}[1+\cos(\theta)]2\pi r'^2\sin(\theta)\mathrm{d}\theta \\ <pre> & \approx -\ \frac{ik\psi_0 r'e^{ik(r'+r)}}{2r}\int_0^{\pi} e^{-ikr'\cos(\theta)}[1+\cos(\theta)]\sin(\theta)\mathrm{d}\theta \\ & \approx -\ \frac{ik\psi_0 r'e^{ik(r'+r)}}{2r}\ \frac{2[\sin(kr')+i\cos(kr')]}{kr'} \\ & \approx \psi_0\frac{e^{ikr}}{r} \\ </pre> \end{align}$。

### 有限尺寸波源

假設波源為有限尺寸，位於曲面$\mathbb{S}$的波擾表達為$\psi(\mathbf{r}')$，則位於點P的波擾為

$\psi(\mathbf{r})=\frac{1}{4\pi}\oint_\mathbb{S} \left[\psi(\mathbf{r}')\frac{\partial}{\partial n'}\left( \frac{e^{ikR}}{R} \right) -\left( \frac{e^{ikR}}{R} \right)\frac{\partial \psi(\mathbf{r}')}{\partial n'}\right] \,\mathrm{d}S'$。

假定$k\gg 1/R$，則

$\psi(\mathbf{r})=-\ \frac{1}{4\pi}\oint_\mathbb{S}\left( \frac{e^{ikR}}{R} \right) \left[ik\psi(\mathbf{r}')\cos{(\hat{\mathbf{R}},\hat{\mathbf{n}}}) +\frac{\partial \psi(\mathbf{r}')}{\partial n'}\right] \,\mathrm{d}S'$。

這是基爾霍夫衍射公式最廣義的形式。解析涉及到有限尺寸波源的問題，必須用體積分來將波源的每一點所給出的貢獻總合在一起。

## 純量理論

光波是傳播於空間的電磁輻射，理當被視為一種電磁場向量現象。但是，基爾霍夫的理論是純量理論，將光波當作純量處理，這可能會造成偏差。因此，物理學者做了很多實驗來檢查結果是否準確。他們發現，只要孔徑尺寸比波長大很多、孔徑與觀察屏之間的距離不很近，則使用純量理論可以得到相當準確的答案。但是對於某些問題，例如高解析度光柵衍射，純量理論就不適用，必須使用向量理論。

## 參閱

- 帕松光斑
- 基爾霍夫衍射公式
- 基爾霍夫積分定理