CN102243065B - 一种基于背向补偿的透明基底薄膜厚度测量系统 - Google Patents

一种基于背向补偿的透明基底薄膜厚度测量系统 Download PDF

Info

Publication number
CN102243065B
CN102243065B CN2011100922833A CN201110092283A CN102243065B CN 102243065 B CN102243065 B CN 102243065B CN 2011100922833 A CN2011100922833 A CN 2011100922833A CN 201110092283 A CN201110092283 A CN 201110092283A CN 102243065 B CN102243065 B CN 102243065B
Authority
CN
China
Prior art keywords
film
reflectivity
lambda
thickness
compensation
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.)
Expired - Fee Related
Application number
CN2011100922833A
Other languages
English (en)
Other versions
CN102243065A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2011100922833A priority Critical patent/CN102243065B/zh
Publication of CN102243065A publication Critical patent/CN102243065A/zh
Application granted granted Critical
Publication of CN102243065B publication Critical patent/CN102243065B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于背向补偿的透明基底薄膜厚度测量系统,包括:光源、光谱仪、反射率光纤探头、测量平台和PC机,反射率光纤探头设于测量平台上方,分别与光源、光谱仪连接,光谱仪与PC机连接;光源发出的光经过反射率光纤探头照射到放置在测量平台上的待测具有透明基底的薄膜上,基底以及薄膜的反射光同时又进入反射率光纤探头后被光谱仪接收,接受的数据经过PC机的补偿处理计算,得到薄膜各层的实际厚度。本发明系统无须在样片背面涂覆消光物质来消除背向反射,无需破坏薄膜基底。本发明系统具有简单灵活,自动适应各种不同材料不同厚度的透明基底,整个系统具有结构简单,成本较低,小型化,无损探测等优点。

Description

