CN105043294B - 适用于纳米级三维形貌测量的远场矢量光学特性建模方法 - Google Patents

适用于纳米级三维形貌测量的远场矢量光学特性建模方法 Download PDF

Info

Publication number
CN105043294B
CN105043294B CN201510364705.6A CN201510364705A CN105043294B CN 105043294 B CN105043294 B CN 105043294B CN 201510364705 A CN201510364705 A CN 201510364705A CN 105043294 B CN105043294 B CN 105043294B
Authority
CN
China
Prior art keywords
mrow
msup
vector
msubsup
msub
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
Application number
CN201510364705.6A
Other languages
English (en)
Other versions
CN105043294A (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 CN201510364705.6A priority Critical patent/CN105043294B/zh
Publication of CN105043294A publication Critical patent/CN105043294A/zh
Application granted granted Critical
Publication of CN105043294B publication Critical patent/CN105043294B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

本发明公开了一种适用于纳米级三维形貌测量的远场矢量光学特性建模方法,包括:获得探测样品近场处的矢量电磁场分布,并经由采用高NA物镜的偏振光学系统执行传播,由此依次执行近场至入瞳、入瞳至出瞳,以及出瞳至探测平面的传播过程;将近场处的矢量电磁场分布转换为入瞳处的矢量电场分布,接着将入瞳处的矢量电场分布相应转变为出瞳处的矢量电场分布;计算得出最终探测平面的矢量电场分布。以上过程还可包括根据琼斯矩阵理论,获得探测样品在不同入射角下的穆勒矩阵分布的步骤。通过本发明,能够以便于操控、高灵敏度和高测量精度的方式实现对纳米级三维形貌特征的测量,并尤其适用于微电子集成电路或微机电系统之类的应用场合。

Description

适用于纳米级三维形貌测量的远场矢量光学特性建模方法
技术领域
本发明属于散射光学测量技术领域,更具体地,涉及一种适用于纳米级三维形貌测量的远场矢量光学特性建模方法。
背景技术
近年来,传统的微电子集成电路(IC)与微机电系统(MEMS)加工从微米量级突破到纳米量级。随着加工尺寸的不断减小,其三维形貌参数对器件最终性能的影响也越来越显著。这些三维形貌参数不仅包括特征线宽(即关键尺寸)、周期间距、高度、侧壁角等轮廓参数,而且包含线宽粗糙度(LWR)、线边粗糙度(LER)等重要特征。由于三维形貌参数是IC制造中影响器件性能的主要特征参数,因此对三维形貌参数的准确测量及获得成为了IC制造中的关键环节。相应地,光学散射仪目前已成为IC制造工艺线上不可或缺的一种测量设备,可以实现小至22nm技术节点的关键尺寸测量。
对于基于远场的纳米结构三维形貌测量技术而言,其成功与否主要取决于正向光学特性建模和逆向参数提取这两个方面的有效性。其中由于正向光学特性建模在逆求中会被多次调用,所以快速准确的正向光学特性建模对于纳米结构三维形貌参数的准确重构起着至关重要的作用。对于穆勒矩阵椭偏仪和光学衍射层析仪这类典型的远场光学散射测量仪器而言,现有技术中已经对其正向光学特性建模方式作出了一些研究:例如,美国托莱多大学的柯林斯等人(R.W.Collins,Joohyun Koh.J.Opt.Soc.Am.A,Vol.16(8),pp.1997-2006,1999)公开了一种通过数值计算方法获得偏振光源与探测样品相互作用的0级衍射光谱,并在此基础上实现勒矩阵椭偏仪的正向光学特性建模;法国的菲涅尔研究所的YI RUAN等人(YI RUAN.3D digital imaging with tomographic diffractive microscopy,PHDthesis,Institute Fresnel UMR 7249,2012),指出可通过计算偏振光源与探测样品相互作用之后到高NA投影物镜入瞳处的散射电磁场分布,由此实现对光学衍射层析仪的正向光学特性建模。
然而,进一步的研究表面,上述现有技术仍然存在以下的缺陷或不足:首先,其并未对探测样品近场处的电磁场分布传播到最终的探测平面的整个过程及其传播机理作出更为深入的研究,尤其是缺乏更为严格和精确的矢量光学远场建模系统的整体设计,并导致最终探测样品的三维形貌重构信息不足及丢失;其次,并未采用穆勒矩阵这一包含样品更多的偏振信息,上述现有技术中采用的设备往往是只获得0级衍射光的信息或者只获得到高NA物镜入瞳处的电场分布而没有考虑高NA物镜的一系列因素等,相应导致可获得的测量信息相对偏少;基于以上分析,对纳米级三维表面形貌检测领域而言,有必要作出进一步的研究和改进,以便更好地适用于微电子集成电路和微机电系统之类的高精度应用场合。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种适用于纳米级三维形貌测量的远场矢量光学特性建模方法,其中通过结合纳米产品三维形貌测量与穆勒矩阵椭偏仪自身的特点,基于散射光学原理将整个建模处理过程划分为三个阶段并对其处理算法进行改进和设计,测试表明能够构建适用性更强的远场矢量光学特性建模体系,以便于操控、高灵敏度和高测量精度的方式实现对纳米级三维形貌特征的测量,并尤其适用于微电子集成电路或微机电系统之类的测量及重构应用场合。
为实现上述目的,按照本发明,提供了一种适用于纳米级三维形貌测量的远场矢量光学特性建模方法,其特征在于,该建模方法包括下列步骤:
(a)首先获得探测样品近场处的矢量电场分布E(r′)和矢量磁场分布H(r′);将上述矢量电场分布E(r′)和矢量磁场分布H(r′)经由具备高NA物镜的偏振光学系统实现传播,由此依次执行探测样品近场处至偏振光学系统的入瞳、该入瞳至偏振光学系统的出瞳,以及该出瞳至最终的探测平面的传播过程;其中所述高NA物镜的数值孔径NA被设定为0.9~1.0;
(b)基于以下转换公式计算得出所述矢量电场分布E(r′)和矢量磁场分布H(r′)传播到所述入瞳处的矢量电场分布E(r):
其中,S表示探测样品的积分表面;i、j均表示虚数单位;T表示电磁波在所述偏振光学系统中的传播周期,则表示相应的角频率;ε表示探测样品的介电常数;μ表示探测样品的磁导率;n表示探测样品表面上的外法线矢量;r′表示探测样品表面上任意一点所对应的矢径,r则表示入瞳空间中任意一点到探测表面的矢径;e表示自然指数;表示真空中的波数;R=|r-r′|也即矢径差的模,也即单位矢径差,其等于矢径差r-r′与矢径差的模|r-r′|之间的比值;I表示3x3的单位矩阵;表示与所述单位矢径差相对应的矩阵Kronecter积;·表示数学上的点乘法且所获得的结果为标量的形式,×则表示数学上的差乘且所获得的结果为矢量的形式;▽表示散度算子;
(c)将步骤(b)所计算得出的所述入瞳处的矢量电场分布E(r)相应转变为出瞳处的矢量电场分布,该出瞳处的矢量电场分布被表示为的矩阵形式,并且其中分别表示出瞳处的电场分布沿着笛卡尔坐标系X轴、Y轴和Z轴方向的分量;Pexit表示所述高NA物镜的出瞳的坐标,并且其在X轴、Y轴和Z轴上的坐标分量分别用来表示;kCCD表示所述高NA物镜的出瞳的波矢;
(d)基于以下Debye的矢量积分表达式,进一步计算得出所述出瞳处的矢量电场分布传播到最终的探测平面的矢量电场分布Edetector((x,y,z);kCCD):
其中,λ表示出瞳处的波长,表示真空中的波数;S1表示所述出瞳的表面;表示所述出瞳处的矢量电场分布在X轴和Y轴方向上的分量;e表示自然指数;i表示虚数单位;分别表示所述高NA物镜的出瞳在X轴、Y轴和Z轴上的坐标分量;x、y和z则分别表示笛卡尔坐标系在X轴、Y轴和Z轴方向的矢量单位。
作为进一步优选地,在步骤(d)之后,还可以包括如下的步骤:
(e1)选择特定的入射角执行以上步骤(a)至(d),获得该特定角下的矢量电场分布且采用琼斯矢量JCCD来表示,同时记录对应的光源电场分布,并且同样采用琼斯矢量Jsource来表示;
(e2)根据琼斯矢量理论,对上述探测样品近场处至偏振光学系统的入瞳、该入瞳至偏振光学系统的出瞳,以及该出瞳至最终的探测平面的整个传播过程采用琼斯矩阵JL予以模块化处理;接着,通过以下表达式计算得出琼斯矩阵JS
JCCD=JL·JS·Jsource
(e3)根据步骤(e2)所计算得出的琼斯矩阵JS,相应获得上述特定入射角所对应的仿真穆勒矩阵;
(e4)多次改变入射角,并重复上述步骤(e1)至步骤(e3)以获得各个入射角分别所对应的仿真穆勒矩阵;结合所有这些仿真穆勒矩阵与实际测量所获得的穆勒矩阵进行逆求算法处理,并将作为逆求算法处理的最优值用于探测样品三维形貌特征参数的最终检测结果。
作为进一步优选地,对于步骤(c)而言,优选还结合以下因素对所述高NA物镜从入瞳处的矢量电场分布到出瞳处的矢量电场分布执行统一处理:高NA物镜的衍射效应H(Pentrance)、出瞳球面电场与入瞳球面电场之间的映射关系E(Pentrance)、入瞳的矢量光线通过高NA物镜之后的方向变化ψ(Pentrance),以及高NA物镜自身的琼斯光瞳对入射偏振光所产生的偏振改变量J(Pentrance);并且优选采用以下表达式来执行上述统一处理的操作:
其中,表示所述出瞳处的矢量电场分布,并且分别表示出瞳处的矢量电场分布沿着笛卡尔坐标系X轴、Y轴和Z轴方向的分量,kCCD表示所述高NA物镜的出瞳的波矢;分别表示入瞳处的矢量电场分布沿着笛卡尔坐标系X轴和Y轴方向的分量,并且ks表示所述高NA物镜的入瞳的波矢。
作为进一步优选地,所述探测样品优选为微电子集成电路或者微机电系统。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,主要具备以下的技术优点:
1、通过结合微电子集成电路之类探测对象自身的特点来设计其正向光学特性建模工艺并对其关键处理步骤进行改进,相应能够构建更为严格和精确的整体矢量光学远场建模体系,并与现有技术相比显著改善其信息不足和丢失现象,同时具备便于操控、精度高和适用性强等优点,因而尤其适用于微电子集成电路或微机电系统之类的应用场合;
2、通过采用琼斯矢量理论,并基于入射角的改变来获得多个仿真穆勒矩阵,测试表明能够获得更多的测量信息,显著提高探测样品三维形貌特征参数检测结果的精度;
3、本发明中充分考虑了对高NA物镜产生影响的一系列因素,并对其相互关系及其统一处理算法进行设计,相应能够进一步提高检测的最终精度,同时具备流程实现简便的优点。
附图说明
图1是按照本发明所构建的远场矢量光学特性建模方法的整体工艺流程图;
图2是按照本发明采用高NA物镜的偏振光学系统执行整个传播过程的示意图;
图3示范性显示了采用穆勒矩阵椭偏仪实现等效经典控制理论的模型示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
图1是按照本发明所构建的远场矢量光学特性建模方法的整体工艺流程图。如图1中所示,该工艺流程主要包括以下的处理步骤:
首先,在步骤一中,可以采用各种常规测量仪器,获得探测样品近场处的矢量电场分布E(r′)和矢量磁场分布H(r′)。如同背景技术部分所指出地,对于孤立纳米结构的正向光学特性建模,现有技术中除了可以采用传统的数值计算方法,如有限元法、有限差分法、有限体积法、边界元法等方法进行求解之外,还譬如发展处了适用于光学衍射层析仪的耦合双极子法、表面积分法等方法。
在本发明中,如图2中所示,选择将上述矢量电场分布E(r′)和矢量磁场分布H(r′)经由采用高NA物镜(本领域中一般将数值孔径NA为0.9~1.0的称之为高NA物镜)的偏振光学系统实现传播,同时将探测样品近场处的电磁场分布传播到远场的探测平面的远场过程进一步分为3个部分:探测样品近场处的矢量电磁场分布至偏振光学系统也即高NA物镜的入瞳处的矢量电场分布的传播、高NA物镜将入瞳处的矢量电场分布转变成出瞳处的矢量电场分布,以及出瞳处的矢量电场分布传播到最终的探测平面的矢量电场分布的传播,由此依次执行近场至入瞳、入瞳至出瞳,以及出瞳至探测平面的整个传播过程。
接着,在步骤二中,可以利用基于麦克斯韦方程组和矢量格林定理的弗朗兹(Franz)公式(1),来初步反映探测样品近场处的矢量电磁场分布传播到所述入瞳处的矢量电场分布:
考虑到在采用高NA物镜的偏振光学系统中,矢量电磁场的远场传播会发生耦合作用,为此,我们将公式(1)中的标量格林函数用并矢格林函数代替,并矢格林函数如公式(2)示:
结合以上公式,即可获得更为准确地反映探测样品近场处的矢量电磁场分布传播到所述入瞳处的矢量电场分布;换而言之,本发明中优选基于以下转换公式计算得出所述矢量电场分布E(r′)和矢量磁场分布H(r′)传播到所述入瞳处的矢量电场分布E(r):
其中,S表示探测样品的积分表面;i、j均表示虚数单位;T表示电磁波在所述偏振光学系统中的传播周期,则表示相应的角频率;ε表示探测样品的介电常数;μ表示探测样品的磁导率;n表示探测样品表面上的外法线矢量;r′表示探测样品表面上任意一点所对应的矢径,r则表示入瞳空间中任意一点到探测表面的矢径;e表示自然指数;表示真空中的波数;R=|r-r′|也即矢径差的模,也即单位矢径差,其等于矢径差r-r′与矢径差的模|r-r′|之间的比值;I表示3x3的单位矩阵;表示与所述单位矢径差相对应的矩阵的Kronecter积;·表示数学上的点乘法且所获得的结果为标量的形式,×则表示数学上的差乘且所获得的结果为矢量的形式;▽表示散度算子。此外,考虑到实际场合下磁场作用会远远小于电场作用,对最终测量结果并无明显影响,因此在本发明的以上处理步骤和转换公式中仅考虑光学偏振系统也即高NA物镜的入瞳处的矢量电场分布E(r),同时忽略对其矢量磁场分布的计算。
接着,在步骤三中,将以上步骤计算得出的所述入瞳处的矢量电场分布E(r)相应转变为出瞳处的矢量电场分布,该出瞳处的矢量电场分布被表示为的矩阵形式,并且其中分别表示出瞳处的电场分布沿着笛卡尔坐标系X轴、Y轴和Z轴方向的分量;Pexit表示所述高NA物镜的出瞳的坐标,并且其在X轴、Y轴和Z轴上的坐标分量分别用来表示;kCCD表示所述高NA物镜的出瞳的波矢。该步骤的具体处理过程和计算为本领域的技术人员所熟知,因此在此不再赘述。
与此同时,除了常规的入瞳矢量电场分布至出瞳矢量电场分布的转变计算方式之外,按照本发明的一个优选实施例,还可以进一步考虑高NA物镜对矢量电场分布改变的一系列因素并对其进行适当处理,以便适当提高测量精度。
具体而言,由于高NA物镜实质上是一低通滤波器,因此可以考虑将物镜的衍射效应纳入影响因素之一。关于物镜的衍射效应,可以用关于入瞳光瞳坐标的透过率函数来表征,如公式(4):
其次,对于采用高NA物镜的偏振光学系统,出瞳球面电场和入瞳球面电场之间的映射关系需要满足能量守恒和正弦条件,对于两者之间的映射关系可以采用入瞳光瞳坐标的倾斜因子函数来表征,如式(5):
其次,由于高NA物镜会对入射的偏振光的方向发生改变,所以可以采用关于入瞳光瞳坐标为函数的转换矩阵来表征矢量光线通过高NA物镜的变化情况。
再次,还可以进一步考虑高NA物镜自身的琼斯光瞳对入射偏振光所产生的偏振改变量J(Pentrance)。
在以上分析的基础上,按照本发明的优选实施例,可以将入瞳处的矢量电磁场分布通过高NA物镜转变成出瞳的矢量电磁场分布,并且用数学乘积的形式联系起来,由此既充分考虑了以上影响因素,并且针对这些影响因素对高NA物镜从入瞳处的矢量电场分布到出瞳处的矢量电场分布执行统一处理:
其中,表示所述出瞳处的矢量电场分布,并且分别表示出瞳处的矢量电场分布沿着笛卡尔坐标系X轴、Y轴和Z轴方向的分量,kCCD表示所述高NA物镜的出瞳的波矢;分别表示入瞳处的矢量电场分布沿着笛卡尔坐标系X轴和Y轴方向的分量,并且ks表示所述高NA物镜的入瞳的波矢。
最后,在步骤四中,基于以下Debye的矢量积分表达式,进一步计算得出所述出瞳处的矢量电场分布传播到最终的探测平面的矢量电场分布Edetector((x,y,z);kCCD):
其中,λ表示出瞳处的波长,表示真空中的波数;S1表示所述出瞳的表面;表示所述出瞳处的矢量电场分布在X轴和Y轴方向上的分量;e表示自然指数;i表示虚数单位;分别表示所述高NA物镜的出瞳在X轴、Y轴和Z轴上的坐标分量;x、y和z则分别表示笛卡尔坐标系在X轴、Y轴和Z轴方向的矢量单位;以此方式,能够以便于操控、高灵敏度和高测量精度的方式实现远场矢量光学特性的建模过程,并据此完成对纳米级三维形貌特征的高精度测量。
此外,在上述建模体系完成之后,可以获得最终探测平面的严格矢量电场分布。为了能够与实际测量所获得的穆勒矩阵分布更好地进行匹配,按照本发明的另一优选实施例,还可以继续执行如下的步骤:
(1)选择特定的入射角执行以上光学特性建模过程,获得该特定角下的矢量电场分布且采用琼斯矢量JCCD来表示,同时记录对应的光源电场分布,并且同样采用琼斯矢量Jsource来表示;
(2)根据琼斯矢量理论,对上述探测样品近场处至偏振光学系统的入瞳、该入瞳至偏振光学系统的出瞳,以及该出瞳至最终的探测平面的整个传播过程采用琼斯矩阵JL予以模块化处理;接着,构建以下表达式计算得出琼斯矩阵JS
JCCD=JL·JS·Jsource (9)
(3)根据以上所计算得出的琼斯矩阵JS,相应获得上述特定入射角所对应的仿真穆勒矩阵;
(4)多次改变入射角,并重复上述步骤以获得各个入射角分别所对应的仿真穆勒矩阵;结合所有这些仿真穆勒矩阵与实际测量所获得的穆勒矩阵进行逆求算法处理,并将作为逆求算法处理的最优值用于探测样品三维形貌特征参数的最终检测结果。
更具体而言,参看图3,示范性显示了采用穆勒矩阵椭偏仪实现等效经典控制理论的模型示意图。如图3中所示,主要包括激光光源1、扫描振镜2、扩束器3,起偏器4、第一增量式编码器5、第一伺服电机6、第一补偿器7、样品台8、高NA物镜9、第一透镜10、分束器11、第二透镜12、第二补偿器13、第二伺服电机14、第二增量式编码器15、检偏器16、图像采集器17。借助经典控制理论模型,我们将散射场层析穆勒矩阵椭偏仪分为4个部分:第一部分:从光源到探测样品表面包括激光光源1、扫描振镜2、扩束器3,起偏器4、第一增量式编码器5、第一伺服电机6、第一补偿器7、高NA物镜9、第一透镜10;第二部分:探测样品包括样品台8;第三部分:从探测样品表面到最终的探测平面包括高NA物镜9、第一透镜10、分束器11、第二透镜12、第二补偿器13、第二伺服电机14、第二增量式编码器15、检偏器16;第四部分:电荷耦合器件(CCD)或互补金属氧化物半导体(CMOS)图像传感器包括图像采集器17。
在此基础上,可以借助琼斯矢量理论,将第一部分的光源和第四部分的图像传感器处的电场均用一琼斯矢量来表示:Jsource,JCCD,第二部分和第三部分同样均用一琼斯矩阵来表示:Js,JL,这4个部分之间的关系可以用下面的公式(10)来表示:
JCCD=JL·JS·Jsource (10)
为了求得第三部分的Jones矩阵JL,我们将采用一个已知样品的琼斯矩阵Jsample,去获取从探测样品表面到CCD平面之间的光学系统的琼斯矩阵JL,由于要获得琼斯矩阵的每个元素,我们需要两次的入射光源的琼斯矢量电磁场以及两次的CCD的琼斯矢量电磁场。因此有公式:
对于获得的第三部分的Jones矩阵JL(只与测量配置条件中的入射角,方位角有关),为了求得待测样品的Jones矩阵J,可以分别在入射光电场矢量Einc与入射面之间的TE偏振和TM偏振情况下计算。当TM偏振时,电场分量ES=0,相应可得:
同理,当TE偏振时,电场分量EP=0,因此有公式:
对于非退偏样品,其Jones矩阵J和Mueller矩阵M之间存在如下的关系:
其中,表示矩阵的Kronecter积,J*为Jones矩阵J的共轭矩阵,矩阵U定义为以下公式(14):
在实际测量过程中,通常用Mueller矩阵中第一行和第一列的M11元素对其进行归一化处理。归一化的Mueller矩阵定义为以下公式(15):
综上所述,本发明中通过结合纳米产品三维形貌测量与穆勒矩阵椭偏仪自身的特点,基于散射光学原理将整个建模处理过程划分为三个阶段并对其处理算法进行改进和设计,测试表明能够构建适用性更强的远场矢量光学特性建模体系;相应地,与现有技术相比显著改善其信息不足和丢失现象,同时具备便于操控、精度高和适用性强等优点,因而尤其适用于微电子集成电路或微机电系统之类的应用场合。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种适用于纳米级三维形貌测量的远场矢量光学特性建模方法,其特征在于,该建模方法包括下列步骤:
(a)首先获得探测样品近场处的矢量电场分布E(r′)和矢量磁场分布H(r′);将上述矢量电场分布E(r′)和矢量磁场分布H(r′)经由采用高NA物镜的偏振光学系统实现传播,由此依次执行探测样品近场处至偏振光学系统的入瞳、该入瞳至偏振光学系统的出瞳,以及该出瞳至最终的探测平面的传播过程;其中所述高NA物镜的数值孔径NA被设定为0.9~1.0;
(b)基于以下转换公式计算得出所述矢量电场分布E(r′)和矢量磁场分布H(r′)传播到所述入瞳处的矢量电场分布E(r):
<mrow> <mi>E</mi> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <mi>S</mi> </munder> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mi>j</mi> <mi>&amp;omega;</mi> <mi>&amp;epsiv;</mi> <mi>&amp;mu;</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>&amp;times;</mo> <mi>H</mi> <mrow> <mo>(</mo> <msup> <mi>r</mi> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <mfrac> <msup> <mi>e</mi> <mrow> <msub> <mi>ik</mi> <mn>0</mn> </msub> <mi>R</mi> </mrow> </msup> <mrow> <mn>4</mn> <msubsup> <mi>&amp;pi;k</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mrow> <mo>(</mo> <mrow> <mn>3</mn> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>&amp;CircleTimes;</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>-</mo> <mi>I</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <mfrac> <mn>1</mn> <msup> <mi>R</mi> <mn>3</mn> </msup> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>ik</mi> <mn>0</mn> </msub> </mrow> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <mi>I</mi> <mo>-</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>&amp;CircleTimes;</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> </mrow> <mo>)</mo> </mrow> <mfrac> <msubsup> <mi>k</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> <mo>)</mo> </mrow> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <mrow> <mi>j</mi> <mi>&amp;omega;</mi> <mi>&amp;epsiv;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>&amp;times;</mo> <mi>H</mi> <mrow> <mo>(</mo> <msup> <mi>r</mi> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mo>&amp;dtri;</mo> </mrow> <mo>)</mo> </mrow> <mo>&amp;dtri;</mo> <mrow> <mo>(</mo> <mrow> <mfrac> <msup> <mi>e</mi> <mrow> <msub> <mi>ik</mi> <mn>0</mn> </msub> <mi>R</mi> </mrow> </msup> <mrow> <mn>4</mn> <msubsup> <mi>&amp;pi;k</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mrow> <mo>(</mo> <mrow> <mn>3</mn> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>&amp;CircleTimes;</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>-</mo> <mi>I</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <mfrac> <mn>1</mn> <msup> <mi>R</mi> <mn>3</mn> </msup> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>ik</mi> <mn>0</mn> </msub> </mrow> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <mi>I</mi> <mo>-</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>&amp;CircleTimes;</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> </mrow> <mo>)</mo> </mrow> <mfrac> <msubsup> <mi>k</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>&amp;times;</mo> <mi>E</mi> <mrow> <mo>(</mo> <msup> <mi>r</mi> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mo>&amp;dtri;</mo> <mrow> <mo>(</mo> <mrow> <mfrac> <msup> <mi>e</mi> <mrow> <msub> <mi>ik</mi> <mn>0</mn> </msub> <mi>R</mi> </mrow> </msup> <mrow> <mn>4</mn> <msubsup> <mi>&amp;pi;k</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mrow> <mo>(</mo> <mrow> <mn>3</mn> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>&amp;CircleTimes;</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>-</mo> <mi>I</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <mfrac> <mn>1</mn> <msup> <mi>R</mi> <mn>3</mn> </msup> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>ik</mi> <mn>0</mn> </msub> </mrow> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <mi>I</mi> <mo>-</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>&amp;CircleTimes;</mo> <mover> <mi>R</mi> <mo>^</mo> </mover> </mrow> <mo>)</mo> </mrow> <mfrac> <msubsup> <mi>k</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mi>d</mi> <mi>S</mi> </mrow>
其中,S表示探测样品的积分表面;i、j均表示虚数单位;T表示电磁波在所述偏振光学系统中的传播周期,则表示相应的角频率;ε表示探测样品的介电常数;μ表示探测样品的磁导率;n表示探测样品表面上的外法线矢量;r′表示探测样品表面上任意一点所对应的矢径,r则表示入瞳空间中任意一点到探测表面的矢径;e表示自然指数;表示真空中的波数;R=|r-r′|也即矢径差的模,也即单位矢径差,其等于矢径差r-r′与矢径差的模|r-r′|之间的比值;I表示3x3的单位矩阵;表示与所述单位矢径差相对应的矩阵Kronecter积;·表示数学上的点乘法且所获得的结果为标量的形式,×则表示数学上的差乘且所获得的结果为矢量的形式;▽表示散度算子;
(c)将步骤(b)所计算得出的所述入瞳处的矢量电场分布E(r)相应转变为出瞳处的矢量电场分布,该出瞳处的矢量电场分布被表示为的矩阵形式,并且其中分别表示出瞳处的电场分布沿着笛卡尔坐标系X轴、Y轴和Z轴方向的分量;Pexit表示所述高NA物镜的出瞳的坐标,并且其在X轴、Y轴和Z轴上的坐标分量分别用来表示;kCCD表示所述高NA物镜的出瞳的波矢;
(d)基于以下Debye的矢量积分表达式,进一步计算得出所述出瞳处的矢量电场分布传播到最终的探测平面的矢量电场分布Edetector((x,y,z);kCCD):
<mrow> <msup> <mi>E</mi> <mrow> <mi>det</mi> <mi>e</mi> <mi>c</mi> <mi>t</mi> <mi>o</mi> <mi>r</mi> </mrow> </msup> <mrow> <mo>(</mo> <mrow> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>z</mi> </mrow> <mo>)</mo> </mrow> <mo>;</mo> <msub> <mi>k</mi> <mrow> <mi>C</mi> <mi>C</mi> <mi>D</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msup> <mi>&amp;lambda;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <msub> <mi>S</mi> <mn>1</mn> </msub> </munder> <mfrac> <mrow> <msup> <mi>E</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>p</mi> <mi>x</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>p</mi> <mi>y</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> </mrow> <mo>)</mo> </mrow> </mrow> <msubsup> <mi>p</mi> <mi>z</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> </mfrac> <msup> <mi>e</mi> <mrow> <msub> <mi>ik</mi> <mn>0</mn> </msub> <mrow> <mo>{</mo> <mrow> <msubsup> <mi>p</mi> <mi>x</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> <mi>x</mi> <mo>+</mo> <msubsup> <mi>p</mi> <mi>y</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> <mi>y</mi> <mo>+</mo> <msubsup> <mi>p</mi> <mi>z</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> <mi>z</mi> </mrow> <mo>}</mo> </mrow> </mrow> </msup> <msub> <mi>d</mi> <msubsup> <mi>p</mi> <mi>x</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> </msub> <msub> <mi>d</mi> <msubsup> <mi>p</mi> <mi>y</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> </msub> </mrow>
其中,λ表示出瞳处的波长,表示真空中的波数;S1表示所述出瞳的表面;表示所述出瞳处的矢量电场分布在X轴和Y轴方向上的分量;e表示自然指数;i表示虚数单位;分别表示所述高NA物镜的出瞳在X轴、Y轴和Z轴上的坐标分量;x、y和z则分别表示笛卡尔坐标系在X轴、Y轴和Z轴方向的矢量单位。
2.如权利要求1所述的远场矢量光学特性建模方法,其特征在于,在步骤(d)之后,还包括如下的步骤:
(e1)选择特定的入射角执行以上步骤(a)至(d),获得该特定角下的矢量电场分布且采用琼斯矢量JCCD来表示,同时记录对应的光源电场分布,并且同样采用琼斯矢量Jsource来表示;
(e2)根据琼斯矢量理论,对上述探测样品近场处至偏振光学系统的入瞳、该入瞳至偏振光学系统的出瞳,以及该出瞳至最终的探测平面的整个传播过程采用琼斯矩阵JL予以模块化处理;接着,通过以下表达式计算得出琼斯矩阵JS
JCCD=JL·JS·Jsource
(e3)根据步骤(e2)所计算得出的琼斯矩阵JS,相应获得上述特定入射角所对应的仿真穆勒矩阵;
(e4)多次改变入射角,并重复上述步骤(e1)至步骤(e3)以获得各个入射角分别所对应的仿真穆勒矩阵;结合所有这些仿真穆勒矩阵与实际测量所获得的穆勒矩阵进行逆求算法处理,并将作为逆求算法处理的最优值用于探测样品三维形貌特征参数的最终检测结果。
3.如权利要求1或2所述的远场矢量光学特性建模方法,其特征在于,对于步骤(c)而言,还结合以下因素对所述高NA物镜从入瞳处的矢量电场分布到出瞳处的矢量电场分布执行统一处理:高NA物镜的衍射效应出瞳球面电场与入瞳球面电场之间的映射关系E(Pentrance)、入瞳的矢量光线通过高NA物镜之后的方向变化ψ(Pentrance),以及高NA物镜自身的琼斯光瞳对入射偏振光所产生的偏振改变量J(Pentrance);并且采用以下表达式来执行上述统一处理的操作:
<mrow> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>E</mi> <mi>x</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mrow> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msup> <mo>;</mo> <msub> <mi>k</mi> <mrow> <mi>C</mi> <mi>C</mi> <mi>D</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>E</mi> <mi>y</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mrow> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msup> <mo>;</mo> <msub> <mi>k</mi> <mrow> <mi>C</mi> <mi>C</mi> <mi>D</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>E</mi> <mi>z</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mrow> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>x</mi> <mi>i</mi> <mi>t</mi> </mrow> </msup> <mo>;</mo> <msub> <mi>k</mi> <mrow> <mi>C</mi> <mi>C</mi> <mi>D</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mi>&amp;psi;</mi> <mrow> <mo>(</mo> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>c</mi> <mi>e</mi> </mrow> </msup> <mo>)</mo> </mrow> <mi>J</mi> <mrow> <mo>(</mo> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>c</mi> <mi>e</mi> </mrow> </msup> <mo>)</mo> </mrow> <mi>H</mi> <mrow> <mo>(</mo> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>c</mi> <mi>e</mi> </mrow> </msup> <mo>)</mo> </mrow> <mi>E</mi> <mrow> <mo>(</mo> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>c</mi> <mi>e</mi> </mrow> </msup> <mo>)</mo> </mrow> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>E</mi> <mi>x</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>c</mi> <mi>e</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mrow> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>c</mi> <mi>e</mi> </mrow> </msup> <mo>;</mo> <msub> <mi>k</mi> <mi>s</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>E</mi> <mi>y</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>c</mi> <mi>e</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mrow> <msup> <mi>P</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>c</mi> <mi>e</mi> </mrow> </msup> <mo>;</mo> <msub> <mi>k</mi> <mi>s</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,表示所述出瞳处的矢量电场分布,并且分别表示出瞳处的矢量电场分布沿着笛卡尔坐标系X轴、Y轴和Z轴方向的分量,kCCD表示所述高NA物镜的出瞳的波矢;分别表示入瞳处的矢量电场分布沿着笛卡尔坐标系X轴和Y轴方向的分量,并且ks表示所述高NA物镜的入瞳的波矢。
4.如权利要求1或2所述的远场矢量光学特性建模方法,其特征在于,所述探测样品为微电子集成电路或者微机电系统。
CN201510364705.6A 2015-06-26 2015-06-26 适用于纳米级三维形貌测量的远场矢量光学特性建模方法 Active CN105043294B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510364705.6A CN105043294B (zh) 2015-06-26 2015-06-26 适用于纳米级三维形貌测量的远场矢量光学特性建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510364705.6A CN105043294B (zh) 2015-06-26 2015-06-26 适用于纳米级三维形貌测量的远场矢量光学特性建模方法

Publications (2)

Publication Number Publication Date
CN105043294A CN105043294A (zh) 2015-11-11
CN105043294B true CN105043294B (zh) 2017-09-01

Family

ID=54450033

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510364705.6A Active CN105043294B (zh) 2015-06-26 2015-06-26 适用于纳米级三维形貌测量的远场矢量光学特性建模方法

Country Status (1)

Country Link
CN (1) CN105043294B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI624648B (zh) * 2017-04-18 2018-05-21 光寶電子(廣州)有限公司 光學建模方法及其電子裝置
CN111328370B (zh) * 2019-04-23 2023-05-02 合刃科技(深圳)有限公司 一种表面缺陷检测系统及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104090364A (zh) * 2014-05-16 2014-10-08 北京理工大学 一种远场光场全矢量计算的方法
CN104482878A (zh) * 2014-12-04 2015-04-01 华中科技大学 一种纳米结构几何参数大面积在线测量装置及其测量方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3670787B2 (ja) * 1997-02-28 2005-07-13 オリンパス株式会社 三次元形状計測装置
JP2007240344A (ja) * 2006-03-09 2007-09-20 Fujitsu Ltd 動的形状計測方法及び動的形状計測装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104090364A (zh) * 2014-05-16 2014-10-08 北京理工大学 一种远场光场全矢量计算的方法
CN104482878A (zh) * 2014-12-04 2015-04-01 华中科技大学 一种纳米结构几何参数大面积在线测量装置及其测量方法

Also Published As

Publication number Publication date
CN105043294A (zh) 2015-11-11

Similar Documents

Publication Publication Date Title
Marchesini et al. Alternating projection, ptychographic imaging and phase synchronization
Wallace et al. Phase-shifting Zernike interferometer wavefront sensor
Leach et al. Open questions in surface topography measurement: a roadmap
CN102735184B (zh) 一种光学面形的检测装置及检测方法
CN101881599A (zh) 一种纳米结构三维形貌测量方法及装置
CN102261985A (zh) 光学系统波像差标定装置及该装置测试误差的标定方法
CN104482878A (zh) 一种纳米结构几何参数大面积在线测量装置及其测量方法
CN102445280B (zh) 一种用于小孔标定的检测装置及方法
Jiang et al. Lithography alignment techniques based on moiré fringe
CN103335731A (zh) 光栅剪切波像差检测干涉仪及检测方法
KR102424819B1 (ko) 물체 형상 복원 방법
CN103674493A (zh) 光栅剪切干涉仪波像差检测的系统误差的消除方法
Roy et al. Interferometric coherent Fourier scatterometry: a method for obtaining high sensitivity in the optical inverse-grating problem
CN104165755A (zh) 光栅剪切波像差检测干涉仪及其检测方法
WO2016183874A1 (zh) 数字相移点衍射干涉仪及光学系统波像差测量方法
CN105043294B (zh) 适用于纳米级三维形貌测量的远场矢量光学特性建模方法
CN102620680A (zh) 一种三平面绝对测量光学面形的检测装置及方法
CN103759655B (zh) 基于光学倍程法的二自由度零差光栅干涉仪位移测量系统
CN104281011A (zh) 一种高数值孔径成像系统偏振像差的检测方法
Legrand et al. Three-dimensional imaging of single atoms in an optical lattice via helical point-spread-function engineering
CN105527027A (zh) 亚微米孔径光纤点衍射波前测量的结构误差校正方法
CN102768474B (zh) 基于二阶像差模型的投影物镜波像差检测方法
CN112229604A (zh) 光栅剪切干涉高精度波像差检测方法
Lin et al. Optimization and fabrication of chromium grating in self-traceable interferometer
CN113432731B (zh) 一种光栅横向剪切干涉波前重建过程中的补偿方法

Legal Events

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