CN112378861A - 快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法 - Google Patents

快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法 Download PDF

Info

Publication number
CN112378861A
CN112378861A CN202011163917.5A CN202011163917A CN112378861A CN 112378861 A CN112378861 A CN 112378861A CN 202011163917 A CN202011163917 A CN 202011163917A CN 112378861 A CN112378861 A CN 112378861A
Authority
CN
China
Prior art keywords
frequency
phase
phase delay
snapshot
mueller matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011163917.5A
Other languages
English (en)
Other versions
CN112378861B (zh
Inventor
陈修国
王鹏
张劲松
石雅婷
刘世元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202011163917.5A priority Critical patent/CN112378861B/zh
Publication of CN112378861A publication Critical patent/CN112378861A/zh
Application granted granted Critical
Publication of CN112378861B publication Critical patent/CN112378861B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/21Polarisation-affecting properties
    • G01N21/211Ellipsometry
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明属于精密光学测量仪器系统参数校准领域,并具体公开了一种快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其包括如下步骤:S1、按预设相位延迟器厚度比搭建快照式穆勒矩阵椭偏仪,然后对三组标准样品进行测量,得到三组测量光谱;S2、分别对三组测量光谱进行频域分析得到频域信号,对频域信号进行分通道处理得到多阶频率通道,再根据相位延迟器厚度比对频率通道进行选择;S3、对选择的频率通道进行波数域分析计算,获得各频率通道对应的测量光谱三角函数展开式的实频、虚频系数,进而得到各相位延迟器的相位延迟量误差。本发明能准确校准快照式穆勒矩阵椭偏仪全测量光谱内的相位延迟量误差,适用于不同的相位延迟器厚度比。

Description

快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法
技术领域
本发明属于精密光学测量仪器系统参数校准领域,更具体地,涉及一种快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法。
背景技术
自Matthieu Dubreuil等人在文章“Snapshot Mueller matrix polarimeter bywavelength polarization coding”中提出快照式穆勒矩阵椭偏仪以来,快照式穆勒矩阵椭偏仪受到了越来越广泛的关注,其具有测量时间极短、测量过程中不受运动器件干扰等优点。快照式穆勒矩阵椭偏仪通常由光源、起偏器、四片具有一定厚度比与方位角的相位延迟器、检偏器与探测器等部分组成,如图2所示,其核心器件为四片相位延迟器,并以此实现对待测样品穆勒矩阵的调制。该仪器模型可描述为:
Sout=P22)R444)R333)MR222)R111)P11)Sin (1)
式中Sin与Sout分别表示由光源出射的光与探测器接收到的光的斯托克斯矢量;M、P1与P2分别表示待测样品、起偏器与检偏器的穆勒矩阵,α1、α2分别表示起偏器与检偏器的方位角;Ri、θi、δi(i=1,2,3,4)分别表示第i个相位延迟器的穆勒矩阵、相对起偏器的方位角以及其相位延迟量。
记各相位延迟量的基础厚度为d0,各相位延迟量厚度为di d0(di为整数, i=1,2,3,4),则其相位延迟量如式(2-a)所示,为使载波为正弦波以方便后续计算,通常在建立模型时对其相位延迟量δi做如式(2-b)所示的近似:
δi(σ)=2πΔn(σ)did0σ,i=1,2,3,4 (2-a)
Figure BDA0002745134640000011
式中f0为调制基频,σ为波数,Δn为相位延迟器所用材料的双折射率。若光源发出的光的斯托克斯矢量为[1 0 0 0]T,则由探测器接收到的光谱I(σ) 如式(3)所示:
Figure BDA0002745134640000021
式中mo为系统调制的最高频次,mo=d1+d2+d3+d4;pj为光谱做傅里叶展开后所得的各阶频次,其取值为d1、d2、d3、d4的线性组合;gj、hj分别为光谱做傅里叶展开后实频与虚频的系数,其取值为待测样品穆勒矩阵的线性组合。公式(3)给出了测量光谱的三角函数展开式,因此对测量光谱进行频域分析即可求解gj、hj,进而求解待测样品的穆勒矩阵。
然而在实际测量中,由于材料的色散效应、加工误差与热胀冷缩效应导致的相位延迟器厚度误差Δdi、受温度影响材料双折射率产生变化等多方面原因,其相位延迟量的真实值存在如式(4)所示的误差Δδi
Δδi(σ)=2πΔn′(σ)Δdiσ+2π(Δn′(σ)-Δn(σ0))did0σ,i=1,2,3,4 (4)
式中Δn'表示材料受上述影响后的真实双折射率。
相位延迟器是快照式穆勒矩阵椭偏仪完成待测样品穆勒矩阵调制的关键器件,其相位延迟量误差将会对样品穆勒矩阵测量结果产生直接影响,因此必须对相位延迟量误差进行分析校正。Matthieu Dubreuil等在文章“Systematic errors specific to asnapshot Mueller matrix polarimeter”中给出一种近似的校准方法,将全测量波段的相位延迟量误差近似为中心波数处的误差,校准得出的相位延迟量误差为不随波数变化的定值,校准结果仍存在较大误差,且与实际情况不符;Hu Dai等在文章“Theoreticalresearch of retarder phase deviation channeled Mueller matrixspectropolarimeters”中提出一种针对厚度比为1:1:5:5的校准方法,该方法校准效果较好,但这种方法由于使用了低通滤波器,因此存在一定的延迟,且难以适用于其他厚度比配置时相位延迟量的校准。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其目的在于,通过对标准样品测量光谱进行频域分析,并选择频率通道进行波数域分析计算,得到各个相位延迟器的相位延迟量误差,能够准确校准快照式穆勒矩阵椭偏仪的全测量光谱内的相位延迟量误差,且适用于不同的相位延迟器厚度比。
为实现上述目的,本发明提出了一种快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,包括如下步骤:
S1、按预设相位延迟器厚度比搭建快照式穆勒矩阵椭偏仪,通过该快照式穆勒矩阵椭偏仪分别对三组标准样品进行测量,得到三组测量光谱;
S2、分别对三组测量光谱进行频域分析得到频域信号,然后对频域信号进行分通道处理得到多阶频率通道,再根据相位延迟器厚度比对频率通道进行选择;
S3、对选择的频率通道进行波数域分析计算,获得各频率通道对应的测量光谱三角函数展开式的实频系数与虚频系数,进而得到各相位延迟器的相位延迟量误差,完成快照式穆勒矩阵椭偏仪相位延迟量误差校准。
作为进一步优选的,还包括如下步骤:
S4、选择校准后相位延迟量误差的中心波段,对其进行二阶多项式拟合,进一步校准相位延迟量误差。
作为进一步优选的,所述快照式穆勒矩阵椭偏仪中,以起偏器方位角作为参考,检偏器方位角为90°,四片相位延迟器的方位角依次分别为45°、 0°、0°、45°。
作为进一步优选的,所述三组标准样品中,第一组标准样品为方位角为0°的偏振片,第二组标准样品为方位角为0°与45°的两个偏振片组合,第三组标准样品为方位角为45°与0°的两个偏振片组合。
作为进一步优选的,所述步骤S2中,先通过傅里叶变换对测量光谱进行频域分析得到频域信号,然后按由相位延迟器基础厚度决定的基础频率 f0将频域信号划分为多阶频率通道,其中,第j阶频率通道中频率范围为(j– 1/2)f0~(j+1/2)f0
作为进一步优选的,所述相位延迟器厚度的最简整数比为d1:d2:d3: d4;所述步骤S2中,频率通道的选择原则为:对于第一组标准样品测量光谱的频域信号,先选择d1+d4阶频率通道,然后,若d1≠|d1–d4|,则选择d1阶频率通道,否则选择d4阶频率通道;对于第二组标准样品测量光谱的频域信号,选择d1+d3+d4阶频率通道;对于第三组标准样品测量光谱的频域信号,选择d1+d2+d4阶频率通道,以此完成四个频率通道的选择。
作为进一步优选的,所述步骤S3中,对选择的频率通道进行波数域分析计算具体包括:对选择的频率通道进行傅里叶逆变换,然后根据傅里叶逆变换结果的相位信息,求取各频率通道对应的测量光谱三角函数展开式的实频系数与虚频系数。
作为进一步优选的,所述步骤S3中,根据求得的实频系数与虚频系数求解相位延迟量误差的线性组合,进而通过线性运算得到四片相位延迟器的相位延迟量误差。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,主要具备以下的技术优点:
1.本发明通过对标准样品测量光谱进行频域分析,并选择频率通道进行波数域分析计算,得到各个相位延迟器的相位延迟量误差,可以简单、快速实现各相位延迟量误差的准确校准,且进一步将校准所得相位延迟量误差用于待测样品的穆勒矩阵测量,可以极大地提高穆勒矩阵的测量精度,从而较为准确地获得待测样品地穆勒矩阵。
2.本发明对标准样品和频率通道的选择原则的进行设置,能够实现全测量波段内的相位延迟量误差的准确校准,并且适用于各种不同的相位延迟器厚度比。
3.本发明对校准所得的相位延迟量误差进一步处理,取中心波段的校正结果进行二阶多项式拟合,可以排除边缘影响,获得更好的误差校准结果。
附图说明
图1为本发明实施例快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法流程图;
图2为本发明实施例快照式穆勒矩阵椭偏仪结构示意图;
图3为本发明实施例一中相位延迟器在测量波段内的双折射率;
图4中(a)~(d)为本发明实施例一中添加厚度误差后各相位延迟器的相位延迟量误差的理论值与本发明获取的校准值对比图;
图5为本发明实施例一中一个仿真测量实例的穆勒矩阵(横坐标单位σ/mm-1),其中(a)~(p)为穆勒矩阵各参数;
图6为本发明实施例二中相位延迟器在测量波段内的双折射率;
图7为本发明实施例二中仿真测量所用样品的形貌参数;
图8中(a)~(d)为本发明实施例二中添加厚度误差后各相位延迟器的相位延迟量误差的理论值与本发明获取的校准值对比图;
图9为本发明实施例二中一个仿真测量实例的穆勒矩阵(横坐标单位σ/mm-1),其中(a)~(p)为穆勒矩阵各参数。
在所有附图中,相同的附图标记用来表示相同的元件或结构,其中: 201-光源,202-起偏器,203-相位延迟器一,204-相位延迟器二,206-相位延迟器三,207-相位延迟器四,205-测量样品,208-检偏器,209-探测器。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明实施例提供的一种快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,如图1所示,包括如下步骤:
S1、根据设计方案选择相位延迟器厚度比,搭建快照式穆勒矩阵椭偏仪,并建立测量系统模型;
具体的,如图2所示,快照式穆勒矩阵椭偏仪包括依次设置的光源201,起偏器202、相位延迟器一203、相位延迟器二204、相位延迟器三206、相位延迟器四207,检偏器208和探测器209;以起偏器方位角作为参考,检偏器方位角为90°,四片相位延迟器的方位角依次分别为45°、0°、0°、 45°;
S2、通过该快照式穆勒矩阵椭偏仪分别对三组标准样品进行测量,得到三组测量光谱,测量时,测量样品205置于相位延迟器二204和相位延迟器三206之间;
具体的,所述三组标准样品中,第一组标准样品为方位角为0°的偏振片,第二组标准样品为方位角为0°与45°的两个偏振片组合,第三组标准样品为方位角为45°与0°的两个偏振片组合。
S3、分别对三组测量光谱进行频域分析得到频域信号,然后对频域信号进行分通道处理得到多阶频率通道,再根据相位延迟器厚度比对频率通道进行选择;
进一步的,先通过傅里叶变换对测量光谱进行频域分析得到频域信号,然后按由相位延迟器基础厚度决定的基础频率f0将频域信号划分为多阶频率通道,其中,第j阶频率通道中频率范围为(j–1/2)f0~(j+1/2)f0;所述相位延迟器厚度的最简整数比为d1:d2:d3:d4(d1、d2、d3、d4为整数),则频率通道的选择原则为:对于第一组标准样品测量光谱的频域信号,首先选择第一组频率d1+d4阶频率通道;接着选择第二组频率(通常d1≠ d4),若d1≠|d1–d4|则选择第二组频率为d1阶频率通道,否则选择第二组频率为d4阶频率通道;对于第二组标准样品测量光谱的频域信号,选择第三组频率为d1+d3+d4阶频率通道;对于第三组标准样品测量光谱的频域信号,选择第四组频率为d1+d2+d4阶频率通道。
S4、对选择的频率通道进行波数域分析计算,即对选择的频率通道进行傅里叶逆变换,然后根据傅里叶逆变换结果的相位信息,获得各频率通道对应的测量光谱三角函数展开式的实频系数与虚频系数。
S5、根据实频系数与虚频系数求解相位延迟量误差的线性组合,进而通过线性运算得到四片相位延迟器的相位延迟量误差。
S6、对校准所得的相位延迟量误差进一步处理,取中心波段的校正结果进行二阶多项式拟合,可以获得更好的校准结果。
下面分别以厚度比为1:1:5:5和1:4:2:9的两种测量配置为例,详细介绍本发明所述快照时穆勒矩阵椭偏仪相位延迟量误差的校准方法,并说明其在全测量波段相位延迟量校准的有效性和通用性,其能够适用于任何可以求解全穆勒矩阵的快照式穆勒矩阵椭偏仪的相位延迟量的校准。
实施例一
以相位延迟器厚度比为1:1:5:5的测量配置为例,介绍本发明所述相位延迟量校准方法,包括如下步骤:
S1、首先建立快照式穆勒矩阵椭偏仪的测量模型:以起偏器201方位角作为参考,设置检偏器208方位角α2=90°,相位延迟器一203方位角θ1=45°,相位延迟器二204方位角θ2=0°,相位延迟器三206方位角θ3=0°,相位延迟器四207方位角θ4=45°;对于厚度比为1:1:5:5的测量配置,有 d1=d2=1,d3=d4=5;则系统的模型可表述为:
Sout=P2(90)R4(45,δ4)R3(0,δ3)MR2(0,δ2)R1(45,δ1)P1(0)Sin (5)
其中,Sin与Sout分别表示由光源出射的光与探测器接收到的光的斯托克斯矢量;M、P1与P2分别表示待测样品、起偏器与检偏器的穆勒矩阵, Ri、δi(i=1,2,3,4)分别表示第i个相位延迟器的穆勒矩阵及其相位延迟量。
各个相位延迟器的理论相位延迟量δ′i(σ)为:
δ′i(σ)=δi(σ)+Δδi(σ)=2πdif0σ+Δδi(σ),i=1,2,3,4 (6)
其中,f0为调制基频,σ为波数,Δδi(σ)为相位延迟量误差。
通常光源发出的光的斯托克斯矢量Sin=[1 0 0 0]T,代入至系统模型可得测量光谱为:
Figure BDA0002745134640000081
上述测量光谱的三角函数展开式为:
Figure BDA0002745134640000082
其中,mo为系统的最高调制频率,mo=1+1+5+5=12;pj为测量光谱三角函数展开式的频率级次;gj′与hj′分别表示展开式的实频与虚频系数,二者具有如下形式:
Figure BDA0002745134640000091
其中,αk,l为常数系数,mk,l为待测样品穆勒矩阵元素。
对于本实施例,设定测量波段为λ=739.6nm~750nm(波数σ= 1333.3mm-1~1352.1mm-1),选择材料为冰洲石,其在测量波段内的双折射率如图3所示;选择相位延迟器基础厚度d0=1.46mm,对各相位延迟器添加厚度误差,分别是Δd1=0.0073mm、Δd2=–0.0146mm、Δd3=–0.0292mm、Δd4=0.0146mm。
S2、分别以三组标准样品为测量样品,获取三组测量光谱;
三组样品分别是方位角为0°的偏振片、方位角为0°与45°的两个偏振片组合、方位角为45°与0°的两个偏振片组合,分别用下标“(0)”、“(0,45)”与“(45,0)”标识上述标准样品,其穆勒矩阵为:
Figure BDA0002745134640000092
Figure BDA0002745134640000093
Figure BDA0002745134640000094
上述三个校准样品的测量光谱分别为:
Figure BDA0002745134640000095
Figure BDA0002745134640000101
Figure BDA0002745134640000102
S3、对测量光谱进行频域分析与分通道处理,选择相应频率通道;
首先对测量光谱进行频域分析,各阶频率通道与测量光谱三角展开式各阶频率之间的关系如式(12)所示:
Figure BDA0002745134640000103
式中,ν表示光程差,pj表示频率阶次,
Figure BDA0002745134640000104
表示傅里叶运算,角标“*”表示共轭运算,Cpj(ν)表示第pj阶频率通道,Cpj *(ν)表示第–pj阶频率通道。
然后选择通道:首先,对于标准样品1(方位角为0度的偏振片)的测量光谱的频域信号,选择第d1+d4=1+5=6阶频率通道作为第一组频率;接着选择d1=1阶频率通道作为第二组频率;然后对于标准样品2(方位角为0°与45°的两个偏振片组合)的测量光谱的频域信号,选择第d1+d3+d4=1+5+5=11阶频率通道作为第三组频率;最后对于标准样品3(方位角为45°与0°的两个偏振片组合)的测量光谱的频域信号,选择第d1+d2+d4=1+1+5=7阶频率通道作为第四组频率。
S4、将步骤S3所挑选出的频率通道进行波数域分析计算,获得测量光谱三角函数展开式对应频率的实频系数g′sec与虚频系数h′sec
求解g′sec与h′sec的过程如式(13-a)至式(13-f)所示:
Figure BDA0002745134640000105
Figure BDA0002745134640000106
Figure BDA0002745134640000111
Figure BDA0002745134640000112
Figure BDA0002745134640000113
Figure BDA0002745134640000114
其中,下标“sec”表示步骤S3选定的频率,“(σ)”表示变量为波数σ的函数,
Figure 1
表示傅里叶逆运算,“Re”与“Im”分别表示取复数的实部与虚部,c表示通道C做傅里叶逆变换后的结果,amp表示c的幅值,ζ表示c的相位角,
Figure BDA0002745134640000116
表示g′与h′的相位角,
Figure BDA0002745134640000117
S5、根据测量系统模型,由步骤S4所得的实频与虚频系数求解相位延迟量的线性组合,进一步通过线性运算求解四片相位延迟器的相位延迟量;
由测量系统模型可知,步骤S3所选定的频率的实频系数g′sec与虚频系数h′sec与相位延迟量误差的关系为:
Figure BDA0002745134640000118
Figure BDA0002745134640000119
Figure BDA00027451346400001110
Figure BDA00027451346400001111
Figure BDA00027451346400001112
Figure BDA00027451346400001113
Figure BDA00027451346400001114
Figure BDA00027451346400001115
求解时,先计算相位延迟量误差的线性组合Δe1、Δe2、Δe3、Δe4
Figure BDA0002745134640000121
Figure BDA0002745134640000122
Figure BDA0002745134640000123
Figure BDA0002745134640000124
进一步由:
Figure BDA0002745134640000125
可得相位延迟量误差Δδ1、Δδ2、Δδ3、Δδ4为:
Figure BDA0002745134640000126
下面给出通过上述步骤对系统的相位延迟量误差的仿真校正结果。
在考虑上述色散效应以及厚度误差后,各个相位延迟器的相位延迟量误差如图4中实线所示。经过上述校准方法得到的相位延迟量误差如图4 中点划线所示;进一步地,选择校准结果的中心的1/3波段,将相位延迟量的校准结果做波数的二阶多项式拟合,进一步校准相位延迟量误差,如图4 中虚线所示,可见本方法所述相位延迟器校准方法对快照式穆勒矩阵椭偏仪相位延迟量校准的有效性。
为了进一步说明校准结果的有效性,对上述相位延迟量误差以及校准结果进行了仿真测量实验。所用的仿真测量样品为方位角15°放置的偏振片,其理论穆勒矩阵如图5中实线所示。图5中点划线展示了存在相位延迟量误差,且未经校准的测量结果,其穆勒矩阵测量结果完全偏离真实值;图5中虚线给出了校准相位延迟量误差后的仿真测量结果,从中可以明显看出,在校准相位延迟量误差后,快照式穆勒矩阵椭偏仪的测量结果非常准确,与待测样品的理论值差别极小(小于1%),进一步验证了所述校准方法的有效性。
实施例二
为了进一步说明本发明所述快照式穆勒矩阵椭偏仪相位延迟量误差校准方法的通用性,即其适用于任何可以求解全穆勒矩阵的厚度比配置,此处给出另一组具体实施例,即相位延迟器厚度比为1:4:2:9的测量配置。
鉴于不同厚度比的测量配置,其测量系统基本模型与相位延迟量校准基本流程并未发生改变,因此在本实施例中并不详细给出各个步骤的处理与计算过程,仅给出因厚度比不同而产生的不同选择结果。
在步骤S1中,其厚度比改为d1=1,d2=4,d3=2,d4=9;系统最高调制频率mo为mo=1+4+2+9=16;
对于本实施例,设定测量波段为λ=505.4nm~547.5nm(波数σ= 1826.6mm-1~1978.6mm-1),选择材料为冰洲石,其测量波段内的双折射率如图6所示。相位延迟器基础厚度d0=0.67mm。对各相位延迟器添加厚度误差,分别是Δd1=–0.005mm、Δd2=0.01mm、Δd3=0.02mm、Δd4=0.01mm。
在步骤S3中,当厚度比为1:4:2:9时,通道选择结果为,首先,对于标准样品1(方位角为0度的偏振片)的测量光谱的频域信号,选择第d1+d4=1+9=10阶频率通道作为第一组频率;接着选择d1=1阶频率通道作为第二组频率;然后对于标准样品2(方位角为0°与45°的两个偏振片组合) 的测量光谱的频域信号,选择第d1+d3+d4=1+2+9=12阶频率通道作为第三组频率;最后对于标准样品3(方位角为45°与0°的两个偏振片组合) 的测量光谱的频域信号,选择第d1+d2+d4=1+4+9=14阶频率通道作为第四组频率。
在步骤S5中,由于步骤S3的通道选择结果发生改变,因此在步骤S5 中对应的频率也会发生改变,如式(18-a)至式(18-h)所示。
Figure BDA0002745134640000141
Figure BDA0002745134640000142
Figure BDA0002745134640000143
Figure BDA0002745134640000144
Figure BDA0002745134640000145
Figure BDA0002745134640000146
Figure BDA0002745134640000147
Figure BDA0002745134640000148
后续计算过程与实施例一中相同,下面给出一组待测样品仿真测量实例。
在本实施例中,在所述色散及相位延迟器厚度误差的影响下,各个相位延迟器的相位延迟量误差如图8中实线所示。经过本发明所述方法校准所得的相位延迟量误差如图8中点划线所示;然后同样选择校准结果的中心的1/3波段,将相位延迟量的校准结果做波数的二阶多项式拟合,进一步校准相位延迟量误差,如图8中虚线所示,可见本方法可以准确校准快照式穆勒矩阵椭偏仪在全测量波段内的相位延迟量误差。
进一步的,对所述测量配置进行仿真测量实验,仿真测量样品形貌结构如图7所示,为Si基底上周期Λ=600nm的Si光栅,其线宽CD=300nm,光栅高度H=100nm,以入射角为60°、方位角90°(沿光栅线方向),利用严格耦合波算法仿真获得其穆勒矩阵,如图9中实线所示。图9中点划线展示了存在相位延迟量误差,且未经校准的测量结果,此时穆勒矩阵测量结果完全偏离真实值;图9中虚线给出了校准相位延迟量误差后的仿真测量结果,从中可以明显看出,在校准相位延迟量误差后,快照式穆勒矩阵椭偏仪的测量结果非常准确,与待测样品的理论值差别极小,进一步验证了所述校准方法的有效性,也说明了本发明所述校准方法的通用性,可以用于任意相位延迟器厚度比的快照式穆勒矩阵椭偏仪的相位延迟量误差校准。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其特征在于,包括如下步骤:
S1、按预设相位延迟器厚度比搭建快照式穆勒矩阵椭偏仪,通过该快照式穆勒矩阵椭偏仪分别对三组标准样品进行测量,得到三组测量光谱;
S2、分别对三组测量光谱进行频域分析得到频域信号,然后对频域信号进行分通道处理得到多阶频率通道,再根据相位延迟器厚度比对频率通道进行选择;
S3、对选择的频率通道进行波数域分析计算,获得各频率通道对应的测量光谱三角函数展开式的实频系数与虚频系数,进而得到各相位延迟器的相位延迟量误差,完成快照式穆勒矩阵椭偏仪相位延迟量误差校准。
2.如权利要求1所述的快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其特征在于,还包括如下步骤:
S4、选择校准后相位延迟量误差的中心波段,对其进行二阶多项式拟合,进一步校准相位延迟量误差。
3.如权利要求1所述的快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其特征在于,所述快照式穆勒矩阵椭偏仪中,以起偏器方位角作为参考,检偏器方位角为90°,四片相位延迟器的方位角依次分别为45°、0°、0°、45°。
4.如权利要求3所述的快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其特征在于,所述三组标准样品中,第一组标准样品为方位角为0°的偏振片,第二组标准样品为方位角为0°与45°的两个偏振片组合,第三组标准样品为方位角为45°与0°的两个偏振片组合。
5.如权利要求1所述的快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其特征在于,所述步骤S2中,先通过傅里叶变换对测量光谱进行频域分析得到频域信号,然后按由相位延迟器基础厚度决定的基础频率f0将频域信号划分为多阶频率通道,其中,第j阶频率通道中频率范围为(j–1/2)f0~(j+1/2)f0
6.如权利要求4所述的快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其特征在于,所述相位延迟器厚度的最简整数比为d1:d2:d3:d4;所述步骤S2中,频率通道的选择原则为:对于第一组标准样品测量光谱的频域信号,先选择d1+d4阶频率通道,然后,若d1≠|d1–d4|,则选择d1阶频率通道,否则选择d4阶频率通道;对于第二组标准样品测量光谱的频域信号,选择d1+d3+d4阶频率通道;对于第三组标准样品测量光谱的频域信号,选择d1+d2+d4阶频率通道,以此完成四个频率通道的选择。
7.如权利要求1所述的快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法,其特征在于,所述步骤S3中,对选择的频率通道进行波数域分析计算具体包括:对选择的频率通道进行傅里叶逆变换,然后根据傅里叶逆变换结果的相位信息,求取各频率通道对应的测量光谱三角函数展开式的实频系数与虚频系数。
8.如权利要求1-7任一项所述的快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准系统,其特征在于,所述步骤S3中,根据求得的实频系数与虚频系数求解相位延迟量误差的线性组合,进而通过线性运算得到四片相位延迟器的相位延迟量误差。
CN202011163917.5A 2020-10-27 2020-10-27 快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法 Active CN112378861B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011163917.5A CN112378861B (zh) 2020-10-27 2020-10-27 快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011163917.5A CN112378861B (zh) 2020-10-27 2020-10-27 快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法