一种基于背向补偿的透明基底薄膜厚度测量系统
技术领域
本发明涉及薄膜测量技术领域,尤其涉及一种基于背向补偿的透明基底薄膜厚度测量系统。
背景技术
薄膜光学是应用光学学科的一个重要分支。虽然起初发展受到了科学技术条件的限制,然而随着光谱干涉技术、激光技术及空间光学等一系列学科的快速发展,薄膜光学也得到了迅速的发展。光学薄膜可以通过特殊的结构控制光束,和产生特殊的光束,在很多领域都有着独立的应用比如镜头制造,随着科学技术的进步,光学薄膜在光电子,光通讯方面也开始发挥其重要的作用。为了获得较好的薄膜光学性能,需要精确的控制薄膜的厚度以及其它光学参数,从而使实际参数和设计参数尽可能一致。于是对于薄膜参数的控制和测量,一直是光学领域的热点问题。
薄膜厚度的测量分为两大类:无损测量和有损测量。由于有损测量会对薄膜表面产生破坏效果,所以利用光学方法的无损测量是薄膜测厚中最常采用的。两种典型的光学方法是椭偏法和光谱反射法。
椭圆偏振法,精度高,但是需要利用椭偏仪步骤繁琐且成本较高,有一定的局限性。
光谱反射法是利用光谱仪测量样品在一定波长范围内的反射率,通过改变多层薄膜各层的厚度在一定范围内进行迭代,然后将计算的反射率与实际测量的反射率进行比对,从而计算实际薄膜的厚度。理论上多层薄膜的反射率可以根据以下公式计算:
B C = { Π j = 1 K cos δ j i n j sin δ j in j sin δ j cos δ j } 1 n g - - - ( 1 )
δ j = 2 π λ n j d j - - - ( 2 )
Y = C B - - - ( 3 )
r = n 0 - Y n 0 + Y - - - ( 4 )
R=r×r*    (5)
其中λ代表波长,j代表层数,ng代表基底折射率,dj代表多层薄膜第j层的厚度,nj代表多层薄膜第j层的折射率,n0代表薄膜所在环境介质的折射率(一般为情况下为空气折射率即1.0),R代表薄膜的反射率。r*代表r的共轭值。
光谱反射法通过改变式2中的dj,根据式1-5计算理论反射率,当所计算的理论反射率和利用光谱仪测得的实际反射率一致或者基本拟合,即认为当前的dj为薄膜的实际厚度。随着薄膜层数的增加,各层厚度的组合情况增多,线性的改变各层厚度然后计算各种组合的反射率需要耗费大量时间,所以需要借助于全局优化算法在给定的厚度范围内搜索,从而获得多层薄膜的各层实际厚度。
在利用全局优化算法之前,需要建立评价函数,在薄膜厚度计算中,我们一般将实际薄膜反射率和计算薄膜反射率的改革波长均方差之和作为评价函数,该加和越小代表计算的薄膜厚度越接近实际薄膜厚度。
全局优化算法种类很多,常用的有模拟退火算法,遗传算法,神经网络算法,共轭梯度算法等,也有将各种算法进行混合,获得更高的效率。各种算法原理不同,但是基本思想一致,都是寻找到评价函数在给定范围内的最优值,并且记录获得最优值时变量的大小。各种算法的步骤也基本一致,即首先确定评价函数,然后确定迭代变量的上下限(在薄膜计算中,迭代参数为薄膜的各层厚度值),接着根据不同的原理迭代变量,每迭代一次便计算一次评价函数,然后根据评价函数的大小确定下面的迭代过程中的迭代参数。如此反复,直到满足迭代终止条件。具体过程可以参考“宫兴致,陈燕平,刘玉玲,余飞鸿,遗传算法在薄膜特性参数测量中的应用[J].光子学报,2007,36(11),2053-2056”。
综上一般的光谱反射法测量薄膜厚度的步骤为:
1.利用光谱仪测量样品反射率。
2.确定样品多层薄膜各层的厚度的上下限(不需要很精确)。
3.利用全局优化算法搜索上下限内的最优的各层薄膜厚度值(在计算反射率时利用式1-5)。
在薄膜测量过程中,有很多薄膜的基底是透明的,比如玻璃基底薄膜,但是在薄膜计算中,基底被当做不透明的,从而忽略基底背面和空气层之间所产生的反射。这些反射会和普通的薄膜反射混合在一起,被反射率探头接收。由于这种背向反射率通常很小,当薄膜反射率较高时对薄膜厚度测量影响较小。但是当薄膜的反射率本身较小的时候,这种背向反射的影响使得无法将其忽略,否则会造成极大的测量误差。通常解决方法是通过对基底背面涂覆消光材料来消除背向反射,这种方法会破坏薄膜的基底;另外有些方法对混合反射率乘以一个固定的系数从而消除背向反射,但是背向反射是一个复杂的变量,和薄膜的透明基底材料以及基底厚度都有关系,固定系数补偿缺乏通用性和准确性。
发明内容
本发明提供一种基于背向补偿的透明基底薄膜厚度测量系统,该系统结构简单灵活,自适应于各种不同材料不同厚度的透明基底,且成本低,无需破坏基底。
一种基于背向补偿的透明基底薄膜厚度测量系统,包括:光源、光谱仪、反射率光纤探头、测量平台和PC机,所述的反射率光纤探头设于测量平台上方,分别与光源、光谱仪连接,光谱仪与PC机连接;光源发出的光经过反射率光纤探头的出射端面,照射到放置在测量平台上的待测具有透明基底的薄膜上,基底以及薄膜的反射光同时又进入反射率光纤探头的接收端面,被光谱仪接收,接收的数据经过PC机的补偿处理计算,得出待测薄膜各层的厚度。
所述的光谱仪工作在可见光波段,所述的光源在整个可见光波段拥有较高的效率,光谱仪的中心波长分辨率不大于5nm。所述的反射率探头采用标准反射率探头,具有7芯,其中6芯为发射,1芯为接收。
在所述的在反射率测量过程中,基底的背向反射会影响测量结果,在一般计算,基底被当作不透明或者无限长度,从而不考虑基底的背面和空气形成的背向反射光。当薄膜反射率较高的情况下,由于背向反射较小对整体结果影响不大。当样品的反射率较低且是透明基底,背向反射光将会极大的影响反射率的测量,从而使测量反射率的和理论反射率出入较大,最后无法根据理论公式拟合反射率光谱,从而无法测量薄膜的厚度。
当光束直接照射到具有透明基底的薄膜上的时候,发生了多次反射,当基底透明,透明基底的背面和空气层之间发生反射,其示意图如图2所示。考虑相干以及非相干反射,薄膜的总反射率R可以表示为:
R = R f + T f ′ R 3 T f + T f ′ R 3 2 R ′ f T f . . . . . . . . . = R f + T f 2 R 3 1 - R 3 R ′ f ≈ R f + ( 1 - R f ) 2 R 3 - - - ( 6 )
其中Rf表示从正面到背面的相干的反射率,Rf是由薄膜的各层厚度以及折射率计算得到;R′f表示背面到正面的相干反射率,R3表示基底背面的背向反射率,Tf表示从正面到背面的相干透过率,T′f表示从背面到正面的相干透过率。
从正面到背面的相干反射率Rf为计算反射率用于计算薄膜的厚度,当薄膜基底是不透明时:R=Rf,式(6)与式(1)表达式相同。
Figure BDA0000055122220000042
称为背向反射率。如果薄膜基地透明且自身反射率较低,忽略式(6)中的背向反射率,最终厚度计算结果会受到极大的影响。
实际测量中,接收到的背向反射率与式(6)中的背向反射率差别很大。因为在理论分析中,反射率测量总是垂直于薄膜表面,在实际测量中入射光并不是完全垂直与样品,即准垂直测量。因为薄膜通常比较薄,并不会影响相干反射率Rf,大部分的反射的相干光都可以被接收器接收到。但是当基底是透明的时候,一部分有角度的背向反射无法被垂直于薄膜表面的接收器所接收到(相对于薄膜,基底通常比较厚,所以反射光偏离探头的距离较大),如图3所示。这部分没有接收到的背向反射就是薄膜的背向反射的损失。背向反射损失难以计算,和基底材料以及基底厚度都有直接的关系。也就是说,因为有了背向反射损失的存在,导致测量的反射率和理论上计算的反射率有很大的差别,从而无法按照一般方法拟合计算出实际的薄膜厚度。
如图3所示,背向反射没有完全被反射率光纤探头接收到,所以产生了背向反射率的损失。这种背向反射可以被分为两部分:
(1)对于某些角度所有的波长,接收器接收不到,如图4a中所示,这部分损失和波长无关。
(2)对于某些角度,有一部分波长的背向反射接收不到,而一部分可以接收得到,如图4b所示,这部分损失和波长有关。
通过改变式(6)中的R3(如图2中所示)的大小来补偿背向反射的损失。根据上面的分析,两种不同的损失需要两种不同的补偿,即波长无关的补偿和波长相关的补偿。波长无关的补偿比较简单,可以写成:
f1=C1
其中C1是一个常数。
对于波长相关的部分,假设最小可以接收到的角度为θmin,最大的可以接收到的角度为θreceive(λ),Rreceive(λ)代表不同波长的反射率,λ1,λ2代表两个不同的波长,n(λ)代表薄膜第一层材料的折射率,这两个波长所接收到的能量可以表示成:
R receive ( λ 1 ) R receive ( λ 2 ) = θ receive ( λ 1 ) - θ min θ receceive ( λ 2 ) - θ min ≈ θ receive ( λ 1 ) θ receive ( λ 2 ) = 1 / n ( λ 1 ) 1 / n ( λ 2 ) - - - ( 7 )
那么不同波长损失的部分的能量比可以表示为,其中Loss(λ)表示不同波长损失的能量,C为一个常数:
Loss ( λ 1 ) Loss ( λ 2 ) = C - 1 / n ( λ 1 ) C - 1 / n ( λ 2 ) - - - ( 8 )
那么波长相关的部分的损失可以表示成:
f 2 ( λ ) = C ′ × ( C - 1 n ( λ ) ) = C 3 - C 2 n ( λ ) - - - ( 9 )
其中C2和C3是常数,C1′=C1+C3,ng(λ)为基底折射率。这样最终的经过补偿修正的R3可以表示为:
R 3 ( λ ) = ( 1 - n g ( λ ) 1 + n g ( λ ) ) 2 + C 1 - C 2 n ( λ ) + C 3 = ( 1 - n g ( λ ) 1 + n g ( λ ) ) 2 + C 1 ′ - C 2 n ( λ ) - - - ( 10 )
将式(10)带入式(6),R3×R′f非常小可以忽略,对于无吸收的薄膜Tf=1-Rtf,最终的经过补偿的薄膜总反射率R′可表示为:
R ′ = R f + ( 1 - R f ) 2 [ ( 1 - n g ( λ ) 1 + n g ( λ ) ) 2 + C 1 ′ - C 2 n ( λ ) ] - - - ( 11 )
其中R′可以根据式1-5计算,在计算过程中式(2)的薄膜各层厚度dj为未知数。因为C1′和C2是推导的参数,无法直接计算,可以将这两个参数同薄膜各层厚度dj,利用全局优化算法在一定范围内进行搜索寻找最优值。每搜索一次就利用dj以及C1′、C2计算一次总反射率R,比较总反射率R以及测量的反射率。最终使最优的各层薄膜厚度dj以及最优的C1′、C2计算出来的总反射率R和测量的反射率最接近。此时的认为薄膜各层厚度dj为实际薄膜厚度,此时的C1′、C2为实际的补偿系数。计算过程如下:
1.利用光谱仪测量透明基底薄膜的反射率。
2确定待测的各层薄膜的厚度dj以及C1′和C2的上下限。
3.利用全局优化算法搜索上下限内的最优的各层薄膜厚度值dj以及C1′、C2,在计算总反射率时利用经过背向反射率补偿的式(11)。
本发明系统无须在样片背面涂覆消光物质来消除背向反射,无需破坏薄膜基底。本发明系统具有简单灵活,自动适应各种不同材料不同厚度的透明基底,整个系统具有结构简单,成本较低,小型化,无损探测等优点。
附图说明
图1是基于背向补偿的透明基底薄膜厚度测量系统的结构示意图。
图2是薄膜反射率叠加示意图。
图3是薄膜基底背向反射损失示意图。
图4不同类型背向反射损失示意图。
图5为本发明实施例中样品测试的波长-反射率数据曲线图。
图6为本发明实施例未经过补偿计算拟合的波长-反射率数据曲线图。
图7为本发明实施例经过补偿计算后拟合的波长-反射率数据曲线图。
具体实施方式
下面结合附图详细说明本发明的具体实施方式。
如图1所示,基于背向补偿的透明基底薄膜厚度测量系统,包括光源1、光谱仪2、PC机3、测量平台4和反射率光纤探头5。反射率光纤探头5设于测量平台4上方,分别与光源1、光谱仪2连接,光谱仪2与PC机3连接;光源1发出的光经过反射率光纤探头5的出射端面,照射到放置在测量平台4上的待测具有透明基底7的薄膜6上,透明基底7以及薄膜6的反射光同时又进入反射率光纤探头5的接收端面,被光谱仪2接收,光谱仪2接受的数据经过PC机3的补偿处理计算,得出待测薄膜的厚度。
测量实例采用卤钨灯光源,根据光谱仪其波长在450nm到800nm。所述光谱仪采用Scopetek_STDFSM3666,探测器像素数目为3666,光谱范围为380nm-980nm,光谱分辨率在中心波长为0.7nm。所述的反射率光纤探头为海洋光学QP230-0.25-XSR标准反射率光纤。样品为BK7玻璃基底的MgF2薄膜,其厚度为192nm±2nm。
对样品进行测试,具体过程为:
1.测量样品的反射率,测量结果如图5所示。
2.该样品为单层薄膜其厚度范围根据设计指标约为0-500nm,一般的薄膜C1′,C2范围为0-1。
3.根据一般方法,不对背向反射进行补偿即利用式(1)-(5)计算反射率,得到结果如图6所示,测量厚度为10nm,和实际薄膜厚度差距较大。
4.利用式(11)以及式(1)-(5)计算反射率,得到的结果如图7所示,测量厚度为192.3nm和实际厚度一致,C1′为0.0078,C2为0.0006。图7中的测量值已经减去背向反射,所以比图6显示的测量值低。

