CN106653559A - 一种具有宽能量聚焦反射器的飞行时间质谱仪 - Google Patents

一种具有宽能量聚焦反射器的飞行时间质谱仪 Download PDF

Info

Publication number
CN106653559A
CN106653559A CN201611028303.XA CN201611028303A CN106653559A CN 106653559 A CN106653559 A CN 106653559A CN 201611028303 A CN201611028303 A CN 201611028303A CN 106653559 A CN106653559 A CN 106653559A
Authority
CN
China
Prior art keywords
reflector
ion
time
electric field
max
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.)
Granted
Application number
CN201611028303.XA
Other languages
English (en)
Other versions
CN106653559B (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.)
Northwest Institute of Nuclear Technology
Original Assignee
Northwest Institute of Nuclear 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 Northwest Institute of Nuclear Technology filed Critical Northwest Institute of Nuclear Technology
Priority to CN201611028303.XA priority Critical patent/CN106653559B/zh
Publication of CN106653559A publication Critical patent/CN106653559A/zh
Application granted granted Critical
Publication of CN106653559B publication Critical patent/CN106653559B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01JELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J49/00Particle spectrometers or separator tubes
    • H01J49/26Mass spectrometers or separator tubes
    • H01J49/34Dynamic spectrometers
    • H01J49/40Time-of-flight spectrometers
    • H01J49/406Time-of-flight spectrometers with multiple reflections

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明提出了一种针对飞行时间质谱的宽能量聚焦反射器设计方法。针对具有较宽能量分布的离子,采用数值计算方法设计了由多级电场组成的反射器的各级电场尺寸与各极板需要加载的电压,实现离子的飞行时间聚焦,从而提高飞行时间质谱的分辨率。该反射器对机械加工和装配精度的要求较低。该反射器不必在每两级电场之间都设计栅网,只需保留前两级栅网即可,这样可以保证反射器具有较高的离子传输效率,此时,可以通过微调电压补偿电场渗透带来的影响。

Description

