Influence of partial structure error of rectangular grating on diffraction efficiency

高光谱遥感是当前遥感技术的前沿领域，集成像和光谱分析技术于一体，在分析了解地物信息方面具有十分重要的作用。衍射光栅作为一种重要的分光器件，在光谱分析中应用广泛。为了提高高光谱遥感技术的分析精度，需要了解光栅更多级次衍射光的特性，对光栅的加工精度也提出了更高的要求。目前制作衍射光栅的方法主要有3种：全息法、复制法和机械刻划法^[1-3]。但无论哪种制作技术，在实际应用过程中因为各种因素：如设备精度、人为操控等问题都会产生一定的误差^[4-6]。衍射效率是衍射光栅一个十分重要的性能指标^[7]。因此，分析加工误差对光栅衍射效率的影响，对于改进衍射光栅的制作技术具有重要的启示作用。

POMMET等人研究了制作误差对矩形光栅衍射效率的影响^[8]。CAO等人研究了制作误差对亚波长介质光栅衍射效率的影响^[9]。GE等人研究了制作误差对锯齿槽型光栅衍射效率的影响^[10], 他们着重于分析光栅整体周期、占空比、槽深和面型等误差对光栅衍射效率的影响，而本文中则侧重于利用标量衍射理论来分析光栅局部周期或缝宽误差对整个光栅各级衍射效率的影响。

光栅的衍射效率^[11]是指第m级的衍射光能量与入射光总能量之比。计算光栅衍射效率的方法一般分为两大类：分别是标量衍射理论^[12]和矢量衍射理论^[13]。在标量衍射理论中，把光波当作标量来处理，只考虑电场或磁场的某一个横向分量，而假设其它有关分量也可以用同样方法来独立处理，完全忽略了电磁场矢量分量之间的耦合特性。基尔霍夫衍射理论^[14]和瑞利-索末菲理论^[15]都属于标量衍射理论。光栅的标量衍射理论是计算光栅衍射效率的一种近似方法，相较于矢量衍射理论其具有数学计算简单、物理概念清晰等特点^[16-19]。下面将利用标量衍射理论对矩形槽光栅衍射光强度分布的一般公式进行推导。

矩形槽衍射光栅工作原理如图 1所示。图中，a为光栅缝宽，d为光栅周期(光栅常数)。本文中将分别讨论因加工造成光栅局部缝宽a和周期d出现的误差对光栅各级衍射效率产生的影响。

Figure 1.  Schematic diagram of rectangular grating diffraction

首先，对光栅衍射强度分布的一般公式进行推导。假设一束单位振幅的单色平行光垂直入射矩形光栅，用积分法求得通过该第n个单缝的衍射光在焦平面上的合成振动公式^[20]为:

式中，y为通过第n个单缝的衍射光在焦平面上的合成振幅；E₀为光波振幅；a_n为第n条狭缝宽度(n=0, 1, 2, …, N-1)；θ为衍射角；λ为单色入射光波长。

令${\beta _n} = {a_n}\frac{{\sin \left({\frac{{{\rm{ \mathsf{ π} }}{a_n}\sin \theta }}{\lambda }} \right)}}{{\frac{{{\rm{ \mathsf{ π} }}{a_n}\sin \theta }}{\lambda }}}$，则有：

式中，由d₀sinθ=λ可得θ=arcsin(λ/d₀)；d₀为起始栅距, β_n为第n条狭缝衍射光合成振幅系数。

通过N条光缝的N束衍射平行光，其相位均不相同，则通过第n条光缝的衍射光，其光波的向量表达式为：

式中，E_n为通过第n条光缝的光波振幅, ω为圆频率, t为时间, r为初相位, k为波数, δ_n为第n条狭缝衍射光相位角, d_n为第n个栅距。

N束相干光合成的向量表达式为：

则(5)式相对应的共轭波表达式为：

则可得N束相干光合成光强度I, 表示为：

式中, I₀=E₀²为入射光强度, F为衍射-干涉复合因子。(7)式为光栅衍射强度分布的一般公式。

下面将分别对以上加工误差进行理论分析，这里主要讨论了两种情况。

