CN109269636B - 一种像散椭圆光束的产生及其轨道角动量的测量方法 - Google Patents

一种像散椭圆光束的产生及其轨道角动量的测量方法 Download PDF

Info

Publication number
CN109269636B
CN109269636B CN201810943246.0A CN201810943246A CN109269636B CN 109269636 B CN109269636 B CN 109269636B CN 201810943246 A CN201810943246 A CN 201810943246A CN 109269636 B CN109269636 B CN 109269636B
Authority
CN
China
Prior art keywords
elliptical
angular momentum
astigmatic
spatial light
light modulator
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
CN201810943246.0A
Other languages
English (en)
Other versions
CN109269636A (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.)
South China Normal University
Original Assignee
South China Normal University
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 South China Normal University filed Critical South China Normal University
Priority to CN201810943246.0A priority Critical patent/CN109269636B/zh
Publication of CN109269636A publication Critical patent/CN109269636A/zh
Application granted granted Critical
Publication of CN109269636B publication Critical patent/CN109269636B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00Photometry, e.g. photographic exposure meter
    • G01J1/42Photometry, e.g. photographic exposure meter using electric radiation detectors
    • G01J1/4257Photometry, e.g. photographic exposure meter using electric radiation detectors applied to monitoring the characteristics of a beam, e.g. laser beam, headlamp beam
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/09Beam shaping, e.g. changing the cross-sectional area, not otherwise provided for
    • G02B27/0927Systems for changing the beam intensity distribution, e.g. Gaussian to top-hat
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/09Beam shaping, e.g. changing the cross-sectional area, not otherwise provided for
    • G02B27/0938Using specific optical elements
    • G02B27/095Refractive optical elements
    • G02B27/0955Lenses
    • G02B27/0966Cylindrical lenses

Abstract

本发明公开了一种像散椭圆光束的产生及其轨道角动量的测量方法。不同于利用倾斜柱透镜产生像散椭圆光束的传统方法,本发明利用空间光调制器在椭圆光束束腰处加载交叉相位因子exp(iΩxy),从而产生像散椭圆光束。根据空间光调制器不加载和加载交叉相位因子时,透镜后焦面处椭圆光斑长轴和短轴的测量值可直接确定像散椭圆光束的轨道角动量。本发明方法新颖,装置简单,为像散椭圆光束的产生及其轨道角动量的测量提供新方案。

Description

