基爾霍夫衍射公式

在光學裏，菲涅耳-基爾霍夫衍射公式（Fresnel-Kirchoff's diffraction formula）可以應用於光波傳播的理論分析模型或數值分析模型。從菲涅耳-基爾霍夫衍射公式，可以推導出惠更斯－菲涅耳原理，並且解釋一些惠更斯－菲涅耳原理無法解釋的物理現象與結果。菲涅耳-基爾霍夫衍射公式常被稱為「基爾霍夫衍射公式」（Kirchoff's diffraction formula）。 從基爾霍夫積分定理，在假定一些近似之後，可以推導出菲涅耳-基爾霍夫衍射公式。 惠更斯原理是克里斯蒂安·惠更斯于1678年提出的关于波传播的理论。惠更斯原理表明，假设在时间 t = t 0 {\displaystyle t=t_{0}} 由主波源Q0发射出的球面波，在时间 t = t 1 {\displaystyle t=t_{1}} 传播到波前 S {\displaystyle \mathbb {S} } ，那么位於波前 S {\displaystyle \mathbb {S} } 的每一个面元素向量 d S {\displaystyle \mathrm {d} \mathbf {S} } 都可以被视为一个次波源，所有从这些次波源发射出的次波，在之后时间 t = t 2 {\displaystyle t=t_{2}} 波前的包络面就是主波源Q0所发射出的球面波在时间 t = t 2 {\displaystyle t=t_{2}} 的波前。 波动有两个基本属性： 它是波扰的传播 它具有时空周期性，能够相干叠加。 惠更斯原理只阐述了前一条属性，奧古斯丁·菲涅耳将惠更斯提出的次波的概念加以延伸，提出用“次波相干叠加”的点子来解释衍射现象，这就是惠更斯-菲涅耳原理。这原理表明，波前 S {\displaystyle \mathbb {S} } 的每个面元素向量 d S ′ {\displaystyle \mathrm {d} \mathbf {S} '} 都可以视为次波源，它们会发射出次波，在空间任意一点P的波扰是所有这些次波在该点P的相干叠加。设定位於波前 S {\displaystyle \mathbb {S} } 的任意一点Q，它在点P贡献的複振幅为 d ψ ( r , r ′ ) {\displaystyle \mathrm {d} \psi (\mathbf {r} ,\mathbf {r} ')} ；其中， r {\displaystyle \mathbf {r} } 、 r ′ {\displaystyle \mathbf {r} '} 分别为点P、点Q的位置。在点P的总波扰为 ψ ( r ) = ∮ S d ψ ( r , r ′ ) {\displaystyle \psi (\mathbf {r} )=\oint _{\mathbb {S} }\,\mathrm {d} \psi (\mathbf {r} ,\mathbf {r} ')} 。 为了将这公式具体化，菲涅耳凭借直觉对 d ψ ( r , r ′ ) {\displaystyle \mathrm {d} \psi (\mathbf {r} ,\mathbf {r} ')} 作出了如下假设： 它应当正比于面元素的面积： d ψ ( r , r ′ ) ∝ d S ′ {\displaystyle \mathrm {d} \psi (\mathbf {r} ,\mathbf {r} ')\propto \mathrm {d} S'} 。 它应当正比于次波源的複振幅： d ψ ( r , r ′ ) ∝ ψ ( r ′ ) {\displaystyle \mathrm {d} \psi (\mathbf {r} ,\mathbf {r} ')\propto \psi (\mathbf {r} ')} 。 次波源发射出的次波应是球面波，其中 k {\displaystyle k} 是波数： d ψ ( r , r ′ ) ∝ e i k R R {\displaystyle \mathrm {d} \psi (\mathbf {r} ,\mathbf {r} ')\propto {\frac {e^{ikR}}{R}}} ；其中， R = r − r ′ {\displaystyle \mathbf {R} =\mathbf {r} -\mathbf {r} '} 是从点Q到点P的位移向量。 次波源发射出的次波是各向异性的。假设 n ^ {\displaystyle {\hat {\mathbf {n} }}} 是与面元素向量 d S ′ {\displaystyle \mathrm {d} \mathbf {S} '} 同方向的单位向量， χ {\displaystyle \chi } 是 n ^ {\displaystyle {\hat {\mathbf {n} }}} 与 R ^ {\displaystyle {\hat {\mathbf {R} }}} 之間的夾角，则倾斜因子 K ( χ ) {\displaystyle K(\chi )} 与 d ψ ( r , r ′ ) {\displaystyle \mathrm {d} \psi (\mathbf {r} ,\mathbf {r} ')} 的关系为 d ψ ( r , r ′ ) ∝ K ( χ ) {\displaystyle \mathrm {d} \psi (\mathbf {r} ,\mathbf {r} ')\propto K(\chi )} 。 根据以上假设可以得到如下菲涅耳衍射积分公式 ψ ( r ) = c ∮ S ψ ( r ′ ) K ( χ ) e i k R R d S ′ {\displaystyle \psi (\mathbf {r} )=c\oint _{\mathbb {S} }\,\psi (\mathbf {r} ')K(\chi ){\frac {e^{ikR}}{R}}\,\mathrm {d} S'} ； 其中， c {\displaystyle c} 是比例常数。 在菲涅耳衍射积分公式提出六十余年后，古斯塔夫·基爾霍夫用严格的数学理论推导出菲涅耳-基尔霍夫衍射公式： ψ ( r ) = −   i ψ 0 2 λ ∮ S ( e i k ( r ′ + R ) r ′ R ) [ cos ⁡ α + cos ⁡ χ ] d S ′ {\displaystyle \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'} ； 其中， α {\displaystyle \alpha } 、 χ {\displaystyle \chi } 分別是 r ′ ^ {\displaystyle {\hat {\mathbf {r} '}}} 、 R ^ {\displaystyle {\hat {\mathbf {R} }}} 與 n ^ {\displaystyle {\hat {\mathbf {n} }}} 之間的夾角。 推论從點光源Q0發射的單色光波，其波擾的數值大小與傳播距離成反比，在位置 r ′ {\displaystyle \mathbf {r} '} 以方程式表達為 ψ ( r ′ ) = ψ 0 e i k r ′ / r ′ {\displaystyle \psi (\mathbf {r} ')=\psi _{0}e^{ikr'}/r'} 。又在其發射出的球面波的波前任意位置， r ′ ^ {\displaystyle {\hat {\mathbf {r} '}}} 與 n ^ {\displaystyle {\hat {\mathbf {n} }}} 同向，夾角 α = 0 {\displaystyle \alpha =0} 。設定比例常数 c = − i / λ {\displaystyle c=-i/\lambda } ， K ( χ ) = ( 1 + cos ⁡ χ ) / 2 {\displaystyle K(\chi )=(1+\cos \chi )/2} ，則可得到菲涅耳衍射积分公式。 基爾霍夫積分定理應用格林第二恆等式來推導出齊次波動方程式的解答，這解答是以波動方程式在任意閉合曲面 S {\displaystyle \mathbb {S} } 的每一個點的解答和其一階導數來表達。 對於單頻率波，解答為 ψ ( r ) = 1 4 π ∮ S [ ψ ( r ′ ) ∇ ′ ( e i k R R ) − ( e i k R R ) ∇ ′ ψ ( r ′ ) ] ⋅ d S ′ {\displaystyle \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} '} ， 或者 ψ ( r ) = 1 4 π ∮ S [ ψ ( r ′ ) ∂ ∂ n ′ ( e i k R R ) − ( e i k R R ) ∂ ψ ( r ′ ) ∂ n ′ ] d S ′ {\displaystyle \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'} ； 其中， r {\displaystyle \mathbf {r} } 、 r ′ {\displaystyle \mathbf {r} '} 分別是從點Q0到點P、點Q的位移向量， ψ ( r ) {\displaystyle \psi (\mathbf {r} )} 是在點P的波擾， R = r − r ′ {\displaystyle \mathbf {R} =\mathbf {r} -\mathbf {r} '} 是從點Q到點P的位移向量， R {\displaystyle R} 是其數值大小， k {\displaystyle k} 是波數， ∇ ′ {\displaystyle \nabla '} 是對於源位置 r ′ {\displaystyle \mathbf {r} '} 的梯度， d S ′ {\displaystyle \mathrm {d} \mathbf {S} '} 是從閉合曲面 S {\displaystyle \mathbb {S} } 向外指出的微小面元素向量， ∂ ∂ n ′ {\displaystyle {\frac {\partial }{\partial n'}}} 是閉合曲面 S {\displaystyle \mathbb {S} } 的法向導數。 在推導基爾霍夫衍射公式的過程中，基爾霍夫做了以下假定： 點波源與孔隙之間的距離 r ′ {\displaystyle r'} 超大於波長 λ = 2 π / k {\displaystyle \lambda =2\pi /k} 。 R {\displaystyle R} 超大於波長 λ {\displaystyle \lambda } 。 從點波源Q0發射的單頻率波，其能量與傳播距離平方成反比，波擾的數值大小與傳播距離成反比，在點Q的波擾以方程式表達為 ψ ( r ′ ) = ψ 0 e i k r ′ r ′ {\displaystyle \psi (\mathbf {r} ')=\psi _{0}{\frac {e^{ikr'}}{r'}}} ； 其中， ψ 0 {\displaystyle \psi _{0}} 是複值波幅。 假設點P在閉合曲面 S {\displaystyle \mathbb {S} } 之外，應用基爾霍夫積分定理的方程式，可以得到在點P的波擾： ψ ( r ) = ψ 0 4 π ∮ S [ ( e i k r ′ r ′ ) ∇ ′ ( e i k R R ) − ( e i k R R ) ∇ ′ ( e i k r ′ r ′ ) ] ⋅ n ^ d S ′ {\displaystyle \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'} ； 其中， n ^ {\displaystyle {\hat {\mathbf {n} }}} 是與 d S ′ {\displaystyle \mathrm {d} \mathbf {S} '} 同方向的單位向量。 注意到球面出射波的梯度為 ∇ ′ ( e i k R R ) = − ( e i k R R ) ( i k −   1 R ) R ^ {\displaystyle \nabla '\left({\frac {e^{ikR}}{R}}\right)=-\left({\frac {e^{ikR}}{R}}\right)\left(ik-\ {\frac {1}{R}}\right){\hat {\mathbf {R} }}} 、 ∇ ′ ( e i k r ′ r ′ ) = ( e i k r ′ r ′ ) ( i k −   1 r ′ ) r ′ ^ {\displaystyle \nabla '\left({\frac {e^{ikr'}}{r'}}\right)=\left({\frac {e^{ikr'}}{r'}}\right)\left(ik-\ {\frac {1}{r'}}\right){\hat {\mathbf {r} '}}} 。 從基爾霍夫所做的假定， k ≫ 1 / R {\displaystyle k\gg 1/R} 、 k ≫ 1 / r ′ {\displaystyle k\gg 1/r'} （例如，假設距離大約為1mm，則對於波長在0.4μm至0.7μm之間的可見光，可以做這假定；但對於波長在1mm至1m之間的微波，這假定不適用），則上述兩個公式近似為 ∇ ′ ( e i k R R ) = − i k ( e i k R R ) R ^ {\displaystyle \nabla '\left({\frac {e^{ikR}}{R}}\right)=-ik\left({\frac {e^{ikR}}{R}}\right){\hat {\mathbf {R} }}} 、 ∇ ′ ( e i k r ′ r ′ ) = i k ( e i k r ′ r ′ ) r ′ ^ {\displaystyle \nabla '\left({\frac {e^{ikr'}}{r'}}\right)=ik\left({\frac {e^{ikr'}}{r'}}\right){\hat {\mathbf {r} '}}} 。 所以，在點P的波擾 ψ ( r ) = −   i ψ 0 2 λ ∮ S ( e i k ( r ′ + R ) r ′ R ) ( r ′ ^ ⋅ n ^ + R ^ ⋅ n ^ ) d S ′ = −   i ψ 0 2 λ ∮ S ( e i k ( r ′ + R ) r ′ R ) ( cos ⁡ α + cos ⁡ χ ) d S ′ {\displaystyle {\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}}} ； 其中， α {\displaystyle \alpha } 、 χ {\displaystyle \chi } 分別是 r ′ ^ {\displaystyle {\hat {\mathbf {r} '}}} 、 R ^ {\displaystyle {\hat {\mathbf {R} }}} 與 n ^ {\displaystyle {\hat {\mathbf {n} }}} 之間的夾角。 這就是菲涅耳-基爾霍夫衍射公式，或基爾霍夫衍射公式。 如右圖所示，假設閉合曲面 S {\displaystyle \mathbb {S} } 是圓球面，點波源Q0與圓球面 S {\displaystyle \mathbb {S} } 的圓心同點。在圓球面 S {\displaystyle \mathbb {S} } 的任意位置， r ′ ^ {\displaystyle {\hat {\mathbf {r} '}}} 與 n ^ {\displaystyle {\hat {\mathbf {n} }}} 同向，所以， cos ⁡ α = 1 {\displaystyle \cos \alpha =1} 。 注意到 r ′ {\displaystyle r'} 是圓球面 S {\displaystyle \mathbb {S} } 的半徑，對於這積分， r ′ {\displaystyle r'} 值不變，可以從積分裏提出。在點P的波擾為 ψ ( r ) = −   i ψ ( r ′ ) λ ∮ S ( e i k R R ) K ( χ ) d S ′ {\displaystyle \psi (\mathbf {r} )=-\ {\frac {i\psi (\mathbf {r} ')}{\lambda }}\oint _{\mathbb {S} }\left({\frac {e^{ikR}}{R}}\right)K(\chi )\,\mathrm {d} S'} ； 其中， K ( χ ) = 1 + cos ⁡ χ 2 {\displaystyle K(\chi )={\frac {1+\cos \chi }{2}}} 為傾斜因子。 應用惠更斯－菲涅耳原理，所得到在點P的波擾的方程式，就是這方程式。但是，惠更斯－菲涅耳原理無法解釋相位差與傾斜因子的物理原因。傾斜因子使得次波的波幅會因為傳播方向而不同；朝著主波方向，波幅較大；逆著主波方向，波幅較小。這解釋了為甚麼波動只會朝著前方傳播的物理現象。 仔細詮釋惠更斯－菲涅耳原理的方程式：從點波源Q0發射的波幅為 ψ 0 {\displaystyle \psi _{0}} 的球面波，在點Q的波擾為 ψ ( r ′ ) = ψ 0 e i k r ′ / r ′ {\displaystyle \psi (\mathbf {r} ')=\psi _{0}e^{ikr'}/r'} ；而從點Q發射的次波，將傾斜因子與相位差納入考量，所貢獻出的波擾，在點P為 −   i ψ ( r ′ ) λ ( e i k R R ) K ( χ ) {\displaystyle -\ {\frac {i\psi (\mathbf {r} ')}{\lambda }}\left({\frac {e^{ikR}}{R}}\right)K(\chi )} 。 總合所有與點Q同波前的點次波源在點P所貢獻出的波擾，就可以得到 ψ ( r ) {\displaystyle \psi (\mathbf {r} )} 。 換另一種直接方法來詮釋，從點波源Q0發射的球面波，在點P的波擾為 ψ ( r ) = ψ 0 e i k r r {\displaystyle \psi (\mathbf {r} )=\psi _{0}{\frac {e^{ikr}}{r}}} 。 假若這兩種詮釋都正確，則從這兩種 ψ ( r ) {\displaystyle \psi (\mathbf {r} )} 的表達式分別計算出的結果，應該可以被核對為相等： ψ 0 e i k r r = −   i λ ∮ S 1 ( ψ 0 e i k r ′ r ′ ) ( e i k R R ) K ( χ ) d S ′ {\displaystyle \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 ≫ r ′ {\displaystyle r\gg r'} ，則以下近似成立： R ≈ r − r ′ cos ⁡ ( θ ) {\displaystyle R\approx r-r'\cos(\theta )} 、 R 2 ≈ r 2 − 2 r r ′ cos ⁡ ( θ ) {\displaystyle R^{2}\approx r^{2}-2rr'\cos(\theta )} 、 K ( χ ) = 1 + cos ⁡ ( χ ) 2 ≈ 1 + cos ⁡ ( θ ) 2 {\displaystyle K(\chi )={\frac {1+\cos(\chi )}{2}}\approx {\frac {1+\cos(\theta )}{2}}} ； 其中， θ {\displaystyle \theta } 為 r ′ {\displaystyle \mathbf {r} '} 與 r {\displaystyle \mathbf {r} } 之間的夾角。 所以，在點P的波擾可以近似為 ψ ( r ) ≈ −   i ψ 0 2 λ e i k ( r ′ + r ) r ′ r ∫ 0 π e − i k r ′ cos ⁡ ( θ ) [ 1 + cos ⁡ ( θ ) ] 2 π r ′ 2 sin ⁡ ( θ ) d θ ≈ −   i k ψ 0 r ′ e i k ( r ′ + r ) 2 r ∫ 0 π e − i k r ′ cos ⁡ ( θ ) [ 1 + cos ⁡ ( θ ) ] sin ⁡ ( θ ) d θ ≈ −   i k ψ 0 r ′ e i k ( r ′ + r ) 2 r   2 [ sin ⁡ ( k r ′ ) + i cos ⁡ ( k r ′ ) ] k r ′ ≈ ψ 0 e i k r r {\displaystyle {\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}}} 。 假設波源為有限尺寸，位於曲面 S {\displaystyle \mathbb {S} } 的波擾表達為 ψ ( r ′ ) {\displaystyle \psi (\mathbf {r} ')} ，則位於點P的波擾為 ψ ( r ) = 1 4 π ∮ S [ ψ ( r ′ ) ∂ ∂ n ′ ( e i k R R ) − ( e i k R R ) ∂ ψ ( r ′ ) ∂ n ′ ] d S ′ {\displaystyle \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 ≫ 1 / R {\displaystyle k\gg 1/R} ，則 ψ ( r ) = −   1 4 π ∮ S ( e i k R R ) [ i k ψ ( r ′ ) cos ⁡ ( R ^ , n ^ ) + ∂ ψ ( r ′ ) ∂ n ′ ] d S ′ {\displaystyle \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'} 。 這是基爾霍夫衍射公式最廣義的形式。解析涉及到有限尺寸波源的問題，必須用體積分來將波源的每一點所給出的貢獻總合在一起。 光波是傳播於空間的電磁輻射，理當被視為一種電磁場向量現象。但是，基爾霍夫的理論是純量理論，將光波當作純量處理，這可能會造成偏差。因此，物理學者做了很多實驗來檢查結果是否準確。他們發現，只要孔徑尺寸比波長大很多、孔徑與觀察屏之間的距離不很近，則使用純量理論可以得到相當準確的答案。但是對於某些問題，例如高解析度光柵衍射，純量理論就不適用，必須使用向量理論。 帕松光斑 基爾霍夫衍射公式 基爾霍夫積分定理 M. Born and E. Wolf, Principles of Optics, 1999, Cambridge University Press, Cambridge RS Longhurst, Gemoetrical and Physical Optics, 1969, Longmans, London Hecht, Eugene. Optics 4th. United States of America: Addison Wesley. 2002: pp. 510–512. ISBN 0-8053-8566-5 （英语）.  引文格式1维护：冗余文本 (link) G. Kirchhoff, Ann. d. Physik. 1883, 2, 18, p663 Goodman, Joseph. Introduction to Fourier Optics 3rd. Roberts and Company Publishers. 2004: pp. 35. ISBN 978-0974707723.  引文格式1维护：冗余文本 (link)