菲涅耳衍射

在光學裏，菲涅耳衍射（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{align}R & =\sqrt{(x-x')^2+(y-y')^2+z^2}=\sqrt{\rho^2+z^2} \\ <pre> & = z \sqrt{ 1 + \frac{\rho^2}{z^2} } \\ & = z \left[ 1 + \frac{\rho^2}{2 z^2} - \frac{1}{8} \left( \frac{\rho^2}{z^2} \right)^2 + \cdots \right] \\ & = z + \frac{\rho^2}{2 z} - \frac{\rho^4}{8z^3} + \cdots \\ </pre> \end{align}$。

假若保留所有項目，則這級數式為精確解。

$e^{i k r}/r$

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

$\frac{ k \rho^4}{8 z^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}{2 z} = z + \frac{(x-x')^2 +(y-y')^2}{2 z}$。

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

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

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

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

### 菲涅耳衍射積分式

菲涅耳數 $F=a^2/L\lambda$ 菲涅耳衍射區域：$F \ge 1$ 夫朗和斐繞射區域：$F \ll 1$ $a$ － 孔徑或狹縫的尺寸 $\lambda$ － 波長 $L$ － 離開孔徑或狹縫的距離

假設孔徑尺寸超小於傳播路徑長度，則 $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{align}\psi(0,0,z) & =-\ \frac{ie^{ikz} \psi_0 }{\lambda z} \int_0^a e^{ik\rho'^2/2z}\ \rho'\mathrm{d}\rho' \\ <pre> & =-\psi_0 e^{ikz} (e^{ika^2/2z}-1) \\ </pre> \end{align}$。

輻照度 $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\cdot10^{-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}{2 z} (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 -2 x x'$。

$(y-y')^2 = y^2 + y'^2 -2 y y'$。

則菲涅耳衍射積分式可以以二維傅立葉變換來表達。設定函數 $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{align}\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}\} \\ <pre> & =-\ \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) \\ </pre> \end{align}$。

其中，函數 $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）
- 基爾霍夫積分定理
- 基爾霍夫衍射公式