一种具有宽能量聚焦反射器的飞行时间质谱仪
技术领域
本发明属于质谱仪器技术领域,可应用于反射式飞行时间质谱仪设计,尤其涉及一种飞行时间质谱仪中的反射器,该反射器能够实现宽能量分布的离子的飞行时间聚焦。
背景技术
飞行时间质谱的分辨率R=t/2Δt,其中,t是离子飞行的平均时间,Δt是离子飞行的时间差。离子的飞行时间差会随着离子的能量色散的增大而增大,因此离子的能量色散是影响飞行时间质谱仪分辨的一个主要因素。
在飞行时间质谱仪中采用反射器可以实现能量时间聚焦,从而提高飞行时间质谱仪的分辨率。目前,飞行时间质谱仪中反射器的设计方法主要是通过使飞行时间对离子能量色散的一阶导数和二阶导数甚至更高阶导数为零,从而设计反射器的尺寸、电压等。当飞行时间对离子能量色散的一阶导数为零时,就可称为一阶聚焦;当飞行时间对离子能量色散的一阶和二阶导数同时为零时,就可称为二阶聚焦。如图4和图5分别是一阶聚焦和二阶聚焦时,离子飞行时间与离子能量的关系示意图。其特点是在该中心离子能量附近,飞行时间相对于离子能量的变化率非常小,因此在较小能量范围内,离子的飞行时间差很小,因此能够实现很高的分辨率。但由于在该中心能量之前或之后,飞行时间相对于能量的变化是单调的,因此,在离子能量进一步偏离中心能量的时候,总的飞行时间差将迅速增大,致使分辨率迅速降低。
发明内容
为了解决离子在较大能量色散下的飞行时间聚焦,提高飞行时间质谱的分辨率,本发明提供了一种具有宽能量聚焦反射器的飞行时间质谱仪。该飞行时间质谱仪可以在相对能量色散很高的情况下,实现较高的分辨率。
本发明的技术解决方案是提供一种具有宽能量聚焦反射器的飞行时间质谱仪,包括离子源、加速区、自由漂移区、反射器及探测器,其特征在于:所述反射器包括依次间隔距离放置的多片电极,每两片电极组成一级电场;
通过下述步骤计算每一级电场的长度和电压:
1)确定基本参数
确定自由漂移区长度为L,确定离子能量范围为[Uminq,Umaxq],确定要实现的分辨率为R;其中q为元电荷,Umin=U+Uini_min,Umax=U+Uini_max,U为加速场的加速电压,Uini_min为离子的最小初始能量,Uini_max为离子的最大初始能量;
2)设反射器的初始级数i=2;
3)通过使离子飞行时间对离子能量色散的二阶导数为零,确定反射器前两级电场的长度D1和D2、电势差U1和U2、电场强度E1和E2以及能量为Umin的离子的飞行时间t0
具体通过下述步骤计算:
3.1)确定反射器第一级长度(即第一级电场的长度)D1=αL,其中α<0.25,本方法取0.01;
3.2)根据公式(1)计算反射器的第一级电势差
3.3)根据公式(2)计算反射器的第二级电场强度
3.4)计算反射器的第二级电场电势差U2=Umin-U1,计算反射器的第二级电场长度
3.5)确定能量为Umin的离子的飞行时间t0
4)计算离子的最大飞行时间差
5)当时,结束计算;否则,增加反射器级数i=i+1,通过下述步骤计算该级电场长度Di、电势差Ui和电场强度Ei
6)令U′max=2U′min,Emin=Ei-1,Emax=2Ei-1
7)
计算Ei=Emiddle时,能量介于[U'minq,U'max q]的离子的最小飞行时间tmin
当tmin>t0-Δt时,取Emin=Emiddle,否则,取Emax=Emiddle
当Emax-Emin<0.0001时,进行步骤8),否则返回步骤7);
8)
在(U'minq,U'max q)中寻找飞行时间为t0的离子能量值Uxq;
Ui=Ux
返回步骤5),当时,结束计算;否则继续增加反射器的级数。
上述所有电极的圆孔上均设置有栅网。
依据本方法设计的反射器实际上只需保留前两级栅网即可。
在计算过程中,第三级及之后的电场长度采用向下取整的原则保留至mm量级。
本发明的有益效果是:
1、本发明反射器采用多级电场设计,每一级电场通过相邻两片电极加载一定电压差,间隔一定距离放置实现,通过设计每一级电场的尺寸(长度)和电压差,使在要聚焦的能量范围内,离子的飞行时间不再是随着能量的增加或减小而单调变化,而是如图6和图7所示,始终控制在一定范围内,从而可以在很宽的能量范围内实现较小的飞行时间差,进而实现较高的分辨率。
2、对电场尺寸的精度要求较低,除前两级实现二阶聚焦,需精确到0.1mm外,其余尺寸均只需精确到mm量级即可。
附图说明
图1为反射器各级电场长度和电势差计算总流程图;
图2为反射器前两级电场长度和电势差计算流程图;
图3为反射器第二级及之后电场长度和电势差计算流程图;
图4为传统聚焦方式,奇数阶聚焦示例;
图5为传统聚焦方式,偶数数阶聚焦示例;
图6为采用宽能量聚焦反射器,在自由漂移区长1m,针对离子[2000eV,3000eV]的能量范围,实现分辨率4000时,离子飞行时间相对于离子能量的关系;
图7为降低装配精度至mm量级(除第二级电场以外),在自由漂移区长1m,针对离子[2000eV,3000eV]的能量范围,实现分辨率4000时,离子飞行时间相对于离子能量的关系;
图8为宽能量聚焦多级反射器结构示意图;
图9为反射器在飞行时间质量分析器中安装示意图;
图10为只保留前两级栅网的宽能量反射器实现示意图;
图11针只保留前两级栅网,经微调电压改善电场渗透后,离子飞行时间相对于离子能量的关系。
具体实施方式
以下结合附图及实施例对本发明做进一步的描述。
如图8所示,本发明宽能量聚焦多级反射器的结构由多级电场组成,每级电场由两片带栅网的电极控制。图9为该多级反射器在飞行时间质谱仪中的安装示意图,根据图1、2和3流程图设计各级电场长度及电势差,
一、确定基本参数,包括自由漂移区长度为L,离子能量范围[Uminq,Umaxq],及要实现的分辨率为R;
二、假设反射器的级数为两级,即i=2;
通过使离子飞行时间对离子能量色散的二阶导数为零,确定反射器前两级电场的长度D1和D2、电势差U1和U2、电场强度E1和E2以及能量为Umin的离子的飞行时间t0
具体为:
1)确定反射器第一级长度D1=αL,其中α<0.25,本方法取0.01;
2)根据公式(1)计算反射器的第一级电势差
3)根据公式(2)计算反射器的第二级电场强度
4)计算反射器的第二级电场电势差U2=Umin-U1,计算反射器的第二级电场长度
5)根据运动过程中匀加速、匀减速、匀速运动规律确定能量为Uminq的离子的飞行时间t0
t0计算步骤:
取离子质量数M=150(注:以下质量数取相同值,该值任意取,最终设计结果与该值无关),则可以计算离子在自由漂移区的飞行速度v1
其中mu为一个质子或一个中子的质量,取为1.67e-27Kg。
则离子在自由漂移区中的飞行时间t1=L/v1
离子在第一级反射场中的加速度
则离子飞出该场时的速度v2可通过下式计算
则离子在该场中的飞行时间t2=2(v1-v2)/a1
离子从第二级反射场反射回去,在该场中离子的加速度为
在该场中的飞行时间为t3=2v2/a2
t0=t1+t2+t3
三、在确定了前两级电场的各参数时,接着通过公式(4)计算满足分辨率要求时,离子的最大飞行时间差;
四、然后计算结束计算;否则,增加反射器级数i=i+1,通过下述步骤计算该级电场长度Di、电势差Ui和电场强度Ei
a)令U′max=2U′min,Emin=Ei-1,Emax=2Ei-1
b)
计算Ei=Emiddle时,能量介于[U'minq,U'max q]的离子的最小飞行时间tmin
计算过程中每次增加0.01eV,计算飞行时间t',并从中得知最小飞行时间;
飞行时间t'的计算过程如下:
先计算离子在自由漂移区的飞行速度v1
其中mu为一个质子或一个中子的质量,取为1.67e-27Kg。U'q为计算过程中离子的能量。
则离子在自由漂移区中的飞行时间t1=L/v1
离子在第j(j=1,2,…,i-1)级反射场中的加速度
则离子飞出该场时的速度vj+1
可通过下式计算
则离子在该场中的飞行时间tj+1=2(vj-vj+1)/aj
离子从第i级反射场反射回去,在该场中离子的加速度为
在该场中的飞行时间为ti+1=2vi/ai
当tmin>t0-Δt时,取Emin=Emiddle,否则,取Emax=Emiddle
当Emax-Emin<0.0001时,进行步骤c),否则返回步骤b);
c)
在(U'minq,U'max q)中寻找飞行时间为t0的离子能量值Uxq;
Ui=Ux
五、返回步骤四。
其离子飞行时间与离子能量的关系不再如传统的聚焦方式那样,随着能量色散的进一步增大,离子的飞行时间差迅速增大,而是始终控制在一定范围内,从而针对较宽能量分布的离子实现较佳的能量时间聚焦,进而实现较高的分辨率。
实施例一
本实施例针对分辨率为3000、4000、5000、6000和7000,利用上述计算方法分别计算了反射器各级电场长度相对于自由漂移区长度的比例关系和反射器各级电场电势差(V)相对于要聚焦的离子的最低能量(eV)的比例关系如表1所示。在计算中,计算了15阶电场的长度与电势差,实际阶数可根据能量范围的大小进行选取。
表1针对不同的分辨率,宽能量聚焦反射器各级电场长度相对于自由漂移区长度的关系
以及各级电场电势差相对于聚焦的离子最小能量(单位:eV)的关系
假设需要聚焦的离子能量范围为[2000eV,3000eV],需要实现的分辨率为4000,自由漂移区长度为1m。则根据表1分辨率为4000时的各级电场电压差相对于聚焦范围内离子最小能量(此处为2000eV)的百分比和各级电场长度相对于自由漂移区(此处为1m)的百分比,可以计算得到各级电场长度、电势差以及极板电压如表2所示。且可以计算得到需要13级电场即可(因第13级电极电压已经大于3000V,此时3000eV的离子被该级电场反射了回去,因此后续电场不起作用)。根据带电粒子在电场中运动规律,可以计算得到离子飞行时间和能量的关系如图6所示。
表2针对离子能量范围[2000eV,3000eV],需要实现的分辨率为4000,自由漂移区长度为1m的情况,计算得到的各级电场长度、电势差及各极板电压。
实施例二
如图10所示,本实施例二与实施例一不同的是该反射器只保留了前两级栅网,可以显著提高离子的传输效率,后续电场间微弱的电场渗透可以通过微调电压予以补偿。此时,针对分辨率为3000、4000、5000、6000和7000,依据图1、2和3计算流程可计算得到,反射器各级电场长度相对于自由漂移区长度的比例关系和反射器各级电场电势差(V)相对于要聚焦的离子的最低能量(eV)的比例关系如表3所示。由于去掉了一些栅网,导致相邻电场之间存在微弱渗透,因此,该计算电压非最终电压,还需经过微调来补偿电场渗透带来的影响。
表3降低装配精度为mm级,针对不同的分辨率,宽能量聚焦反射器各级电场长度相对于自由漂移区长度的关系以及各级电场电势差相对于聚焦的离子最小能量(单位:eV)的关
假设需要聚焦的离子能量范围为[2000eV,3000eV],需要实现的分辨率为4000,自由漂移区长度为1m。则根据表3分辨率为4000时的各级电场电压差相对于聚焦范围内离子最小能量(此处为2000eV)的百分比和各级电场长度相对于自由漂移区(此处为1m)的百分比,可以计算得到各级电场长度、电势差以及极板电压如表4所示。且可以计算得到需要13级电场即可。根据带电粒子在电场中运动规律,可以计算得到离子飞行时间和能量的关系如图7所示。
表4电场长度精确到mm量级(除第二级电场长度精确到0.1mm)时,宽能量聚焦反射器尺寸及应加电压实例(自由漂移区长度L=1m,能量聚焦范围[2000eV,3000eV],实现分辨率为4000)
电场序号 长度/cm 电势差/V 极板电压/V
1 1 1360.00 1360.00
2 4.73 640.00 2000.00
3 1.7 230.37 2230.37
4 0.8 108.79 2339.16
5 0.7 95.53 2434.69
6 0.6 82.24 2516.93
7 0.5 68.81 2585.74
8 0.5 69.06 2654.80
9 0.5 69.37 2724.17
10 0.5 69.67 2793.84
11 0.5 70 2863.84
12 0.5 70.35 2934.19
13 0.5 70.7 3004.89
由于第三级电场之后已经去掉了栅网,后续电场间存在微弱电场渗透,这一状况可以通过依次微调后续电极电压予以改善。在调节之后,离子的飞行时间和离子的能量关系如图11所示。此时实现分辨率约为

