搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

激光照明条件对超振荡平面透镜聚焦性能的影响

刘康 何韬 刘涛 李国卿 田博 王佳怡 杨树明

引用本文:
Citation:

激光照明条件对超振荡平面透镜聚焦性能的影响

刘康, 何韬, 刘涛, 李国卿, 田博, 王佳怡, 杨树明

Effect of laser illumination conditions on focusing performance of super-oscillatory lens

Liu Kang, He Tao, Liu Tao, Li Guo-Qing, Tian Bo, Wang Jia-Yi, Yang Shu-Ming
PDF
HTML
导出引用
  • 超振荡平面透镜(super-oscillatory lens, SOL)是近几年出现的新型平面光学元件, 基于矢量角谱理论设计了振幅型和相位型SOL, 采用时域有限差分法对衍射聚焦光场进行严格电磁仿真计算, 研究发现: 当照明激光束腰半径w0小于SOL半径a时, 超衍射极限聚焦能力明显下降, 聚焦光斑强度衰减超过50%; 束腰半径w0对相位型SOL影响更剧烈, 且会发生显著正向焦移; 当w0不小于2a时可获得接近理想的聚焦特性. 倾斜照明条件下, 大数值孔径振幅型SOL一般允许的倾斜角度小于10°, 而相位型SOL具有宽广的倾角适应性(可超过40°), 聚焦光斑会发生横向展宽, 且强度急剧下降. 大数值孔径SOL对无限远点物成像会产生显著的负畸变和波动变化的场曲, 小数值孔径SOL在宽视场范围内则无畸变. 本文研究结果为SOL在超衍射极限聚焦、超分辨显微成像、飞秒激光直写微纳加工等领域的实际应用提供重要理论支撑.
    Super-oscillatory lens (SOL), a new type of planar optical element developed in recent years, may play an important role in the integrated optics, microscopy, advanced sensor, and astronomy. Based on the vectorial angular spectrum theory and genetic algorithm, both binary amplitude-type and phase-type SOLs are designed. Various sub-diffraction focusing properties can be realized by optimizing the design procedure. In order to investigate the focusing characteristics of SOLs under different illumination conditions, rigorous electromagnetic simulation calculations of the diffracted focusing light field are implemented by the finite-difference time-domain method. The results show that when the beam waist radius w0 of the illuminating laser is less than the SOL radius a, not only the capability of super-diffraction limit focusing will decrease significantly, but also the intensity of the focal spot will attenuate by more than 50%. Comparing with the amplitude-type SOL, the waist radius w0 has a strong effect on the phase-type SOL and causes a significant focus to shift in the positive direction. However, if w0 is larger than 2a, the ideal focusing characteristics of SOL can be maintained. Under the condition of oblique illumination, the high numerical aperture amplitude-type SOL generally only allows a small inclination angle of less than 10°, while the phase-type SOL has a wide inclination adaptability (can exceed 40°) regardless of the numerical aperture. For the latter, the focal spot will expand laterally and the intensity will decrease sharply with the increase of inclination angle. As for low numerical aperture phase-type SOL, the focusing characteristics, including focal spot size, focusing intensity and the angular position of the focus, can keep stable within an inclination angle of 18°. For imaging infinitely distant objects, the oblique illumination will produce a fluctuating field curvature and significant negative distortion for high numerical aperture SOLs, while for the low numerical aperture SOLs, the field curvature increases with inclination angle increasing and the distortion disappears almost. The research results of this paper provide an important theoretical basis for practical applications of super-oscillatory lens in the fields of sub-diffraction light focusing, super-resolution microscopic imaging, and micro-nano processing of femtosecond laser direct writing.
      通信作者: 刘涛, liu8483@xjtu.edu.cn ; 杨树明, shuming.yang@mail.xjtu.edu.cn
    • 基金项目: 国家自然科学基金优秀青年科学基金(批准号: 51722509)、国家重点研发计划(批准号: 2017YFB1104700)和陕西省自然科学基础研究计划(批准号: 2020SF-170)资助的课题
      Corresponding author: Liu Tao, liu8483@xjtu.edu.cn ; Yang Shu-Ming, shuming.yang@mail.xjtu.edu.cn
    • Funds: Project supported by the Excellent Young Scientists Fund of the National Natural Science Foundation of China (Grant No. 51722509), the National Key R&D Program of China (Grant No. 2017YFB1104700), and the Natural Science Foundation for Basic Research of Shaanxi Province, China (Grant No. 2020SF-170)
    [1]

    Rogers E T F, Lindberg J, Roy T, Savo S, Chad J E, Dennis M R, Zhelidev N I 2012 Nat. Mater. 11 432Google Scholar

    [2]

    刘涛 2013 博士学位论文 (哈尔滨: 哈尔滨工业大学)

    Liu T 2013 Ph. D. Dissertation (Harbin: Harbin Institute of Technology) (in Chinese)

    [3]

    Francia G T D 1952 IL Nuovo Cimento 9 426Google Scholar

    [4]

    邱丽荣 2005 博士学位论文 (哈尔滨: 哈尔滨工业大学)

    Qiu L R 2005 Ph. D. Dissertation (Harbin: Harbin Institute of Technology) (in Chinese)

    [5]

    Diao J S, Yuan W Z, Yu Y T, Zhu Y C, Wu Y 2016 Opt. Express 24 1924Google Scholar

    [6]

    Wu J, Wu Z X, He Y G, Yu A P, Zhang Z H, Wen Z Q, Chen G 2017 Opt. Express 25 6274Google Scholar

    [7]

    Liu T, Tan J B, Liu J, Wang H T 2013 Opt. Lett. 38 2742Google Scholar

    [8]

    Liu T, Liu Q, Yang S M, Jiang Z D, Wang T, Zhang G F 2017 Appl. Opt. 56 3725Google Scholar

    [9]

    Liu T, Tan J B, Liu J, Wang H T 2013 Opt. Express 21 15090Google Scholar

    [10]

    Liu T, Shen T, Yang S M, Jiang Z D 2015 J. Opt. 17 035610Google Scholar

    [11]

    Liu T, Liu Q, Yang S M, Jiang Z D, Wang T, Yang X K 2017 Opt. Commun. 393 72Google Scholar

    [12]

    Yuan G H, Rogers Edward T F, Zheludev N I 2017 Light: Sci. Appl. 6 e17036Google Scholar

    [13]

    Qin F, Huang K, Wu J F, Teng J H, Qiu C W, Hong M H 2017 Adv. Mater. 29 1602721Google Scholar

    [14]

    武靖 2018 硕士学位论文 (重庆: 重庆大学)

    Wu J 2018 M. S. Thesis (Chongqing: Chongqing university) (in Chinese)

    [15]

    Nagarajan A, Stoevelaar L P, Silvestri F, et al. 2019 Opt. Express 27 20012Google Scholar

    [16]

    Liu T, Wang T, Yang S M, Sun L, Jiang Z D 2015 Opt. Express 23 32139Google Scholar

    [17]

    Yang S M, Wang T, Liu T, Jiang Z D 2016 Opt. Commun. 372 166Google Scholar

    [18]

    Yu Y T, Li W L, Li H Y, Li M Y, Yuan W Z 2018 Nanomaterials 8 185Google Scholar

    [19]

    Ni H B, Yuan G H, Sun L D, Chang N, Zhang D, Chen R P, Jiang L Y, Chen H Y, Gu Z Z, Zhao X W 2018 RSC Adv. 8 20117Google Scholar

    [20]

    Luneburg R K 1966 Mathematical Theory of Optics (Berkeley: University of California Press) p305

    [21]

    Liu T, Yang S M, Jiang Z D 2016 Opt. Express 24 16297Google Scholar

    [22]

    Chen G, Wu Z X, Yu A P, et al. 2016 Sci. Rep. 6 37776Google Scholar

    [23]

    张蕾, 蔡阳健, 陆漩辉 2004 物理学报 53 1777Google Scholar

    Zhang L, Cai Y J, Lu X H 2004 Acta Phys. Sin. 53 1777Google Scholar

    [24]

    Rogers E T F, Savo S, Lindberg J, Roy T, Dennis M R, Zheludev N 2013 Appl. Phys. Lett. 102 031108Google Scholar

    [25]

    郁道银, 谈恒英 2011 工程光学 (北京: 机械工业出版社) 第118页

    Yu D Y, Tan H Y 2011 Engineering Optics (Beijing: China Machine Press) p118 (in Chinese)

  • 图 1  SOL衍射聚焦示意图

    Fig. 1.  Schematic diagram of diffraction focusing of SOL.

    图 2  SOL1的聚焦光场FDTD计算结果 (a) 轴向光强分布; (b) 焦平面内横向光强分布(Y方向)

    Fig. 2.  Focused light field of SOL1 by FDTD simulation: (a) On-axis intensity distribution; (b) transverse intensity distribution in the focal plane (Y-direction).

    图 3  SOL聚焦性能变化(FDTD) (a) 光斑横向尺寸; (b) 光斑中心强度

    Fig. 3.  Focusing performance of SOL (FDTD): (a) Transverse size; (b) central intensity.

    图 4  SOL2的聚焦光场FDTD计算结果 (a) 轴向光强分布; (b) 焦平面内横向光强分布(Y方向)

    Fig. 4.  Focused light field of SOL2 by FDTD simulation: (a) On-axis intensity distribution; (b) transverse intensity distribution in the focal plane (Y-direction).

    图 5  SOL1聚焦光斑电场分布的FDTD结果 (a) 均匀平面波照明; (b) w0 = a/2的高斯光束照明

    Fig. 5.  Electric field distribution of SOL1 by FDTD simulation: (a) Uniform plane beam illumination; (b) w0 = a/2 Gaussian beam illumination.

    图 6  BG径向偏振光束照明条件下SOL聚焦光场的VAS计算结果 (a), (b) SOL1; (c), (d) SOL2

    Fig. 6.  Focused light intensity of SOL by VAS calculation under BG radially polarized illumination: (a), (b) SOL1; (c), (d) SOL2.

    图 7  SOL对倾斜平行光的衍射聚焦成像示意图

    Fig. 7.  Schematic diagram of diffraction focused imaging by SOL under oblique illumination.

    图 8  倾斜平行光照明的SOL衍射聚焦示意图

    Fig. 8.  Schematic diagram of diffraction focusing by SOL under oblique illumination.

    图 9  振幅型SOL的倾斜照明FDTD聚焦计算结果 (a), (c), (e) SOL5; (b), (d), (f) SOL6

    Fig. 9.  Focusing results of amplitude-type SOL under oblique illumination: (a), (c), (e) SOL5; (b), (d), (f) SOL6.

    图 10  倾斜照明SOL的聚焦光斑尺寸 (a) FWHMx ; (b) FWHMz

    Fig. 10.  Focal spot size of SOL under oblique illumination: (a) FWHMx ; (b) FWHMz .

    图 11  不同照明角度下SOL的聚焦光斑强度

    Fig. 11.  Focal intensity of SOL varying with oblique illumination angle.

    图 12  聚焦光斑角位置变化 (a) $\omega $-$\omega '$关系曲线; (b)$\omega $-$\varGamma $关系曲线

    Fig. 12.  Angle position of focusing spot: (a) $\omega $-$\omega '$; (b) $\omega $-$\varGamma $.

    图 13  SOL的场曲 (a) SOL3; (b) SOL4

    Fig. 13.  Field curvature of SOL: (a) SOL3; (b) SOL4.

    表 1  优化设计的SOL结构参数

    Table 1.  Structural parameters of optimized SOL.

    SOL类型r/μmD/μmfsol/μmNAri /μmti
    SOL1振幅0.2163.50.92[0.2, 0.4, 1.0, 1.2, 2.0, 2.2, 2.4, 2.6, 3.2,
    3.6, 3.8, 4.0, 4.2, 4.6, 4.8, 5.2, 5.4,
    6.0, 6.4, 6.8, 7.0, 7.4, 7.6, 8.0]
    [0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1,
    0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1]
    SOL2相位0.2163.50.92[0.2, 0.6, 1.0, 1.2, 1.6, 1.8, 2.0, 2.4,
    3.0, 3.2, 4.8, 5.0, 5.4, 5.8, 6.2,
    6.6, 6.8, 7.2, 7.4, 8.0]
    [–1, 1, –1, 1, –1, 1, –1, 1,
    –1, 1, –1, 1, –1, 1,
    –1, 1, –1, 1, –1, 1]
    SOL3相位0.2102.00.93[0.473, 0.885, 1.458, 1.969, 2.382,
    2.842, 3.044, 4.239, 4.745, 5.000]
    [1, j, 1, j, 1, j, 1, j, 1, j]
    SOL4相位0.2109.90.45[0.268, 0.468, 1.165, 1.365, 1.734,
    3.189, 3.399, 4.278, 4.792, 5.000]
    [j, 1, j, 1, j, 1, j, 1, j, 1]
    SOL5振幅0.2101.90.93[0.600, 0.899, 1.915, 2.190, 2.440,
    3.756, 4.076, 4.357, 4.769, 5.000]
    [0, 1, 0, 1, 0, 1, 0, 1, 0, 1]
    SOL6振幅0.21010.10.44[0.300, 0.506, 1.248, 1.460, 1.660, 2.885,
    3.085, 3.335, 3.755, 4.094, 4.294, 5.000]
    [0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1]
    下载: 导出CSV
  • [1]

    Rogers E T F, Lindberg J, Roy T, Savo S, Chad J E, Dennis M R, Zhelidev N I 2012 Nat. Mater. 11 432Google Scholar

    [2]

    刘涛 2013 博士学位论文 (哈尔滨: 哈尔滨工业大学)

    Liu T 2013 Ph. D. Dissertation (Harbin: Harbin Institute of Technology) (in Chinese)

    [3]

    Francia G T D 1952 IL Nuovo Cimento 9 426Google Scholar

    [4]

    邱丽荣 2005 博士学位论文 (哈尔滨: 哈尔滨工业大学)

    Qiu L R 2005 Ph. D. Dissertation (Harbin: Harbin Institute of Technology) (in Chinese)

    [5]

    Diao J S, Yuan W Z, Yu Y T, Zhu Y C, Wu Y 2016 Opt. Express 24 1924Google Scholar

    [6]

    Wu J, Wu Z X, He Y G, Yu A P, Zhang Z H, Wen Z Q, Chen G 2017 Opt. Express 25 6274Google Scholar

    [7]

    Liu T, Tan J B, Liu J, Wang H T 2013 Opt. Lett. 38 2742Google Scholar

    [8]

    Liu T, Liu Q, Yang S M, Jiang Z D, Wang T, Zhang G F 2017 Appl. Opt. 56 3725Google Scholar

    [9]

    Liu T, Tan J B, Liu J, Wang H T 2013 Opt. Express 21 15090Google Scholar

    [10]

    Liu T, Shen T, Yang S M, Jiang Z D 2015 J. Opt. 17 035610Google Scholar

    [11]

    Liu T, Liu Q, Yang S M, Jiang Z D, Wang T, Yang X K 2017 Opt. Commun. 393 72Google Scholar

    [12]

    Yuan G H, Rogers Edward T F, Zheludev N I 2017 Light: Sci. Appl. 6 e17036Google Scholar

    [13]

    Qin F, Huang K, Wu J F, Teng J H, Qiu C W, Hong M H 2017 Adv. Mater. 29 1602721Google Scholar

    [14]

    武靖 2018 硕士学位论文 (重庆: 重庆大学)

    Wu J 2018 M. S. Thesis (Chongqing: Chongqing university) (in Chinese)

    [15]

    Nagarajan A, Stoevelaar L P, Silvestri F, et al. 2019 Opt. Express 27 20012Google Scholar

    [16]

    Liu T, Wang T, Yang S M, Sun L, Jiang Z D 2015 Opt. Express 23 32139Google Scholar

    [17]

    Yang S M, Wang T, Liu T, Jiang Z D 2016 Opt. Commun. 372 166Google Scholar

    [18]

    Yu Y T, Li W L, Li H Y, Li M Y, Yuan W Z 2018 Nanomaterials 8 185Google Scholar

    [19]

    Ni H B, Yuan G H, Sun L D, Chang N, Zhang D, Chen R P, Jiang L Y, Chen H Y, Gu Z Z, Zhao X W 2018 RSC Adv. 8 20117Google Scholar

    [20]

    Luneburg R K 1966 Mathematical Theory of Optics (Berkeley: University of California Press) p305

    [21]

    Liu T, Yang S M, Jiang Z D 2016 Opt. Express 24 16297Google Scholar

    [22]

    Chen G, Wu Z X, Yu A P, et al. 2016 Sci. Rep. 6 37776Google Scholar

    [23]

    张蕾, 蔡阳健, 陆漩辉 2004 物理学报 53 1777Google Scholar

    Zhang L, Cai Y J, Lu X H 2004 Acta Phys. Sin. 53 1777Google Scholar

    [24]

    Rogers E T F, Savo S, Lindberg J, Roy T, Dennis M R, Zheludev N 2013 Appl. Phys. Lett. 102 031108Google Scholar

    [25]

    郁道银, 谈恒英 2011 工程光学 (北京: 机械工业出版社) 第118页

    Yu D Y, Tan H Y 2011 Engineering Optics (Beijing: China Machine Press) p118 (in Chinese)