一种像散椭圆光束的产生及其轨道角动量的测量方法
技术领域
本发明涉及的技术领域是光场调控(optical field manipulation)及其应用,具体涉及光束相位控制(control of the phase of optical beams)。
背景技术
一般来说,携带轨道角动量的光束有两类:空心光束(hollow beams)和像散(astigmatic)椭圆光束。空心光束是包括拉盖尔-高斯(Laguerre-Gaussian)光束和贝塞尔(Bessel)光束在内的光斑中心光强为零的环状光束。这些光束的表达式中均含有exp(ilφ)的相位因子(l为整数,φ为方位角),轨道角动量就来源于该因子。空心光束携带的单光子轨道角动量为
Figure GDA0002619137550000011
(
Figure GDA0002619137550000012
为约化普朗克常数),其轨道角动量往往较低。像散椭圆光束是光斑为椭圆且具有像散(astigmatism)的光束,其轨道角动量来源于交叉相位因子exp(iΘxy)(Θ为交叉相位系数,x、y为与椭圆光斑主轴重合的坐标系),其单光子轨道角动量为
Figure GDA0002619137550000016
(c、b为椭圆长短半轴)[Opt.Commun.,1997,144(4):210]。比较两者的轨道角动量可见,要实现对轨道角动量的调控,无论是调控的维度还是调控的范围,像散椭圆光束均比空心光束更大。像散椭圆光束的单光子轨道角动量在理论上甚至可以高达
Figure GDA0002619137550000014
[Opt.Commun.,1997,144(4):210],而最新实验得到的量值为
Figure GDA0002619137550000015
[Opt.Express,2018,26(1):141]。携带轨道角动量的光束可以实现光学扳手(optical spanner)技术,在光场调控领域有重要的应用,不仅可以实现对介观(mesoscale)尺度微粒的操控,还可以用来操控微观尺度的原子。其基本原理是将光束携带的轨道角动量转移给微粒并使之旋转,实现可操控的微粒旋转,即微纳尺度的光学马达。而要将携带更高轨道角动量的像散椭圆光束用于实际的光操控,则需要探索这类光束更简便易控的产生方法及其轨道角动量更精确的测量方法。
1997年,J.Courtial等人首次发现像散椭圆光束具有轨道角动量[Opt.Commun.,1997,144(4):210]。实验上产生这类光束的传统方法是在椭圆光束束腰处放置对称轴与光束对称轴成一定倾斜角的柱透镜,其后表面的光场相位因子为exp[ik(x2sin2α+y2cos2α+xysin2α)/(2f)](k为波数,α为倾斜角,f为透镜焦距)[Opt.Commun.,1997,144(4):210;物理学报,2006,55(11):5755;Opt.Express,2018,26(1):141],而实际对轨道角动量有贡献的仅仅是最后一项,即exp[ikxysin2α/(2f)]。2006年,董一鸣等人提出一种直接测量像散椭圆光束轨道角动量的方法[物理学报,2006,55(11):5755]。其原理是把像散椭圆光束照射到金属靶上,使其在光束轨道角动量作用下发生转动,通过测量靶转动的角度计算光束的轨道角动量值,然而此方法的测量误差很大。随后Kotlyar等人利用π相移干涉法测出这类像散椭圆光束的轨道角动量,实验测量值与理论值有50%的误差[Opt.Express,2018,26(1):141]。总之,以往的像散椭圆光束都是在椭圆光束束腰处放置倾斜的柱透镜而得到,其轨道角动量随α(椭圆光束对称轴与柱透镜对称轴的夹角)的改变而大范围地变动,不容易精细调节,且实验测量得到的轨道角动量值均有较大的误差。
发明内容
有鉴于此,为了解决现有技术中的上述问题,本发明提出一种像散椭圆光束的产生及其轨道角动量的测量方法。不同于利用倾斜柱透镜产生像散椭圆光束的传统方法,本发明利用空间光调制器在椭圆光束束腰处加载“纯净的”交叉相位因子exp(iΩxy),从而获得像散椭圆光束。该新方法的优点是:产生的像散椭圆光束没有传统方法(椭圆光束经过倾斜柱透镜)必然附加的、使光束发生汇聚或发散而对轨道角动量却无贡献的相位因子exp[ik(x2sin2α+y2cos2α)/2f]。通过改变交叉相位系数Ω的大小即可对像散椭圆光束的轨道角动量进行精细调控,弥补了产生像散椭圆光束的传统技术中的缺陷。通过测量不加载和加载交叉相位因子两种情况下透镜后焦面处椭圆光斑的长轴和短轴,即可确定该像散光束所携带的轨道角动量。本发明方法新颖,装置简单,轨道角动量的调控和测量精度高,为像散椭圆光束的产生及其轨道角动量的测量提供一种新方案。
本发明的技术方案:激光器输出的基模高斯光束经过平行放置的柱透镜对(CL1、CL2,4f系统)产生准直的椭圆高斯光束,空间光调制器(SLM)放置于4f系统的后焦面及透镜F的前焦面,从空间光调制器反射的光束经透镜F聚焦于后焦面(输出面,如图1中虚线所示)。光束通过显微物镜放大后在电荷耦合器上成像,进行相关测量。
一方面,本发明提供一种像散椭圆光束的产生方法,包括如下步骤:
激光器输出的高斯光束经过柱透镜对产生椭圆高斯光束,通过空间光调制器在椭圆高斯光束束腰处加载交叉相位因子,调制入射光束的相位,从空间光调制器反射的像散椭圆高斯光束为:
Figure GDA0002619137550000031
其中π为圆周率,wx、wy为空间光调制器处椭圆光斑长短半轴,x、y为与椭圆光斑主轴重合的坐标系,i为虚数单位,Ω为相息图中的交叉相位系数。
进一步地,所述空间光调制器为反射式纯相位液晶空间光调制器。
另一方面,本发明还提供一种像散椭圆光束轨道角动量的测量方法,包括如下步骤:
激光器输出的高斯光束经过柱透镜对产生椭圆高斯光束,通过空间光调制器在椭圆高斯光束束腰处加载交叉相位因子,调制入射光束的相位,从空间光调制器反射的像散椭圆高斯光束为:
Figure GDA0002619137550000041
其中π为圆周率,wx、wy为空间光调制器处椭圆光斑长短半轴,x、y为与椭圆光斑主轴重合的坐标系,i为虚数单位,Ω为相息图中的交叉相位系数;
从空间光调制器反射的像散椭圆光束经透镜聚焦于输出面并通过显微物镜放大后在电荷耦合器上成像,进行相关测量,从空间光调制器反射的光束经过一次傅里叶变换,可得输出面的光束:
Figure GDA0002619137550000042
其中b、c、Θ为输出光束的短半轴、长半轴及交叉相位系数,分别为:
Figure GDA0002619137550000043
Figure GDA0002619137550000044
Figure GDA0002619137550000045
上式中λ为波长,f为透镜焦距,由(3)-(5)式可得输出光束的交叉相位系数为:
Figure GDA0002619137550000046
其中b0,c0分别为Ω=0时由(3)、(4)式确定的输出面处椭圆光斑短半轴和长半轴,由此可得输出面处像散椭圆光束的单光子轨道角动量为:
Figure GDA0002619137550000051
至此,利用空间光调制器加载交叉相位因子产生了像散椭圆光束,并根据空间光调制器不加载(Ω=0)和加载(Ω≠0)交叉相位因子时输出面处椭圆光斑长短半轴的测量值,由公式(6)即可计算出该像散椭圆光束的单光子轨道角动量。
进一步地,所述空间光调制器为反射式纯相位液晶空间光调制器。
进一步地,所述反射式纯相位液晶空间光调制器设置在所述透镜的前焦面,所述输出面位于所述透镜的后焦面,通过所述透镜构成傅里叶变换关系。
与现有技术相比,本发明的有益效果至少包括:
(1)通过改变空间光调制器所加载的相息图中的交叉相位系数,可精细简便地调节产生的像散椭圆光束的轨道角动量;
(2)根据不加载和加载交叉相位因子时,透镜后焦面处椭圆光斑长轴和短轴的测量值,可对产生的像散椭圆光束的轨道角动量进行精细测量;
(3)本发明集像散椭圆光束的产生及其轨道角动量的测量于一体,装置简单,轨道角动量的调控及测量精度高。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明方法采用的光路示意图。
图2为本发明方法采用的相息图(a)及相应的调制后的光束等相位面示意图(b)。
图3为输出面处不同束宽的椭圆光斑。(a)Ω=0时输出面的椭圆光斑,(b)加载Ω=0.5mm-2相息图时输出面的椭圆光斑,(c)加载Ω=1.0mm-2相息图时输出面的椭圆光斑,(d)加载Ω=2.0mm-2相息图时输出面的椭圆光斑。其中光斑椭圆度c/b=2.0。
图1中:
L:Verdi12固体激光器; CL1,2:柱透镜;
LC-SLM:液晶空间光调制器; CGH:相息图;
M:平面反射镜; F:透镜;
O:输出面; MO:显微物镜;
CCD:电荷耦合器。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面将结合附图和具体的实施例对本发明的技术方案进行详细说明。需要指出的是,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
如图1,2所示,本发明提供一种像散椭圆光束的产生方法,包括如下步骤:
激光器(Verdi 12固体激光器,工作波长532nm)输出的基模高斯光束经过平行竖直放置的柱透镜对产生准直的椭圆高斯光束。通过空间光调制器(1920×1080分辨率)在椭圆高斯光束束腰处加载包含横向坐标交叉相位信息的相息图[如图2(a)所示],调制入射光束的相位,使反射后的椭圆高斯光束携带交叉相位因子。从空间光调制器反射的像散椭圆高斯光束为:
Figure GDA0002619137550000071
其中π为圆周率,wx、wy为空间光调制器处椭圆光斑长短半轴,x、y为与椭圆光斑主轴重合的坐标系,i为虚数单位,Ω为相息图中的交叉相位系数。在我们的实验中,束宽和交叉相位系数的典型数据为:wx=2.0mm,wy=1.0mm,Ω1=0,Ω2=0.5mm-2,Ω3=1.0mm-2,Ω4=2.0mm-2
实施例2
如图1-3所示,本发明还提供一种像散椭圆光束轨道角动量的测量方法,包括如下步骤:
(1)激光器(Verdi 12固体激光器,工作波长532nm)输出的基模高斯光束经过平行竖直放置的柱透镜对产生准直的椭圆高斯光束。通过空间光调制器(1920×1080分辨率)在椭圆高斯光束束腰处加载包含横向坐标交叉相位信息的相息图[如图2(a)所示],调制入射光束的相位,使反射后的椭圆高斯光束携带交叉相位因子。从空间光调制器反射的像散椭圆高斯光束为:
Figure GDA0002619137550000072
其中π为圆周率,wx、wy为空间光调制器处椭圆光斑长短半轴,x、y为与椭圆光斑主轴重合的坐标系,i为虚数单位,Ω为相息图中的交叉相位系数。在我们的实验中,束宽和交叉相位系数的典型数据为:wx=2.0mm,wy=1.0mm,Ω1=0,Ω2=0.5mm-2,Ω3=1.0mm-2,Ω4=2.0mm-2
(2)从空间光调制器反射的像散椭圆高斯光束经透镜聚焦于输出面并通过显微物镜放大后在电荷耦合器上成像,进行相关测量。经过一次傅里叶变换(空间光调制器和显微物镜分别置于透镜的前后焦平面),可得输出面的光束:
Figure GDA0002619137550000081
其中b、c、Θ为输出光束的短半轴、长半轴及交叉相位系数,分别为:
Figure GDA0002619137550000082
Figure GDA0002619137550000083
Figure GDA0002619137550000084
由(3)、(4)式可知输出面处像散椭圆光斑短半轴b,长半轴c随加载的交叉相位系数Ω的变化而变化(如图3所示)。由(3)-(5)式可得输出光束的交叉相位系数为:
Figure GDA0002619137550000085
其中b0,c0分别为Ω=0时由(3)、(4)式确定的输出面处椭圆光斑短半轴和长半轴。由此可得输出面处像散椭圆光束的单光子轨道角动量为:
Figure GDA0002619137550000086
实验上,所用激光波长λ=532nm,透镜焦距f=400mm。当Ω1=0时,实验测得b0=18μm,c0=36μm;当Ω2=0.5mm-2时,实验测得b=25μm,c=51μm,单光子轨道角动量为
Figure GDA0002619137550000091
当Ω3=1.0mm-2时,实验测得b=40μm,c=81μm,单光子轨道角动量为
Figure GDA0002619137550000092
当Ω4=2.0mm-2时,实验测得b=72μm,c=142μm,单光子轨道角动量为
Figure GDA0002619137550000093
至此,本发明利用空间光调制器加载交叉相位因子产生了像散椭圆光束,并根据空间光调制器不加载(Ω=0)和加载(Ω≠0)交叉相位因子时输出面处椭圆光斑长轴和短轴的测量值,由公式(6)计算出该像散椭圆光束的单光子轨道角动量。
发明的有益效果至少包括:
(1)通过改变空间光调制器所加载的相息图中的交叉相位系数,可精细简便地调节产生的像散椭圆光束的轨道角动量;
(2)根据不加载和加载交叉相位因子时,透镜后焦面处椭圆光斑长轴和短轴的测量值,可对产生的像散椭圆光束的轨道角动量进行精细测量;
(3)本发明集像散椭圆光束的产生及其轨道角动量的测量于一体,装置简单,轨道角动量的调控及测量精度高。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (3)

