广角衍射光学元件的优化设计

杨亮亮，刘成林，赵勇兵，沈法华，赵 琪，刘建利

〈系统与设计〉

广角衍射光学元件的优化设计

杨亮亮，刘成林，赵勇兵，沈法华，赵 琪，刘建利

（盐城师范学院 物理与电子工程学院，江苏省智能光电器件与测控工程研究中心，江苏 盐城 224007）

为进一步研究入射角度的增大对衍射光学元件（diffractive optical element, DOE）衍射效率及微结构高度等参数的影响，分析了入射角度和周期宽度对带宽积分平均衍射效率的影响。基于扩展标量衍射理论，建立了DOE的微结构高度与入射角度和周期宽度的数学模型，提出了工作在一定入射角度范围内，基于复合带宽积分平均衍射效率（comprehensive polychromatic integral diffraction efficiency, CPIDE）最大化实现设计波长和微结构高度等结构参数的优化设计方法。以工作在红外波段的DOE为例进行分析。结果表明：当相对周期宽度为20，入射角度范围为0°～40°时，该DOE的CPIDE为94.15%，微结构高度为1.3396mm。该设计方法可以实现广角DOE的优化设计。

衍射光学元件；衍射效率；微结构高度；入射角度

0 引言

近年来，鉴于衍射光学元件能够为系统提供一定的设计自由度，被广泛地应用于各种光学系统中，如成像系统、光束整形系统、复眼和3D显示等[1-5]。衍射效率和带宽积分平均衍射效率（polychromatic integral diffraction efficiency，PIDE）是决定DOE工作波段的重要参数。DOE的衍射效率对入射角度具有一定的依赖性。入射角度的增大会降低其衍射效率，进一步影响折衍混合光学系统的成像质量。

目前，成像光学系统中的DOE都是利用标量衍射理论（scalar diffraction theory，SDT）进行设计的，但该理论没有考虑入射角度对微结构高度的影响[6]。随着微结构表面入射角的增大，DOE的衍射效率会不断下降[7-8]。当微结构高度和波长处于同一数量级时，采用SDT计算衍射效率的准确度会大幅度下降，此时可以利用矢量衍射理论（vector diffraction theory，VDT）进行分析计算[9-12]。但是，VDT很难通过优化设计微结构高度等参数实现衍射效率的最大化。扩展标量衍射理论（extended scalar diffraction theory，ESDT）考虑了入射角度这一参数，其计算结果要比SDT更加精确，能够简化VDT的计算时间，并实现DOE的优化设计[13]。对于工作在可见光波段的双层DOE，文献[14]讨论了基于SDT和VDT计算的入射角度对衍射效率的影响。文献[15]基于ESDT讨论了周期宽度对微结构高度和PIDE的影响，但并没有给出DOE结构参数的优化设计。对于工作在一定入射角度范围内的DOE，基于ESDT对PIDE的优化设计未见报道。

本文基于ESDT，提出了工作在一定入射角度范围内，基于复合带宽积分平均衍射效率（comprehensive PIDE，CPIDE）最大化实现DOE设计波长、微结构高度等结构参数的优化设计方法。以工作在红外波段的DOE为例进行了分析与讨论。该方法可以实现工作在一定入射角度范围内的DOE结构参数的优化设计，特别是在相对周期宽度不是很大的情况下。

1 理论模型

当光线以入射角传播到DOE的微结构上时，如图1所示，根据衍射光栅公式，得到DOE的衍射光栅方程为：

(rsind－isin)＝(1)

式中：为光栅周期；i与r分别为相应介质材料在波长时的折射率；d为衍射角；为衍射级次。依据Snell折射定律，考虑衍射微结构对光线传播的影响，有：

isin(＋)＝rsin(r＋) (2)

图1 光线经过DOE的传播模型

式中：为微结构表面的倾角，tan＝/，r为折射角。当衍射角等于折射角，即r＝d时，第衍射级次的衍射效率最大。利用公式(1)和(2)，得到DOE的表面微结构高度为：