Claims (3)

1.一种基于背向补偿的透明基底薄膜厚度测量系统,包括:光源、光谱仪、反射率光纤探头、测量平台和PC机,所述的反射率光纤探头设于测量平台上方,分别与光源、光谱仪连接,光谱仪与PC机连接;光源发出的光经过反射率光纤探头的出射端面,照射到放置在测量平台上的待测具有透明基底的薄膜上,基底以及薄膜的反射光同时又进入反射率光纤探头的接收端面,被光谱仪接收,接受的数据经过PC机的补偿处理计算,得到薄膜各层的实际厚度,其特征在于:所述的补偿处理计算方法为:
(1)首先测量待测样品反射率;
(2)根据薄膜理论计算薄膜的总反射率R:
R = R f + T f ′ R 3 T f + T f ′ R 3 2 R ′ f T f . . . . . . . . . = R f + T f 2 R 3 1 - R 3 R ′ f ≈ R f + ( 1 - R f ) 2 R 3
其中Rf表示从正面到背面的相干的反射率,Rf由各层厚度和折射率计算得到;R3表示基底背面的背向反射率,由基底材料折射率计算得到;R′f表示背面到正面的相干反射率;Tf表示从正面到背面的相干透过率;T′f表示从背面到正面的相干透过率;
(3)对基底背面的背向反射率R3进行补偿修正,补偿修正的R3(λ)表示为:
R 3 ( λ ) = ( 1 - n g ( λ ) 1 + n g ( λ ) ) 2 + C 1 - C 2 n ( λ ) + C 3 = ( 1 - n g ( λ ) 1 + n g ( λ ) ) 2 + C 1 ′ - C 2 n ( λ )
其中:C1、C2、C3是常数,ng(λ)为基底折射率,n(λ)是薄膜的折射率;C1′=C1+C3,补偿计算后的总反射率R′表示为:
R ′ = R f + ( 1 - R f ) 2 [ ( 1 - n g ( λ ) 1 + n g ( λ ) ) 2 + C 1 ′ - C 2 n ( λ ) ]
(4)在一定范围内改变薄膜各层的厚度以及C1′、C2,每改变一次计算对应的总反射率R′,与测量的反射率比较;利用全局优化算法在设定范围内搜索出最优的各层薄膜厚度以及C1′、C2,使补偿计算后的总反射率R′与实际反射率最接近;
(5)当经过补偿计算的总反射率R′与测量的反射率误差小于阈值,则此时的薄膜各层厚度为实际厚度,此时的C1′、C2为实际的补偿值。
2.根据权利要求1所述的基于背向补偿的透明基底薄膜厚度测量系统,其特征在于:所述的光源采用卤钨灯光源,波长在450nm-800nm。
3.根据权利要求1所述的基于背向补偿的透明基底薄膜厚度测量系统,其特征在于:所述的反射率光纤探头为七芯反射率光纤,其中六芯用于出光,一芯用于接收。
CN2011100922833A 2011-04-13 2011-04-13 一种基于背向补偿的透明基底薄膜厚度测量系统 Expired - Fee Related CN102243065B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100922833A CN102243065B (zh) 2011-04-13 2011-04-13 一种基于背向补偿的透明基底薄膜厚度测量系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100922833A CN102243065B (zh) 2011-04-13 2011-04-13 一种基于背向补偿的透明基底薄膜厚度测量系统

