CN112179851A - 一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法 - Google Patents
一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法 Download PDFInfo
- Publication number
- CN112179851A CN112179851A CN202011044645.7A CN202011044645A CN112179851A CN 112179851 A CN112179851 A CN 112179851A CN 202011044645 A CN202011044645 A CN 202011044645A CN 112179851 A CN112179851 A CN 112179851A
- Authority
- CN
- China
- Prior art keywords
- error
- azimuth
- mueller matrix
- azimuth angle
- sample
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/21—Polarisation-affecting properties
- G01N21/211—Ellipsometry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/21—Polarisation-affecting properties
- G01N21/211—Ellipsometry
- G01N2021/213—Spectrometric ellipsometry
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Chemical & Material Sciences (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明属于精密光学仪器的系统参数校准相关技术领域,其公开了一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法,所述方法首先建立包含器件方位角误差的系统模型,然后对模型进行合理约减,通过标准样品获得足够的线性无关的方程组以求解器件方位角误差,并将求解结果拟合为不随波长变动的定值作为最终校准结果,此外,将校准结果代入至系统模型中可以进一步校准待测样品穆勒矩阵。本发明公开的方法能够准确地校准快照式穆勒矩阵椭偏仪各个器件的方位角误差,能够同时适用于透射式与反射式测量配置;同时,所述方法校准出方位角误差后,可以代入至包含方位角误差的系统模型,能够有效提升穆勒矩阵的测量精度。
Description
技术领域
本发明属于精密光学器件的系统参数校准相关技术领域,更具体地,涉及一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法。
背景技术
穆勒矩阵椭偏仪可以实现待测样品穆勒矩阵的测量,穆勒矩阵可以表征样品的偏振响应特性,在样品光学特性与物理特性表征、偏振成像、遥感等领域具有广泛应用。传统商用椭偏仪采用基于时间调制的方法求解穆勒矩阵,其测量时间通常为数秒,在偏振成像等需要快速测量穆勒矩阵的领域面临时间过长的困境。Kazuhiko Oka等在文章“Spectroscopic polarimetry with a channeled spectrum”中提出的基于波长调制的通道式斯托克斯偏振仪,单次测量即可获得待测光束偏振态,而后Matthieu Dubreuil等在文章“Snapshot Mueller matrix polarimeter by wavelength polarization coding”中将其发展为快照式穆勒矩阵椭偏仪,能够快速获得待测样品的穆勒矩阵。
快照式穆勒矩阵椭偏仪的仪器组成如图2所示,通常由光源201、起偏器202、四片具有一定厚度比的相位延迟器203~206、检偏器207、探测器208等部分组成。为了实现对待测样品的调制,通常要求上述器件按照一定方位角放置。在透射式测量系统中,若以起偏器方位角为参考,其余各个器件相对于参考位置的方位角误差通常设置为:第1片至第4片相位延迟器的方位角为45°、0°、0°、45°,检偏器方位角为0°或90°;若考虑快照式穆勒矩阵椭偏仪的反射式模型,则需要以入射面作为参考,记入射面内方位角为0°,则需要设置起偏器方位角为0°,第1片至第4片相位延迟器的方位角为45°、0、0°、45°,检偏器方位角为0°或90°。
但是在实际测量过程中,各个器件无法完美地安装至设计位置,总是存在一定的安装方位角误差。在仿真过程中发现,当各个器件存在方位角误差时,即使较小的方位角误差也会对待测样品的穆勒矩阵测量结果产生较大影响。Matthieu Dubreuil等在文章“Systematic errors specific to a snapshot Mueller matrix polarimeter”中初步评估了相位延迟器方位角误差对穆勒矩阵测量结果的影响,但是由于方位角误差模型太过繁琐而未给出校准方法。Tingkui Mu等在文章“Alignment and retardance errors,andcompensation of a channeled spectropolarimeter”中给出了针对通道式偏振仪的校准方法,且效果较好,但是该方法仅考虑了相位延迟器的方位角误差,较为繁琐,而且是适用于通道式斯托克斯偏振仪的校准方法,难以用于快照式穆勒矩阵椭偏仪的器件方位角误差校准。因此有必要对快照式穆勒矩阵椭偏仪器件方位角进行校准。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法,所述方法能够实现快照式穆勒矩阵椭偏仪器件方位角的准确校准,将校准结果代入至包含方位角误差的系统模型中可以有效降低因仪器方位角误差引起的穆勒矩阵测量误差。
为实现上述目的,本发明提供了一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法,所述方法包括以下步骤:
步骤S1,建立包含器件方位角误差的快照式穆勒矩阵椭偏仪的系统模型;
步骤S2,根据器件方位角误差的特点对系统模型进行约减;
步骤S3,通过测量标准样品以建立方位角误差与测量结果之间的关系,并且选出预定数量的线性无关的方程组;
步骤S4,由步骤S3选出的线性无关的方程组,初步求解各个器件的方位角误差;
步骤S5,将步骤S4所得的方位角误差拟合为不随测量波长变化的定值,以作为方位角误差的最终校准结果;
步骤S6,将步骤S5所求得的最终方位角误差代入至步骤S2所得的约减后的系统模型,以校正由方位角误差引起的穆勒矩阵测量误差。
进一步地,对于第i个器件的穆勒矩阵Mi,将其旋转一个微小的误差角度εi,则考虑方位角误差后的器件实际穆勒矩阵Mi′为:
M′i=RM(-εi)MiRM(εi),i=1,2,3,4,5,6 (1-a)
其中,RM表示坐标旋转矩阵。
进一步地,包含各个器件方位角误差的系统模型的表达式为:
式中,P1(α1)、P2(α2)分别表示方位角为α1的起偏器与方位角为α2的检偏器的穆勒矩阵;Rt(θt,δt)(t=1,2,3,4)分别表示第t片相位延迟器的穆勒矩阵,其方位角为θt,相位延迟量为δt;M为待测样品的穆勒矩阵;Sin与Sout分别表示入射光与出射光的斯托克斯矢量,通常Sin=I0[1 0 0 0]。
进一步地,经步骤S2进行约减后的包含各个器件方位角误差的系统模型的表达式为:
S′out=RM(-ε6)P2(90°)RM(ε6)RM(-ε5)R4(45°,δ4)RM(ε5)RM(-ε4)R3(0°,δ3)RM(ε4)MRM(-ε3)R2(0°,δ2)RM(ε3)RM(-ε2)R1(45°,δ1)RM(ε2)RM(-ε1)P1(0°)RM(ε1)Sin (3)。
进一步地,步骤S2中,安装时将方位角控制在一定误差范围内,故能得到近似公式(4-a)与公式(4-b):
sin2εi=2εi,i=1,2,3,4,5,6 (4-a)
εi小于0.01时,二阶误差小于10-4,将二阶及以上误差项忽略,得到公式(4-c):
εiεn=0,i,n=1,2,3,4,5,6 (4-c)
进而测量光谱的表达式为,
其中,fmax表示系统最高频率,fmax=k1+k2+k3+k4,k1:k2:k3:k4为相位延迟器的最简整数厚度比;gj′与hj′分别表示展开式的实频与虚频系数,二者具有如下形式:
式中,ml表示待测样品穆勒矩阵的第l个元素;aj,l,0与aj,l,i表示常数项系数。
进一步地,步骤S3中,标准样品为穆勒矩阵已知的样品。
进一步地,通过对校准样品的测量光谱进行频域分析与分通道处理即可获得测量光谱的三角展开式的实频系数gj′与虚频系数hj′,若记,
ε=[ε1 ε2 ε3 ε4 ε5 ε6]T (6-a)
则有,
g′j,h′j=[qj,1 qj,2 qj,3 qj,4 qj,5 qj,6]ε+bj (6-d)
将校准样品的实频系数gj′与虚频系数hj′代入至公式(6-a)至公式(6-d)所表达的系统模型中,从中挑选出预定数量的线性无关的方程组。
进一步地,步骤S6中,将公式(5-b)的实频系数gj′与虚频系数hj′进行变形得到:
将校准所得最终误差代入至公式(7)中以求解系数pj,l,将公式(7)按频率排列为矩阵形式,得到公式(8):
式中,M表示将待测样品穆勒矩阵重新排列的列向量;P为由方位角误差及系统设置共同决定的系数矩阵;由此求得校正方位角误差后的待测样品穆勒矩阵为:
M=P+V (9)
式中P+表示矩阵P的广义逆矩阵。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,本发明提供的快照式穆勒矩阵椭偏仪器件的方位角误差校准方法主要具有以下有益效果:
1.本发明提供了一种快照式穆勒矩阵椭偏仪器件方位角误差的校准方法,能够有效校准出各器件安装方位角误差,同时适用于透射式与反射式测量配置。
2.本发明所述方法对方位角误差进行了有效的约减,建立了包含器件方位角误差的快照式穆勒矩阵椭偏仪的系统模型。
3.本发明所述方法校准出方位角误差后,可以代入至包含方位角误差的系统模型,能够有效提升穆勒矩阵的测量精度。
4.本发明所述校准方法能够准确地校准器件方位角误差与待测样品穆勒矩阵,即使仪器存在较大的器件方位角误差,仍能准确获得待测样品穆勒矩阵,从而极大地降低了器件安装精度要求。
附图说明
图1是本发明提供的快照式穆勒矩阵椭偏仪器件的方位角误差校准方法的流程示意图100;
图2是快照式穆勒矩阵椭偏仪仪器组成示意图;
图3是本发明较佳实施例中校准样品的穆勒矩阵Mcs,图中的mij表示Mcs的第i行、第j列(i=1,2,3,4,j=1,2,3,4)个元素在测量波段内的变化曲线;
图4是本发明较佳实施例中存在方位角误差且按照本发明所述方法校准所得的方位角误差的初步结果,图中的ε1c、ε2c、ε3c、ε4c、ε5c及ε6c分别是起偏器、四片相位延迟器、检偏器的方位角误差的初步校准结果;
图5是本发明较佳实施例中仿真测量样品的穆勒矩阵Ms,图中的mij表示Ms的第i行、第j列(i=1,2,3,4,j=1,2,3,4)个元素在测量波段内的变化曲线;
图6是本发明较佳实施例中存在方位角误差且未校准时待测样品穆勒矩阵仿真测量误差ΔM,图中的Δmij表示ΔM的第i行、第j列(i=1,2,3,4,j=1,2,3,4)个元素在测量波段内的变化曲线;
图7是本发明较佳实施例中存在方位角误差且按照本发明所述方法进行校准后的待测样品穆勒矩阵仿真测量误差ΔM',图中的Δm'ij表示ΔM'的第i行、第j列(i=1,2,3,4,j=1,2,3,4)个元素在测量波段内的变化曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
请参阅图1,本发明提供的快照式穆勒矩阵椭偏仪器件的方位角误差校准方法,所述校准方法主要包括以下步骤:
步骤S1,建立包含器件方位角误差的快照式穆勒矩阵椭偏仪的系统模型。
具体地,包含器件方位角误差的系统模型的建立方法为:对于第i个器件的穆勒矩阵Mi,将其旋转一个微小的误差角度εi,则考虑方位角误差后器件实际穆勒矩阵Mi′为,
M′i=RM(-εi)MiRM(εi),i=1,2,3,4,5,6 (1-a)
其中,RM表示坐标旋转矩阵。因此,包含各个器件方位角误差的系统模型可以描述为:
式中,P1(α1)、P2(α2)分别表示方位角为α1的起偏器与方位角为α2的检偏器的穆勒矩阵,Rt(θt,δt)(t=1,2,3,4)分别表示第t片相位延迟器的穆勒矩阵,其方位角为θt,相位延迟量为δt,M为待测样品的穆勒矩阵,Sin与Sout分别表示入射光与出射光的斯托克斯矢量,通常Sin=I0[1 0 0 0]。当相位延迟器基础厚度为d0,厚度比为k1:k2:k3:k4时,相位延迟量可以近似为:
δi(σ)=2πΔn(σ0)kid0σ,i=1,2,3,4 (3)
其中,Δn表示材料的双折射率,σ表示波数,σ0表示中心波数。
本实施方式选用相位延迟器厚度比为k1:k2:k3:k4=1:1:5:5,系统最高调制频率fmax=1+1+5+5=12;以起偏器方位角为参考,检偏器方位角为90°,相位延迟器方位角依次为45°、0°、0°、45°。记起偏器、第一相位延迟器至第四相位延迟器、检偏器的方位角误差分别为ε1~ε6,此时系统模型为:
此外,本实施例的其他测量配置为,相位延迟器材料为冰洲石,测量波段为σ=1696.4mm-1~1751.3mm-1,基础厚度d0=0.89mm;其在测量波段中心处双折射率为Δn(σ0)=Δn(1723.9mm-1)=0.1724,此时各个相位延迟器的相位延迟量为:
δt(σ)=2πΔn(σ0)ktd0σ=2πktf0σ=2π×0.1534×ktσ,t=1,2,3,4 (11)。
步骤S2,根据器件方位角误差的特点,对系统模型进行约减。
具体地,所述的约减包括,由于安装时可以将方位角控制在一定误差范围内,如0.5°之内,因此可以做如公式(4-a)与公式(4-b)所示近似,
sin2εi=2εi,i=1,2,3,4,5,6 (4-a)
由于方位角误差角小时,εi小于0.01,因此二阶误差小于10-4,可以将二阶及以上误差项忽略,如公式(4-c)所示,
εiεn=0,i,n=1,2,3,4,5,6 (4-c)
经过上述约减后,测量光谱可以表述为,
其中,fmax表示系统最高频率,fmax=k1+k2+k3+k4,gj′与hj′分别表示展开式的实频与虚频系数,二者具有如下形式:
式中,ml表示待测样品穆勒矩阵的第l个元素,aj,l,0与aj,l,i表示常数项系数。
本实施方式中,设入射光的斯托克斯矢量Sin=I0[1 0 0 0],将式(10)进行三角展开,根据式(4-a)至式(4-c)所示方法对步骤S1所得的系统模型进行约减(也可以根据式(4-a)与式(4-b)首先对坐标旋转矩阵RM进行约减,在将系统模型三角展开后进一步约减),可得如式(7)所示展开光谱,其中各项实频系数gj′与虚频系数hj′分别为,
g′4=1/16(2ε2m23-2ε1m23-m22-2ε3m23-2ε4m32+2ε5m32-2ε6m32) (12-e)
g′6=1/16(2ε2m23-2ε1m23-m22-2ε3m23-2ε4m32+2ε5m32-2ε6m32) (12-g)
h1′=1/4(-2ε1m14+ε1m34+2ε2m14-ε2m34)(12-n)
h′4=1/8(-ε1m24+ε2m24+ε5m42-ε6m42) (12-q)
h′5=1/8(2ε5m41+ε5m43-2ε6m41-ε6m43) (12-r)
h′6=1/8(ε1m24-ε2m24+ε5m42-ε6m42) (12-s)
步骤S3,通过测量标准样品以建立方位角误差与测量结果之间的关系,并且选出预定数量的线性无关的方程组。
具体地,所述标准样品为穆勒矩阵已知的样品,对于透射式仪器,优选穆勒矩阵已知且不随波长变化的样品,如空气、具有一定方位角的偏振片等;对于反射式仪器,优选性质稳定且已知的薄膜样品,如具有一定厚度的SiO2薄膜等;
建立方位角误差与测量结果之间关系的方法为,通过对校准样品的测量光谱进行频域分析与分通道处理即可获得测量光谱的三角展开式的实频系数gj′与虚频系数hj′,若记,
ε=[ε1 ε2 ε3 ε4 ε5 ε6]T (6-a)
则有,
g′j,h′j=[qj,1 qj,2 qj,3 qj,4 qj,5 qj,6]ε+bj (6-d)
将校准样品的实频系数gj′与虚频系数hj′代入至式(6-a)至式(6-d)所示的系统模型中,从中挑选出足够的线性无关的方程组。
本实施方式选择标准样品为400nm厚的SiO2薄膜,校准时入射角设置为45°,其穆勒矩阵Mcs如图3所示。获得其测量光谱后,经过频域分析,分通道计算等,可得到各项实频系数gj′与虚频系数hj′。
将待测样品穆勒矩阵与各项实频系数gj′与虚频系数hj′代入至式(6-a)至式(6-d),可得到不同的关于器件方位角误差ε的方程,本实施例选择全部的方程用于求解ε,此时求解方法为,
ε=Q+Vε (13-b)
步骤S4,由步骤S3选出的线性无关的方程组,初步求解各个器件的方位角误差。
具体地,由式(13-a)至式(13-b)可以初步获得各器件的方位角误差。在本实施例中给出一组仿真案例,设置起偏器、四片相位延迟器、检偏器的方位角误差依次为ε1=0.5°=0.0087rad,ε2=0.5°=0.0087rad,ε3=-0.5°=-0.0087rad,ε4=-0.5°=-0.0087rad,ε5=0.5°=0.0087rad,ε6=0.5°=0.0087rad,经过上述案例仿真获得的结果如图4所示。
步骤S5,将步骤S4所得的方位角误差拟合为不随测量波长变化的定值,以作为方位角误差的最终校准结果。
具体地,步骤S5采用的拟合方法可以为Levenberg-Marquardt算法、遗传算法等常见的拟合算法。由步骤S4所得的方位角误差是随波长变化曲线,但实际方位角误差为一定值,因此对所得结果进一步拟合,获取一组不随波长变化的定值,作为最终器件安装方位角误差的最终校准结果。
本仿真案例采用Levenberg-Marquardt算法进行拟合,由于边缘波段数据受截断与加窗等的影响较大,因此在拟合时丢弃边缘区域,所得最终方位角误差如表一所示,可见本发明所述方法能够较好的校准器件方位角误差。
表一:方位角误差最终校准结果表
步骤S6,将步骤S5所求得的最终方位角误差代入至步骤S2所得的约减系统模型,校正由方位角误差引起的穆勒矩阵测量误差。
具体地,步骤S6所述校准方法为,将式(5-b)所示的实频系数gj′与虚频系数hj′进行变形,可得,
将校准所得最终误差代入至式(7)中,可以求解系数pj,l,将式(7)按频率排列为矩阵形式,如式(8)所示,
式中,M表示将待测样品穆勒矩阵重新排列为列向量,P为由方位角误差及系统设置共同决定的系数矩阵;由此可求得校正方位角误差后的待测样品穆勒矩阵为:
M=P+V (9)
式中P+表示矩阵P的广义逆矩阵。
若存在上述方位角误差而不对其进行校准,则会对待测样品穆勒矩阵产生较大影响。本实施例以Si基底上200nm厚的SiO2薄膜为测量样品,入射角设置为50°,其穆勒矩阵Ms(图中各元素具有相同的坐标)如图5所示。
如果测量过程中存在上述方位角误差而不对其进行校准,那么测量结果将会产生较大的误差,图6展示了在存在方位角误差而不校准的情况下,通过快照式穆勒矩阵椭偏仪仿真测量所得的穆勒矩阵与理论穆勒矩阵的误差ΔM=Ms-Mc(图中各元素具有相同的坐标),从中可以看出,即使较小的方位角误差,也会产生较大的穆勒矩阵测量误差,在本实施例中,穆勒矩阵测量误差最大已经超过0.06,这在穆勒矩阵测量中是不允许的。
可见,必须对方位角误差进行校准,并且将其考虑至系统模型中。将步骤S5校准所得的方位角误差代入至步骤S2中的包含了方位角误差的系统模型,用以求解待测样品穆勒矩阵。图7展示了校准后仿真测量所得的穆勒矩阵与理论穆勒矩阵的误差ΔM'=Ms-Mc'(图中各元素具有相同的坐标),从中可以看出,除了边缘波段受到截断与加窗等计算的影响外,其余测量波段的测量误差基本为0,最大的穆勒矩阵测量误差也小于0.004,穆勒矩阵测精度提高了一个数量级以上,且测量结果能够满足测量精度需求。
从本实施例中可以明显看出,本发明所述方位角误差校准方法能够有效地校准各个器件的方位角误差,此外将校准结果代入至本发明提出的包含了方位角误差的系统模型,能够极大地提升待测样品穆勒矩阵的测量精度,降低了快照式穆勒矩阵椭偏仪器件方位角的安装精度要求。
对于透射式快照式穆勒矩阵椭偏仪而言,其校准流程与反射式相同,可以将起偏器方位角为参考,此时起偏器方位角误差ε1=0,可以忽略该项。
本发明同时建立了包含器件方位角误差的系统模型,能够有效校准器件方位角误差以及由此引起的穆勒矩阵测量误差。所述方法首先建立包含器件方位角误差的系统模型,然后对模型进行合理约减,通过标准样品获得足够的线性无关的方程组以求解器件方位角误差,并将求解结果拟合为不随波长变动的定值作为最终校准结果,此外,将校准结果代入至系统模型中可以进一步校准待测样品穆勒矩阵。本发明公开的方法能够准确地校准快照式穆勒矩阵椭偏仪各个器件的方位角误差,能够同时适用于透射式与反射式测量配置;同时,本发明建立了包含器件方位角误差的快照式穆勒矩阵椭偏仪的理论模型,所述方法校准出方位角误差后,可以代入至包含方位角误差的系统模型,能够有效提升穆勒矩阵的测量精度;此外,所述方法可以在较大范围内准确校准器件方位角,从而极大地降低了器件安装精度要求。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法,其特征在于:
步骤S1,建立包含器件方位角误差的快照式穆勒矩阵椭偏仪的系统模型;
步骤S2,根据器件方位角误差的特点对系统模型进行约减;
步骤S3,通过测量标准样品以建立方位角误差与测量结果之间的关系,并且选出预定数量的线性无关的方程组;
步骤S4,由步骤S3选出的线性无关的方程组,初步求解各个器件的方位角误差;
步骤S5,将步骤S4所得的方位角误差拟合为不随测量波长变化的定值,以作为方位角误差的最终校准结果;
步骤S6,将步骤S5所求得的最终方位角误差代入至步骤S2所得的约减后的系统模型,以校正由方位角误差引起的穆勒矩阵测量误差。
4.如权利要求2所述的快照式穆勒矩阵椭偏仪器件的方位角误差校准方法,其特征在于:经步骤S2进行约减后的包含各个器件方位角误差的系统模型的表达式为:
S′out=RM(-ε6)P2(90°)RM(ε6)RM(-ε5)R4(45°,δ4)RM(ε5)RM(-ε4)R3(0°,δ3)RM(ε4)MRM(-ε3)R2(0°,δ2)RM(ε3)RM(-ε2)R1(45°,δ1)RM(ε2)RM(-ε1)P1(0°)RM(ε1)Sin (10)。
5.如权利要求1所述的快照式穆勒矩阵椭偏仪器件的方位角误差校准方法,其特征在于:步骤S2中,安装时将方位角控制在一定误差范围内,故能得到近似公式(4-a)与公式(4-b):
sin2εi=2εi,i=1,2,3,4,5,6 (4-a)
εi小于0.01时,二阶误差小于10-4,将二阶及以上误差项忽略,得到公式(4-c):
εiεn=0,i,n=1,2,3,4,5,6 (4-c)
进而测量光谱的表达式为,
其中,fmax表示系统最高频率,fmax=k1+k2+k3+k4,k1:k2:k3:k4为相位延迟器的最简整数厚度比;gj′与hj′分别表示展开式的实频与虚频系数,二者具有如下形式:
式中,ml表示待测样品穆勒矩阵的第l个元素;aj,l,0与aj,l,i表示常数项系数。
6.如权利要求1-5任一项所述的快照式穆勒矩阵椭偏仪器件的方位角误差校准方法,其特征在于:步骤S3中,标准样品为穆勒矩阵已知的样品。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011044645.7A CN112179851B (zh) | 2020-09-29 | 2020-09-29 | 一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011044645.7A CN112179851B (zh) | 2020-09-29 | 2020-09-29 | 一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112179851A true CN112179851A (zh) | 2021-01-05 |
CN112179851B CN112179851B (zh) | 2021-08-03 |
Family
ID=73946551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011044645.7A Active CN112179851B (zh) | 2020-09-29 | 2020-09-29 | 一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112179851B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220350261A1 (en) * | 2021-04-29 | 2022-11-03 | Huazhong University Of Science And Technology | Snapshot type overlay error measuring device and measuring method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT506599A4 (de) * | 2008-06-05 | 2009-10-15 | Aziz Abdulrasagh Dipl Ing Dr | Multiobjektidentifikation |
CN101936774A (zh) * | 2010-08-27 | 2011-01-05 | 中国科学院上海光学精密机械研究所 | 偏振检测装置中器件误差的测量方法 |
CN103134592A (zh) * | 2013-01-31 | 2013-06-05 | 华中科技大学 | 一种透射式全穆勒矩阵光谱椭偏仪及其测量方法 |
CN103217385A (zh) * | 2013-03-22 | 2013-07-24 | 华中科技大学 | 椭偏仪测量系统的随机误差评估方法 |
CN103411890A (zh) * | 2013-07-18 | 2013-11-27 | 华中科技大学 | 一种旋转补偿器型椭偏仪的系统误差评估及消除方法 |
CN104677837A (zh) * | 2013-11-26 | 2015-06-03 | 北京智朗芯光科技有限公司 | 一种全穆勒矩阵椭圆偏振仪的校准方法 |
-
2020
- 2020-09-29 CN CN202011044645.7A patent/CN112179851B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT506599A4 (de) * | 2008-06-05 | 2009-10-15 | Aziz Abdulrasagh Dipl Ing Dr | Multiobjektidentifikation |
CN101936774A (zh) * | 2010-08-27 | 2011-01-05 | 中国科学院上海光学精密机械研究所 | 偏振检测装置中器件误差的测量方法 |
CN103134592A (zh) * | 2013-01-31 | 2013-06-05 | 华中科技大学 | 一种透射式全穆勒矩阵光谱椭偏仪及其测量方法 |
CN103217385A (zh) * | 2013-03-22 | 2013-07-24 | 华中科技大学 | 椭偏仪测量系统的随机误差评估方法 |
CN103411890A (zh) * | 2013-07-18 | 2013-11-27 | 华中科技大学 | 一种旋转补偿器型椭偏仪的系统误差评估及消除方法 |
CN104677837A (zh) * | 2013-11-26 | 2015-06-03 | 北京智朗芯光科技有限公司 | 一种全穆勒矩阵椭圆偏振仪的校准方法 |
Non-Patent Citations (3)
Title |
---|
BHATTACHARYYA, K等: ""Mueller Matrix Polarimeter with Diattenuation Error Calibration Approach"", 《ADVANCES IN OPTICAL SCIENCE AND ENGINEERING》 * |
DENNIS H.GOLDSTEIN等: ""Error analysis of a Mueller matrix polarimeter"", 《J.OPT.SOC.AM.A》 * |
陈修国等: ""Study of the retardance of a birefringent waveplate at tilt incidence by Mueller matrix ellipsometer"", 《JOURNAL OF OPTICS》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220350261A1 (en) * | 2021-04-29 | 2022-11-03 | Huazhong University Of Science And Technology | Snapshot type overlay error measuring device and measuring method |
US11619883B2 (en) * | 2021-04-29 | 2023-04-04 | Huazhong University Of Science And Technology | Snapshot type overlay error measuring device and measuring method |
Also Published As
Publication number | Publication date |
---|---|
CN112179851B (zh) | 2021-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103163077B (zh) | 旋转器件型光谱椭偏仪系统参数校准方法 | |
CN102878940B (zh) | 一种包含相位补偿器的椭圆偏振仪的校准方法 | |
US11187649B2 (en) | Method for conducting optical measurement usingfull Mueller matrix ellipsometer | |
CN104677838A (zh) | 一种自校准的全穆勒矩阵椭偏仪测量系统 | |
CN112378860B (zh) | 旋转器件型穆勒矩阵椭偏仪系统参数的校准方法 | |
CN109901281A (zh) | 一种穆勒显微偏振像差校准方法及装置 | |
CN104864815B (zh) | 校准椭偏测量中应力元件带来的误差影响的方法 | |
CN112179851B (zh) | 一种快照式穆勒矩阵椭偏仪器件的方位角误差校准方法 | |
CN107677370B (zh) | 一种通道型偏振光谱成像仪的偏振辐射标定系统以及方法 | |
CN112345464A (zh) | 一种椭偏仪优化校准方法 | |
CN104677833A (zh) | 一种利用全穆勒矩阵椭偏仪进行光学测量的方法 | |
CN104677837A (zh) | 一种全穆勒矩阵椭圆偏振仪的校准方法 | |
CN106197949A (zh) | 宽波段消色差复合波片的定标方法 | |
CN104677835A (zh) | 一种全穆勒矩阵椭圆偏振仪的校准方法 | |
CN106547132A (zh) | 液晶盒制造装置及方法 | |
CN109459138A (zh) | 基于四次光强测量的Mueller型椭偏仪椭偏参数测量方法及装置 | |
CN112284541B (zh) | 一种绝对辐射与偏振联合的定标方法及装置 | |
CN114384017A (zh) | 一种基于椭偏仪的光谱匹配校准方法 | |
CN107314839A (zh) | 基于穆勒矩阵的应力检测装置及方法 | |
CN112378861B (zh) | 快照式穆勒矩阵椭偏仪相位延迟量误差的通用校准方法 | |
CN104677836A (zh) | 一种自校准的全穆勒矩阵椭圆偏振仪测量系统 | |
JPH1038694A (ja) | エリプソメーター | |
US8570514B2 (en) | Optical system polarizer calibration | |
Riedling | Error effects in the ellipsometric investigation of thin films | |
CN116817768A (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 |