由公式(3)可知，当衍射面两端介质材料确定后，DOE的微结构高度与波长、入射角度和周期宽度有关。当DOE工作在正入射的状态下，微结构高度可以表示为：

当周期宽度远大于波长并且光线正入射时，得到基于SDT的DOE的微结构高度为：

可见，基于SDT，当介质材料和衍射级次确定后，微结构高度仅由设计波长0决定。

基于ESDT，斜入射时DOE的衍射效率为：

利用公式(3)和(6)可以计算周期宽度和入射角度对DOE衍射效率的影响。若把公式(6)中的换成SDT0即得到SDT的计算结果。

若DOE工作在min～max波段范围时，其PIDE为：

要基于PIDE的最大化确定微结构高度的大小，需要利用公式(7)，在周期宽度已知的前提下，确定公式(3)中的设计波长0和设计入射角度0，即可计算得到微结构高度0，即：

若工作在成像光学系统中的DOE，其入射角度范围为min～max，则DOE在整个工作入射角度范围内的CPIDE为：

若DOE工作在几个分离的入射角度范围内，如变焦系统，则其CPIDE为：

式中：min和max分别表示第个入射角范围的最小和最大入射角；为第个入射角范围的权重；表示总的入射角范围数量。

2 基于SDT的仿真与分析

以工作在红外波段1.4～2.2mm的DOE为例，基底材料采用硫化锌，衍射级次＝1。假设光束从空气介质入射到衍射基底，如图2所示为基于SDT计算得到的正入射时DOE的PIDE与波长的关系。在整个波段范围内，PIDE最高为94.47%，此时对应的峰值波长为设计波长，即1.7410mm，利用公式(5)计算得到微结构高度SDT0为1.3726mm。

图2 正入射时的PIDE与波长的关系

当入射角度分别为20°和40°时，DOE的衍射效率与波长的关系如图3所示。正入射时，在设计波长处的衍射效率为100%，最低衍射效率为81.07%。随着入射角度的增大，100%衍射效率对应的设计波长向长波方向移动；而且在工作波段范围内的最低衍射效率呈现快速下降的变化趋势。

图4给出了PIDE与入射角度的关系。当入射角度分别为20°、40°和60°时，DOE在整个波段的PIDE分别为94.27%、91.13%%和78.05%。假设各个视场的权重相同，当DOE分别工作在0°～20°、0°～40°和0°～60°入射角度范围内时，利用公式(10)计算得到CPIDE分别为94.43%、93.81%和91.13%。可见，基于正入射计算得到的微结构高度，随着入射角度或入射角度范围的增大，对应的PIDE或CPIDE逐渐减小。

图3 衍射效率与波长的关系

图4 PIDE与入射角度的关系

3 基于ESDT的仿真与分析

基于ESDT，首先分析周期宽度和入射角度对DOE微结构高度的影响；然后，分析一定的周期宽度和不同入射角度时DOE的衍射效率；最后，基于CPIDE最大化实现微结构高度等结构参数的优化设计。

3.1 微结构高度的分析

由公式(3)可知，DOE的微结构高度的大小与周期宽度和入射角度有关。当入射角度分别为0°、20°、40°和60°时，DOE的微结构高度与相对周期宽度（周期宽度与波长1.7410mm的比值）的关系如图5所示。正入射时，当相对周期宽度为5时，衍射微结构高度为1.3758mm，与SDT0相比，增大了0.0032mm。如表1所示，当入射角度为40°时，相对周期宽度分别为10或无穷大时，对应的微结构高度分别为1.2511mm和1.2355mm，与SDT0相比，分别减小了8.85%和9.99%。如图5和表1所示，当入射角度或者相对周期宽度改变时，基于ESDT计算得到的微结构高度与SDT的偏差不同。所以，当入射角度较大时，需要考虑入射角度和相对周期宽度对微结构高度的影响。

图5 微结构高度与周期宽度的关系

