CN103983367B - 基于光强分析的分数阶涡旋光束拓扑荷值测量方法 - Google Patents
基于光强分析的分数阶涡旋光束拓扑荷值测量方法 Download PDFInfo
- Publication number
- CN103983367B CN103983367B CN201410181175.7A CN201410181175A CN103983367B CN 103983367 B CN103983367 B CN 103983367B CN 201410181175 A CN201410181175 A CN 201410181175A CN 103983367 B CN103983367 B CN 103983367B
- Authority
- CN
- China
- Prior art keywords
- topological charge
- measured
- charge values
- light beam
- speck
- 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
Links
Landscapes
- Instruments For Measurement Of Length By Optical Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
基于光强分析的分数阶涡旋光束拓扑荷值测量方法,利用光束分束装置将拓扑荷值为m的待测涡旋光束分为两束光,将其中一束光倒置,变为待测光束的共轭光束;待测涡旋光束与其共轭涡旋光束相互干涉,形成的干涉图案在CCD相机中成像;然后,储存进计算机;利用计算机对干涉条纹图进行处理,将待测涡旋光束与其共轭光束相干涉,并记录干涉条纹图,后根据干涉条纹图亮条纹分布以及各个亮斑的强度比较计算得到拓扑荷值,本发明方法能实现涡旋光束任意阶(0.1阶)精度拓扑荷值的测量,具有实质性特点和显著进步,可广泛应用于玻色‑爱因斯坦凝聚、量子通信、信息编码与传输、粒子囚禁、光镊、光扳手等领域的拓扑荷值测量。
Description
技术领域
本发明涉及一种测量分数阶涡旋光束拓扑荷值的方法,具体的说是涉及一种利用光强分析测量分数阶涡旋光束拓扑荷值测量方法。
背景技术
由于涡旋光束在玻色-爱因斯坦凝聚、量子信息编码、粒子囚禁、光镊及光扳手等领域具有重要的应用前景,成为近年来信息光学领域一个非常重要的研究热点。由于涡旋光束的拓扑荷值(即光子轨道角动量)是携带了光束的信息和能量,因此,精确测量拓扑荷值是涡旋光束研究时首先要解决的问题。
从目前研究看,涡旋光束拓扑荷值的测量主要分为干涉测量和衍射测量。其中,P.VaⅠty利用倾斜双凸透镜的迈克尔逊干涉光路测量了整数阶涡旋[Optics Letters,37(2012)1301-1303],通过数干涉亮条纹数可测量拓扑荷值为14以内的涡旋光束;而典型的衍射测量方法有三角孔衍射法,该方法可测量7以内的拓扑荷值[Optics Letters,36(2011)787-789]。但这两种方法都是通过数干涉/衍射条纹数测量来实现,仅能达到半整数阶(0.5阶)精度[A.Mourka et al.,Optics Express 19(2011)5760]的拓扑荷值测量。
因此,如何实现任意阶(0.1阶)精度的拓扑荷值的测量是该技术领域面临的一个亟待解决的技术难题。
发明内容
本发明要解决的技术问题:提供一种能实现任意阶(0.1阶)精度的基于光强分析的分数阶涡旋光束拓扑荷值测量方法。
本发明所采用的技术方案为:基于光强分析的分数阶涡旋光束拓扑荷值测量方法,该方法利用由准直扩束器、分束镜Ⅰ、分束镜Ⅱ、道威棱镜、反射镜Ⅰ、反射镜Ⅱ、成像透镜、CCD相机和计算机构成的装置进行检测,待测涡旋光束经准直扩束器扩束后照射在分束镜Ⅰ上,分为透射光束和反射光束;透射光束照射在一道威棱镜上,经过道威棱镜后,照射在反射镜Ⅰ上,反射后照射在分束镜Ⅱ上;反射光束照射在反射镜Ⅱ上,反射后照射在分束镜Ⅱ上;透射光束和反射光束经分束镜Ⅱ合束后,经成像透镜后在CCD相机中干涉成像,CCD相机记录的干涉条纹图传输进计算机进行处理。
基于光强分析的分数阶涡旋光束拓扑荷值测量装置的测量方法,主要为:
将待测涡旋光束与其共轭光束相干涉,并记录干涉条纹图;
若干涉条纹图亮条纹分布具有圆对称性,则利用公式m=n0/2,求得涡旋光束的拓扑荷值,其中,m为拓扑荷值,n0为亮斑个数;
若干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各个亮斑的强度相同,则通过公式m=n1/2+ε1,求得涡旋光束的拓扑荷值,其中,n1为亮斑个数,ε1的值通过干涉图样对称轴上两个光强极大值的比值α1来确定;
若干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各亮斑的强度不同时,则通过公式m=n2/2+ε2,求得涡旋光束的拓扑荷值,其中,n2为强度相同的亮斑个数,ε2的值通过干涉图样对称轴上两个光强极大值的比值α2来确定。
基于光强分析的分数阶涡旋光束拓扑荷值测量装置的测量方法,具体的步骤如下:
步骤一、利用光束分束装置将拓扑荷值为m的待测涡旋光束分为两束光,将其中一束光倒置,变为待测光束的共轭光束;
步骤二、待测涡旋光束与其共轭涡旋光束相互干涉,形成的干涉图案在CCD相机中成像;然后,储存进计算机;
步骤三、利用计算机对干涉条纹图进行处理,若干涉条纹图亮条纹分布具有圆对称性,则利用公式m=n0/2,求得涡旋光束的拓扑荷值,其中,m为拓扑荷值,n0为亮斑个数;
步骤四、若干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各个亮斑的强度相同,则通过公式m=n1/2+ε1,求得涡旋光束的拓扑荷值,其中,n1为亮斑个数,ε1的值通过干涉图样对称轴上两个光强极大值的比值α1来确定;
步骤五、在步骤四条件下,若α1=0.9,则ε1=0.1,此时待测光束的拓扑荷值为m=n1/2+0.1;
步骤六、在步骤四条件下,若α1=0.655,则ε1=-0.3,此时待测光束的拓扑荷值为m=n1/2-0.3;
步骤七、在步骤四条件下,若α1=0.346,则ε1=-0.2,此时待测光束的拓扑荷值为m=n1/2-0.2;
步骤八、在步骤四条件下,若α1=0.1,则ε1=-0.1,此时待测光束的拓扑荷值为m=n1/2-0.1;
步骤九、若干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各亮斑的强度不同时,则通过公式m=n2/2+ε2,求得涡旋光束的拓扑荷值,其中,n2为强度相同的亮斑个数,ε2的值通过干涉图样对称轴上两个光强极大值的比值α2来确定;
步骤十、在步骤九条件下,若α2=0.9,则ε2=0.4,此时待测光束的拓扑荷值为m=n2/2+0.4;
步骤十一、若α2=0.655,则ε2=0.3,此时待测光束的拓扑荷值为m=n2/2+0.3;
步骤十二、若α2=0.346,则ε2=0.2,此时待测光束的拓扑荷值为m=n2/2+0.2;
步骤十三、若α2=0.1,则ε2=0.1,此时待测光束的拓扑荷值为m=n2/2+0.1;
步骤十四、最终,利用干涉条纹图的光强分析可实现任意阶精度涡旋光束拓扑荷值的测量。
本发明的工作原理是:
假设有一束待测拉盖尔-高斯涡旋光束[Phys.Rev.Lett.92,143905(2004)],其在观察面的复振幅为:
其中,w0为光束的束腰,m为涡旋光束的拓扑荷值。其共轭光束为,
两束光在成像平面叠加的复振幅为,
E=E1+E2(3)
由(1)、(2)、(3)式得两束光在成像平面相干叠加的光强为
分析(4)式可知,该光强图的干涉亮条纹分布在一个圆上。通过分析亮条纹的特点,我们就可以得到涡旋光束的拓扑荷值。最终,该方法可实现任意阶(0.1阶)精度的拓扑荷值的测量。
与以往技术相比,本发明方法能实现涡旋光束任意阶(0.1阶)精度拓扑荷值的测量,具有实质性特点和显著进步,可广泛应用于玻色-爱因斯坦凝聚、量子通信、信息编码与传输、粒子囚禁、光镊、光扳手等领域的拓扑荷值测量。
附图说明
图1为发明中测量装置的结构示意图。
图2为利用图1测量装置测得的拓扑荷值分别为m=2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3.0的10束待测涡旋光束的干涉条纹图。其中,白色虚线为干涉条纹图的对称轴,白色实线为对称轴上的归一化光强值曲线。
附图标记:100、待测涡旋光束,110、准直扩束器,121、分束镜Ⅰ,122、分束镜Ⅱ,130、道威棱镜,141、反射镜Ⅰ,142、反射镜Ⅱ,150、成像透镜,200、CCD相机,300、计算机。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步详细的说明。
附图1为测量装置的结构示意图,该装置采用了改进的马赫-曾德尔干涉光路,待测涡旋光束100经准直扩束器110扩束后照射在分束镜Ⅰ121上,分为透射光束和反射光束;透射光束照射在一道威棱镜130上,经过道威棱镜130后,照射在反射镜Ⅰ141上,反射后照射在分束镜Ⅱ122上;反射光束照射在反射镜Ⅱ142上,反射后照射也在分束镜Ⅱ122上;透射光束和反射光束经分束镜Ⅱ122合束后,经成像透镜150后在CCD相机200中干涉成像;CCD相机200记录的干涉条纹图传输进计算机300进行处理;
所述的透射光束经道威棱镜130后,变成待测涡旋光束的共轭光束;
所述的反射光束与待测涡旋光束的特性相同(即拓扑荷值不变)。
在该实施例中,利用本发明的方法,我们分别对m=2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3.0的10束待测涡旋光束的拓扑荷值进行了测量,测量结果见附图2。下面结合附图2,对本实施例的测量过程进一步说明。
一种基于光强分析的分数阶涡旋光束拓扑荷值测量方法,其主要是:
其一、将待测涡旋光束与其共轭光束相干涉,记录干涉条纹图;
其二、若干涉条纹图亮条纹分布具有圆对称性,则利用公式m=n0/2,求得涡旋光束的拓扑荷值,其中,m为拓扑荷值,n0为亮斑个数;
其三、若干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各个亮斑的强度相同,则通过公式m=n1/2+ε1,求得涡旋光束的拓扑荷值,其中,n1为亮斑个数,ε1的值通过干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α1来确定;
其四、若干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各亮斑的强度不同时,则通过公式m=n2/2+ε2,求得涡旋光束的拓扑荷值,其中,n2为强度相同的亮斑个数,ε2的值通过干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α2来确定。
所述的基于光强分析的分数阶涡旋光束拓扑荷值测量方法,其具体步骤如下:
步骤一、利用图1中的分束镜Ⅰ121将拓扑荷值为m的待测涡旋光束分为两束光,利用道威棱镜130将其中的透射光束倒置,变为待测光束的共轭光束(其拓扑荷值变为-m);
步骤二、待测涡旋光束与其共轭涡旋光束相互干涉,形成的干涉图案在CCD相机200中成像;然后,储存进计算机300;
步骤三、利用计算机300对干涉条纹图进行处理,图2(e)中干涉条纹图亮条纹分布具有圆对称性,且亮斑个数n0=5,求得涡旋光束的拓扑荷值为m=n0/2=2.5;图2(j)中干涉条纹图亮条纹分布具有圆对称性,且亮斑个数n0=6,求得涡旋光束的拓扑荷值为m=n0/2=3.0;
步骤四、图2(a)中干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各个亮斑的强度相同,亮斑个数n1=4,计算可得干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α1=0.9045(相对误差<10%),则ε1=0.1;则涡旋光束的拓扑荷值为m=n1/2+ε1=2.1;
步骤五、图2(b)中干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各个亮斑的强度相同,亮斑个数n1=5,计算可得干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α1=0.6545,所以ε1=-0.3;则涡旋光束的拓扑荷值为m=n1/2+ε1=2.2;
步骤六、图2(c)中干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各个亮斑的强度相同,亮斑个数n1=5,计算可得干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α1=0.3455,所以ε1=-0.2;则涡旋光束的拓扑荷值为m=n1/2+ε1=2.3;
步骤七、图2(d)中干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各个亮斑的强度相同,亮斑个数n1=5,计算可得干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α1=0.3455,所以ε1=-0.1;则涡旋光束的拓扑荷值为m=n1/2+ε1=2.4;
步骤八、图2(f)中干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各亮斑的强度不同,亮斑个数n2=5,计算可得干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α2=0.0955,所以ε2=0.1;则涡旋光束的拓扑荷值为m=n2/2+ε2=2.6;
步骤九、图2(g)中干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各亮斑的强度不同,亮斑个数n2=5,计算可得干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α2=0.3455,所以ε2=0.2;则涡旋光束的拓扑荷值为m=n2/2+ε2=2.7;
步骤十、图2(h)中干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各亮斑的强度不同,亮斑个数n2=5,计算可得干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α2=0.6545,所以ε2=0.3;则涡旋光束的拓扑荷值为m=n2/2+ε2=2.8;
步骤十一、图2(i)中干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各亮斑的强度不同,亮斑个数n2=5,计算可得干涉图样对称轴上两个光强极大值的比值(两者中的小值除以大值)α2=0.9045,所以ε2=0.4;则涡旋光束的拓扑荷值为m=n2/2+ε2=2.9;
步骤十二、与这10束涡旋光束生成时的拓扑荷值比较,该方法利用干涉条纹图的光强分析准确实现了这10束分数阶涡旋光束拓扑荷值的测量。
经实验表明:本发明方法能实现任意阶(0.1阶)涡旋光束拓扑荷值的测量,与现有测量方法相比,测试精度提高了一个数量级;并且具有方法简单、快速、准确的特点。
Claims (1)
1.基于光强分析的分数阶涡旋光束拓扑荷值测量方法,该方法利用由准直扩束器(110)、分束镜Ⅰ(121)、分束镜Ⅱ(122)、道威棱镜(130)、反射镜Ⅰ(141)、反射镜Ⅱ(142)、成像透镜(150)、CCD相机(200)和计算机(300)构成的装置进行检测,待测涡旋光束(100)经准直扩束器(110)扩束后照射在分束镜Ⅰ(121)上,分为透射光束和反射光束;透射光束照射在一道威棱镜(130)上,经过道威棱镜(130)后,照射在反射镜Ⅰ(141)上,反射后照射在分束镜Ⅱ(122)上; 反射光束照射在反射镜Ⅱ(142)上,反射后照射在分束镜Ⅱ(122)上;透射光束和反射光束经分束镜Ⅱ(122)合束后,经成像透镜(150)后在CCD相机(200)中干涉成像,CCD相机(200)记录的干涉条纹图传输进计算机(300)进行处理;其特征在于:测量步骤如下:
步骤一、利用光束分束装置将拓扑荷值为m的待测涡旋光束分为两束光,将其中一束光倒置,变为待测光束的共轭光束;
步骤二、待测涡旋光束与其共轭涡旋光束相互干涉,形成的干涉图案在CCD相机中成像;然后,储存进计算机;
步骤三、利用计算机对干涉条纹图进行处理,若干涉条纹图亮条纹分布具有圆对称性,则利用公式,求得涡旋光束的拓扑荷值,其中,m为拓扑荷值,n 0为亮斑个数;
步骤四、若干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各个亮斑的强度相同,则通过公式,求得涡旋光束的拓扑荷值,其中,n 1为亮斑个数,的值通过干涉图样对称轴上两个光强极大值的比值来确定;
步骤五、在步骤四条件下,若,则,此时待测光束的拓扑荷值为m=n 1/2+0.1;
步骤六、在步骤四条件下,若,则,此时待测光束的拓扑荷值为;
步骤七、在步骤四条件下,若,则,此时待测光束的拓扑荷值为;
步骤八、在步骤四条件下,若,则,此时待测光束的拓扑荷值为;
步骤九、若干涉条纹亮斑分布不具有圆对称性但具有轴对称性且各亮斑的强度不同时,则通过公式,求得涡旋光束的拓扑荷值,其中,n 2为强度相同的亮斑个数,的值通过干涉图样对称轴上两个光强极大值的比值来确定;
步骤十、在步骤九条件下,若,则,此时待测光束的拓扑荷值为;
步骤十一、若,则,此时待测光束的拓扑荷值为;
步骤十二、若,则,此时待测光束的拓扑荷值为;
步骤十三、若,则,此时待测光束的拓扑荷值为;
步骤十四、最终,利用干涉条纹图的光强分析可实现任意阶精度涡旋光束拓扑荷值的测量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410181175.7A CN103983367B (zh) | 2014-04-30 | 2014-04-30 | 基于光强分析的分数阶涡旋光束拓扑荷值测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410181175.7A CN103983367B (zh) | 2014-04-30 | 2014-04-30 | 基于光强分析的分数阶涡旋光束拓扑荷值测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103983367A CN103983367A (zh) | 2014-08-13 |
CN103983367B true CN103983367B (zh) | 2017-01-25 |
Family
ID=51275427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410181175.7A Expired - Fee Related CN103983367B (zh) | 2014-04-30 | 2014-04-30 | 基于光强分析的分数阶涡旋光束拓扑荷值测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103983367B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104833314B (zh) * | 2015-06-03 | 2018-07-03 | 中国科学技术大学 | 一种高分辨率光纤应变传感器及测量方法 |
CN105353463B (zh) * | 2015-12-04 | 2018-07-03 | 东南大学 | 一种检测和接收涡旋光场的装置及方法 |
CN105445943B (zh) * | 2015-12-24 | 2018-04-13 | 河南科技大学 | 一种分数阶完美涡旋光束的产生装置及产生方法 |
CN106767428B (zh) * | 2016-11-24 | 2019-01-22 | 李达成 | 基于全息共轭光补偿空气扰动的激光准直、位移测量系统 |
CN110402375A (zh) * | 2018-11-05 | 2019-11-01 | 东莞理工学院 | 基于模式变换的oam模式解调系统 |
CN109297930B (zh) * | 2018-11-14 | 2020-11-06 | 深圳大学 | 一种基于涡旋光束共轭干涉的三阶非线性测量装置及方法 |
CN111175969B (zh) * | 2020-01-03 | 2020-12-11 | 浙江大学 | 一种基于涡旋对光束的光镊系统 |
CN114964527B (zh) * | 2022-05-10 | 2023-06-02 | 苏州大学 | 一种部分相干分数阶涡旋光束拓扑荷数测量方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706444A (zh) * | 2012-06-25 | 2012-10-03 | 杭州电子科技大学 | 一种光旋涡拓扑数测量方法 |
-
2014
- 2014-04-30 CN CN201410181175.7A patent/CN103983367B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706444A (zh) * | 2012-06-25 | 2012-10-03 | 杭州电子科技大学 | 一种光旋涡拓扑数测量方法 |
Non-Patent Citations (3)
Title |
---|
Spatial Correlation Singularity of a Vortex Field;D.M.Palacios et al.;《PHYSICAL REVIEW LETTERS》;20040409;第92卷(第14期);第143905-1到143905-4页 * |
Visualization of the birth of an optical vortex using diffraction from a triangular aperture;A Mourka et al.;《Optics Express》;20110314;第19卷(第7期);第5760-5771页 * |
太阳能硅电池表面特性的激光散斑研究;李新忠 等;《河南科技大学学报》;20120831;第33卷(第4期);第28-31页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103983367A (zh) | 2014-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103983367B (zh) | 基于光强分析的分数阶涡旋光束拓扑荷值测量方法 | |
CN103954367B (zh) | 一种测量分数阶光学涡旋拓扑荷值的装置及其测量方法 | |
Bouchet et al. | Maximum information states for coherent scattering measurements | |
CN105466577B (zh) | 基于光强分析的完美涡旋光束拓扑荷值的测量装置及方法 | |
CN102818786B (zh) | 正弦相位调制的并行复频域光学相干层析成像系统和方法 | |
Lue et al. | Single-shot quantitative dispersion phase microscopy | |
CN102735184B (zh) | 一种光学面形的检测装置及检测方法 | |
Claus et al. | Dual wavelength optical metrology using ptychography | |
CN103424190B (zh) | 双楔板色散剪切干涉超光谱成像装置及方法 | |
CN102230827B (zh) | 马赫-曾德点衍射干涉仪及激光复振幅重建方法 | |
CN104748672A (zh) | 干涉量分离单频激光干涉仪非线性误差修正方法及装置 | |
CN208060025U (zh) | 实现测量任意光束的偏振态和相位的两幅干涉图的光路 | |
CN107339944B (zh) | 一种稀土钕铁硼表面镀膜厚度在线测量系统及测量方法 | |
CN106949842A (zh) | 二维位移测量装置及测量方法 | |
CN104748671A (zh) | 角位移单频激光干涉仪非线性误差修正方法及装置 | |
Villiger et al. | Image formation and tomogram reconstruction in optical coherence microscopy | |
CN105181141B (zh) | 一种外差式偏振干涉光谱成像方法及光谱仪 | |
Wang et al. | Quantitative interferometric microscopy with two dimensional Hilbert transform based phase retrieval method | |
Zhang et al. | Ultrahigh-accuracy measurement of refractive index curves of optical materials using interferometry technology | |
CN102539381B (zh) | 基于微离轴显微干涉投影的折射率层析装置 | |
Marthy et al. | Cross-grating phase microscopy (CGM): In silico experiment (insilex) algorithm, noise and accuracy | |
Sokolenko et al. | Phase shifting profilometry with optical vortices | |
Yang et al. | Realizing depth measurement and edge detection based on a single metasurface | |
CN203837820U (zh) | 基于光强分析的分数阶涡旋光束拓扑荷值测量装置 | |
CN109540017A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170125 Termination date: 20170430 |