Claims (5)

1.一种具有宽能量聚焦反射器的飞行时间质谱仪,包括离子源、加速区、自由漂移区、反射器及探测器,其特征在于:所述反射器包括依次间隔距离放置的多片电极,所述电极中心设置有圆孔,每两片电极组成一级电场;
通过下述步骤计算每一级电场的长度和电压:
1)确定基本参数
确定自由漂移区长度为L,确定离子能量范围为[Uminq,Umaxq],确定要实现的分辨率为R;其中q为元电荷,Umin=U+Uini_min,Umax=U+Uini_max,U为加速场的加速电压,Uini_min为离子的最小初始能量,Uini_max为离子的最大初始能量;
2)设反射器的初始级数i=2;
3)通过使离子飞行时间对离子能量色散的二阶导数为零,确定反射器前两级电场的长度D1和D2、电势差U1和U2、电场强度E1和E2以及能量为Umin的离子的飞行时间t0
4)计算离子的最大飞行时间差
Δ t = t 0 2 R
5)当时,结束计算;否则,增加反射器级数i=i+1,通过下述步骤计算该级电场长度Di、电势差Ui和电场强度Ei
6)令U′max=2U′min,Emin=Ei-1,Emax=2Ei-1
7)则
E m i d d l e = E m a x + E min 2
计算Ei=Emiddle时,能量介于[U'minq,U'maxq]的离子的最小飞行时间tmin
当tmin>t0-Δt时,取Emin=Emiddle,否则,取Emax=Emiddle
当Emax-Emin<0.0001时,进行步骤8),否则返回步骤7);
8)
E i = E m a x + E m i n 2
在(U'minq,U'maxq)中寻找飞行时间为t0的离子能量值Uxq;
Ui=Ux
D i = U i E i
9)返回步骤5)。
2.根据权利要求1所述的具有宽能量聚焦反射器的飞行时间质谱仪,其特征在于:步骤3)具体为:
1)确定反射器第一级长度D1=αL,其中α<0.25,本方法取0.01;
2)根据公式(1)计算反射器的第一级电势差
U 1 = 4 3 ( 1 2 + D 1 L ) U min - - - ( 1 )
3)根据公式(2)计算反射器的第二级电场强度
E 2 = U m i n L 2 ( 1 - U 1 U min ) - 1 / 2 1 2 - 2 D 1 L U min U 1 ( 1 - ( 1 - U 1 U min ) - 1 / 2 ) - - - ( 2 )
4)计算反射器的第二级电场电势差U2=Umin-U1,计算反射器的第二级电场长度
D 2 = U 2 E 2 - - - ( 3 )
5)确定能量为Umin的离子的飞行时间t0
3.根据权利要求1所述的具有宽能量聚焦反射器的飞行时间质谱仪,其特征在于:所述所有电极的圆孔上均设置有栅网。
4.根据权利要求1所述的具有宽能量聚焦反射器的飞行时间质谱仪,其特征在于:所述前两片电极上设置有栅网。
5.根据权利要求1所述的具有宽能量聚焦反射器的飞行时间质谱仪,其特征在于:步骤6)和步骤7)中的计算结果保留至mm量级。
CN201611028303.XA 2016-11-18 2016-11-18 一种具有宽能量聚焦反射器的飞行时间质谱仪 Active CN106653559B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611028303.XA CN106653559B (zh) 2016-11-18 2016-11-18 一种具有宽能量聚焦反射器的飞行时间质谱仪

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611028303.XA CN106653559B (zh) 2016-11-18 2016-11-18 一种具有宽能量聚焦反射器的飞行时间质谱仪

