CN113281256B - 穆勒矩阵测量装置及其测量方法 - Google Patents
穆勒矩阵测量装置及其测量方法 Download PDFInfo
- Publication number
- CN113281256B CN113281256B CN202110605491.2A CN202110605491A CN113281256B CN 113281256 B CN113281256 B CN 113281256B CN 202110605491 A CN202110605491 A CN 202110605491A CN 113281256 B CN113281256 B CN 113281256B
- Authority
- CN
- China
- Prior art keywords
- light
- light path
- optical system
- mueller matrix
- polarization
- 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.)
- Active
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/01—Arrangements or apparatus for facilitating the optical investigation
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种穆勒矩阵的测量装置及其测量方法。穆勒矩阵的测量装置,包括起偏光路、检偏光路和传递光路,传递光路包括分光组件、第一透镜组、小孔光阑和平面反射镜。待测光学系统放置在传递光路中的小孔光阑和平面反射镜之间,并使得待测光学系统的焦点与第一透镜组的焦点重合,平面反射镜的口径大于待测光学系统的有效通光口径。穆勒矩阵的测量装置通过增加传递光路来实现对大口径、非平面的复杂光学系统偏振特性的测量。在该测量装置及其测量方法中,采用普通口径的偏振光学元件即可实现大口径复杂光学系统偏振特性的测量。本发明的自准直光路设计确保了光线均正入射到平面反射镜,减小了平面反射镜引入的测量误差,提高了测量精度。
Description
技术领域
本发明属于光学技术领域,具体涉及一种穆勒矩阵的测量装置及其测量方法。
背景技术
穆勒矩阵椭偏仪能够测量得到目标的穆勒矩阵。穆勒矩阵是一个4×4矩阵,它可以表征样本自身的全部偏振性质,包括二向衰减、相位延迟和退偏。近年来越来越多的研究人员通过测量样本的穆勒矩阵来获取薄膜的折射率和厚度、微纳结构的特征尺寸和光学常数、生物组织的散射特性等重要信息,使得穆勒矩阵椭偏仪广泛应用于通信光电子、光伏太阳能、光学薄膜、微纳器件、集成电路、生物组织等领域。
随着对地遥感、天文观测等水平的提高,望远镜的非理想偏振特性逐渐成为一个重要的系统误差源,限制了类地行星观测、太阳望远镜、空间目标探测、偏振激光雷达和偏振光谱遥感等技术领域的进一步发展。因此,上述应用中均迫切需要精确掌握望远镜的偏振特性。然而,望远镜的大口径和非平面光学系统属性大幅度增加了精确获取其穆勒矩阵的难度。
现有的穆勒矩阵椭偏仪中起偏臂和检偏臂关于待测样品轴线对称布置,起偏臂出射光经过待测目标后直接进入检偏臂,导致可检测的样品只能是平面。即便存在能够测量曲面的椭偏仪,但是其原理是通过减小起偏臂出射光斑的口径,以多次小口径平面检测的方式覆盖整个待测曲面样品,本质上依然属于平面测量。更重要的是,现有穆勒矩阵椭偏仪中待测样品的口径与起偏臂和检偏臂中偏振光学元件的口径相当。然而,目前技术条件下,高精度、大口径偏振光学元件很难制造,因此现有的穆勒矩阵椭偏仪均只能针对小口径的样品进行测量,无法实现大口径、复杂光学系统的偏振特性测量。
发明内容
针对现有技术的不足,本发明提供一种穆勒矩阵的测量装置;本发明还提供了使用该穆勒矩阵的测量装置进行测量的方法,以解决现有技术缺陷。
为实现上述目的,本发明提供了一种穆勒矩阵的测量装置,包括起偏光路和检偏光路,还包括位于所述起偏光路光束出射方向的传递光路;所述传递光路包括依次沿光路方向设置的分光组件、第一透镜组、位于所述第一透镜组焦点处的小孔光阑、平面反射镜;待测光学系统放置于所述小孔光阑和所述平面反射镜之间,所述待测光学系统的焦点与所述第一透镜组的焦点重合;所述起偏光路发出偏振态被调制的光束,经过所述分光组件分光,透射光束经第一透镜组汇聚于所述待测光学系统的焦点处,再经所述小孔光阑入射至所述待测光学系统,变为平行光束后经所述平面反射镜自准直反射,沿原光路返回,依次通过所述小孔光阑和第一透镜组后,经所述分光组件反射进入所述检偏光路,通过计算获得待测光学系统的穆勒矩阵。
进一步地,平面反射镜的口径大于待测光学系统的有效通光口径。
进一步地,所述分光组件包括分光面相互正交放置的第一分光棱镜和第二分光棱镜;经所述平面反射镜反射沿原光路返回的光束经所述第一分光棱镜反射折转至所述第二分光棱镜,再经所述第二分光棱镜的反射折转进入所述检偏光路。
进一步地,所述第一透镜组经消偏振处理,以降低所述第一透镜组对光束偏振态的改变。
进一步地,所述起偏光路包括依次沿光路方向设置的光源、第二透镜组、起偏器、第一相位补偿器;所述光源发出光强、偏振态、波长稳定的光束,经所述第二透镜组准直后,经所述起偏器变为线偏振光束,所述线偏振光束进入所述第一相位补偿器被调制后,入射至所述传递光路。
进一步地,所述检偏光路包括依次沿光路方向设置的第二相位补偿器、检偏器、第三透镜组、接收光强信号的光电探测器;经所述平面反射镜反射沿原光路返回的光束,经所述分光组件折转进入所述第二相位补偿器,其偏振态被调制后,进入所述检偏器,经所述第三透镜组汇聚至所述光电探测器进行光强探测。
本发明还提供一种穆勒矩阵测量的方法,包括如下步骤:
S1、起偏光路发出的偏振光束进入传递光路;所述偏振光束经过所述传递光路的分光组件后被分光,其中透射光束依次经第一透镜组、小孔光阑、待测光学系统、平面反射镜,经平面反射镜自准直反射后再沿原光路返回,进入检偏光路,其偏振态被解调,最后被光电探测器接收;
S2、旋转所述起偏光路中第一相位补偿器和所述检偏光路中第二相位补偿器,使第一相位补偿器和第二相位补偿器保持固定的转速比;
S3、重复步骤S2,使用光电探测器进行数据采集,获取频率在0倍到2×(u+v)倍范围内的光强信号,重复次数不少于16次,计算获得经待测光学系统在自准直光路中的穆勒矩阵M;
S4、根据经待测光学系统在自准直光路中的穆勒矩阵M,计算获得待测光学系统的穆勒矩阵M0。
进一步地,步骤S3中计算获得经待测光学系统在自准直光路中的穆勒矩阵M的具体方法如下:
起偏光路的矩阵公式为:g=R1(δ1,θ1(t))P1(α1);检偏光路的矩阵公式为:a=P2(α2)R2(δ2,θ2(t));
式中,P1、P2分别表示起偏器、检偏器的穆勒矩阵;R1、R2分别表示第一相位补偿器、第二相位补偿器的穆勒矩阵;α1,θ1,θ2和α2分别表示起偏器、第一相位补偿器、第二相位补偿器和检偏器的方位角;δ1和δ2分别是第一相位补偿器、第二相位补偿器的相位延迟量,
光电探测器采集的Q次光强信号IQ表示为:
式中,aq,j表示起偏光路穆勒矩阵第一行的各个元素,mj,k(j,k=1,2,3,4)表示待测光学系统在自准直光路中的穆勒矩阵M的各个元素,gq,k表示检偏光路穆勒矩阵第一列的各个元素,
WQ=(wq,1,1,wq,1,2,wq,1,3,…,wq,4,4)=(aq,1gq,1,aq,1gq,2,aq,1gq,3,…,aq,4gq,4);
M=(m11,m12,m13,m14,…,m44)T;
其中WQ表示第Q次测量时的仪器矩阵,T表示矩阵的转置运算,
待测光学系统在自准直光路中表现出的穆勒矩阵通过伪逆求解:
M=(WT·W)-1·WT·I。
进一步地,步骤S4中计算获得待测光学系统的穆勒矩阵M0的计算公式如下所示:
本发明所提供的一种穆勒矩阵的测量装置及其测量方法,通过增加传递光路来实现对大口径、非平面的复杂光学系统偏振特性的测量。在该测量装置及其测量方法中,由于平面反射镜与待测光学系统构成了自准直光路,光线均为正入射到平面反射镜,减小了平面反射镜引入的测量误差,提高了穆勒矩阵测量精度。
附图说明
图1是本发明提供的穆勒矩阵的测量装置的第一结构示意图;
图2是本发明提供的穆勒矩阵的测量装置的第二结构示意图;
图3是本发明提供的测量穆勒矩阵的方法的流程示意图。
具体附图标记包括:
分光组件1、第一分光棱镜101、第二分光棱镜102、第一透镜组2、小孔光阑3、平面反射镜4、光源5、第二透镜组6、起偏器7、第一相位补偿器8、第二相位补偿器9、检偏器10、第三透镜组11、光电探测器12、待测光学系统13。
具体实施方式
下面结合附图和实施例对本发明的实施方式作进一步详细描述。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互结合。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一方面,本发明提供一种穆勒矩阵的测量装置。
图1示出了本发明提供的穆勒矩阵测量装置的第一结构示意图,图2示出了本发明提供的穆勒矩阵测量装置的第二结构示意图。
如图1所示,穆勒矩阵测量装置包括起偏光路、检偏光路、传递光路,传递光路位于起偏光路的光束出射方向,包括依次沿光路方向设置的分光组件1、第一透镜组2、位于第一透镜组2焦点处的小孔光阑3,位于光路最末端的有效通光口径为1000mm的平面反射镜4。有效通光口径为500mm的待测光学系统13放置于小孔光阑3和平面反射镜4之间,待测光学系统13的焦点与第一透镜组2的焦点重合。
起偏光路发出偏振态被调制的偏振光束,透射经过分光组件1分光,透射光束经第一透镜组2后汇聚于待测光学系统13的焦点处的小孔光阑3。光束经小孔光阑3入射至待测光学系统13变为平行光束,并入射至平面反射镜4后,经平面反射镜4自准直反射再沿原光路返回,入射至分光组件1,经反射进入检偏光路,其偏振态被解调后最后由光电探测器接收。由于平面反射镜4与待测光学系统13构成了自准直光路,光束均是正入射,因此平面反射镜4对光束偏振态的改变极小,可以忽略。
在起偏光路、检偏光路之间加入传递光路,一方面突破了传统的穆勒矩阵椭偏仪的起偏臂出射光束经过待测光学系统13后直接进入检偏臂,导致可检测的样品只能是平面,而本发明提供的穆勒矩阵的测量装置可实现非平面的复杂光学系统偏振特性的测量。另一方面,传统的穆勒矩阵椭偏仪只能测量与起偏臂和检偏臂中偏振光学元件的口径相当的待测光学系统13,而本发明提供的穆勒矩阵的测量装置的传递光路通过分光组件1、第一透镜组2、小孔光阑3和平面反射镜4构成的光路,只需要将待测光学系统13的焦点与第一透镜组2焦点重合,并且平面反射镜4的口径大于待测光学系统13的有效通光口径,采用普通口径的偏振光学元件即可实现大口径复杂光学系统偏振特性的测量。
在本发明的一个优选实施例中,平面反射镜4的口径大于待测光学系统13的有效通光口径。
在本发明的一个优选实施例中,小孔光阑3的通光口径大于待测光学系统13的艾利斑直径其中R、f和D分别为待测光学系统(13)的艾利斑直径、有效焦距和有效通光口径,λ为有效光谱范围内的最大波长,在具体实施中光源5后方可以添加滤光片等装置,因此光源5的光谱范围并不代表是实际有效光谱范围。本实施例中小孔光阑3的通光口径为30μm。
在本发明的一个优选实施例中,如图2所示,分光组件1还包括分光面相互正交放置的第一分光棱镜101和第二分光棱镜102,经平面反射镜4自准直后原路返回的光束入射到第一分光棱镜101,经反射折转至第二分光棱镜102,再经第二分光棱镜102反射折转进入检偏光路,用于偏振态的解调和光强探测。第一分光棱镜101和第二分光棱镜102的型号和规格应当相同,第一分光棱镜101和第二分光棱镜102的分光面相互正交放置,使得第二分光棱镜102能够补偿第一分光棱镜101对光线偏振态的改变量,从而大幅度降低整个分光组件对光线偏振态的影响,提高系统穆勒矩阵的测量精度。
在本发明的一个优选实施例中,第一透镜组2进行消偏振设计,降低第一透镜组2对光束偏振态的改变。
在本发明的一个优选实施例中,起偏光路包括依次沿光路方向设置的光强、偏振态、波长稳定的光源5、第二透镜组6、起偏器7为格兰泰勒棱镜消光比为100000:1、第一相位补偿器8。氦氖激光器发出光强、偏振态稳定的光束,经第二透镜组6准直后,进入起偏器7变为线偏振光束,线偏振光束进入第一相位补偿器8被调制后,入射至传递光路。常用的光源5有激光、卤素灯、氙灯、SiC棒热发光光源等。根据需要测得的光谱范围,选择相应的光源5,本实施例中采用氦氖激光器作为光源5,其波长为632.8nm。
在本发明的一个优选实施例中,检偏光路包括依次沿光路方向设置的第二相位补偿器9、检偏器10、第三透镜组11、光电探测器12。如图1、如图2所示,分光组件1分光的透射光束沿传递光路原路返回后,又经分光组件1分光反射折转进入第二相位补偿器9,其偏振态被调制后进入检偏器10,其偏振态被解调,最后被第三透镜组11汇聚后由光电探测器12接收。
本发明实施例中第一相位补偿器8和第二相位补偿器9的转速比为1:5,其相位延迟量均为128°。本发明实施例中第二透镜组6、第三透镜组11均要求在设定的光谱范围内消色差,即保证透镜组在全光谱范围内有一致的焦距。
起偏光路、检偏光路可参考现有的穆勒矩阵椭偏仪技术方案进行设计,本发明对此不进行限定。起偏光路、检偏光路中的起偏器7和检偏器10的均能将任意光束变为线偏振光。目前常用的偏振器有格兰泰勒棱镜、汤普森棱镜、金属线栅偏振器等。不同材料不同类型的偏振器都有其适用的光谱范围,因此符合光谱范围的偏振片均适用于本发明实施例,本发明实施例中起偏器和检偏器均选择格兰泰勒棱镜,其消光比高达100000:1。起偏光路、检偏光路中的第一相位补偿器8和第二相位补偿器9可以是用精密旋转电机驱动的波片,也可以是液晶、光弹调制器等元件。本发明实施例中选择石英波片,相位延迟量均为128°,通过精密电机实现高精度旋转。
另一方面,本发明还提供一种测量穆勒矩阵的方法,如图3所示,包括如下步骤:
S1、起偏光路发出的偏振光束进入传递光路,经传递光路中分光组件1后被分光,透射光束依次经第一透镜组2、小孔光阑3、待测光学系统13、平面反射镜4,由平面反射镜4自准直反射后再沿原光路返回,进入检偏光路,最后被探测器接收。
S2、通过精密电机,旋转第一相位补偿器8和第二相位补偿器9,并使得二者以恒定的转速比u:v同步旋转。
S3、重复步骤S2使用光电探测器12进行数据采集,获取频率在0倍到2×(u+v)倍范围内的光强信号,重复次数不少于16次,计算获得经待测光学系统13在自准直光路中的穆勒矩阵M。
根据图1中的光路,将两块相位延迟器旋转到不同的方位角,调制不同偏振态测得多组光强,其中获取的t时刻的第Q组光强IQ,按照穆勒矩阵-斯托克斯矢量理论,整个光路可以用穆勒矩阵描述为:
IQ=Sout(t)=P2(α2)·R2(δ2,θ2(t))·M·R1(δ1,θ1(t))·P1(α1)·Sin;
式中IQ表示第Q组光强,Sin表示光源5发出光束的斯托克斯矢量,Sout表示经检偏光路出射光的斯托克斯矢量,R1(δ1,θ1(t))·P1(α1)、P2(α2)·R2(δ2,θ2(t))分别表示起偏光路、检偏光路的穆勒矩阵,P1、P2分别表示起偏器7、检偏器10的穆勒矩阵;R1、R2分别表示第一相位补偿器8、第二相位补偿器9的穆勒矩阵;α1,θ1,θ2和α2分别表示起偏器7、第一相位补偿器8、第二相位补偿器9和检偏器10的方位角,δ1和δ2分别表示第一相位补偿器8、第二相位补偿器9的相位延迟量,M表示经待测光学系统13在自准直光路中的穆勒矩阵。
起偏光路的矩阵公式为:a=P2(α2)R2(δ2,θ2(t));检偏光路的矩阵公式为:g=R1(δ1,θ1(t))P1(α1);;
为获得待测光学系统13穆勒矩阵全部16个元素,测量时要至少调制16组独立的偏振态进行测量。光电探测器12采集的第Q次光强信号IQ表示为:
式中,aq,j表示起偏光路穆勒矩阵第一行的各个元素,mj,k表示待测光学系统13在自准直光路中的穆勒矩阵M的各个元素,gq,k表示检偏光路穆勒矩阵第一列的各个元素。
WQ=(wq,1,1,wq,1,2,wq,1,3,…,wq,4,4)=(aq,1gq,1,aq,1gq,2,aq,1gq,3,…,aq,4gq,4);
M=(m11,m12,m13,m14,…,m44)T;
式中,WQ表示第Q次测量时的仪器矩阵,T表示矩阵的转置运算。
待测光学系统在自准直光路中表现出的穆勒矩阵通过伪逆求解:
M=(WT·W)-1·WT·I。
步骤S4中计算获得待测光学系统13的穆勒矩阵M0的计算公式如下所示:
其中mj,k(j,k=1,2,3,4)表示待测光学系统13在自准直光路中的穆勒矩阵M的各个元素,M0是指光束单次经过待测光学系统13的穆勒矩阵,也是最终需要得到的测量结果。
本发明所提供的一种穆勒矩阵的测量装置及其测量方法,通过增加传递光路来实现对大口径、非平面的复杂光学系统偏振特性的测量。在该测量装置及其测量方法中,由于平面反射镜与待测光学系统构成了自准直光路,光线均为正入射到平面反射镜,减小了平面反射镜引入的测量误差,提高了穆勒矩阵测量精度。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制。本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
以上本发明的具体实施方式,并不构成对本发明保护范围的限定。任何根据本发明的技术构思所做出的各种其他相应的改变与变形,均应包含在本发明权利要求的保护范围内。
Claims (8)
1.一种穆勒矩阵的测量装置,包括起偏光路和检偏光路,其特征在于,还包括位于所述起偏光路光束出射方向的传递光路;
所述传递光路包括依次沿光路方向设置的分光组件(1)、第一透镜组(2)、位于所述第一透镜组(2)焦点处的小孔光阑(3)、位于待测光学系统(13)的光束出射端的平面反射镜(4);
所述待测光学系统(13)放置于所述小孔光阑(3)和所述平面反射镜(4)之间,所述待测光学系统(13)的焦点与所述第一透镜组(2)的焦点重合;
所述起偏光路包括依次沿光路方向设置的光源(5)、第二透镜组(6)、起偏器(7)、第一相位补偿器(8);所述光源(5)发出光强、偏振态、波长稳定的光束,经所述第二透镜组(6)准直后,经所述起偏器(7)变为线偏振光束,所述线偏振光束进入所述第一相位补偿器(8)被调制后,入射至所述传递光路;
所述检偏光路包括依次沿光路方向设置的第二相位补偿器(9)、检偏器(10)、第三透镜组(11)、接收光强信号的光电探测器(12);经所述平面反射镜(4)反射沿原光路返回的光束,经所述分光组件(1)反射折转进入所述第二相位补偿器(9),其偏振态被调制后,进入所述检偏器(10),经所述第三透镜组(11)汇聚至所述光电探测器(12)进行光强探测;
所述起偏光路发出偏振态被调制的光束,经过所述分光组件(1)分光,透射光束经第一透镜组(2)汇聚于所述待测光学系统(13)的焦点处,再经所述小孔光阑(3)入射至所述待测光学系统(13),变为平行光束后经所述平面反射镜(4)自准直反射,并沿原光路返回,依次通过所述小孔光阑(3)和所述第一透镜组(2)后,经所述分光组件(1)反射进入所述检偏光路,通过计算获得所述待测光学系统(13)的穆勒矩阵。
2.根据权利要求1所述的穆勒矩阵的测量装置,其特征在于,平面反射镜(4)的口径大于待测光学系统(13)的有效通光口径。
4.根据权利要求1所述的穆勒矩阵的测量装置,其特征在于,所述分光组件(1)包括分光面相互正交放置的第一分光棱镜(101)和第二分光棱镜(102);
经所述平面反射镜(4)反射沿原光路返回的光束经所述第一分光棱镜(101)反射折转至所述第二分光棱镜(102),再经所述第二分光棱镜(102)的反射折转进入所述检偏光路。
5.根据权利要求1所述的穆勒矩阵的测量装置,其特征在于,所述第一透镜组(2)经消偏振处理,以降低所述第一透镜组(2)对光束偏振态的改变。
6.一种穆勒矩阵测量的方法,其特征在于,包括如下步骤:
S1、起偏光路发出的偏振光束进入传递光路;所述偏振光束经过所述传递光路的分光组件(1)后被分光,其中透射光束依次经第一透镜组(2)、小孔光阑(3)、待测光学系统(13)、平面反射镜(4),经平面反射镜(4)自准直反射后再沿原光路返回,进入检偏光路,其偏振态被解调,最后被光电探测器接收;
S2、旋转所述起偏光路中第一相位补偿器(8)和所述检偏光路中第二相位补偿器(9),使第一相位补偿器(8)和第二相位补偿器(9)保持固定的转速比;
S3、重复步骤S2,使用光电探测器(12)进行数据采集,获取频率在0倍到2×(u+v)倍范围内的光强信号,重复次数不少于16次,计算获得经待测光学系统(13)在自准直光路中的穆勒矩阵M;
S4、根据经待测光学系统(13)在自准直光路中的穆勒矩阵M,计算获得待测光学系统(13)的穆勒矩阵M0。
7.根据权利要求6所述的穆勒矩阵测量的方法,其特征在于,
步骤S3中计算获得经待测光学系统(13)在自准直光路中的穆勒矩阵M的具体方法如下:
起偏光路的矩阵公式为:g=R1(δ1,θ1(t))P1(α1);检偏光路的矩阵公式为:a=P2(α2)R2(δ2,θ2(t));
式中,P1、P2分别表示起偏器(7)、检偏器(10)的穆勒矩阵;R1、R2分别表示第一相位补偿器(8)、第二相位补偿器(9)的穆勒矩阵;α1,θ1,θ2和α2分别表示起偏器(7)、第一相位补偿器(8)、第二相位补偿器(9)和检偏器(10)的方位角;δ1和δ2分别是第一相位补偿器(8)、第二相位补偿器(9)的相位延迟量,
光电探测器(12)采集的Q次光强信号IQ表示为:
式中,aq,j表示起偏光路穆勒矩阵第一行的各个元素,mj,k(j,k=1,2,3,4)表示待测光学系统(13)在自准直光路中的穆勒矩阵M的各个元素,gq,k表示检偏光路穆勒矩阵第一列的各个元素,
WQ=(wq,1,1,wq,1,2,wq,1,3,…,wq,4,4)=(aq,1gq,1,aq,1gq,2,aq,1gq,3,…,aq,4gq,4);
M=(m11,m12,m13,m14,…,m44)T;
其中WQ表示第Q次测量时的仪器矩阵,T表示矩阵的转置运算,
待测光学系统在自准直光路中表现出的穆勒矩阵通过伪逆求解:
M=(WT·W)-1·WT·I。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110605491.2A CN113281256B (zh) | 2021-05-31 | 2021-05-31 | 穆勒矩阵测量装置及其测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110605491.2A CN113281256B (zh) | 2021-05-31 | 2021-05-31 | 穆勒矩阵测量装置及其测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113281256A CN113281256A (zh) | 2021-08-20 |
CN113281256B true CN113281256B (zh) | 2022-06-03 |
Family
ID=77282888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110605491.2A Active CN113281256B (zh) | 2021-05-31 | 2021-05-31 | 穆勒矩阵测量装置及其测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113281256B (zh) |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1345024A1 (en) * | 2002-03-11 | 2003-09-17 | Centre National De La Recherche Scientifique | Method and device for polarimetric measurement of the Mueller matrix coefficients of a sample in the far ultraviolet to visible spectral range |
JP2014035257A (ja) * | 2012-08-08 | 2014-02-24 | National Institute Of Advanced Industrial & Technology | ミューラー行列顕微エリプソメータ |
CN104297744A (zh) * | 2014-10-16 | 2015-01-21 | 西安理工大学 | 偏振激光雷达的偏振标定与补偿装置及标定与补偿方法 |
CN104501738A (zh) * | 2014-12-31 | 2015-04-08 | 华中科技大学 | 纳米尺度下大面积散射场的快速测量方法及装置 |
CN104536146A (zh) * | 2014-12-17 | 2015-04-22 | 中国科学院长春光学精密机械与物理研究所 | 使用在光纤点衍射干涉仪波面参考源中的偏振控制方法 |
CN104781650A (zh) * | 2012-10-24 | 2015-07-15 | 科磊股份有限公司 | 用于高纵横比及大横向尺寸结构的度量系统及方法 |
WO2015125997A1 (ko) * | 2014-02-20 | 2015-08-27 | 코스맥스 주식회사 | 뮬러 매트릭스 측정에 의한 액정에멀젼의 다층상 구조 확인 방법 |
CN105157618A (zh) * | 2015-06-01 | 2015-12-16 | 北京航天控制仪器研究所 | 一种计算强度关联成像自准直仪及测量方法 |
CN205562384U (zh) * | 2016-03-25 | 2016-09-07 | 中国人民解放军理工大学 | 气溶胶穆勒矩阵测量仪的散射系数标定装置 |
CN106597414A (zh) * | 2016-10-18 | 2017-04-26 | 浙江大学 | 一种定标偏振激光雷达增益比的方法 |
CN108918425A (zh) * | 2018-06-01 | 2018-11-30 | 中国科学院西安光学精密机械研究所 | 一种穆勒矩阵测量系统及方法 |
CN109029932A (zh) * | 2018-08-03 | 2018-12-18 | 长春理工大学 | 一种折反式望远系统的全口径全视场光线追迹方法 |
CN109164048A (zh) * | 2018-09-18 | 2019-01-08 | 天津大学 | 一种对导管偏振敏感光学相干层析成像偏振解调方法 |
CN109901281A (zh) * | 2019-03-12 | 2019-06-18 | 清华大学深圳研究生院 | 一种穆勒显微偏振像差校准方法及装置 |
CN112229799A (zh) * | 2020-09-21 | 2021-01-15 | 华南师范大学 | 偏振光声穆勒矩阵成像方法和系统 |
US10901241B1 (en) * | 2018-03-14 | 2021-01-26 | Onto Innovation Inc. | Optical metrology system using infrared wavelengths |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7202950B2 (en) * | 2003-07-08 | 2007-04-10 | Marine Biological Laboratory | Retardance measurement system and method |
TWI352878B (en) * | 2005-06-13 | 2011-11-21 | Asml Netherlands Bv | Lithographic device, and method |
WO2013042679A1 (ja) * | 2011-09-19 | 2013-03-28 | 株式会社ニコン | 照明光学装置、光学系ユニット、照明方法、並びに露光方法及び装置 |
US9952140B2 (en) * | 2012-05-29 | 2018-04-24 | Kla-Tencor Corporation | Small spot size spectroscopic ellipsometer |
IL258195B (en) * | 2018-03-18 | 2019-03-31 | Technion Res & Development Found Ltd | A system and method for monitoring the condition of an object |
-
2021
- 2021-05-31 CN CN202110605491.2A patent/CN113281256B/zh active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1345024A1 (en) * | 2002-03-11 | 2003-09-17 | Centre National De La Recherche Scientifique | Method and device for polarimetric measurement of the Mueller matrix coefficients of a sample in the far ultraviolet to visible spectral range |
JP2014035257A (ja) * | 2012-08-08 | 2014-02-24 | National Institute Of Advanced Industrial & Technology | ミューラー行列顕微エリプソメータ |
CN104781650A (zh) * | 2012-10-24 | 2015-07-15 | 科磊股份有限公司 | 用于高纵横比及大横向尺寸结构的度量系统及方法 |
WO2015125997A1 (ko) * | 2014-02-20 | 2015-08-27 | 코스맥스 주식회사 | 뮬러 매트릭스 측정에 의한 액정에멀젼의 다층상 구조 확인 방법 |
CN104297744A (zh) * | 2014-10-16 | 2015-01-21 | 西安理工大学 | 偏振激光雷达的偏振标定与补偿装置及标定与补偿方法 |
CN104536146A (zh) * | 2014-12-17 | 2015-04-22 | 中国科学院长春光学精密机械与物理研究所 | 使用在光纤点衍射干涉仪波面参考源中的偏振控制方法 |
CN104501738A (zh) * | 2014-12-31 | 2015-04-08 | 华中科技大学 | 纳米尺度下大面积散射场的快速测量方法及装置 |
CN105157618A (zh) * | 2015-06-01 | 2015-12-16 | 北京航天控制仪器研究所 | 一种计算强度关联成像自准直仪及测量方法 |
CN205562384U (zh) * | 2016-03-25 | 2016-09-07 | 中国人民解放军理工大学 | 气溶胶穆勒矩阵测量仪的散射系数标定装置 |
CN106597414A (zh) * | 2016-10-18 | 2017-04-26 | 浙江大学 | 一种定标偏振激光雷达增益比的方法 |
US10901241B1 (en) * | 2018-03-14 | 2021-01-26 | Onto Innovation Inc. | Optical metrology system using infrared wavelengths |
CN108918425A (zh) * | 2018-06-01 | 2018-11-30 | 中国科学院西安光学精密机械研究所 | 一种穆勒矩阵测量系统及方法 |
CN109029932A (zh) * | 2018-08-03 | 2018-12-18 | 长春理工大学 | 一种折反式望远系统的全口径全视场光线追迹方法 |
CN109164048A (zh) * | 2018-09-18 | 2019-01-08 | 天津大学 | 一种对导管偏振敏感光学相干层析成像偏振解调方法 |
CN109901281A (zh) * | 2019-03-12 | 2019-06-18 | 清华大学深圳研究生院 | 一种穆勒显微偏振像差校准方法及装置 |
CN112229799A (zh) * | 2020-09-21 | 2021-01-15 | 华南师范大学 | 偏振光声穆勒矩阵成像方法和系统 |
Non-Patent Citations (3)
Title |
---|
Mueller-matrix ellipsometry:a review;Rasheed M.A.Azzam;《PROCEEDINGS OF SPIE》;19971003;第399-409页 * |
基于穆勒矩阵的偏振显微镜及其在生物医学领域的应用;王晔等;《世界复合医学》;20150315(第01期);第76-80页 * |
穆勒矩阵图像的获取及处理;杨蔚等;《红外与激光工程》;20151225(第12期);第331-336页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113281256A (zh) | 2021-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5581350A (en) | Method and system for calibrating an ellipsometer | |
CN105929382B (zh) | 一种主动光电系统的收发同轴辅助光校装置及方法 | |
CN109580177B (zh) | 机载三光轴一致性测试组件、系统及测试方法 | |
CN107727008B (zh) | 一种测量主动光电系统收发同轴的装置及方法 | |
CN108548658B (zh) | 一种单层膜光学元件应力和光学损耗同时测量的方法 | |
CN110146259A (zh) | 一种大口径离轴反射式多光轴一致性定量测试和校准装置 | |
CN103162831B (zh) | 宽带偏振光谱仪及光学测量系统 | |
CN103048047B (zh) | 包含相位元件的垂直入射宽带偏振光谱仪和光学测量系统 | |
CN110806266A (zh) | 一种偏振态检测系统中偏振态分析器的选择方法 | |
EP1666871A2 (en) | Swept-angle measurement system | |
CN109115690A (zh) | 实时偏振敏感的太赫兹时域椭偏仪及光学常数测量方法 | |
CN107941477B (zh) | 一种能精确控制入射角的分光镜测量方法及装置 | |
CN107782697B (zh) | 宽波段共焦红外透镜元件折射率测量方法与装置 | |
CN107727368B (zh) | 一种标定平行光管焦面位置的装置及方法 | |
CN209927419U (zh) | 一种大口径离轴反射式多光轴一致性定量测试和校准装置 | |
CN108731650A (zh) | 一种标定具备自检功能的激光发射系统光轴的装置及方法 | |
CN110488252B (zh) | 一种地基气溶胶激光雷达系统的重叠因子定标装置和标定方法 | |
CN105758625A (zh) | 一种测量遥感仪器的线偏振灵敏度的装置及方法 | |
CN114812889A (zh) | 一种大口径光学元件应力检测装置及其检测方法 | |
CN209043571U (zh) | 机载三光轴一致性测试组件及测试系统 | |
CN108801465B (zh) | 一种激光偏振态测量装置及其测量方法 | |
US8314992B1 (en) | Field widening lens | |
CN113281256B (zh) | 穆勒矩阵测量装置及其测量方法 | |
CN106352985B (zh) | 一种非对称空间外差光谱仪结构 | |
US4171910A (en) | Retroreflectance measurement system |
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 |