1.一种像散椭圆光束轨道角动量的测量方法,其特征在于,包括如下步骤:
激光器输出的基模高斯光束经过平行放置的柱透镜对产生椭圆高斯光束,通过空间光调制器在椭圆高斯光束束腰处加载交叉相位因子,调制入射光束的相位,从空间光调制器反射的像散椭圆高斯光束为:
Figure FDA0002619137540000011
其中π为圆周率,wx、wy为空间光调制器处椭圆光斑长短半轴,x、y为与椭圆光斑主轴重合的坐标系,i为虚数单位,Ω为相息图中的交叉相位系数;
从空间光调制器反射的像散椭圆光束经透镜聚焦于输出面并通过显微物镜放大后在电荷耦合器上成像,进行相关测量,从空间光调制器反射的光束经过一次傅里叶变换,可得输出面的光束:
Figure FDA0002619137540000012
其中b、c、Θ为输出光束的短半轴、长半轴及交叉相位系数,分别为:
Figure FDA0002619137540000013
Figure FDA0002619137540000014
Figure FDA0002619137540000015
上式中λ为波长,f为透镜焦距,由(3)-(5)式可得输出光束的交叉相位系数为:
Figure FDA0002619137540000021
其中b0,c0分别为Ω=0时由(3)、(4)式确定的输出面处椭圆光斑短半轴和长半轴,由此可得输出面处像散椭圆光束的单光子轨道角动量为:
Figure FDA0002619137540000022
其中,
Figure FDA0002619137540000023
为普朗克常数,至此,利用空间光调制器加载交叉相位因子产生了像散椭圆光束,并根据空间光调制器不加载(Ω=0)和加载(Ω≠0)交叉相位因子时输出面处椭圆光斑长轴和短轴的测量值,由公式(6)即可计算出该像散椭圆光束的单光子轨道角动量。
2.根据权利要求1所述的像散椭圆光束轨道角动量的测量方法,其特征在于,所述空间光调制器为反射式纯相位液晶空间光调制器。
3.根据权利要求2所述的像散椭圆光束轨道角动量的测量方法,其特征在于,所述反射式纯相位液晶空间光调制器设置在所述透镜的前焦面,所述输出面位于所述透镜的后焦面,通过所述透镜构成傅里叶变换关系。
CN201810943246.0A 2018-08-17 2018-08-17 一种像散椭圆光束的产生及其轨道角动量的测量方法 Expired - Fee Related CN109269636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810943246.0A CN109269636B (zh) 2018-08-17 2018-08-17 一种像散椭圆光束的产生及其轨道角动量的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810943246.0A CN109269636B (zh) 2018-08-17 2018-08-17 一种像散椭圆光束的产生及其轨道角动量的测量方法