计量
  • 文章访问数:  2812
  • PDF下载量:  50
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-19
  • 修回日期:  2020-05-30
  • 上网日期:  2020-06-05
  • 刊出日期:  2020-09-20

激光照明条件对超振荡平面透镜聚焦性能的影响

    基金项目: 国家自然科学基金优秀青年科学基金(批准号: 51722509)、国家重点研发计划(批准号: 2017YFB1104700)和陕西省自然科学基础研究计划(批准号: 2020SF-170)资助的课题

摘要: 超振荡平面透镜(super-oscillatory lens, SOL)是近几年出现的新型平面光学元件, 基于矢量角谱理论设计了振幅型和相位型SOL, 采用时域有限差分法对衍射聚焦光场进行严格电磁仿真计算, 研究发现: 当照明激光束腰半径w0小于SOL半径a时, 超衍射极限聚焦能力明显下降, 聚焦光斑强度衰减超过50%; 束腰半径w0对相位型SOL影响更剧烈, 且会发生显著正向焦移; 当w0不小于2a时可获得接近理想的聚焦特性. 倾斜照明条件下, 大数值孔径振幅型SOL一般允许的倾斜角度小于10°, 而相位型SOL具有宽广的倾角适应性(可超过40°), 聚焦光斑会发生横向展宽, 且强度急剧下降. 大数值孔径SOL对无限远点物成像会产生显著的负畸变和波动变化的场曲, 小数值孔径SOL在宽视场范围内则无畸变. 本文研究结果为SOL在超衍射极限聚焦、超分辨显微成像、飞秒激光直写微纳加工等领域的实际应用提供重要理论支撑.

