菲涅尔衍射

在光學裏，**菲涅耳繞射**（Fresnel diffraction），又稱**近場繞射**（Near-field diffraction），指的是光波在近場區域的繞射。菲涅耳繞射積分式可以用來計算光波在近場區域的傳播，因法國物理學者奧古斯丁・菲涅耳而命名，是克希荷夫繞射公式的近似。

從每一個光學系統特徵的菲涅耳數，可以辨別光波傳播的區域是近場還是遠場。設想光波入射於任意孔徑，對於這光學系統，菲涅耳數定義為

　
: $F\ {\stackrel {def}{=}}\ {\frac {a^{2}}{L\lambda }}$ ；

其中， $a$ 是孔徑的尺寸， $L$ 是孔徑與觀察屏之間的距離， $\lambda $ 是入射波的波長。

假若 $F\gtrsim 1$ ，則繞射波是處於近場，可以使用菲涅耳繞射積分式來計算其物理性質。

## 菲涅耳繞射

假設照射光波於開有孔徑的不透明擋板，則會有繞射圖樣出現於觀察屏。根據惠更斯 - 菲涅耳原理，從孔徑內部任意點次波源 Q 發射出的圓球面次波，在觀察屏點 P 的波擾 $\psi (x,y,z)$ 為

　
: $\psi (x,y,z)=-\ {\frac {i}{\lambda }}\int _{\mathbb {S} }\psi (x',y',0){\frac {e^{ikR}}{R}}K(\chi )\ \mathrm {d} x'\mathrm {d} y'$ ；

其中， $\mathbf {r} =(x,y,z)$ 是點 P 的直角坐標， $\mathbf {r} '=(x',y',0)$ 是點 Q 的直角坐標， $\lambda $ 是波長， $\mathbb {S} $ 是積分平面（孔徑）， $\psi (x',y',0)$ 是位於點次波源 Q 的波擾， $\mathbf {R} $ 是從點 Q 到點 P 的位移向量， $R$ 是 $\mathbf {R} $ 的數值大小， $K(\chi )$ 是傾斜因子， $\chi $ 是垂直於孔徑平面的法向量與 $\mathbf {R} $ 之間的夾角。

古斯塔夫・克希荷夫給出了傾斜因子 $K(\chi )$ 的表達式：

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

除了最簡單的繞射案例以外，幾乎不可能找到這積分式的解析解。通常，必須使用數值分析方法來解析這積分式。

### 菲涅耳近似

為了要計算這積分式的解答，必須先使積分項目更簡單化。設定

　
: $\rho ={\sqrt {(x-x')^{2}+(y-y')^{2}}}$ 。

 $(x',y',0)$ 與 $(x,y,z)$ 之間的距離 $R$ 可以以泰勒級數表示為

　
: ${\begin{aligned}R&={\sqrt {(x-x')^{2}+(y-y')^{2}+z^{2}}}={\sqrt {\rho ^{2}+z^{2}}}\\&=z{\sqrt {1+{\frac {\rho ^{2}}{z^{2}}}}}\\&=z\left[1+{\frac {\rho ^{2}}{2z^{2}}}-{\frac {1}{8}}\left({\frac {\rho ^{2}}{z^{2}}}\right)^{2}+\cdots \right]\\&=z+{\frac {\rho ^{2}}{2z}}-{\frac {\rho ^{4}}{8z^{3}}}+\cdots \\\end{aligned}}$ 。

假若保留所有項目，則這級數式為精確解。將這 $R$ 的級數式代入被積函數的相位。菲涅耳近似的要點是在假定級數式的第三個項目非常微小，可以被忽略。為了達到這目的，第三個項目必須超小於相位的週期 $2\pi $ ：

　
: ${\frac {k\rho ^{4}}{8z^{3}}}\ll 2\pi $ 。

改換以波長 $\lambda =2\pi /k$ 來表達，

　
: ${\frac {\rho ^{4}}{8z^{3}\lambda }}\ll 1$ 。

將先前 $\rho $ 的表達式代入，

　
: ${\frac {[(x-x')^{2}+(y-y')^{2}]^{2}}{8z^{3}\lambda }}\ll 1$ 。

假若，對於所有 $(x',y',0)$ 與 $(x,y,z)$ 的可能值，這條件成立，則泰勒級數式的第三個項目和更高階項目都可以忽略。

從這些論述， $R$ 可以近似為

　
: $R\approx z+{\frac {\rho ^{2}}{2z}}=z+{\frac {(x-x')^{2}+(y-y')^{2}}{2z}}$ 。

這方程式稱為「菲涅耳近似」。這近似成立的條件是上述不等式。

