高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

宽光谱穆勒矩阵椭偏技术的拓展应用

方慧雯 杨锦宏 章美娟 贺胜男 汪卫华

引用本文:
Citation:

宽光谱穆勒矩阵椭偏技术的拓展应用

    通讯作者: 杨锦宏, 271778116@qq.com
  • 基金项目:

    国家重点研发计划资助项目 2019YFE03080400

  • 中图分类号: O436.3

Extended application of wide spectrum Mueller matrix ellipsometry

    Corresponding author: YANG Jinhong, 271778116@qq.com ;
  • CLC number: O436.3

  • 摘要: 为了解决现有光谱穆勒矩阵椭偏检测仪的测量功能固定、不能根据实验需求满足更多物理量测量的问题, 提出了结合多物理场光学仿真与宽光谱穆勒矩阵椭偏测量数据, 以实现更多参数测量的新方案。以膜厚测量为实例, 通过比较硅基底上不同厚度、不同入射角下, 二氧化硅薄膜仿真值与椭偏仪实验测量值所得穆勒矩阵的匹配度, 得到均方误差值相对最小时的二氧化硅薄膜厚度值。结果表明, 所得膜厚结果与实际值符合得较好。该研究验证了光谱穆勒矩阵椭偏测量与仿真模拟相结合方法的可行性和有效性。
  • 图 1  光在单层薄膜系统的反射与折射示意图

    Figure 1.  Schematic diagram of reflection and refraction of light single-layer film system

    图 2  椭偏仪结构示意图

    Figure 2.  Structure diagram of ellipsometer

    图 3  Si基底上100.0 nm SiO2薄膜实验值及两种SiO2薄膜仿真穆勒矩阵参数对比图

    Figure 3.  Comparison of experimental values of 100.0 nm SiO2 thin film on Si substrate and parameters of two types of SiO2 thin film simulation Mueller matrix

    图 4  入射角55°下M12M33M34实验与仿真数据对比图

    Figure 4.  Comparison of M12, M33 and M34 experimental and simulation data at 55° incident angle

    图 5  入射角65°下M12M33M34实验与仿真数据对比图

    Figure 5.  Comparison of M12, M33 and M34 experimental and simulation data at 65°incident angle

    图 6  入射角75°下M12M33M34的实验与仿真数据对比图

    Figure 6.  Comparison of M12, M33 and M34 experimental and simulation data at 75° incident angle

    图 7  不同入射角下121.8 nm SiO2薄膜的M12M33M34实验数据图

    Figure 7.  M12, M33 and M34 experimental data of 121.8 nm SiO2 film at different incident angles

    图 8  不同入射角下121.8 nm SiO2薄膜的M12M33M34仿真数据图

    Figure 8.  M12, M33 and M34 simulation data of 121.8 nm SiO2 film at different incident angles

    图 9  实验测量Si基底上15.0 nm SiO2与仿真的16.5 nm SiO2M12M33M34对比图

    Figure 9.  Comparison between M12, M33 and M34 of 15.0 nm SiO2 on Si substrate and the simulated 16.5 nm SiO2

    图 10  实验测量Si基底上100.0 nm SiO2与仿真的105.2 nm SiO2M12M33M34对比图

    Figure 10.  Comparison between M12, M33 and M34 of 100.0 nm SiO2 on Si substrate and the simulated 105.2 nm SiO2

    图 11  实验测量Si基底上340.0 nm SiO2与仿真的340.0 nm SiO2M12M33M34对比图

    Figure 11.  Comparison between M12, M33 and M34 of 340.0 nm SiO2 on Si substrate and the simulated 340.0 nm SiO2

    表 1  不同入射角下M12M33M34的实验与仿真数据的均方差值

    Table 1.  eMSE of experimental and simulation data of M12, M33 and M34 at different incident angles

    incident angle M12 M33 M34
    55° 1.61×10-4 4.45×10-5 9.32×10-5
    65° 8.91×10-4 1.47×10-4 1.79×10-5
    75° 1.11×10-4 7.54×10-4 4.14×10-5
    下载: 导出CSV

    表 2  Si基底上15.0 nm SiO2薄膜实验测量与不同厚度仿真的M12M33M34eMSE

    Table 2.  Experimental measurements of eMSE values of 15.0 nm SiO2 films on Si substrate with different thickness simulations of M12, M33, M34

    thickness M12 M33 M34
    15.0 nm SiO2 2.62×10-5 2.11×10-4 4.25×10-4
    16.4 nm SiO2 3.26×10-5 7.29×10-5 3.79×10-5
    16.5 nm SiO2 3.31×10-5 7.07×10-5 4.26×10-5
    16.6 nm SiO2 3.41×10-5 6.46×10-5 5.16×10-5
    下载: 导出CSV

    表 3  Si基底上100.0 nm SiO2薄膜实验测量与不同厚度仿真的M12M33M34eMSE

    Table 3.  Experimental measurements of eMSE values of 100.0 nm SiO2 films on Si substrate with different thickness simulations of M12, M33, M34

    thickness M12 M33 M34
    100.0 nm SiO2 1.16×10-2 2.12×10-3 1.81×10-2
    105.1 nm SiO2 2.89×10-4 1.10×10-4 3.96×10-5
    105.2 nm SiO2 2.83×10-4 9.87×10-5 4.59×10-5
    105.3 nm SiO2 2.85×10-5 8.80×10-5 6.62×10-5
    下载: 导出CSV

    表 4  Si基底上340.0 nm SiO2薄膜实验测量与不同厚度仿真的M12M33M34eMSE

    Table 4.  Experimental measurements of eMSE values of 340.0 nm SiO2 films on Si substrate with different thickness simulations of M12, M33, M34

    thickness M12 M33 M34
    340.1 nm SiO2 3.23×10-4 5.62×10-4 1.40×10-3
    340.0 nm SiO2 3.13×10-4 5.60×10-4 1.40×10-3
    339.9 nm SiO2 3.20×10-4 5.66×10-4 1.40×10-3
    340.6 nm SiO2 3.67×10-4 5.59×10-4 1.46×10-3
    下载: 导出CSV
  • [1] 吴曼玉, 黄水平. Zr基薄膜金属玻璃的制备及其光学常数测定[J]. 光子学报, 2020, 49(10): 1031001.

    WU M Y, HUANG Sh P. Preparation of Zr-based thin-film metallic glass and determination of optical constants [J]. Acta Photonica Sinica, 2020, 49(10): 1031001(in Chinese).
    [2]

    LI J, HWU J J, LIU Y, et al. Mueller matrix measurement of asymmetric gratings[J]. Journal of Micro/Nanolithography, MEMS, and MOEMS, 2010, 9(4): 041305. doi: 10.1117/1.3514708
    [3]

    CHEN X, ZHANG C, LIU S, et al. Mueller matrix ellipsometric detection of profile asymmetry in nanoimprinted grating structures[J]. Journal of Applied Physics, 2014, 116(19): 194305. doi: 10.1063/1.4902154
    [4]

    KU S, MAHATO K K, MAZUMDER N. Polarization-resolved Stokes-Mueller imaging: A review of technology and applications[J]. Lasers in Medical Science, 2019, 34(7): 1283-1293. doi: 10.1007/s10103-019-02752-1
    [5]

    NOVIKOVA T, de MARTINO A, OSSIKOVSKI R, et al. Metrological applications of Mueller polarimetry in conical diffraction for overlay characterization in microelectronics[J]. The European Physical Journal-Applied Physics, 2005, 31(1): 63-69. doi: 10.1051/epjap:2005034
    [6] 陈修国, 刘世元, 张传维, 等. 基于Mueller矩阵椭偏仪的纳米压印模板与光刻胶光栅结构准确测量[J]. 物理学报, 2014, 63(18): 180701.

    CHEN X G, LIU Sh Y, ZHANG Ch W, et al. Accurate measurement of nanoimprint template and photoresist grating structure based on Mueller matrix ellipsometer [J]. Acta Physica Sinica, 2014, 63(18): 180701(in Chinese).
    [7] 陈超. 高分辨层析成像穆勒矩阵椭偏仪研制与应用研究[D]. 武汉: 华中科技大学, 2021: 11-14.

    CHEN Ch. Research on the development and application of high-resolution tomography Mueller matrix ellipsometer [D]. Wuhan: Huazhong University of Science and Technology, 2021: 11-14(in Chinese).
    [8] 郑州, 侯俊峰. 单波长中红外穆勒矩阵椭偏仪的设计、定标与测试[J]. 光学学报, 2022, 42(18): 1812004.

    ZHENG Zh, HOU J F. Design, calibration and test of single-wavelength mid-infrared Mueller matrix ellipsometer [J]. Journal of Optics, 2022, 42 (18): 1812004(in Chinese).
    [9] 李伟奇. 高精度宽光谱穆勒矩阵椭偏仪研制与应用研究[D]. 武汉: 华中科技大学, 2016: 10-12.

    LI W Q. Research on the development and application of high-precision wide-spectrum Mueller matrix ellipsometer [D]. Wuhan: Huazhong University of Science and Technology, 2016: 10-12(in Chinese).
    [10]

    FUJIWARA H. Spectroscopic ellipsometry principles and applications[M]. New York, USA: John Wiley & Sons Inc., 2007: 23-25.
    [11] 石玉君. 超导薄膜正常态下光电性质的椭偏研究[D]. 济南: 山东大学, 2020: 44-45.

    SHI Y J. Ellipsometry study on photoelectric properties of superconducting thin films in normal state [D]. Ji'nan: Shandong University, 2020: 44-45(in Chinese).
    [12]

    AZZAM R M A, BASHARA N M. Ellipsometry and polarized light [M]. Beijing: Science Press, 1986: 103-113.
    [13]

    GARCIA-CAUREL E, de MARTINO A, GASTON J P, et al. Application of spectroscopic ellipsometry and Mueller ellipsometry to optical characterization[J]. Applied Spectroscopy, 2013, 67(1): 1-21. doi: 10.1366/12-06883
    [14]

    ZHAO Y, LI Y, HE W, et al. Polarization scattering imaging experiment based on Mueller matrix[J]. Optics Communications, 2021, 490: 126892. doi: 10.1016/j.optcom.2021.126892
    [15]

    OSSIKOVSKI R, BUGAMI B A L, GARCIA-CAUREL E, et al. Polarizer calibration method for Mueller matrix polarimeters[J]. Applied Optics, 2020, 59(33): 10389-10395. doi: 10.1364/AO.409799
    [16]

    YADAV H N S, KUMAR M, KUMAR A, et al. COMSOL simulation of microwave plasma polishing on different surfaces[J]. Materials Today, 2021, 45: 4803-4809.
    [17]

    CHEN C, AN I, FERREIRA G M. Multichannel Mueller matrix ellipsometer based on the dual rotating compensator principle [J]. Thin Solid Films, 2004, 455(1): 14-23.
    [18]

    FENG R, QIU J, CAO Y Y. Omnidirectional and polarization insensitive nearly perfect absorber in one dimensional meta-structure[J]. Applied Physics Letters, 2014, 105(18): 207402.
    [19]

    MIAO Z, TANG Y, WEI K, et al. Random error analysis of normalized Fourier coefficient in dual-rotating compensator Mueller matrix ellipsometer[J]. Measurement Science and Technology, 2021, 32(12): 125602. doi: 10.1088/1361-6501/ac1a80
  • [1] 栗开婷吴福全刘前彭敦云李丁丁 . 三元波片复合退偏器退偏性能的穆勒矩阵分析. 激光技术, 2015, 39(4): 549-551. doi: 10.7510/jgjs.issn.1001-3806.2015.04.027
    [2] 陈国赵长明纪荣祎李鲲罗雄白羽 . 基于ZEMAX的半导体激光准直仿真设计. 激光技术, 2012, 36(3): 318-321,325.
    [3] 王顺程高峰李强杨建昌闫宗群 . 激光发射天线仿真平台设计. 激光技术, 2019, 43(1): 131-136. doi: 10.7510/jgjs.issn.1001-3806.2019.01.026
    [4] 杨仁付朱孝立陈军宁 . Ronchi光栅误差的数值仿真. 激光技术, 2012, 36(1): 37-41. doi: 10.3969/j.issn.1001-3806.2012.01.011
    [5] 陈煌飞陈勇李怡勇王志 . 1.06μm激光的大气传输仿真研究. 激光技术, 2014, 38(2): 266-269. doi: 10.7510/jgjs.issn.1001-3806.2014.02.025
    [6] 陈永庆张陈涛张建寰 . 激光化学气相沉积石墨烯的基底温度场仿真. 激光技术, 2015, 39(5): 648-653. doi: 10.7510/jgjs.issn.1001-3806.2015.05.013
    [7] 刘亚萍彭绪金赵刚陶刚高恒白杨唐伟罗杰平 . 重频50Hz风冷YAG固体激光器热设计及仿真分析. 激光技术, 2021, 45(6): 735-739. doi: 10.7510/jgjs.issn.1001-3806.2021.06.010
    [8] 葛卫龙张晓晖徐翔 . 反射式LED二次光学系统研究. 激光技术, 2011, 35(6): 857-860,864. doi: 10.3969/j.issn.1001-3806.2011.06.03
    [9] 谢洪波祝晓晨吕二阳李勇 . 一种可用于激光周视探测的发射光学系统设计. 激光技术, 2013, 37(2): 191-194. doi: 10.7510/jgjs.issn.1001-3806.2013.02.013
    [10] 买买提艾力巴克艾尔肯扎克尔沙依甫加马力达吾来提 . 激光-等离子体相互作用过程中光子加速的研究. 激光技术, 2011, 35(2): 282-284. doi: 10.3969/j.issn.1001-3806.2011.02.038
    [11] 许国侯春萍沈丽丽 . 狭缝光栅立体显示光栅倾斜角度理论分析. 激光技术, 2016, 40(3): 388-391. doi: 10.7510/jgjs.issn.1001-3806.2016.03.018
    [12] 屈直张伯虎 . 一种改进的小波阈值算法在激光侦听中的应用. 激光技术, 2014, 38(2): 218-224. doi: 10.7510/jgjs.issn.1001-3806.2014.02.016
    [13] 赵麒白忠臣周骅陆安江张正平刘桥 . 拉盖尔-高斯光束作用下熔石英温度及应力研究. 激光技术, 2018, 42(1): 121-126. doi: 10.7510/jgjs.issn.1001-3806.2018.01.024
    [14] 乔迁赵虎谢兰强 . 全固态风冷式激光器散热设计及优化. 激光技术, 2024, 48(2): 235-239. doi: 10.7510/jgjs.issn.1001-3806.2024.02.014
    [15] 王可宁刘允雷陈海滨郭子龙 . 移频延时自外差法的DFB激光器线宽测量. 激光技术, 2018, 42(5): 633-637. doi: 10.7510/jgjs.issn.1001-3806.2018.05.010
    [16] 张俊张为国 . Micro-LED蓝宝石衬底AlN上GaN激光剥离研究. 激光技术, 2023, 47(1): 25-31. doi: 10.7510/jgjs.issn.1001-3806.2023.01.004
    [17] 任亚辉林菊平童勇张强刘攀夏春燕王国帅 . 高能脉冲氙灯驱动电路特性分析. 激光技术, 2022, 46(3): 397-401. doi: 10.7510/jgjs.issn.1001-3806.2022.03.016
    [18] 张新陆程元丽王骐 . X射线激光在柱状等离子体中的传播与放大. 激光技术, 2002, 26(3): 221-222.
    [19] 张新陆王骐程元丽 . 柱状等离子体侧漏光线对光强分布的影响. 激光技术, 2003, 27(2): 122-123,125.
    [20] 张新陆王骐程元丽 . 毛细管放电软X射线激光传播的解析研究. 激光技术, 2002, 26(5): 392-393,397.
  • 加载中
