衍射导论

当光波遇到障碍物的时候，或多或少偏离几何光学的直线传播，这种现象统称为光的衍射。光的波长和孔径线度相当时，衍射现象尤为明显，且关系满足：

$\rho \Delta \theta\approx\lambda$

衍射现象可以使用惠更斯原理来解释，菲涅尔推导了衍射的数学表达式，称为惠更斯-菲涅尔原理：

$\tilde{U}(P)=\oiint d\tilde{U}(P)$

菲涅尔衍射积分

从物理的角度考虑，可以给出$d\tilde{U}(P)$的决定因素：

面元大小$dS$；
次波源的振幅$\tilde{U_0}$；
次波源发射球面波到场点$\dfrac{e^{ikr}}{r}$；
倾角因子$f(\theta_0,\theta)=\frac12(\cos{\theta_0}+\cos{\theta})$。

引入比例常数$K$，可以得到：

$\tilde{U}(P)=K\oiint \tilde{U_0} \frac{e^{ikr}}{r} f(\theta_0,\theta)dS$

此即为菲涅尔衍射积分表达式。

基尔霍夫衍射积分

基尔霍夫在$kr\gg 1$和$r\gg \lambda$的情况下，导出了无源空间边值定解的表达式：

$\tilde{U}(P)=\frac{-i}{\lambda}\oiint_{\Sigma} \frac12(\cos{\theta_0}+\cos{\theta})\tilde{U_0} \frac{e^{ikr}}{r}dS$

表面上看，这只是明确了比例系数和倾斜因子；实际上，基尔霍夫衍射积分允许任意隔离光源和场点的闭合面，而不要求等相面。

基尔霍夫边界条件与傍轴衍射积分公式

取衍射屏面为$\Sigma_1$，无穷远半球面为$\Sigma_2$，衍射孔面为$\Sigma_0$，则有：

$\Sigma=\Sigma_1+\Sigma_2+\Sigma_0$

将基尔霍夫衍射积分公式代入，得到：

$\tilde{U}(P)=\oiint_{\Sigma_1} \cdots+\oiint_{\Sigma_2} \cdots+\oiint_{\Sigma_0} \cdots$

可以做基尔霍夫边界条件的假设，$\Sigma_1$和$\Sigma_2$的积分可以忽略，得到：

$\tilde{U}(P)=\frac{-i}{\lambda}\oiint_{\Sigma_0} \frac12(\cos{\theta_0}+\cos{\theta})\tilde{U_0} \frac{e^{ikr}}{r}dS$

在傍轴条件下，还可以进一步简化：

$\tilde{U}(P)=\frac{-i}{\lambda r_0}\oiint_{\Sigma_0} \tilde{U_0} e^{ikr}dS$

衍射的分类

衍射现象可以分为两类：

菲涅尔衍射：近场衍射，光源和屏幕距离较近；
夫琅禾费衍射：远场衍射，光源和屏幕距离无限远（或者衍射孔线度无线小）。

巴比涅原理

如果$\Sigma_a$和$\Sigma_b$是两个不同的互补衍射屏面，且它们的衍射场分别为$\tilde{U_a}$和$\tilde{U_b}$，则有：

$\tilde{U_a}+\tilde{U_b}=\tilde{U_0}$

对于夫琅禾费衍射尤其方便，因为光线聚集在后焦点，轴外的自由光场为零，那么：

$\tilde{U_a}=-\tilde{U_b},I_a(P)=I_b(P)$

这说明两者的衍射图样除了像点完全相同。

菲涅尔圆孔衍射

半波带法

处理圆孔衍射很自然的一个想法是将衍射孔的波前分割为一个个波带，每个波带的相位差为$\lambda/2$，这样，只需要确定最后一个波带的相位，就能定性确定接受屏的强度分布。

不过，我们需要考证每个波带的振幅是否相同。写出每个波带的振幅：

$A_k\propto f(\theta_k)\frac{dS_k}{r_k}$

alt

通过简单的数学运算，可以得到：

$\begin{aligned} dS_k&=d[2\pi R^2(1-\cos{\alpha})]\\ &=2\pi R^2\sin{\alpha}d\alpha\\ &=2\pi R^2\frac{rdr}{R(R+b)}\\ &\propto r \end{aligned}$ $A_k=\frac{f(\theta_k)}{R(R+b)}$

这说明每个半波带的振幅和$r$是无关的，只和$f(\theta_k(r_k))$有关。我们可以写出这一泛函的表达式：

$\begin{aligned} f(\theta_k(r_k))&=\frac12(1+\cos{\theta_k})\\ &=\frac12(1+\frac{(R+b)^2-R^2-(b+\frac{k\lambda}{2})}{2R(b+\frac{k\lambda}{2})})\\ &\approx \frac12[1-k\lambda(\frac1R+\frac1b)] \end{aligned}$

代入$\lambda\sim 100nm,b\sim 1m,R\sim 1m$可知，振幅随着波带数的变化及其缓慢。

做题的时候，只需要知道以下公式：

最大波带半径：$r_{max}=\sqrt{b^2+R^2\sin^2\alpha}$
最大波带数：$k_{max}=\frac{r_{max}}{\lambda}$
第k级波带半径：$\rho_k=\sqrt{k}\rho_1,\rho_1=\sqrt{\frac{Rb\lambda}{R+b}}$

以$\rho$为变量，最大半径和最小半径之间的差值为：
$\begin{aligned}\Delta r&=\sqrt{b^2+\rho^2}-(a+b-\sqrt{a^2+\rho^2})\\&\approx \frac{\rho^2}{2a}+\frac{\rho^2}{2b}\\&\approx \frac{\rho^2}{2R}+\frac{\rho^2}{2b}\end{aligned}$
所以：
$\rho=\sqrt{\frac{2Rb\Delta r}{R+b}}$

利用以上特性，可以绘制菲涅尔波带片，在奇数或偶数的波带上涂黑，即可增强光强。菲涅尔波带片就像透镜，可以将入射光汇聚起来，可以将上式改写为：

$\frac1R+\frac1b=\frac{k\lambda}{\rho_k^2}$

等效焦距为$f=\frac{\rho_k^2}{k\lambda}=\frac{\rho_1^2}{\lambda}$。

单缝夫琅禾费衍射

利用傍轴的基尔霍夫衍射公式：

$\tilde{U}(P)=\frac{-i}{\lambda r_0}\oiint_{\Sigma} \tilde{U_0}e^{ikr}dS$

假定缝宽为$b$，长度为$a$，坐标原点沿$\theta$方向到场点P的距离为$L_0$。那么可以计算出：

$\begin{aligned} \tilde{U}(P)&=\frac{-i}{\lambda f}Abe^{ikL_0}\int_{-\frac{a}{2}}^{\frac{a}{2}} e^{-ik\sin\theta x}dx\\ &=\frac{-i}{\lambda f}Abe^{ikL_0}\times a\text{sinc}{(\frac{\pi a \sin \theta }{\lambda})} \end{aligned}$

光强分布为：

$I=I_0(\frac{\sin \alpha}{\alpha})^2,\alpha=\frac{\pi a \sin \theta }{\lambda}$