English Abstract

    • 为了实现远场亚波长聚焦和超分辨成像, 英国南安普顿大学于2012年首次提出了超振荡透镜(super-oscillatory lens, SOL)[1], 作为一种新的平面光学元件, 其基本原理是利用微纳环带结构实现多光束衍射, 并由相干光场的精细干涉叠加在光学远场区域获得亚波长尺寸的聚焦光斑[1,2]. SOL是一种微纳同心圆环结构, 与常规的同心圆环光瞳滤波器(或称Toraldo光瞳滤波器)[3,4], 工作原理截然不同, 前者是利用光的衍射对光场直接聚焦, 后者必须工作于透镜或物镜系统中. SOL是一种随机优化产生的平面光学元件[1,2,5,6], 与熟知的菲涅耳波带片[7,8]也存在明显区别, 后者是利用波带法直接构造的疏密有致的同心圆环结构. 根据光场调制参数的不同, SOL可分为振幅型、相位型和复振幅型, 振幅和相位调制分别通过刻蚀金属膜和透明介质材料实现. 由于SOL具有结构轻薄、易于实现亚波长聚焦和衍射光场定制设计的突出优点, 近几年受到广泛的关注[5,6,9-12]. 主要包括广泛深入的优化设计方法研究[1,5,9-11], 和基于SOL的超分辨成像应用研究[1,13-15]. 对于引入加工误差后SOL聚焦性能变化的研究较少, 西安交通大学分析了金属膜的材料和厚度以及照明激光波长对SOL聚焦特性的影响[16,17]; 西北工业大学初步研究了SOL环带宽度误差和环带侧壁倾斜度对聚焦性能的影响[18]; 东南大学基于激光直写技术利用视场拼接的办法加工了大口径SOL并分析了拼接误差对聚焦性能的影响[19].

      目前为止, 激光照明条件对SOL聚焦性能的影响尚未见相关研究报道, 尚不清楚SOL实际应用中照明光源振幅分布、倾斜、安装偏差等客观条件下SOL聚焦性能如何变化, 这些情形都偏离了SOL优化设计的基本假设条件[1,2,9,10], 目前研究也尚未揭示SOL对无限远物点成像的像差分布. 本文结合矢量角谱(vectorial angular spectrum, VAS)理论和遗传算法(genetic algorithm, GA), 在激光均匀照明、垂直入射条件下优化设计典型SOL, 重点基于时域有限差分法(finite-difference time-domain, FDTD)对衍射光场进行严格电磁仿真计算, 通过大量的FDTD定量计算研究揭示实际应用中激光照明条件对SOL聚焦特性的影响规律, 并研究揭示离轴倾斜照明(或无限远点物成像)条件下的畸变和场曲像差分布特性, 研究结果将为SOL在当前超分辨聚焦、显微成像、微纳光刻等领域的广泛工程应用奠定坚实的理论基础.

    • 当矢量偏振光束垂直照射SOL表面时, 其衍射聚焦过程如图1所示. 根据VAS理论, 如果已知SOL后表面处的横向电场分布, 则SOL之后空间中任意一点的电场分布均可由VAS积分表达式计算得到[20], 而SOL后表面的电场分布近似等于入射场与SOL衍射孔径标量透过率函数的乘积[1,2,5,6,9]. 根据VAS理论对任意典型偏振态激光光束的衍射传播过程进行积分描述, 再设定需求目标对衍射光场进行单目标或多目标约束, 建立数学优化模型, 进而通过优化算法(如GA算法[2]或粒子群优化算法[1])反复迭代计算, 可以优化设计出性能满足(或接近)设计需求的SOL结构[1,2,5,6,9].

      图  1  SOL衍射聚焦示意图

      Figure 1.  Schematic diagram of diffraction focusing of SOL.

      针对最典型的线偏振光(偏振方向设为沿X方向), 基于VAS理论导出的积分公式[2,9], 利用GA算法和快速汉克尔变换算法[9,16], 分别优化设计出三组振幅型和相位型的SOL结构, 参数见表1. 其中, 照明激光波长λ = 632.8 nm, 光源垂直入射, 振幅设为均匀分布. SOL工作介质为空气, SOL的最大外径为D, 设计焦距为fsol, 环带最小径向宽为∆r, NA是SOL的数值孔径($NA = \sin \alpha $, $\alpha $是最大会聚半角). SOL的几何结构用每个环带的外环半径值ri表示(i = 1, 2, ···, N, N是最大环带数), 每个环带的透过率用ti表示, 其中振幅型SOL的透过率“1”和“0”分别表示该环带是透光和不透光, 而相位型SOL的相位调制量为ψ = 0, π/2或π, 对应ti = exp(jψ)等于“1”, “j”和“-1”(j是虚数单位). 第3节将基于典型SOL结构 (SOL1—6)采用严格电磁仿真方法——FDTD进行计算分析.

      SOL类型r/μmD/μmfsol/μmNAri /μmti
      SOL1振幅0.2163.50.92[0.2, 0.4, 1.0, 1.2, 2.0, 2.2, 2.4, 2.6, 3.2,
      3.6, 3.8, 4.0, 4.2, 4.6, 4.8, 5.2, 5.4,
      6.0, 6.4, 6.8, 7.0, 7.4, 7.6, 8.0]
      [0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1,
      0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1]
      SOL2相位0.2163.50.92[0.2, 0.6, 1.0, 1.2, 1.6, 1.8, 2.0, 2.4,
      3.0, 3.2, 4.8, 5.0, 5.4, 5.8, 6.2,
      6.6, 6.8, 7.2, 7.4, 8.0]
      [–1, 1, –1, 1, –1, 1, –1, 1,
      –1, 1, –1, 1, –1, 1,
      –1, 1, –1, 1, –1, 1]
      SOL3相位0.2102.00.93[0.473, 0.885, 1.458, 1.969, 2.382,
      2.842, 3.044, 4.239, 4.745, 5.000]
      [1, j, 1, j, 1, j, 1, j, 1, j]
      SOL4相位0.2109.90.45[0.268, 0.468, 1.165, 1.365, 1.734,
      3.189, 3.399, 4.278, 4.792, 5.000]
      [j, 1, j, 1, j, 1, j, 1, j, 1]
      SOL5振幅0.2101.90.93[0.600, 0.899, 1.915, 2.190, 2.440,
      3.756, 4.076, 4.357, 4.769, 5.000]
      [0, 1, 0, 1, 0, 1, 0, 1, 0, 1]
      SOL6振幅0.21010.10.44[0.300, 0.506, 1.248, 1.460, 1.660, 2.885,
      3.085, 3.335, 3.755, 4.094, 4.294, 5.000]
      [0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1]

      表 1  优化设计的SOL结构参数

      Table 1.  Structural parameters of optimized SOL.

    • 针对振幅型SOL1和相位型SOL2 (见表1), 分别采用束腰半径w0 = 2a, a, 2a/3, a/2的高斯光束进行照明, 其中$a$为SOL的半径(a = D/2), 高斯光束振幅分布(设在SOL衍射孔径平面内)表示为

      $g(r) = {A_0}\exp \left( { - {{{r^2}}}/{{w_0^2}}} \right), $

      其中A0是光束中心处的幅度. 利用FDTD三维严格电磁仿真, 振幅型SOL1采用100 nm厚的铝膜作为光的遮挡介质[1,16,21], 相位型SOL2的相位调制介质采用Si3N4材料, 在λ = 632.8 nm波长下Si3N4的折射率为1.91[22], 通过对Si3N4层的刻蚀产生深度为348 nm的沟槽进行相位π的调制. FDTD采用完全匹配层(perfectly matched layer, PML)作为吸收边界, 设置总场散射场(total-field scattered-field, TFSF)光源边界, FDTD仿真区域网格划分为20 nm × 20 nm × 20 nm (x, y, z), 更精细的网格划分可提高计算精度, 但同时对计算机存储空间的需求迅速增大.

      振幅型SOL1仿真计算得到的衍射场沿轴光强分布和焦平面内垂直于偏振方向的光强分布如图2所示, 其中束腰半径w0 = $ \infty $对应均匀平面波. 从图3(a)可见, 光源束腰半径w0$ \infty $减小至a/2时, 聚焦光斑的横向半高全宽(full width at half maximum, FWHM)逐渐展宽. 当高斯光的束腰半径w0大于SOL的半径a时, 聚焦光斑的FWHM仅展宽了1.4%, 光斑强度衰减为等幅度均匀平面波时的43.3%; 而当束腰半径w0小于SOL衍射孔半径a时, 聚焦光斑由287 nm (0.42λ/NA, w0 = a)迅速展宽为341 nm (0.50λ/NA, w0 = a/2), 光斑尺寸增大了19%, SOL的超衍射极限聚焦的能力逐渐下降(标量衍射理论框架下, 艾里斑的FWHM为0.51λ/NA), 同时聚焦光斑的光场强度急剧下降, 仅为原来的18%, 如图3(b)所示.

      图  2  SOL1的聚焦光场FDTD计算结果 (a) 轴向光强分布; (b) 焦平面内横向光强分布(Y方向)

      Figure 2.  Focused light field of SOL1 by FDTD simulation: (a) On-axis intensity distribution; (b) transverse intensity distribution in the focal plane (Y-direction).

      图  3  SOL聚焦性能变化(FDTD) (a) 光斑横向尺寸; (b) 光斑中心强度

      Figure 3.  Focusing performance of SOL (FDTD): (a) Transverse size; (b) central intensity.

      类似于振幅型SOL1, 相位型SOL2的聚焦光斑同样随着照明高斯光束的束腰半径的减小而逐渐展宽, 光斑强度迅速减弱, 如图3图4所示. 相位型SOL2的光强分布变化更为剧烈, 照明光束的振幅分布对其影响更大, 特别是当束腰半径w0a减小到a/2时, 轴向光斑的旁瓣显著增强, 并且伴随着显著的轴向正向焦移(160 nm), 横向光斑的FWHM由297 nm展宽为365 nm, 尺寸增大了23%, 光斑强度也减弱为原来的25%. 而当高斯光束的束腰半径w0 = a时, 横向光斑的FWHM仅展宽了7.2%, 光斑强度衰减为均匀照明时的49%.

      图  4  SOL2的聚焦光场FDTD计算结果 (a) 轴向光强分布; (b) 焦平面内横向光强分布(Y方向)

      Figure 4.  Focused light field of SOL2 by FDTD simulation: (a) On-axis intensity distribution; (b) transverse intensity distribution in the focal plane (Y-direction).

      从SOL1和SOL2的FDTD仿真计算结果不难发现, 对于均匀照明条件下优化设计的SOL, 若采用高斯光束进行照明使用(实际情形), 无论是振幅型还是相位型SOL, 聚焦光斑的分布总体仍然保持了一致的单焦点特性, 但激光的束腰半径会对SOL聚焦性能产生显著影响, 强度和光斑尺寸在${w_0} \geqslant D$(即${w_0} \geqslant 2 a$)时接近理想设计情形. 当束腰半径w0减小至SOL半径a时, 聚焦光斑尺寸明显展宽, 而光斑强度衰减超过50%, SOL聚焦光斑的分布影响相对较小. 当束腰半径w0小于SOL半径a后, 不仅聚焦光斑尺寸会发生显著扩展, 超衍射极限聚焦能力明显下降, 特别对于相位型SOL还会出现显著的正向焦移现象, 光斑的轴向旁瓣也会明显退化.

      当前大多数微纳光学研究使用的SOL半径在几十至百微米量级[1,5,9-11], 由于实际激光光束的束腰半径都在亚毫米至毫米(准直扩束后可增大至数毫米至厘米量级), 远大于SOL的半径, 因此实际照明情形可按均匀平面波对待. 但为了与常规的光学元件配合使用以及提高SOL的光能利用率, 需要设计毫米级甚至厘米级的大口径SOL[18,19], 此时SOL的衍射半径已经等于或大于照明激光束腰半径, 激光光源的振幅分布特性对SOL聚焦光斑的影响变得尤为显著. 因此, 如果光源束腰半径与SOL半径相近或更小时, 必须将高斯光束振幅分布代入优化设计过程, 对SOL进行更加精准的优化设计.

      进一步地对X方向线偏振光照明的SOL聚焦光斑矢量特性进行研究, 以SOL1为例, 在均匀平面波和w0 = a/2高斯光照明条件下焦平面各电场分布如图5所示. 可以看到平面波入射条件下聚焦光斑中电场EzEy分量在电场总场中的比重较大, 导致聚焦光斑呈扁椭圆形; 而高斯光照明条件下SOL聚焦光斑的EzEy分量却相对较小, 对合成光斑的分布影响减弱, 光斑最终呈轻微椭圆状. 通过分析不难发现, 由于采用高斯光束进行SOL照明时, 对聚焦光斑EzEy分量贡献最大的外环结构接收到的激光照明强度大大降低, 因此聚焦光斑中的电场EzEy分量减小而Ex成为最主要的分量, 合成光斑将在Ex分布基础上沿着照明光偏振方向轻微展宽.

      图  5  SOL1聚焦光斑电场分布的FDTD结果 (a) 均匀平面波照明; (b) w0 = a/2的高斯光束照明

      Figure 5.  Electric field distribution of SOL1 by FDTD simulation: (a) Uniform plane beam illumination; (b) w0 = a/2 Gaussian beam illumination.

      除了最常见的高斯平面波, 进一步对空心高斯光束[23]和贝塞尔-高斯(Bessel-Gaussian, BG)光束[9]等不同振幅分布照明光源条件下的SOL聚焦特性进行了分析.

      首先基于VAS理论对X方向线偏振空心高斯光束照明条件下的SOL聚焦光场分布进行了计算, 发现随着光束空心半径的增大焦平面光斑的EzEy分量占比急剧上升, 呈明显的“哑铃”形分布; 而光场轴向分布则逐渐展宽呈不均匀的“光针”状, 这样的变化与对SOL进行中心遮挡空间滤波情形类似[24].

      其次采用典型的BG振幅分布的径向偏振态光束, 针对SOL1和SOL2进行矢量衍射计算, 振幅服从[9]

      $g(r) = \exp \left( { - \frac{{{r^2}}}{{w_0^2}}} \right){{\rm{J}}_1}\left( {2\frac{r}{{{w_0}}}} \right),$

      其中J1为第一类一阶贝塞尔函数. 设高斯束腰半径等于SOL1半径(w0 = 8 μm), 基于VAS理论进行聚焦光场计算, 结果如图6所示. 可见, 采用径向偏振BG光束照明的SOL, 聚焦光斑横向尺寸与X线偏振均匀平面光照明条件下的Y向光斑尺寸大小接近且略微偏大, 而光斑轴向分布严重变形, 根本原因是这里使用的SOL是采用线偏振光设计的, 如果改为径向偏振光照明则其聚焦性能不能达到最优, 此时应该在径向偏振光照明条件下重新优化设计[9].

      图  6  BG径向偏振光束照明条件下SOL聚焦光场的VAS计算结果 (a), (b) SOL1; (c), (d) SOL2

      Figure 6.  Focused light intensity of SOL by VAS calculation under BG radially polarized illumination: (a), (b) SOL1; (c), (d) SOL2.

    • SOL对有一定倾斜角度(倾角)的单色相干平行光束直接聚焦(等价于对轴外无穷远物点进行衍射聚焦成像), 或者可将SOL用作管镜进行无限远物点的光学成像, 如图7所示. 此外, 由于装调误差等原因也可能导致入射光束方向与SOL元件表面的法线方向存在一定的夹角, 如图8所示. 这些都偏离了目前SOL的设计基本条件要求, 即均假设单色激光垂直照明, 此时SOL是否依旧能有效聚焦, 或者对倾斜平行光有一定的视场角范围, 以及聚焦光斑分布会如何变化, 将直接影响SOL的实际应用.

      图  7  SOL对倾斜平行光的衍射聚焦成像示意图

      Figure 7.  Schematic diagram of diffraction focused imaging by SOL under oblique illumination.

      图  8  倾斜平行光照明的SOL衍射聚焦示意图

      Figure 8.  Schematic diagram of diffraction focusing by SOL under oblique illumination.

      首先对振幅型SOL在倾斜入射下的聚焦性能进行了初步研究, 发现振幅型SOL对入射光束倾角(偏转角度)十分敏感. 图9为振幅型SOL的倾斜照明FDTD聚焦计算结果, 其中大数值孔径SOL5 (NA = 0.93)的聚焦光斑在倾斜角度$\omega = $ $ {12^ \circ } $时出现了明显的弥散现象, 而在$\omega = {15^ \circ }$时失去聚焦能力; 小数值孔径SOL6 (NA = 0.44)在$\omega = {18^ \circ }$倾斜照明条件下光斑出现显著变形, 在$\omega = {21^ \circ }$时光斑严重退化, 此时SOL失去对该倾角平行光束的良好聚焦. 总体而言, 振幅型SOL聚焦性能严重依赖于照明光束倾角大小, 一般仅允许相对较小角度($\omega \leqslant {10^ \circ }$)的倾斜照明, 小数值孔径SOL具有相对较弱的倾角依赖性.

      图  9  振幅型SOL的倾斜照明FDTD聚焦计算结果 (a), (c), (e) SOL5; (b), (d), (f) SOL6

      Figure 9.  Focusing results of amplitude-type SOL under oblique illumination: (a), (c), (e) SOL5; (b), (d), (f) SOL6.

      相比之下, 相位型SOL则可以在较大角度的倾斜照明条件下进行光束聚焦, FDTD仿真计算结果显示, 在倾斜角度最大达到$\omega = {42^ \circ }$时相位型SOL3和SOL4仍然能够实现良好聚焦, 远大于振幅型SOL允许的入射倾角大小. 除了对入射光倾斜角度具有更高的适应性, 相位型SOL的聚焦效率以及对入射光场调制的灵活性等也远优于振幅型SOL, 有望在未来获得更广泛的使用, 因此, 下面针对相位型SOL进行详细比较分析, 包含大数值孔径和小数值孔径SOL两类.

      对于NA = 0.93的相位型SOL3NA = 0.45的相位型SOL4, 将平行光束进行间隔为3°的等间距倾斜($\omega = {0^ \circ }—\;{42^ \circ }$), 对聚焦光场分别进行FDTD严格电磁仿真计算. FDTD仿真模型中SOL的相位调制层采用SiO2材料, 在632.8 nm波长下折射率为1.457, 通过对SiO2层的刻蚀, 产生深度为346 nm的沟槽进行相位π/2的调制. 仿真采用PML边界条件、TFSF光源, 仿真区域网格划分为20 nm × 20 nm × 20 nm (x, y, z).

      首先以聚焦光斑中强度最高的点作为焦点中心, 定量统计聚焦光斑的横向和轴向尺寸(FWHMx和FWHMz), 如图10所示. 随着入射激光倾斜角度的逐步增加, SOL3聚焦光斑的横向和轴向尺寸波动较大, 光斑的横向尺寸大致呈增大趋势, 轴向尺寸则呈现先增大后减小的变化趋势, 其中光斑横向FWHM最大展宽了52%, 而光斑轴向尺寸波动更为明显, 最大展宽达到了94%. 对于小数值孔径的SOL4, 光斑的横向和轴向FWHM变化则呈现出一定的规律: 当倾斜角度处于18°以内时, 光斑的FWHMx保持着较高的稳定性, 倾斜角度达到18°时聚焦光斑也仅有47 nm的展宽, 光斑的FWHMz则随着倾斜角度的增加逐步减小且整体呈现一定的线性关系; 当倾斜角度超过18°后, 光斑的横、轴向尺寸均开始出现明显波动, 并且在42°的倾斜角度下光斑的横向尺寸出现急剧变化, SOL的聚焦能力退化严重.

      图  10  倾斜照明SOL的聚焦光斑尺寸 (a) FWHMx ; (b) FWHMz

      Figure 10.  Focal spot size of SOL under oblique illumination: (a) FWHMx ; (b) FWHMz .

      其次, 聚焦光斑强度也随着倾斜角度的增大而出现显著下降, 如图11所示. 对于大数值孔径的SOL3, 倾斜角度在0°—9°的变化过程中, 光斑强度急剧衰减了61%, 而在倾斜角度达到18°以后光斑的强度衰减逐渐趋于平缓, 总体处于较低水平. 对于小数值孔径的SOL4, 与前述光斑尺寸一样, 光斑强度在18°以内的倾斜入射角度下能够保持较高稳定性, 最高仅衰减了11%, 而在18°以后光斑强度变化开始变得较为明显, 到42°的倾斜角度时强度衰减已经达到了69%.

      图  11  不同照明角度下SOL的聚焦光斑强度

      Figure 11.  Focal intensity of SOL varying with oblique illumination angle.

    • 对聚焦光斑的位置进行了定量统计分析, 发现在倾斜角度$\omega $小于18°范围内SOL3的聚焦光斑倾斜角度位置$\omega '$与入射激光倾斜角度$\omega $存在较好的线性关系, 如图12(a)所示, 经最小二乘拟合后该线性斜率为2.03, 拟合均方根误差为1.077. 进一步对SOL3的出射、入射角度的正切之比$\varGamma \!=\! \dfrac{\tan \omega '}{\tan \omega}$进行计算, 发现该比值随着倾斜角度的增加而逐渐减小, 如图12(b)所示, 而$\varGamma $表示图7所示成像系统的垂轴放大率的调制系数, 该调制系数随入射角度变化则意味着SOL3不同的垂轴位置具有不同的放大倍数$\beta $, 可利用下式表示:

      图  12  聚焦光斑角位置变化 (a) $\omega $-$\omega '$关系曲线; (b)$\omega $-$\varGamma $关系曲线

      Figure 12.  Angle position of focusing spot: (a) $\omega $-$\omega '$; (b) $\omega $-$\varGamma $.

      $\beta = \dfrac{{{f_{{\rm{sol}}}}\tan \omega '}}{{{f_{\rm{l}}}\tan \omega }} = \varGamma \dfrac{{{f_{{\rm{sol}}}}}}{{{f_{\rm{l}}}}}.$

      因此像面会产生畸变. 对于SOL3, 当倾斜角度$\omega $小于18°时, 利用正切倍角公式, 得到$\varGamma \approx 2/(1 - {\tan ^2}\omega )$, 即可以通过入射角估算出垂轴放大倍数$\beta $. 对于小数值孔径的SOL4, 聚焦光斑的角度位置$\omega '$与入射角$\omega $之间的线性关系更为明显, 倾斜角度$\omega $在0°—27°的角度范围内, 最小二乘拟合的线性斜率为1.15, 拟合均方根误差为0.142 (对应微小畸变), 比值$\varGamma $保持良好的一致性, 显示出远比SOL3好的规律性, 根据(2)式, 此时具有类似透镜成像系统, 垂轴放大倍数为常数(无畸变).

      利用(2)式, 结合折射透镜光学系统中相对畸变的定义[25], 以小角度(本文可取3°)对应的实际垂轴放大倍率为基准, 则相对畸变

      $q' = \frac{{{\beta _{\text{实}}} - {\beta _{\text{理}}}}}{{{\beta _{\text{理}}}}} \times 100\% .$

      显然由图12(b)可知, $q' < 0$(负畸变), 对应桶形畸变.

      由聚焦光斑在横向和轴向位置的变化可以得到SOL聚焦产生的场曲和畸变(针对无限远物体成像), 轴向和横向位置分别用XfZf表示, 类似于单色几何像差折射透镜系统的分析, 以激光垂直入射时焦点所在的平面为基准面, SOL3和SOL4XZ平面(子午平面)的场曲曲线如图13所示, 随着入射激光倾斜角度的增大, 场曲逐渐增大, 在倾斜角度分别达到18°和27°过后, 两者的场曲开始出现明显波动, 场曲的存在意味着聚焦像点(或离轴焦点)的轴向位置会发生明显的移动.

      图  13  SOL的场曲 (a) SOL3; (b) SOL4

      Figure 13.  Field curvature of SOL: (a) SOL3; (b) SOL4.

      根据上述两种数值孔径SOL的仿真结果和分析, 可以发现虽然在0°—42°的倾斜角度范围内两者均能保持聚焦功能, 但聚焦性能完全不同; 大数值孔径的SOL3对照明光的倾斜角度比较敏感, 即便在较小的倾斜角度, 聚焦光斑的横、纵向尺寸、光强等都会出现明显的波动和降低, 小数值孔径的SOL4则在一定的倾斜角度内表现出了良好的稳定性, 在18°的范围内不论是光斑的横纵向尺寸、强度、角度位置还是场曲和畸变等像差都表现出了良好的稳定性.

    • 本文基于VAS理论设计了多组典型的SOL, 重点采用FDTD严格电磁仿真计算分析激光照明条件(振幅分布、倾斜角度等)对SOL聚焦的影响. 仿真和分析结果表明, 在照明高斯光束的束腰半径大于SOL半径(${w_0} \geqslant D$)时, 振幅型和相位型SOL的聚焦性能均接近理想情形, 而当束腰半径小于SOL半径(${w_0} < D/2$)时聚焦光斑迅速展宽, 且相位型SOL还出现了显著焦移和轴向旁瓣增强. 大数值孔径振幅型SOL对照明光束的倾斜敏感, 仅允许几度的倾角; 大数值孔径相位型SOL允许较为宽广的倾斜角度范围, 但聚焦光斑尺寸和强度对入射激光的倾斜角度极为敏感, 即便几度的倾角下聚焦光斑强度也会产生迅速降低, 且存在较大的场曲和负畸变. 相比之下, 小数值孔径相位型SOL聚焦光斑的尺寸、强度、角度位置和场曲、畸变等能够保持较高的稳定性, 这表明小数值孔径SOL对照明激光倾斜误差的容忍性更高, 其安装使用条件相比大数值孔径SOL较为宽松, 同时也可以用于一些需要倾斜照明的场合下.

参考文献 (25)

目录

    /

    返回文章
    返回