(1) 改变光栅的局部周期d，其它周期d保持相同不变，以此来分析这种加工的周期变化对矩形光栅衍射效率能产生的影响。

由(7)式可知，当光栅狭缝不变时，即a为常数，则有：

式中, $\beta = a\frac{{\sin \left({\frac{{{\rm{ \mathsf{ π} }}a\sin \theta }}{\lambda }} \right)}}{{\frac{{{\rm{ \mathsf{ π} }}a\sin \theta }}{\lambda }}}$。

(2) 改变矩形光栅的部分缝宽，但各周期保持不变(即加工周期不存在误差)，以此来分析这种加工光栅时产生的缝宽误差对矩形光栅衍射效率的影响。

由(7)式可知，当光栅周期d(光栅常数)不变时，即d为常数，则有：

当d为常数时，由(4)式可得：

将(10)式代入(9)式中，可得：

为分析方便，假设在讨论加工周期误差时，不讨论其它加工误差的存在。当N=5，入射光强为λ=632.8nm，理想狭缝宽a=0.0020mm，理想周期(光栅常数)d=0.0100mm时，如果有3个周期是理想的，而有两个周期因加工出现误差0.5μm或1.0μm的误差。针对这一情况，则利用(8)式进行计算，得到衍射光分布结果如图 2所示。

Figure 2.  Comparison of diffraction light intensity distribution between ideal grating and error grating with the increase of cylce by formula (8)

图 2中曲线的变化显示了存在周期误差时, ±1, ±2, ±3和±4级次的衍射光强均会发生变化，随着周期误差的增加，各衍射主极大将朝着零级方向靠近。另外衍射光强会有相应降低，具体降低的值如表 1所示。

表 1中给出了光栅周期存在5%和10%的误差时，各级次衍射光强的变化率。相较于理想光栅，随着误差增大，各级次衍射光强会有不同程度的减小，级次越高的衍射光强下降得越快。对所得数据进行归一化处理，如图 3所示。

Figure 3.  Relative diffraction light intensity at ±1 order, ±2 order, ±3 order and ±4 order of each group with the increase of cycle by formula (8)

从图 3中可知，在加工过程中造成的矩形光栅周期误差。对比理想光栅可知，当加工误差使得光栅部分周期变大时，将会造成光栅(除零级外)各级次主极大衍射光强度都变小，但随着周期误差变大，各级次衍射光强的变化程度不同。在该误差范围内，±1级处的衍射光强下降程度相对较小，而其它级次的衍射光强度随着级次和误差的增大，衍射光强的下降程度越来越大。由此可知，该周期误差会降低矩形光栅各级次的衍射效率，而且随着误差增大，级次越高的衍射效率下降得越快。

当加工误差使周期变小时，设其中两个周期因加工出现误差0.5μm或1.0μm的误差。则利用(8)式进行计算，得到的衍射光分布结果如图 4所示。

Figure 4.  Comparison of diffraction light intensity distribution between ideal grating and error grating with the decrease of cylce by formula (8)

图 4中曲线的变化显示了存在周期误差时, ±1, ±2, ±3和±4级次的衍射光强均会发生变化，随着周期误差的增加，各衍射主极大将远离零级方向朝两边移动。另外衍射光强会有相应降低，具体降低的值列于表 2中。

表 2中给出了光栅周期存在-5%和-10%的误差时，各级次衍射光强的变化率。相较于理想光栅，随着误差增大，各级次衍射光强会有不同程度的减小，级次越高的衍射光强下降得越快。对所得数据进行归一化处理，如图 5所示。

Figure 5.  Relative diffraction light intensity at ±1 order, ±2 order, ±3 order and ±4 order of each group with the decrease of cycle by formula (8)

从图 5中可知，对比理想光栅可得，当加工误差使得光栅部分周期变小时，同样会造成(除零级外)的各主极大衍射光强度变小；但随着光栅周期误差变大时，各级次衍射光强度的变化程度不同。在该误差范围内，±1级处的衍射光强下降程度相对较小，而其它级次的衍射光强度随着级次和误差的增大，衍射光强的下降程度越来越大。由此可知，该周期误差会降低矩形光栅各级次的衍射效率，而且随着误差增大，级次越高的衍射效率下降得越快。