图(11) / 表(4)
计量
  • 文章访问数:  2369
  • HTML全文浏览量:  1618
  • PDF下载量:  15
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-12-07
  • 录用日期:  2023-02-24
  • 刊出日期:  2024-01-25

宽光谱穆勒矩阵椭偏技术的拓展应用

    通讯作者: 杨锦宏, 271778116@qq.com
  • 安徽大学 物质科学与信息技术研究院, 合肥 230601, 中国
基金项目:  国家重点研发计划资助项目 2019YFE03080400

摘要: 为了解决现有光谱穆勒矩阵椭偏检测仪的测量功能固定、不能根据实验需求满足更多物理量测量的问题, 提出了结合多物理场光学仿真与宽光谱穆勒矩阵椭偏测量数据, 以实现更多参数测量的新方案。以膜厚测量为实例, 通过比较硅基底上不同厚度、不同入射角下, 二氧化硅薄膜仿真值与椭偏仪实验测量值所得穆勒矩阵的匹配度, 得到均方误差值相对最小时的二氧化硅薄膜厚度值。结果表明, 所得膜厚结果与实际值符合得较好。该研究验证了光谱穆勒矩阵椭偏测量与仿真模拟相结合方法的可行性和有效性。

English Abstract

    • 光谱椭偏技术是一种无损检测技术,通过测量的数据进行反向拟合,得到薄膜的光学常数和厚度[1]。传统的光谱椭偏仪(spectroscopic ellipsometer, SE)只能获得相位差和振幅比这两个测量参数[2]。穆勒矩阵椭偏仪(Mueller matrix ellipsometer,MME)可以获得样品表面的16个穆勒矩阵参数[3],能够对样品的光学性质和结构信息做出一定表述[4]。光谱椭偏技术虽被发现和使用很久,但得益于精密加工、高分辨成像等技术,光谱椭偏在现代科学研究和工业领域仍然应用广泛。光谱椭偏检测技术也在不断发展,如法国国家科学研究中心界面与薄膜物理实验室和设备供应商Nanometrics公司都使用MME对光栅结构的特点进行测量[5]。华中科技大学学者利用MME研究了纳米压印过程中的退偏效应[6],而且还设计了具有高分辨率层析成像性能的MME[7]。ZHENG等人设计搭建了单波长中红外MME, 并且误差测量在0.02[8]。国际国内多个公司研发生产出不同型号宽光谱穆勒矩阵椭偏仪,可进行镀膜工艺过程的检测,被广泛运用到集成电路、太阳能光伏、化学等领域[9]。目前椭偏仪已成为表征材料光学常数及薄膜厚度的标准工具[10]

      现有的SE功能和可测参数已经固定,无法用于扩展测量其它的物理参数。但光谱椭偏仪测量的相关偏振参数数据非常丰富,包含了多种物理量信息,通过采用适当的数据处理,可望拓展光谱椭偏仪的应用范围,获得更多的偏振数据和薄膜参数,如薄膜的反射率、透射率, 并且能够根据实验需求测量多物理场情况下薄膜的光学常数,甚至测量结构复杂薄膜(例如光栅)的光学常数。本文作者将仿真软件与椭偏测量相结合,以膜厚测量为实例,通过比较仿真软件在不同入射角、不同厚度的Si基底上的SiO2薄膜数据与椭偏测量的实验数据两者之间的拟合程度,验证该方法的可行性。

    • 以SiO2薄膜为研究对象,结构示意图如图 1所示。图中,n1n2n3分别为环境/薄膜/基底这三者的折射率; 膜厚为d; θ1为环境与薄膜界面的入射角,θ2为经过薄膜发生折射时的折射角(也是薄膜与基底界面的入射角),θ3为光线经过基底产生折射的折射角。

      图  1  光在单层薄膜系统的反射与折射示意图

      Figure 1.  Schematic diagram of reflection and refraction of light single-layer film system

      光入射在物体表面时会发生反射与折射, 将光波电矢量划分为振动方向平行于入射光面的p光与振动方向垂直于入射光面的s光。由图 1可知,入射光在薄膜内多次反射和折射,回到空气介质中的每束光波的振幅和相位都不同,其振幅的大小由菲涅耳公式决定,相邻两束反射光之间的相位差由光程差乘以折射率决定[11],即:

      $ 2 \delta=\frac{4 {\rm{\mathsf{π}}}}{\lambda} d n_2 \cos \theta_2 $

      (1)

      式中: δ是相邻两束反射光之间的相位差; λ为波长。

      因此根据反射系数的定义,p光和s光的总反射率为:

      $ r_p=\frac{r_{1, p}+r_{2, p} \mathrm{e}^{-\mathrm{i} 2 \delta}}{1+r_{1, p} \mathrm{e}^{-\mathrm{i} 2 \delta}} $

      (2)

      $ r_s=\frac{r_{1, s}+r_{2, s} \mathrm{e}^{-\mathrm{i} 2 \delta}}{1+r_{1, s} \mathrm{e}^{-\mathrm{i} 2 \delta}} $

      (3)

      式中: r1, pr1, sp光和s光在空气和薄膜界面的反射系数; r2, pr2, sp光和s光在薄膜与衬底界面的反射系数。

      引入参量tanφ描述p光和s光振幅比Ep/Es, Δ表示相位差[12],定义如下式:

      $ \begin{aligned} & \rho=\frac{r_p}{r_s}=\tan \varphi \cdot \exp (\mathrm{i} \varDelta)= \\ & f\left(\lambda, \theta_1, n_1, n_2, k_2, n_3, k_3, d\right) \end{aligned} $

      (4)

      式中: k2为薄膜的消光系数; k3为基底的消光系数。由式(4)可见,椭偏参数φ和Δ可以看作f(·)的函数,与薄膜和衬底的负折射率、厚度有关,即可用椭偏参数反应薄膜的光学常数和膜厚等信息[13]。但通过与仿真软件相结合能够得到更多的偏振信息,如薄膜的反射率、透射率,并且能够根据实验需求测量多物理场情况下薄膜的光学常数,甚至测量结构复杂薄膜(例如光栅)的光学常数。

    • 斯托克斯矢量(S0 S1 S2 S3)T是描述光波偏振态的简单数学表达形式[14]。入射光与反射光的偏振态存在联系,如下:

      $ \boldsymbol{S^{\prime}=M S} $

      (5)

      式中: S为入射光的斯托克斯矢量;S′为反射光的斯托克斯矢量;M即为穆勒矩阵[15]

      $ \left[\begin{array}{l} S_0{ }^{\prime} \\ S_1{ }^{\prime} \\ S_2{ }^{\prime} \\ S_3{ }^{\prime} \end{array}\right]=\left[\begin{array}{llll} M_{11} & M_{12} & M_{13} & M_{14} \\ M_{21} & M_{22} & M_{23} & M_{24} \\ M_{31} & M_{32} & M_{33} & M_{34} \\ M_{41} & M_{42} & M_{43} & M_{44} \end{array}\right]\left[\begin{array}{l} S_0 \\ S_1 \\ S_2 \\ S_3 \end{array}\right] $

      (6)
    • 多物理场光学仿真可以在单个域模型中容纳多个偏微分方程并进行多物理场的耦合[16]。本文中仿真使用的是几何光学模块,不仅能对光线轨迹追踪还能通过边界条件的改变,来仿真不同情况。为了更好地模拟椭偏仪的测量过程,在软件中构建了椭偏仪内部结构模型,如图 2所示。从光源发射宽光谱光经过起偏器变成偏振光;然后再经过第一旋转补偿器后打在样品表面,入射光经过样品反射后通过第二旋转补偿器、检偏器,最后到达探测器。根据实际样品的情况对于模型样品进行设置,先设置基底的材料接着在基底上表面添加薄膜材料,再不断地改变薄膜厚度,直到最后与实际测量曲线拟合成功。本文中实例主要测量膜厚,故在仿真中以膜厚为变量,测量不同的物理参数,可以相应地设其为变量。在建模时所使用的入射光源、光路、样品、测量环境等参数与椭偏仪实际测量时一致。

      图  2  椭偏仪结构示意图

      Figure 2.  Structure diagram of ellipsometer

      使用仿真软件建模后,导出不同入射光宽光谱的斯托克斯矢量数据,再根据式(6),通过4组入射光和对应反射光宽光谱的斯托克斯矢量计算样品的宽光谱穆勒矩阵。其中4组入射光初始斯托克斯矢量为:

      $ \left\{\begin{array}{l} \boldsymbol{S}_{0^{\circ}}=I_0\left[\begin{array}{llll} 1 & 1 & 0 & 0 \end{array}\right]^{\mathrm{T}} \\ \boldsymbol{S}_{45^{\circ}}=I_0\left[\begin{array}{llll} 1 & 0 & 1 & 0 \end{array}\right]^{\mathrm{T}} \\ \boldsymbol{S}_{90^{\circ}}=I_0\left[\begin{array}{llll} 1 & -1 & 0 & 0 \end{array}\right]^{\mathrm{T}} \\ \boldsymbol{S}_{\mathrm{c}}=I_0\left[\begin{array}{llll} 1 & 0 & 0 & 1 \end{array}\right]^{\mathrm{T}} \end{array}\right. $

      (7)

      式中: S代表的是水平偏振光; S45°代表的是45°线偏振光; S90°代表的是垂直偏振光; Sc代表的是圆偏振光; I0为光强。

    • 实验中使用MP-L穆勒椭圆偏振仪测量Si基底上不同厚度SiO2薄膜样品在不同入射角下的穆勒矩阵。由于SiO2薄膜是各向同性材料,故穆勒矩阵样式如下[17]

      $ \boldsymbol{M}=\left[\begin{array}{cccc} 1 & -N & 0 & 0 \\ -N & 1 & 0 & 0 \\ 0 & 0 & C & D \\ 0 & 0 & -D & C \end{array}\right] $

      (8)

      式中: 参数NCD的表达式如下:

      $ \left\{\begin{array}{l} N=\cos (2 \varphi) \\ C=\sin (2 \varphi) \cos \varDelta \\ D=\sin (2 \varphi) \sin \varDelta \end{array}\right. $

      (9)

      式中: φ与Δ是式(4)中的椭偏参数; 由式(8)可知,M12M21数值相等,M33M44数值相等,M34M43互为相反数,斜对角的8个参数都为0;且NCD不独立,三者关系如下[18]

      $ N^2+C^2+D^2=1 $

      (10)

      本文中用均方误差(mean squared error, MSE)eMSE来判断实验值与仿真值匹配程度的好坏,均方误差数值越小,则两者之间拟合的越好[19]。其定义式如下:

      $ e_{\mathrm{MSE}}=\frac{1}{t} \sum\limits_{i=1}^t\left(Q_i-q_i\right)^2 $

      (11)

      式中: t是数据个数; Qi是实验值; qi是仿真值。eMSE的值越小,说明实验值与仿真值匹配得越好。通过式(11)进行计算得到eMSE,其最小时所得的厚度即为仿真软件确定的薄膜厚度。

      以Si基底上100.0 nm SiO2薄膜的测量为例,实验中入射光入射角为60°、波长范围380 nm~1000 nm,光路如图 2所示,测量得到样品的穆勒矩阵。然后利用仿真软件仿真Si基底上SiO2薄膜厚度在100.0 nm附近的反射过程,计算得到不同膜厚的穆勒矩阵。由于SiO2薄膜是各向同性薄膜,采用式(8)的穆勒矩阵,因此只需要比较实验与仿真穆勒矩阵参数的M12M33M34,如图 3所示。

      图  3  Si基底上100.0 nm SiO2薄膜实验值及两种SiO2薄膜仿真穆勒矩阵参数对比图

      Figure 3.  Comparison of experimental values of 100.0 nm SiO2 thin film on Si substrate and parameters of two types of SiO2 thin film simulation Mueller matrix

      图 3中横坐标为波长,纵坐标是归一化后的穆勒矩阵参数值。红色实线是实验测得数值,黑色虚线是105.2 nm SiO2薄膜仿真值,蓝色点划线是100.0 nm SiO2薄膜仿真值。由图 3可见, 实验数据与厚度为105.2 nm SiO2薄膜仿真数据拟合较好,均方误差值也很小,并且满足各向同性材料的规律,则105.2 nm为仿真确定的薄膜厚度。图 3中只列出100.0 nm与105.2 nm仿真结果,实际上做了多个膜厚仿真和计算,因为厚度为105.2 nm时仿真值与实验值的eMSE相对最小,所以为了简洁, 图中只列出了两个仿真值进行比较。以下的实验与仿真对比过程中,也是直接列出eMSE相对最小的仿真数据。

    • 以上对入射角为60°时实验与仿真结果进行了对比,下面比较实验与仿真的M12M33M34在不同入射角时匹配情况。实验测量的样件标定厚度是Si基底上121.8 nm SiO2薄膜,椭偏仪测量该薄膜厚度为122.1 nm。入射角分别为55°, 65°, 75°,其它参数同图 3。通过与图 3相同的方法对比不同入射角时实验与仿真的eMSE,如图 4图 5图 6所示。120.5 nm为SiO2薄膜厚度最优结果。

      图  4  入射角55°下M12M33M34实验与仿真数据对比图

      Figure 4.  Comparison of M12, M33 and M34 experimental and simulation data at 55° incident angle

      图  5  入射角65°下M12M33M34实验与仿真数据对比图

      Figure 5.  Comparison of M12, M33 and M34 experimental and simulation data at 65°incident angle

      图  6  入射角75°下M12M33M34的实验与仿真数据对比图

      Figure 6.  Comparison of M12, M33 and M34 experimental and simulation data at 75° incident angle

      图 4图 5图 6中黑色虚线是仿真值,红色实线是实验值,3个图中SiO2薄膜厚度一样,只是在不同入射角下的M12M33M34的实验与仿真数据对比。表 1中将不同入射角下的M12M33M34的实验与仿真数据的均方误差值进行对比,随着入射角的增大,M33的拟合效果变差,这是由于过大的入射角会导致实验误差变大。

      表 1  不同入射角下M12M33M34的实验与仿真数据的均方差值

      Table 1.  eMSE of experimental and simulation data of M12, M33 and M34 at different incident angles

      incident angle M12 M33 M34
      55° 1.61×10-4 4.45×10-5 9.32×10-5
      65° 8.91×10-4 1.47×10-4 1.79×10-5
      75° 1.11×10-4 7.54×10-4 4.14×10-5

      为了更进一步探究入射角对M12M33M34的影响,分别测量了Si基底上121.8 nm SiO2薄膜在50°、52°、55°、60°、65°、70°、75°入射角下的M12M33M34的数据,并使用仿真软件对这些情况同时进行仿真,得到的实验值与仿真值如图 7图 8所示。

      图  7  不同入射角下121.8 nm SiO2薄膜的M12M33M34实验数据图

      Figure 7.  M12, M33 and M34 experimental data of 121.8 nm SiO2 film at different incident angles

      图  8  不同入射角下121.8 nm SiO2薄膜的M12M33M34仿真数据图

      Figure 8.  M12, M33 and M34 simulation data of 121.8 nm SiO2 film at different incident angles

      图 7图 8中不同颜色表示不同入射角度,其它参数同图 3。由图 7图 8可见,不同入射角对M12M33M34影响明显,但相同的入射角实验与仿真结果相符度很好,说明在这些入射角范围内进行实验与仿真拟合效果较好。

    • 下面对Si基底上3种不同厚度的SiO2薄膜进行实验与仿真对比,其它参数同图 3。样品标示SiO2膜厚分别为15.0 nm、100.0 nm、340.0 nm,椭偏仪测量的厚度分别为15.5 nm、101.2 nm、340.6 nm。仿真结果选取eMSE达到最小时,确定的薄膜的厚度分别为16.5 nm、105.2 nm、340.0 nm,在表 2表 3表 4中列出了最佳厚度附近的仿真数据与实测厚度的穆勒矩阵参数的eMSE对比值,如图 9图 10图 11所示。

      表 2  Si基底上15.0 nm SiO2薄膜实验测量与不同厚度仿真的M12M33M34eMSE

      Table 2.  Experimental measurements of eMSE values of 15.0 nm SiO2 films on Si substrate with different thickness simulations of M12, M33, M34

      thickness M12 M33 M34
      15.0 nm SiO2 2.62×10-5 2.11×10-4 4.25×10-4
      16.4 nm SiO2 3.26×10-5 7.29×10-5 3.79×10-5
      16.5 nm SiO2 3.31×10-5 7.07×10-5 4.26×10-5
      16.6 nm SiO2 3.41×10-5 6.46×10-5 5.16×10-5

      表 3  Si基底上100.0 nm SiO2薄膜实验测量与不同厚度仿真的M12M33M34eMSE

      Table 3.  Experimental measurements of eMSE values of 100.0 nm SiO2 films on Si substrate with different thickness simulations of M12, M33, M34

      thickness M12 M33 M34
      100.0 nm SiO2 1.16×10-2 2.12×10-3 1.81×10-2
      105.1 nm SiO2 2.89×10-4 1.10×10-4 3.96×10-5
      105.2 nm SiO2 2.83×10-4 9.87×10-5 4.59×10-5
      105.3 nm SiO2 2.85×10-5 8.80×10-5 6.62×10-5

      表 4  Si基底上340.0 nm SiO2薄膜实验测量与不同厚度仿真的M12M33M34eMSE

      Table 4.  Experimental measurements of eMSE values of 340.0 nm SiO2 films on Si substrate with different thickness simulations of M12, M33, M34

      thickness M12 M33 M34
      340.1 nm SiO2 3.23×10-4 5.62×10-4 1.40×10-3
      340.0 nm SiO2 3.13×10-4 5.60×10-4 1.40×10-3
      339.9 nm SiO2 3.20×10-4 5.66×10-4 1.40×10-3
      340.6 nm SiO2 3.67×10-4 5.59×10-4 1.46×10-3

      图  9  实验测量Si基底上15.0 nm SiO2与仿真的16.5 nm SiO2M12M33M34对比图

      Figure 9.  Comparison between M12, M33 and M34 of 15.0 nm SiO2 on Si substrate and the simulated 16.5 nm SiO2

      图  10  实验测量Si基底上100.0 nm SiO2与仿真的105.2 nm SiO2M12M33M34对比图

      Figure 10.  Comparison between M12, M33 and M34 of 100.0 nm SiO2 on Si substrate and the simulated 105.2 nm SiO2

      图  11  实验测量Si基底上340.0 nm SiO2与仿真的340.0 nm SiO2M12M33M34对比图

      Figure 11.  Comparison between M12, M33 and M34 of 340.0 nm SiO2 on Si substrate and the simulated 340.0 nm SiO2

      图 9图 10图 11是在入射角同为60°时,不同厚度薄膜实验测得数据和仿真数据的对比。黑色虚线是仿真值,红色实线是实验值。由图 9~图 11可见,仿真的厚度与实验测量厚度非常接近,总体上各穆勒矩阵参数与实际测量数据的eMSE值都很小。

    • 为了拓展宽光谱椭偏仪的功能,实现更多薄膜参数的测量,提出了椭偏仪实验数据与多物理场仿真相结合实现更多物理量测量的新方法,并以膜厚测量为例进行了验证。通过对Si基底上SiO2膜厚实验测量与仿真得到的穆勒矩阵进行比对,获得了被测样品薄膜的厚度,所得膜厚结果与实际值符合较好。该方法可以实现多物理场环境中(不限于空气中)薄膜参数的椭偏测量。

参考文献 (19)

目录

    /

    返回文章
    返回