Publications (2)

Publication Number Publication Date
CN106653559A true CN106653559A (zh) 2017-05-10
CN106653559B CN106653559B (zh) 2018-06-26

Family

ID=58807389

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611028303.XA Active CN106653559B (zh) 2016-11-18 2016-11-18 一种具有宽能量聚焦反射器的飞行时间质谱仪

Country Status (1)

Country Link
CN (1) CN106653559B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107219214A (zh) * 2017-07-26 2017-09-29 大连理工大学 一种光谱结合质谱的未知样品中元素的定量分析装置
CN108811293A (zh) * 2018-10-08 2018-11-13 中国科学院上海高等研究院 带电粒子加速电场产生装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101176185A (zh) * 2005-05-11 2008-05-07 埃美格科学仪器公司 反射器
WO2011019457A1 (en) * 2009-08-11 2011-02-17 Regents Of The University Of California Time-of-flight electron energy analyzer
CN105378891A (zh) * 2013-07-10 2016-03-02 布鲁克道尔顿有限公司 带Cassini反射器的飞行时间质谱仪
CN105529240A (zh) * 2014-10-16 2016-04-27 布鲁克道尔顿有限公司 利用宽质量范围空间聚焦的飞行时间质谱仪

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101176185A (zh) * 2005-05-11 2008-05-07 埃美格科学仪器公司 反射器
WO2011019457A1 (en) * 2009-08-11 2011-02-17 Regents Of The University Of California Time-of-flight electron energy analyzer
CN105378891A (zh) * 2013-07-10 2016-03-02 布鲁克道尔顿有限公司 带Cassini反射器的飞行时间质谱仪
CN105529240A (zh) * 2014-10-16 2016-04-27 布鲁克道尔顿有限公司 利用宽质量范围空间聚焦的飞行时间质谱仪

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107219214A (zh) * 2017-07-26 2017-09-29 大连理工大学 一种光谱结合质谱的未知样品中元素的定量分析装置
CN107219214B (zh) * 2017-07-26 2023-10-27 大连理工大学 一种光谱结合质谱的未知样品中元素的定量分析装置
CN108811293A (zh) * 2018-10-08 2018-11-13 中国科学院上海高等研究院 带电粒子加速电场产生装置
CN108811293B (zh) * 2018-10-08 2019-01-01 中国科学院上海高等研究院 带电粒子加速电场产生装置