例如，對於半徑為 1mm 的圓孔，假設觀察屏區域的半徑也是 1mm，入射波的波長為 500nm，則近似成立的條件為

　
: $z\gg \left({\frac {\rho ^{4}}{8\lambda }}\right)^{1/3}=\left[{\frac {0.002^{4}}{8\cdot 500\cdot 10^{-9}}}\right]^{1/3}\approx 0.016[m]$ 。

圓孔與觀察屏之間的距離 $z$ 必須超大於 16mm。實際而言，這條件太過嚴苛，從數值分析的結果，只要圓孔與觀察屏之間的距離 $z$ 大於 16mm 就行了。

### 菲涅耳繞射積分式

假設孔徑尺寸超小於傳播路徑長度，則 $K(\chi )\approx 1$ 。特別是在 z - 軸附近的小範圍區域， $x,y\ll z$ ，分母的 $R$ ，可以近似為 $R\approx z$ ，只取至線性項目。現在，採用菲涅耳近似，則在位置 $(x,y,z)$ 的波擾為

　
: $\psi (x,y,z)=-\ {\frac {ie^{ikz}}{\lambda z}}\int _{\mathbb {S} }\psi (x',y',0)e^{ik[(x-x')^{2}+(y-y')^{2}]/2z}\ \mathrm {d} x'\mathrm {d} y'$ 。

這就是「菲涅耳繞射積分式」。仔細推敲這積分式的含意，假設菲涅耳近似成立，則位於孔徑的次波源發射出的圓球面次波，會沿著 z - 軸方向，傳播到觀察屏。整個積分調製圓球面波的波幅與相位。只有對於少數案例，這方程式存在分析解答。

更進一步近似，將 $e^{ikR}$ 近似為 $e^{ikz}$ ，相位部分僅取至線性項目，這只有當觀察屏與孔徑之間的距離超遠時才成立，請參閱條目夫朗和斐繞射。菲涅耳繞射與夫朗和斐繞射不同的地方，主要是菲涅耳繞射將波前的曲率納入考量，這是為了要精確計算相互干涉的波擾彼此之間的相對相位。

### 圓孔繞射

**軸上照度**

假設孔徑是半徑為 $a$ 的圓孔，首先計算沿著中心軸的波擾， $x,y=0$ ，又假設入射波是波幅為 $\psi _{0}$ 、朝著 z - 軸傳播的平面波。

根據菲涅耳繞射積分式，

　
: $\psi (0,0,z)=-\ {\frac {ie^{ikz}\psi _{0}}{\lambda z}}\int _{\mathbb {S} }e^{ik(x'^{2}+y'^{2})/2z}\ \mathrm {d} x'\mathrm {d} y'$ 。

