克希荷夫繞射公式

在光學裏，**菲涅耳 - 克希荷夫繞射公式**（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{aligned}\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{aligned}}$ ；

其中， $\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{aligned}\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 \\&\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}}\\\end{aligned}}$ 。

### 有限尺寸波源

假設波源為有限尺寸，位於曲面 $\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'$ 。

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

## 純量理論

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

## 參閱

* 帕松光斑
* 克希荷夫繞射公式
* 克希荷夫積分定理

## 參考文獻