表1 微结构高度与周期宽度的关系

3.2 斜入射时的仿真与分析

当相对周期宽度确定为20时，在上述4个入射角度情况下，DOE的PIDE与波长的关系如图6所示。随着入射角度的增大，PIDE最大值对应的峰值设计波长向短波方向移动。利用公式(4)计算得到微结构高度如表2所示。伴随着设计波长的减小，微结构高度随入射角度的增大也减小。

图6 不同入射角度时PIDE与波长的关系

表2 基于带宽积分平均衍射效率最大化确定的结构参数

利用上述计算得到的设计波长和微结构高度，计算DOE的衍射效率如图7所示。入射角度为20°时，采用SDT和ESDT设计时的DOE在整个工作波段范围内的衍射效率图7(a)所示，衍射效率最小值分别为76.18%和81.65%，提高了5.47%。当入射角度分别增大到40°和60°时，如图7(b)、(c)所示，衍射效率最小值分别从59.25%增大到81.17%，从29.39%增大到81.21%，分别提高了21.92%和51.82%。可见，当入射角度偏离正入射时，利用基于ESDT计算得到的设计波长和微结构高度能够显著提高衍射效率。

图7 衍射效率与波长关系的对比

3.3 一定入射角度范围工作时的仿真与分析

同样假设DOE的相对周期宽度确定为20，当DOE工作时的入射角度范围为0°～20°，并假设各个视场的权重因子相同，DOE的PIDE与波长和入射角度的关系如图8所示，基于CPIDE的最大化（94.45%），得到设计波长为1.73mm，设计角度为8.25°，进一步计算得到DOE的微结构高度为1.3615mm。若DOE工作的入射角度范围增大到0°～40°或0°～60°时，其微结构高度和CPIDE等参数如表3所示。可见，随着DOE衍射面入射角度范围的增大，其微结构高度和CPIDE都逐渐减小。

图8 一定周期宽度时的PIDE与入射角度和波长的关系

表3 基于带宽积分平均衍射效率最大化确定的结构参数

4 结论

本文基于ESDT，建立了DOE的微结构高度与入射角度、周期宽度的数学关系模型，提出了工作在一定入射角度范围内，基于CPIDE最大化实现设计波长和微结构高度等结构参数的优化设计方法。对工作在红外波段的DOE进行仿真分析。当入射角度为40°时，对比SDT，基于ESDT计算得到的在工作波段范围内衍射效率的最小值提高了29.39%；当DOE工作在0°～40°范围内时，通过优化设计得到微结构高度为1.3396mm，CPIDE为94.15%。随着DOE在各类光学系统中的广泛应用，该方法为工作在较大入射角度范围内DOE的优化设计提供了理论依据。

[1] Greisukh G I, Ezhov E G, Zakharov O A, Diffractive microstructures of zoom lenses for visible and near-infrared ranges based on novel optical plastics[J]., 2022, 89(3): 127-131.

[2] 田晓航, 薛常喜. 小F数红外双波段无热化折衍摄远物镜设计[J]. 光学学报, 2022, 42(14): 181-187. TIAN Xiaohang, XUE Changxi. Athermalization design of small f-number refractive-diffractive telephoto objective lens in infrared dual-band[J]., 2022, 42(14): 181-187.

[3] 段慧慧, 杨艳芳, 何英,等. 4π聚焦系统中衍射光学元件对聚焦场多光球结构的影响[J]. 光学学报, 2021, 41(20): 174-179. DUAN Huihui, YANG Yanfang, HE Ying, et al. Influence of diffractive optical elements on multiple spherical spots in a 4π focusing system[J]., 2021, 41(20): 174-179.

[4] 张博, 崔庆丰, 朴明旭, 等. 双波段多层衍射光学元件的基底材料选择方法研究及其在变焦系统中的应用[J]. 光学学报, 2020, 40(6): 0605001. ZHANG Bo, CUI Qingfeng, PIAO Mingxu, et al. Substrate material selection method for dual-band multilayer diffractive optical elements and its application in the zoom system[J]., 2020, 40(6): 0605001.