Publications (2)

Publication Number Publication Date
CN102243065A CN102243065A (zh) 2011-11-16
CN102243065B true CN102243065B (zh) 2012-11-14

Family

ID=44961239

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100922833A Expired - Fee Related CN102243065B (zh) 2011-04-13 2011-04-13 一种基于背向补偿的透明基底薄膜厚度测量系统

Country Status (1)

Country Link
CN (1) CN102243065B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107121080A (zh) * 2017-06-16 2017-09-01 东南大学 一种测量有序多孔纳米薄膜厚度的方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435146A (zh) * 2011-12-08 2012-05-02 广州标旗电子科技有限公司 光学透镜中心厚度测量系统及方法
CN103837088B (zh) * 2014-03-25 2016-08-24 北京印刷学院 光纤传感润版液水膜厚度连续测量系统和方法
CN105091762B (zh) * 2015-05-26 2016-05-25 雷艳梅 一种光刀透射式试剂药液涂抹均匀性检测方法及装置
CN105088172B (zh) * 2015-09-21 2017-12-01 京东方科技集团股份有限公司 膜厚控制系统及膜厚控制方法
CN106247958A (zh) * 2016-07-11 2016-12-21 京东方科技集团股份有限公司 测量设备及测量方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6392756B1 (en) * 1999-06-18 2002-05-21 N&K Technology, Inc. Method and apparatus for optically determining physical parameters of thin films deposited on a complex substrate
CN100573036C (zh) * 2008-04-20 2009-12-23 华中科技大学 一种薄膜厚度和折射率的光学测量方法
CN101865641B (zh) * 2010-03-05 2012-05-30 新奥光伏能源有限公司 一种测量半导体薄膜厚度的方法及装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107121080A (zh) * 2017-06-16 2017-09-01 东南大学 一种测量有序多孔纳米薄膜厚度的方法