Also Published As

Publication number Publication date
CN106653559B (zh) 2018-06-26

Similar Documents

Publication Publication Date Title
Dowell Sources of Emittance in RF Photocathode Injectors: Intrinsic emittance, space charge forces due to non-uniformities, RF and solenoid effects
CN106653559A (zh) 一种具有宽能量聚焦反射器的飞行时间质谱仪
Fraser The gain, temporal resolution and magnetic-field immunity of microchannel plates
JP2014049196A (ja) イオン移動度分離装置
CN113606103B (zh) 一种阶梯栅极、栅极结构、阶梯栅极参数确定方法及系统
Dietz et al. Single‐Electron Response of a Porous KCl Transmission Dynode and Application of Polya Statistics to Particle Counting in an Electron Multiplier
Peterson et al. Measurements of energy and angular distributions of secondary electrons produced in electron-impact ionization of helium
CN105702558B (zh) 一种消除飞行时间质谱仪离子初始位置分散的方法
CN203733758U (zh) 一种离子源电参数质量歧视自动校正装置
CN114171367A (zh) 飞行时间质量分析器
CN110390136A (zh) 一种自适应扩展光电倍增管动态范围的方法
Hvelplund et al. Equilibrium charge distributions of light ions in helium, measured with a position-sensitive open electron multiplier
Egberts et al. Detailed experimental characterization of an ionization profile monitor
CN117012608A (zh) 摆线质谱仪及其分辨率调节方法
KR101446083B1 (ko) ExB 탐침의 신호 처리 방법
Freyou et al. Recoil-ion–projectile-ion–cusp-electron triple-coincidence measurements for 1-MeV/u bare and one-electron O and C projectiles on Ne
Belov et al. High responsivity secondary ion energy analyzer
Zarcone et al. Further development with heavy ion sources at Brookhaven National Laboratory's Tandem van de Graaff Facility
CN111965689B (zh) 一种用于中性原子分析的测量装置
Stockli et al. Developments on the KSU-CRYEBIS, a user facility for low energy, highly charged ions
CN117871976A (zh) 一种离子推力器放电室性能等效测试方法
Warwick et al. ACCELERATION, CORRENT AMPLIFICATION AND EMITTANCE IN MBE-4, AN EXPERIMENTAL BEAM INDUCTION LINEAR ACCELERATOR FOR HEAVY IONS
CN118073171A (zh) 解决反应管离子透过率低并提高分子离子反应时间的系统
Bekmirzaeva et al. R. Bekmirzaev
SU1457180A1 (ru) Изохронный циклотрон

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