改採極坐標 $(\rho ',\theta ')$ ，

　
: ${\begin{aligned}\psi (0,0,z)&=-\ {\frac {ie^{ikz}\psi _{0}}{\lambda z}}\int _{0}^{a}e^{ik\rho '^{2}/2z}\ \rho '\mathrm {d} \rho '\\&=-\psi _{0}e^{ikz}(e^{ika^{2}/2z}-1)\\\end{aligned}}$ 。

輻照度 $I(z)$ 為

　
: $I(z)=\psi ^{*}\psi /2=\psi _{0}^{\ 2}\ 2\sin ^{2}(ka^{2}/4z)=I_{0}\sin ^{2}(ka^{2}/4z)$ 。

從這函數繪製的輻照度比率對無量綱距離圖展示出，離孔徑越近，震盪越劇烈。這區域是菲涅耳繞射區域。在這區域裏，輻照度的極值點分別為

* 極大值：當 $z={\frac {a^{2}}{2n\lambda }},\qquad n=1,2,3,\dots $ 。
* 極小值：當 $z={\frac {a^{2}}{(2n-1)\lambda }},\qquad n=1,2,3,\dots $ 。

離孔徑越遠，兩個相鄰極值點之間的間隔越大， $z=a^{2}/\lambda $ 是最後一個極值點。遠於這距離，輻照度呈單調遞減。通常，物理學者規定菲涅耳繞射區域的菲涅耳數大於或等於 1，這對應於 $Z_{F}=a^{2}/\lambda $ 為分界點；超遠於這分界點，是夫朗和斐繞射區域，可以使用夫朗和斐近似，數學計算比較簡單很多。

例如，對於半徑為 1mm 的圓孔，假設入射波的波長為 500nm，則 $Z_{F}$ 為

　
: $Z_{F}={\frac {0.001^{2}}{500\cdot 10^{-9}}}\approx 2[m]$ 。

總結，孔徑與觀察屏之間的距離在 2m 以內是菲涅耳繞射區域，以外是夫朗和斐繞射區域。

**軸側照度**

 
$$
I=\left(V_{0}-\cos \left({\frac {u^{2}+v^{2}}{2u}}\right)\right)^{2}+\left(V_{1}-\sin \left({\frac {u^{2}+v^{2}}{2u}}\right)\right)^{2}
$$

其中 $V_{m}$ 是二元隆梅爾函數 (Lommel function)

$$
V_{m}=\sum _{n=0}^{\infty }*((-1)^{n}*({\frac {v}{u}})^{2*n+m}*J_{2n+m}(v))
$$

 $J_{2n+m}(v)$ 為 第一類 $2n+m$ 階貝索函數

### 圓盤繞射

圓盤繞射在軸上的強度為

$$
I=I_{0}*\lambda ^{2}/4
$$

因此圓盤繞射的軸上強度，和波長的平方成正比，而與圓盤的直徑、與圓盤的距離無關，所以繞射圖形的中心一定是個亮點。這個亮點稱為帕松光斑。菲涅耳圓孔繞射圖形的中心點，根據圓孔直徑和距離之不同，可以是亮點，也可以是暗點。

### 單縫繞射

菲涅耳單縫繞射的強度分布為：.

$$
I=(Cp(Y)-Cq(Y))^{2}+(Sp(Y)-Sq(Y))^{2}
$$

其中 Cp,Cq 為餘弦菲涅耳積分：

$$
Cp(Y):=\int _{0}^{p}(\cos((1/2)*\pi *t^{2})\,dt
$$

 $Cq(Y)=\int _{0}^{q}(\cos((1/2)*\pi *t^{2})\,dt$ ;

 
Sp,Sq 為正弦菲涅耳積分：

$$
Sp(Y):=\int _{0}^{p}(\sin((1/2)*\pi *t^{2})\,dt
$$

 $Sq(Y)=\int _{0}^{q}(\sin((1/2)*\pi *t^{2})\,dt$ ;

菲涅耳單縫繞射圖形與夫琅禾費單縫繞射明顯不同之處在於前者的第一個極小值不等於 0（如圖），而後者為 0。

### 直邊繞射

一個平面波通過與光線傳播方向垂直的不透明直邊，

 
其菲涅耳直邊繞射的強度分布為：.

$$
I=(Cp(Y)+0.5)^{2}+(Sp(Y)+0.5))^{2}
$$

其中 Cq 為餘弦菲涅耳積分：

 
 $Cq(Y)=\int _{0}^{q}(\cos((1/2)*\pi *t^{2})\,dt$ ;

Sp 為正弦菲涅耳積分：

$$
Sp(Y):=\int _{0}^{p}(\sin((1/2)*\pi *t^{2})\,dt
$$

 
菲涅爾直邊繞射圖樣在幾何陰影區附近強度不為零，在明亮區間，光強以阻尼振動形式逐漸衰減至一個穩定數值

## 進階理論

### 卷積

設定函數 $h(x,y,z)$ 為

　
: $h(x,y,z)=-\ {\frac {ie^{ikz}}{\lambda z}}e^{i{\frac {k}{2z}}(x^{2}+y^{2})}$ 。

波擾 $\psi (x,y,z)$ 可以寫為卷積形式：

　
: $\psi (x,y,z)=\iint \limits _{-\infty }^{\ \ \ \infty }\psi (x',y',0)h(x-x',y-y',z)\ \mathrm {d} x'\mathrm {d} y'$ ，

或者，由於 z - 坐標與積分無關，可以將 z - 坐標資訊提出，

　
: $\psi _{z}(x,y)=\iint \limits _{-\infty }^{\ \ \ \infty }\psi _{0}(x',y')h_{z}(x-x',y-y')\ \mathrm {d} x'\mathrm {d} y'$ ，

這卷積又可以標記為

　
: $\psi _{z}(x,y)=\psi _{0}(x,y)*h_{z}(x,y)$ 。

根據卷積定理，函數卷積的傅立葉變換是函數傅立葉變換的乘積，以方程式表達，

　
: ${\mathcal {F}}\{\psi _{z}(x,y)\}={\mathcal {F}}\{\psi _{0}(x,y)*h(x,y)\}={\mathcal {F}}\{\psi _{0}(x,y)\}\cdot {\mathcal {F}}\{h_{z}(x,y)\}$ ；

其中， ${\mathcal {F}}\{f(x,y)\}$ 是函數 $f(x,y)$ 的傅立葉變換。

假設這光學系統是線性系統，滿足空間不變性，即改變波源的位置只會改變繞射圖樣的位置，不會改變繞射圖樣的形狀。這樣，一個有限尺寸波源所產生的繞射圖樣，可以視為是由其每一個點波源所產生的繞射圖樣共同線性疊加而形成。

假設這線性系統的線性算子為 ${\mathcal {L}}$ 、輸入函數為 $f(x,y)$ 、輸出函數為 $G(X,Y)$ ，則這兩個函數之間的關係可以表達為

　
: $G(X,Y)={\mathcal {L}}\{f(x,y)\}$ 。

應用狄拉克 δ 函數的數學性質，

　
: $G(X,Y)={\mathcal {L}}\left\{\iint \limits _{-\infty }^{\ \ \ \infty }f(x',y')\delta (x-x')\delta (y-y')\ \mathrm {d} x'\mathrm {d} y'\right\}$ 。

將 $f(x',y')$ 視為函數 $\delta (x-x')\delta (y-y')$ 權重係數，應用線性系統的性質，可以將積分式寫為

　
: $G(X,Y)=\iint \limits _{-\infty }^{\ \ \ \infty }f(x',y'){\mathcal {L}}\{\delta (x-x')\delta (y-y')\}\ \mathrm {d} x'\mathrm {d} y'$ 。

由此推論，表現觀察屏輻照圖案的函數 $h_{z}(x,y)$ 是線性系統對於在孔徑位置 $(x',y')$ 的狄拉克 δ 函數所做出的響應，因此稱為脈衝響應。

### 傅立葉變換

定義空間頻率 $K_{x},K_{y}$ 為

　
: $K_{x}\ {\stackrel {def}{=}}\ kx/z$ 、

　
: $K_{y}\ {\stackrel {def}{=}}\ ky/z$ 。

將橫向位移的每一個分量展開，

　
: $(x-x')^{2}=x^{2}+x'^{2}-2xx'$ 、

　
: $(y-y')^{2}=y^{2}+y'^{2}-2yy'$ ，

則菲涅耳繞射積分式可以以二維傅立葉變換來表達。設定函數 $G(K_{x},K_{y})$ 為函數 $g(x',y')$ 的傅立葉變換，那麼，根據定義，函數 $G(K_{x},K_{y})$ 為

　
: $G(K_{x},K_{y})\ {\stackrel {def}{=}}\ {\mathcal {F}}\left\{g(x',y')\right\}\ {\stackrel {def}{=}}\ \iint \limits _{-\infty }^{\ \ \ \infty }g(x',y')e^{-i(K_{x}x'+K_{y}y')}\ \mathrm {d} x'\mathrm {d} y'$ ；

設定函數 $g(x',y')$ 為

　
: $g(x',y')=\psi _{0}(x',y')e^{ik(x'^{2}+y'^{2})/2z}$ 。

菲涅耳繞射積分式表達為

　
: ${\begin{aligned}\psi _{z}(x,y)&=-\ {\frac {ie^{ikz}}{\lambda z}}e^{ik(x^{2}+y^{2})/2z}\ {\mathcal {F}}\{\psi _{0}(x',y')e^{ik(x'^{2}+y'^{2})/2z}\}\\&=-\ {\frac {ie^{ikz}}{\lambda z}}e^{ik(x^{2}+y^{2})/2z}\ {\mathcal {F}}\{g(x',y')\}\\&=-\ {\frac {ie^{ikz}}{\lambda z}}e^{ik(x^{2}+y^{2})/2z}\ G(K_{x},K_{y})\\&=h_{z}(x,y)\ G(K_{x},K_{y})\\\end{aligned}}$ ；

其中，函數 $h_{z}(x,y)=-\ {\frac {ie^{ikz}}{\lambda z}}e^{ik(x^{2}+y^{2})/2z}$ 。

在做實例計算時，先計算 $g(x',y')$ 的傅立葉變換，然後將空間頻率 $K_{x},K_{y}$ 替換回原來的變量，最後再乘以 $h_{z}(x,y)$ ，就可以得到 $\psi _{z}(x,y)$ 。假若 $g(x',y')$ 是個常見函數，而且已知道 $g(x',y')$ 的傅立葉變換，則這是一種比較精緻的理論方法。更多相關內容，請參閱條目線性標準轉換。

## 原子反射鏡

由於菲涅耳繞射機制，原子波入射於由物質形成的相互平行凸脊陣列，會被鏡面反射。這效應可以用來實現原子反射鏡。

## 參閱

* 菲涅耳積分
* 菲涅耳波帶（Fresnel zone）
* 克希荷夫積分定理
* 克希荷夫繞射公式

## 參考文獻

* Goodman, Joseph W. Introduction to Fourier optics. New York: McGraw-Hill. 1996. ISBN 0-07-024254-2.