Publications (2)

Publication Number Publication Date
CN112378861A true CN112378861A (zh) 2021-02-19
CN112378861B CN112378861B (zh) 2021-10-15

Family

ID=74576733

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011163917.5A Active CN112378861B (zh) 2020-10-27 2020-10-27 快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法

Country Status (1)

Country Link
CN (1) CN112378861B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101464576A (zh) * 2009-01-14 2009-06-24 北京航空航天大学 液晶相位可变延迟器特性参数定标系统
CN104677837A (zh) * 2013-11-26 2015-06-03 北京智朗芯光科技有限公司 一种全穆勒矩阵椭圆偏振仪的校准方法
CN111307722A (zh) * 2019-11-29 2020-06-19 北京理工大学 一种多功能Stokes-Mueller矩阵成像偏振仪的标定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101464576A (zh) * 2009-01-14 2009-06-24 北京航空航天大学 液晶相位可变延迟器特性参数定标系统
CN104677837A (zh) * 2013-11-26 2015-06-03 北京智朗芯光科技有限公司 一种全穆勒矩阵椭圆偏振仪的校准方法
CN111307722A (zh) * 2019-11-29 2020-06-19 北京理工大学 一种多功能Stokes-Mueller矩阵成像偏振仪的标定方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
HU DAI等: "Theoretical research of retarder phase deviation in channeled Mueller matrix spectropolarimeters", 《OPTICS EXPRESS》 *
MATTHIEU DUBREUIL等: "ystematic errors specific to a snapshot Mueller matrix polarimeter", 《APPLIED OPTICS》 *
XUEPING JU等: "Easily implemented approach for the calibration of alignment and retardation errors in a channeled spectropolarimeter", 《OPTICAL SOCIETY OF AMERICA》 *
李建慧等: "Mueller 矩阵成像偏振仪的误差标定和补偿研究", 《激光与光电子学进展》 *
杨斌: "强度调制型偏振光谱成像系统测量误差及定标研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Also Published As