当N=5、入射光强λ=632.8nm、理想周期(光栅常数)d=0.0100mm、理想狭缝宽a=0.0020mm时，如果有3个狭缝是理想的，而有两个狭缝因加工出现0.5μm或1.0μm的误差。针对这一情况，则利用(11)式进行计算，得到的衍射光分布结果如图 6所示。

Figure 6.  Comparison of diffraction light intensity distribution between ideal grating and error grating with the increase of slit width by formula (11)

图 6中曲线的变化显示了存在缝宽误差时, ±1, ±2, ±3和±4级次的衍射光强均会发生变化，随着缝宽误差的增加，各衍射主极大光强会有相应降低，具体降低的值列于表 3中。

表 3中给出了光栅缝宽存在25%和50%的误差时，各级次衍射光强的变化率。相较于理想光栅，随着误差增大，各级次衍射光强会有不同程度的减小，级次越高的衍射光强下降得越快。对所得数据进行归一化处理，如图 7所示。

Figure 7.  Relative diffraction light intensity at ±1 order, ±2 order, ±3 order and ±4 order of each group with the increase of slit width by formula (11)

从图 7可知，在加工制作矩形光栅过程中产生的缝宽误差，对比理想光栅可得，当加工误差使得光栅部分狭缝宽度变大时，会造成光栅各主极大(除零级外)衍射光强度减小，但随着光栅缝宽误差变大时，各级次衍射光强度的变化程度不同。在该误差范围内，±1级处的衍射光强下降程度相对较小，而其它级次的衍射光强度随着级次和误差的增大，衍射光强的下降程度越来越大。由此可知，该缝宽误差会降低矩形光栅各级次的衍射效率，而且随着误差增大，级次越高的衍射效率下降得越快。

当加工误差使缝宽变小时，设其中两个缝宽因加工出现0.5μm或1μm的误差。则利用(11)式进行计算，得到的衍射光分布结果如图 8所示。

Figure 8.  Comparison of diffraction light intensity distribution between ideal grating and error grating with the decrease of slit width by formula (11)

图 8中曲线的变化显示了存在缝宽误差时, ±1, ±2, ±3和±4级次的衍射光强均会发生变化，随着缝宽误差的增加，各衍射主极大光强会有相应增加，具体增加的值列于表 4中。

表 4中给出了光栅缝宽存在-25%和-50%的误差时，各级次衍射光强的变化率。相较于理想光栅，随着误差增大，各级次衍射光强会有不同程度的减小，级次越高的衍射光强增大得越快。对所得数据进行归一化处理，如图 9所示。

Figure 9.  Relative diffraction light intensity at ±1 order, ±2 order, ±3 order and ±4 order of each group with the decrease of slit width by formula (11)

从图 9可知，对比理想光栅可得，当加工制作光栅出现部分狭缝宽度变小的误差时，会造成光栅各主极大(除零级外)的衍射光强度变大，但随着光栅缝宽误差变大时，各级次衍射光强度的变化程度不同。在该误差范围内，±1级处的衍射光强增大程度相对较小，而其它级次的衍射光强度随着级次和误差的增大，衍射光强的增大程度越来越大。由此可知，该缝宽误差会提高矩形光栅各级次的衍射效率，而且随着误差增大，级次越高的衍射效率增大得越快。

通过对光栅加工误差的计算模拟，分析了光栅局部周期或缝宽误差对光栅各级衍射效率的影响。综合上述实验结果可得:当加工误差使得光栅局部周期变大或是变小，以及光栅局部缝宽变大时，相较于理想光栅而言，都会造成光栅各主极大(除零级)衍射光强变小，在该误差范围内，±1级处的衍射光强变化较小，当误差变大时，级次越高对应的主极大衍射光强下降得越快; 而当加工误差使得光栅局部缝宽变小时，相较于理想光栅而言，会造成光栅各主极大(除零级)衍射光强变大，在该误差范围内，±1级处的衍射光强变化较小，当误差变大时，级次越高对应的主极大衍射光强增大得越快。因此，上述结论对实际加工制作光栅具有一定的启示作用，即尽可能地减小加工制作光栅产生的部分周期误差，而对于光栅部分，缝宽误差则可以适当地保持在合理的范围内。