Publications (2)

Publication Number Publication Date
CN109269636A CN109269636A (zh) 2019-01-25
CN109269636B true CN109269636B (zh) 2020-10-23

Family

ID=65154006

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810943246.0A Expired - Fee Related CN109269636B (zh) 2018-08-17 2018-08-17 一种像散椭圆光束的产生及其轨道角动量的测量方法

Country Status (1)

Country Link
CN (1) CN109269636B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111258163B (zh) * 2020-03-19 2021-04-13 无锡视美乐激光显示科技有限公司 光源装置、光路结构设计方法及投影系统
CN116774450A (zh) * 2023-06-20 2023-09-19 首都师范大学 圆艾里光束诱导空气等离子体调控太赫兹波分布的系统及方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106199981A (zh) * 2016-06-28 2016-12-07 浙江师范大学 一种产生携带轨道角动量的无衍射艾里光束的方法及装置
CN106569369B (zh) * 2016-11-02 2019-05-07 西北大学 一种基于交叉相位调制的贝塞尔光束的获得方法及装置

Also Published As

Publication number Publication date
CN109269636A (zh) 2019-01-25

Similar Documents

Publication Publication Date Title
CN109269636B (zh) 一种像散椭圆光束的产生及其轨道角动量的测量方法
Shmyrov et al. Capillary wave method: An alternative approach to wave excitation and to wave profile reconstruction
US20230398639A1 (en) Method for laser machining inside materials
US20220111469A1 (en) Laser machining inside materials
Seiboth et al. Hard X-ray wavefront correction via refractive phase plates made by additive and subtractive fabrication techniques
Gao et al. Picosecond laser-induced periodic surface structures (LIPSS) on crystalline silicon
Dai et al. Define and measure the dimensional accuracy of two-photon laser lithography based on its instrument transfer function
Šarbort et al. Comparison of three focus sensors for optical topography measurement of rough surfaces
Xu et al. Ultrahigh-aspect-ratio beam generation with super-resolution spot
Das et al. Quantum weak measurement of Goos–Hänchen effect of light in total internal reflection using a Gaussian-mode laser beam
Strauß et al. Evaluation and calibration of LCoS SLM for direct laser structuring with tailored intensity distributions
Ługowski et al. Application of laser-trapping technique for measuring the three-dimensional distribution of viscosity
Nagata et al. The determination of rms roughness and correlation length of rough surface by measuring spatial coherence function
CN214503264U (zh) 一种基于双波长激光调制的细颗粒物测量装置
Liu et al. Robust phase-shifting interferometry resistant to multiple disturbances
Shimizu et al. A PD-edge method associated with the laser autocollimation for measurement of a focused laser beam diameter
Lettieri et al. Certification of SRM1960: nominal 10 μm diameter polystyrene spheres (“space beads”)
Hanus et al. Carrier-envelope phase on-chip scanner and control of laser beams
Gerhard et al. Limits for interferometric measurements on rough surfaces in streaming inhomogeneous media
Zheng et al. Defocusing detection based on asymmetric placement of dual-quadrant detector
Poleshchuk et al. Multibeam laser writing of diffractive optical elements
Miyake et al. Thermal Effusivity Calibration Procedure for Thermal Microscope.
WO2022052197A1 (zh) 皮米尺
Fujita et al. Nanoparticle imaging with polarization interferometric nonlinear confocal microscope
Mercier et al. Ion beam milling fabrication of a small off-axis ellipsoidal mirror, diffraction limited to 1 μm resolution at 14 nm

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
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: 20201023

Termination date: 20210817