Also Published As

Publication number Publication date
CN102243065A (zh) 2011-11-16

Similar Documents

Publication Publication Date Title
CN102243065B (zh) 一种基于背向补偿的透明基底薄膜厚度测量系统
US9176048B2 (en) Normal incidence broadband spectroscopic polarimeter and optical measurement system
KR20160114080A (ko) 막두께 계측 방법 및 막두께 계측 장치
Jäger et al. Angular resolved scattering measurements of nano-textured substrates in a broad wavelength range
CN106885787A (zh) 一种测量玻璃表面粗糙薄膜复折射率的方法和装置
US5526117A (en) Method for the determination of characteristic values of transparent layers with the aid of ellipsometry
KR20140107332A (ko) 막두께 분포 측정 방법
CN105066889A (zh) 一种便携式薄膜测厚仪及其膜厚测量方法
CN107917672B (zh) 一种用于提高超薄金属薄膜测试灵敏度的测试方法
CN107991728A (zh) 宽带及宽视场角补偿器
CN114216559A (zh) 一种星上定标机构的部分孔径因子测量方法和装置
Bai et al. A new method to measure spectral reflectance and film thickness using a modified chromatic confocal sensor
Bae et al. Optical method for simultaneous thickness measurements of two layers with a significant thickness difference
CN107533195A (zh) 用于光学纤维连接的方法和设备
CN103575661A (zh) 包含垂直入射和斜入射的光学测量系统
US8976360B2 (en) Surface plasmon sensor and method of measuring refractive index
Wang et al. Polarized angle-resolved spectral reflectometry for real-time ultra-thin film measurement
Ayupov et al. Searching for the starting approximation when solving inverse problems in ellipsometry and spectrophotometry
US20070174014A1 (en) Method for matching a model spectrum to a measured spectrum
CN109883553B (zh) 一种偏振测量装置及偏振测量方法
US11662197B2 (en) Rapid measurement method for ultra-thin film optical constant
CN109342317B (zh) 一种用于平面黑材料的光学参数测试装置及方法
CN109141260B (zh) 多孔膜厚度与孔隙率二维分布测定方法
JP4136740B2 (ja) 分光エリプソメータを用いた薄膜3層構造の解析方法
CN107525589B (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
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121114

Termination date: 20160413

CF01 Termination of patent right due to non-payment of annual fee