标量衍射理论

标量衍射理论

---

引言

人类对于光的本质的认知：

几何光学是在波长$\lambda \to 0$时的理论。这时光线的传播轨迹为一条直线，这是波动光学的一级近似。
几何光学三大基本定律

• 光的直线传播定律
• 光的独立传播定律
• 光的反射和折射定律

这些基本定律都可以通过费马原理来导出。

当波长不可忽略，这时在光空间中传播可以用惠更斯-菲涅尔原理来解释，即：任意时刻的光场是子波源产生球面波的相干叠加。光波的波前是子波源的包络面。惠更斯-菲涅尔原理可以成功解释光的反射和折射，但无法定量给出强度分布。并且惠更斯原理中的波是带有经典力学性质的机械波，与麦克斯为提出的电磁波理论不同。惠更斯假设是人类对光学现像的一个近似认识。

在历史的早期由于牛顿的权威地位，他提出的光的微粒假说一直受到推崇，对于光本性的探索停滞不前。直到1801杨氏进行了著名的双缝干涉实验，宣告光是一种波动现像。

1818年菲涅尔发展了惠更斯的思想，在杨氏实验证明光是一种波动的前提下，对于次级波源的相位和振幅进行了适当的假设，并且将光的纵波观点改成了横波观点。但菲涅尔理论中存在一个倾斜因子，无法解释其物理含义，理论存在缺陷。

1860年，麦克斯韦建立了其经典的电磁理论，并且预言了光是一种电磁波。而后基尔霍夫在更加坚实的数学基础上建立了比较完善的衍射理论模型。但该理论模型中的边界假设条件后来证明是不相容的，索末菲取消了其中之一的假设，利用格林函数修正了基尔霍夫理论，建立了瑞利-索末菲衍射理论。

这些理论都是线性近似理论，把光当成是标量，忽略了电磁场的适量性，电磁场通过麦克斯韦方程组相互作用。不能假定其相互无关，单独处理。但在满足一下两个条件是，标量衍射理论给出相当精确的结果。

1.衍射孔径比波长大的多
2.在太靠近衍射空的地方观察电磁场

当光在非自由空间传播（衍射受限）会发生明显的衍射现像，这是光波动性的必然结果。

不能用反射或折射解释的光学现像叫做衍射———索末菲

电磁理论

当介质中没有自由电荷存在时，麦克斯韦方程组：
$$\nabla \times E=-\mu \frac{\partial H}{\partial t}\text{\hspace{1em}(1)}$$

$$\nabla \times H=\epsilon \frac{\partial E}{\partial t}\text{\hspace{1em}(2)}$$

$$\nabla \cdot \epsilon E=0\text{\hspace{1em}(3)}$$

$$\nabla \cdot \mu H=0\text{\hspace{1em}(4)}$$
根据矢量计量理论有
$$\nabla \times (\nabla \times E)=\nabla (\nabla \cdot E)-{\nabla }^{2}E$$
带入(1)和(2)得到
$${\nabla }^{2}E-\frac{n^2}{c^2}\frac{{\partial }^{2}E}{\partial t^2}=0$$
其中n表示介质的折射率，定义为
$$n=(\epsilon /{\epsilon }_0{)}^{1/2}$$
${\epsilon }_0$为真空介电常数，$c$表示真空中的光速
$$c=\frac{1}{\sqrt{{\epsilon }_0{\mu }_0}}$$
利用已知的电磁常量${\epsilon }_0=8.8542\times 10^{-12}{C}^2/(N\cdot m^2)$，${\mu }_0=4\pi \times 10^{-7}N\cdot S^2/C^2$，即可算出

$$c=2.99794\times 10^{8}m/s$$

与实验测得的光速相吻合，说明了光是一种电磁波。

同理，得到磁场满足的波动方程
$${\nabla }^{2}H-\frac{n^2}{c^2}\frac{{\partial }^{2}H}{\partial t^2}=0$$
上述结论需假定，介质具有线性性，各向同性，均匀，无色散。

如果介质是非均匀，即$\epsilon$与空间位置有关，这时的波动方程为
$${\nabla }^{2}E+2\nabla (E\cdot \nabla \ln n)-\frac{n^2}{c^2}\frac{{\partial }^{2}E}{\partial t^2}=0$$
空间中任意一点$Q$的惠更斯-菲涅尔原理的数学表达式
$$U(Q)=c{\iint }_{\Sigma }{U}_0(P)K(\theta )\frac{e^{ikr}}{r}\mathrm{d}S$$
其中，$\Sigma$是波面，$U_0(P)$是波面在$P$的振幅，$r$是点$P$和$Q$点的距离，$\theta$是$PQ$连线与波面法线的夹角。$K(\theta )$是倾斜因子，表示与角度有关的函数。$c$表示一个常量。
对于单色光场的任意一点$Q$，满足波动方程
$${\nabla }^{2}u-\frac{1}{c^2}\frac{{\partial }^{2}u}{\partial t^2}=0\text{\hspace{1em}(5)}$$
将复振动表示为
$$u(Q,t)=U(Q)exp(-j2\pi \nu t)$$
带入到波动方程(5)中得到
$$({\nabla }^2+k^2)U(Q)=0\text{\hspace{1em}(6)}$$
该方程称作亥姆霍兹方程，亥姆霍兹方程与时间无关，表示单色场所满足的空间部分关系。

基尔霍夫衍射理论

基尔霍夫公式的推导

格林定理

如果$U$和$P$是任意位置的复值函数，并且具有连续单值的二阶偏导。则有
$${\iiint }_V(U{\nabla }^{2}G-G{\nabla }^{2}U)\mathrm{d}v={\iint }_S(U\frac{\partial G}{\partial n}-G\frac{\partial U}{\partial n})\mathrm{d}s$$
设观察点为$P_0$,在格林公式中，选择自由空间格林函数作为辅助函数
$$G(P_1)=\frac{jk{r}_{01}}{r_{01}}$$
式中的$r_{01}$代表着从点$P_0$到点$P_1$的距离。
为了排除奇异点，在点$P_0$的周围用足够小的面将其包裹，这样和原来的面$S$构成一个复合平面
$$S^{{}^{\prime }}=S+S_{ϵ}$$
在面所包围的体积内，格林函数满足亥姆霍兹方程
$$({\nabla }^2+k^2)G=0$$
代入到格林公式中有
$${\iint }_{S^{{}^{\prime }}}(U\frac{\partial G}{\partial n}-G\frac{\partial U}{\partial n})\mathrm{d}s$$
在面$S^{{}^{\prime }}$上的一般点
$$G(P_1)=\frac{exp(jk{r}_{01})}{r_{01}}$$
以及
$$\frac{\partial G(P_1)}{\partial n}=\cos (n,r_{01})(jk-\frac{1}{r_{01}})\frac{exp(jk{r}_{01})}{r_{01}}$$
最后让$ϵ$任意小，并利用在$P_0$的连续性得到
U ( P 0 ) = 1 4 π ∬ S { ∂ U ∂ n [ e x p ( j k r 01 ) r 01 ] − U ∂ ∂ n [ e x p ( j k r 01 ) r 01 ] } d s U(P_0)=\frac{1}{4\pi}\iint_S \Big\{\frac{\partial U}{\partial n}\left[\frac{exp(jkr_{01})}{r_{01}}\right]-U\frac{\partial}{\partial n}\left[\frac{exp(jkr_{01})}{r_{01}}\right]\Big\}\mathrm{d}s U(P0​)=4π1</