[5] Shima G G, David Fischer, Stefan Sinzinger. Multifocal multi-value phase zone plate for 3D focusing[J].2019, 58(32): 8943-8949.

[6] XUE Changxi, CUI Qingfeng. Design of multilayer diffractive optical elements with polychromatic integral diffraction efficiency[J]., 2010, 35(7): 986-988.

[7] YANG Liangliang, LIU Chenglin, LI Shengqiang. Optimal design of depth-scaling error for multilayer diffractive optical elements with oblique incidence[J]., 2017, 56(15): 4532-4536.

[8] GAO Long, To Suet, YANG Hongfang, et al. Effect of assembling errors on the diffraction efficiency for multilayer diffractive optical elements[J]., 2014, 53: 7341-7347.

[9] Laborde V, Loicq J, Habraken S. Modeling infrared behavior of multilayer diffractive optical elements using Fourier optics[J]., 2021, 60(7): 2037-2045.

[10] YANG Hongfang, XUE Changxi. Sensitivity of diffraction efficiency to period width errors for multilayer diffractive optical elements[J]., 2018, 57(14): 855-860.

[11] PANG H, YIN S Y, DENG Q L, et al. A novel method for the design of diffractive optical elements based on the Rayleigh–Sommerfeld integral[J]., 2015, 70: 38-44.

[12] Noponen E, Turunen J. Binary high-frequency-carrier diffractive optical elements: electromagnetic theory[J]., 1994, 11(3): 1097-1109.

[13] Swanson G J. Binary Optics Technology: Theoretical Limits on the Diffraction Efficiency of Multilevel Diffractive Optical Elements[R/OL]. M.I.T. Technical Report, 1991, https://api.semanticscholar.org/ CorpusID: 26445902.

[14] Greisukh G I, Danilov V A, Ezhov E G, et al. Comparison of electromagnetic and scalar methods for evaluation of efficiency of diffractive lenses for wide spectral bandwidth[J]., 2015, 338: 54-57.

[15] HUO Furong, WANG Wensheng, XUE Changxi. Limits of scalar diffraction theory for multilayer diffractive optical elements[J]., 2016, 127: 5688-5694.

Optimal Design of Wide Angle Diffractive Optical Element

YANG Liangliang，LIU Chenglin，ZHAO Yongbing，SHEN Fahua，ZHAO Qi，LIU Jianli

( Jiangsu Province Intelligent Optoelectronic Devices and Measurement-Control Engineering Research Center, School of Physics and Electronics Engineering, Yancheng Teachers University, Yancheng 224007, China)

The influence of the incident angle on the diffraction efficiency and microstructure height of the diffractive optical element (DOE) was analyzed to further study the influence of the incident angle and period width on the polychromatic integral diffraction efficiency (PIDE). Based on the extended scalar diffraction theory (ESDT), a mathematical model of the relationship among the microstructure height, incident angle, and period width of the DOE was established. An optimal design method for structural parameters, such as the design wavelength and microstructure height, was proposed based on maximizing the comprehensive PIDE (CPIDE) within a certain range of incident angles. A DOE operating within the infrared waveband was considered as an example. The results indicate that when the relative period width is 20 and the incidence angle range is 0° to 40°, the CPIDE of the DOE is 94.15%, and the microstructure height is 1.3396mm. This design method can realize the optimal design of a wide-angle DOE.

diffractive optical element, diffraction efficiency, microstructure height, incident angle

O436

A

1001-8891(2024)03-0256-05

2023-02-15；

2023-03-31.

杨亮亮（1986-），女，博士，副教授，主要从事衍射光学和光学设计方面的研究。E-mail: yang_liangliang@163.com。

国家自然科学基金（61904158）；江苏省智能光电器件与测控工程研究中心开放基金（306054014）；江苏省重点研发计划-社会发展项目（BE2021685）江苏省自然科学基金面上项目（BK20201475）。