Publication number Publication date
CN112378861B (zh) 2021-10-15

Similar Documents

Publication Publication Date Title
CN110261317B (zh) 一种Mueller矩阵光谱的测量系统及方法
US11187649B2 (en) Method for conducting optical measurement usingfull Mueller matrix ellipsometer
CN111122460B (zh) 单旋转补偿器型光谱椭偏仪参数校准方法和装置
CN105628343B (zh) 一种波片检测装置及方法
CN112378860B (zh) 旋转器件型穆勒矩阵椭偏仪系统参数的校准方法
CN113281267B (zh) 一种双旋转补偿器型穆勒矩阵椭偏仪系统参数的校准方法
CN109470154B (zh) 一种适用于光谱椭偏仪的薄膜厚度初值测量方法
CN104677838A (zh) 一种自校准的全穆勒矩阵椭偏仪测量系统
CN104864815B (zh) 校准椭偏测量中应力元件带来的误差影响的方法
CN107677370B (zh) 一种通道型偏振光谱成像仪的偏振辐射标定系统以及方法
TWI615604B (zh) 寬波段消色差複合波片的定標方法
Hinrichs et al. Comparative study of an anisotropic polymer layer by infrared spectroscopic techniques
CN112378861B (zh) 快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法
JP3311497B2 (ja) フーリエ変換分光位相変調偏光解析法
CN104677837A (zh) 一种全穆勒矩阵椭圆偏振仪的校准方法
CN109115695B (zh) 一种各向异性体材料光学常数和欧拉角的提取方法
CN112179851B (zh) 一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法
CN111122458A (zh) 一种单旋转补偿器型椭偏仪随机误差评估方法
JP5361843B2 (ja) 光学的異方性の評価方法及び評価装置
CN113358604B (zh) 一种斜入射式光谱型反射差分测量装置及方法
CN113654996A (zh) 一种复合消色差波片相位延迟量测量装置和方法
CN103968783B (zh) 一种测量双片波片补偿器中光轴偏差角的方法
Schwiecker et al. High-speed ellipsometry for the production of thin metal layers
CN104677836A (zh) 一种自校准的全穆勒矩阵椭圆偏振仪测量系统
CN115727952A (zh) 一种偏振度测量仪的全局校准与测量方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant