CN115754845B - 基于矢量光调制的原子磁强计空间磁场成像装置及方法 - Google Patents

基于矢量光调制的原子磁强计空间磁场成像装置及方法 Download PDF

Info

Publication number
CN115754845B
CN115754845B CN202211560528.5A CN202211560528A CN115754845B CN 115754845 B CN115754845 B CN 115754845B CN 202211560528 A CN202211560528 A CN 202211560528A CN 115754845 B CN115754845 B CN 115754845B
Authority
CN
China
Prior art keywords
space
light
magnetic field
alkali metal
vector
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
CN202211560528.5A
Other languages
English (en)
Other versions
CN115754845A (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 Lab
Original Assignee
Zhejiang Lab
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 Lab filed Critical Zhejiang Lab
Priority to CN202211560528.5A priority Critical patent/CN115754845B/zh
Publication of CN115754845A publication Critical patent/CN115754845A/zh
Application granted granted Critical
Publication of CN115754845B publication Critical patent/CN115754845B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Measuring Magnetic Variables (AREA)

Abstract

本发明公开了一种基于矢量光调制的原子磁强计空间磁场成像装置及方法,该方法使用空间光调制器将通过碱金属气室的抽运激光调制成空间圆锥螺旋线矢量光,并用光束质量分析仪测量抽运激光的光强通过气室前后在平面投影上的衰减信息,利用圆偏振态光强衰减公式,拟合出原子自旋极化率,不断改变圆锥螺旋线的参数,使矢量光扫过气室内所有三维像素点,得到极化率在空间内的分布,进一步利用布洛赫方程的稳态解解算出外部磁场值,得到磁场在三维空间上微米像素量级的成像结果,本发明理论依据充分合理,实验操作简单,可以实现极高空间分辨率的磁场测量,未来有望服务于生物医疗,基础理论研究等多个领域。

Description

基于矢量光调制的原子磁强计空间磁场成像装置及方法
技术领域
本发明属于磁场技术领域,具体是一种基于矢量光调制的原子磁强计空间磁场成像装置及方法。
背景技术
磁场成像技术在天文观测、医疗病灶探查、生物磁场分析、基础物理研究、电子工业等多个领域有着应用,目前常用的磁成像技术包括核磁共振成像、微波磁成像等,但精度有限,且都是二维平面成像,所以在应用上有一定限制,例如核磁共振在对人体进行检查时,要求病人不能带有心脏起搏器和金属异物。因此,需要一种能够实现磁场三维立体成像的方法和装置。在多种磁传感器中,原子磁场计受到了广泛的关注,它利用了原子在外磁场作用下的磁矩旋转效应,因而具有较高的灵敏度,使用时,需要用抽运激光将碱金属原子极化从而实现原子自旋的宏观指向,因此,可以利用圆偏光在碱金属气室内被吸收而导致光强衰减的特性测量出原子的极化率分布,进而测量出外部磁场的大小。但这种方法只能测量出气室外部磁场的平均值,无法实现在较高空间分辨率上的磁场分布。
光场调控技术可以实现对激光振幅、相位、功率密度、偏振态等重要参数的调控,架构出非均匀分布的空间矢量光场。这样,可以更灵活方便的操控光场,生成多类新型空间结构光场,如无衍射光束、矢量涡旋光束等,可以应用到全息成像、光信息存储、激光加工等多个领域。如果将光场调控技术应用到原子磁强计中,对圆偏振抽运光进行矢量调控,使其携带极化率的空间位置信息,就能通过测量极化率在三维空间的分布,实现磁场的三维成像。
发明内容
为了解决现有技术中的上述问题,办发明提出了一种基于矢量光调制的原子磁强计空间磁场成像装置及方法,基于空间光调制器将原子磁强计抽运光的调制成圆锥螺旋矢量光,利用碱金属气室对光的吸收作用,使用高分辨率的光束质量分析仪对光强的变化进行采集,通过获取气室前后圆偏振抽运光光强衰减信息进行极化率测量,用上位机对采集的数据进行拟合和处理,从而实现磁场的三维成像。
本发明提供了一种基于矢量光调制的原子磁强计空间磁场成像装置,包括沿光路方向依次设置的激光器,扩束镜、空间光调制器,分光棱镜、起偏器、四分之一波片、碱金属原子气室、第二光束质量分析仪以及位于分光棱镜另一束光方向的第一光束质量分析仪,所述第一光束质量分析仪和第二光束质量分析仪均电连接上位机数据处理系统;所述的碱金属原子气室外部由内向外依次设置无磁电加热烤箱、隔热保温材料腔、磁补偿线圈、磁屏蔽桶;
所述激光器发射出的抽运激光经过扩束镜将光斑放大,然后进入空间光调制器调控成微米像素级别的空间矢量光,调控后的激光通过分光棱镜按比例分为两束,一束由第一光束质量分析仪接收对入射光光强进行检测,另一束激光经过偏振元件和四分之一波片后确保以圆偏振态入射到经过高温加热的碱金属原子气室中,经过碱金属原子气室的激光被第二光束质量分析仪接收监控出射光光强。
作为优选,所述空间光调制器的分辨率为1024*1024个单元,每个单元的大小为12.5um。
作为优选,所述的碱金属气室内的气体为钾、铷、铯其中的任意一种或两种,且处于SERF态。
作为优选,所述无磁电加热烤箱采用氮化硼材料制成。
本发明还提供了一种基于矢量光调制的原子磁强计空间磁场成像方法,包括以下步骤:
S1:搭建调整光路,开启激光器,开启无磁电加热烤箱使碱金属原子气室内的碱金属原子处于SERF态,令原子磁强计处于测量工作状态,令原子磁强计处于测量工作状态;
S2:用空间光调制器将抽运激光调制成空间螺旋线形状;
S3:用光束质量分析仪测量空间圆锥螺旋线矢量光光强在通过气室前后在二维平面上的投影;
S4:用空间光调制器改变空间圆锥螺旋线矢量光的初始光强,重复步骤S3;
S5:根据圆偏振态光束穿过碱金属原子气室时光强衰减公式和极化率公式,拟合计算出自旋驰豫率和极化率,并计算出抽运率;
S6:改变圆锥螺旋线相关参数,重复步骤S3至S5,直到测量出三维空间所有像素点上的极化率;
S7:根据Bloch方程的稳态解,代入极化率计算出磁场在三维空间的分布。
作为优选,所述空间螺旋线形状包括圆锥螺旋线形状、斐波那契螺旋线形状、边长变化的方形螺旋线在二维平面投影没有重叠点的三维空间曲线形状。
作为优选,所述空间光调制器的分辨率为1024*1024个单元,每个单元的大小为12.5um。
本发明通过使用空间光调制器将通过碱金属气室的抽运激光调制成空间圆锥螺旋线矢量光,并用光束质量分析仪测量抽运激光的光强通过气室前后在平面投影上的衰减信息,利用圆偏振态光强衰减公式,拟合出原子自旋极化率,不断改变圆锥螺旋线的参数,使矢量光扫过气室内所有三维像素点,得到极化率在空间内的分布,进一步利用布洛赫方程的稳态解解算出外部磁场值,得到磁场在三维空间上微米像素量级的成像结果,理论依据充分合理,实验操作简单,可以实现极高空间分辨率的磁场测量,未来有望服务于生物医疗,基础理论研究等多个领域。
附图说明
图1为本发明的装置结构示意图;
附图1标记列示如下:1-抽运光激光器,2-扩束镜,3-空间光调制器, 4-分光棱镜,5、13-光束质量分析仪,6-起偏器,7-四分之一波片,8-磁补偿线圈,9-隔热保温材料腔,10-无磁电加热烤箱,11-碱金属原子气室,12-磁屏蔽桶,14-上位机数据处理系统
图2为本发明的实施流程图;
图3为空间光调制器调制的圆锥螺旋矢量光形状图;
图4为空间光调制器调制的圆锥螺旋矢量光在光束质量分析仪上的投影图;
图中:1-抽运光激光器,2-扩束镜,3-空间光调制器, 4-分光棱镜,5-第一光束质量分析仪,6-起偏器,7-四分之一波片,8-磁补偿线圈,9-隔热保温材料腔,10-无磁电加热烤箱,11-碱金属原子气室,12-磁屏蔽桶,13-第二光束质量分析仪,14-上位机数据处理系统。
具体实施方式
下面通过附图以及具体实施方式进一步说明本发明。
如图1所示,本发明提出一种基于矢量光调制的原子磁强计空间磁场成像装置,在不需要改变光路的前提下,可以测量出磁场在三维空间上的成像。装置包括沿光路方向依次设置的激光器1,扩束镜2、空间光调制器3,分光棱镜4、起偏器6、四分之一波片7、碱金属原子气室11、第二光束质量分析仪13以及位于分光棱镜4另一束光方向的第一光束质量分析仪5,所述第一光束质量分析仪5和第二光束质量分析仪13均电连接上位机数据处理系统14;所述的碱金属原子气室11外部由内向外依次设置无磁电加热烤箱10、隔热保温材料腔9、磁补偿线圈8、磁屏蔽桶12;所述的碱金属气室内的气体为钾、铷、铯其中的任意一种或两种,且处于SERF态。所述无磁电加热烤箱采用氮化硼材料制成。
所述激光器1发射出的抽运激光经过扩束镜2将光斑放大,然后进入空间光调制器3调控成微米像素级别的空间矢量光,调控后的激光通过分光棱镜4按比例分为两束,一束由第一光束质量分析仪5接收对入射光光强进行检测,另一束激光经过偏振元件6和四分之一波片7后确保以圆偏振态入射到经过高温加热的碱金属原子气室11中,经过碱金属原子气室11的激光被第二光束质量分析仪13接收监控出射光光强。
本发明的具体实施步骤包括:
S1:搭建调整光路,开启激光器,开启无磁电加热烤箱使碱金属原子气室内的碱金属原子处于SERF态,令原子磁强计处于测量工作状态,令原子磁强计处于测量工作状态;
S2:用空间光调制器将抽运激光调制成空间螺旋线形状;
S3:用光束质量分析仪测量空间圆锥螺旋线矢量光光强在通过气室前后在二维平面上的投影;
S4:用空间光调制器改变空间圆锥螺旋线矢量光的初始光强,重复步骤S3;
S5:根据圆偏振态光束穿过碱金属原子气室时光强衰减公式和极化率公式,拟合计算出自旋驰豫率和极化率,并计算出抽运率;
S6:改变圆锥螺旋线相关参数,重复步骤S3至S5,直到测量出三维空间所有像素点上的极化率;
S7:根据Bloch方程的稳态解,代入极化率计算出磁场在三维空间的分布。
其中,所述空间螺旋线形状包括圆锥螺旋线形状、斐波那契螺旋线形状、边长变化的方形螺旋线在二维平面投影没有重叠点的三维空间曲线形状。
所述空间光调制器的分辨率为1024*1024个单元,每个单元的大小为12.5um。
实施例1
如图2所示:
(1)调整光路,开启加热烤箱,使原子磁强计处于工作状态;
(2)用空间光调制器3使矢量抽运光在气室内成为空间圆锥螺旋线矢量光(如图3所示):激光器1发射出的抽运激光经过扩束镜2将光斑放大,进入空间光调制器3调控成微米像素级别的圆锥螺旋线矢量光,螺旋线上各点光强一致,幅值为I0,螺旋线上任意一点的坐标为
其中,a,b,c是螺旋线的形状参数,t是位置相关自变量。由公式可知,在螺旋线形状确定的条件下,若已知某点的任意一个方向的坐标,可以计算出另外两个方向的坐标值。如图3所示,定义z轴的方向是光传播的方向(即穿过气室的方向)。
(3)用光束质量分析仪检测入射光和出射光/>光强的二维平面分布:
圆偏振抽运光在碱金属气室中被吸收,沿着其传播方向光强会衰减。则光束质量分析仪13接收到的圆锥螺旋线在xy平面上的投影(如图4所示),其螺旋线上各点的光强值为:
是自旋弛豫率,是未知量;/>为光吸收截面,是已知量;n是原子数密度,是未知量。
(4)改变圆锥螺旋线矢量光的光强,重复步骤(3):得到一系列和出射光/>,带入步骤(3)中公式,得到光强衰减在三维空间的信息/>并拟合出/>
(5)根据公式拟合出圆锥螺旋线所有像素点上的极化率和自旋弛豫率,并计算出抽运率:
光传播方向极化率、自旋驰豫率和光强的关系式:
根据拟合得到的,可以计算出/>的值。
计算出抽运率:
h 是普朗克常量, P 为抽运激光功率,为光传播速度,/>为经典电子半径,/>为抽运光频率,/>为抽运激光光斑面积,/>为碱金属原子 D1 线的振荡强度,/>为原子 D1 线共振频率,/>为碱金属原子 D1 线压力展宽的线宽。
(6)改变圆锥螺旋线的a,b,c参数,重复测量,直到测量出气室内三维空间上所有像素点上的
(7)根据Bloch方程的稳态解,抽运方向的极化率还可以表示为:
是旋磁比参数,由上式可知,在已知/>,/>,/>的情况下,可以计算出,即y方向的磁场在三维空间的分布,实现空间磁场成像。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围内。

Claims (7)

1.一种基于矢量光调制的原子磁强计空间磁场成像装置,其特征在于:包括沿光路方向依次设置的激光器,扩束镜、空间光调制器,分光棱镜、起偏器、四分之一波片、碱金属原子气室、第二光束质量分析仪以及位于分光棱镜另一束光方向的第一光束质量分析仪,所述第一光束质量分析仪和第二光束质量分析仪均电连接上位机数据处理系统;所述的碱金属原子气室外部由内向外依次设置无磁电加热烤箱、隔热保温材料腔、磁补偿线圈、磁屏蔽桶;
所述激光器发射出的抽运激光经过扩束镜将光斑放大,然后进入空间光调制器调控成微米像素级别的空间矢量光,调控后的激光通过分光棱镜按比例分为两束,一束由第一光束质量分析仪接收对入射光光强进行检测,另一束激光经过偏振元件和四分之一波片后确保以圆偏振态入射到经过高温加热的碱金属原子气室中,经过碱金属原子气室的激光被第二光束质量分析仪接收监控出射光光强;
所述成像装置采用的成像方法,具体步骤如下:
S1:搭建调整光路,开启激光器,开启无磁电加热烤箱使碱金属原子气室内的碱金属原子处于SERF态,令原子磁强计处于测量工作状态;
S2:用空间光调制器将抽运激光调制成空间螺旋线形状;
S3:用光束质量分析仪测量空间圆锥螺旋线矢量光光强在通过气室前后在二维平面上的投影;
S4:用空间光调制器改变空间圆锥螺旋线矢量光的初始光强,重复步骤S3;
S5:根据圆偏振态光束穿过碱金属原子气室时光强衰减公式和极化率公式,拟合计算出自旋驰豫率和极化率,并计算出抽运率;
S6:改变圆锥螺旋线相关参数,重复步骤S3至S5,直到测量出三维空间所有像素点上的极化率;
S7:根据Bloch方程的稳态解,代入极化率计算出磁场在三维空间的分布。
2.根据权利要求1所述的一种基于矢量光调制的原子磁强计空间磁场成像装置,其特征在于:所述空间光调制器的分辨率为1024*1024个单元,每个单元的大小为12.5um。
3.根据权利要求1所述的一种基于矢量光调制的原子磁强计空间磁场成像装置,其特征在于:所述的碱金属原子气室内的气体为钾、铷、铯其中的任意一种或两种。
4.根据权利要求2所述的一种基于矢量光调制的原子磁强计空间磁场成像装置,其特征在于:所述碱金属原子气室内的气体处于SERF态。
5.根据权利要求1所述的一种基于矢量光调制的原子磁强计空间磁场成像装置,其特征在于:所述无磁电加热烤箱采用氮化硼材料制成。
6.根据权利要求1所述的一种基于矢量光调制的原子磁强计空间磁场成像装置,其特征在于:所述空间螺旋线形状包括圆锥螺旋线形状、斐波那契螺旋线形状、边长变化的方形螺旋线在二维平面投影没有重叠点的三维空间曲线形状。
7.根据权利要求1所述的一种基于矢量光调制的原子磁强计空间磁场成像装置,其特征在于:所述空间光调制器的分辨率为1024*1024个单元,每个单元的大小为12.5um。
CN202211560528.5A 2022-12-07 2022-12-07 基于矢量光调制的原子磁强计空间磁场成像装置及方法 Active CN115754845B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211560528.5A CN115754845B (zh) 2022-12-07 2022-12-07 基于矢量光调制的原子磁强计空间磁场成像装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211560528.5A CN115754845B (zh) 2022-12-07 2022-12-07 基于矢量光调制的原子磁强计空间磁场成像装置及方法

Publications (2)

Publication Number Publication Date
CN115754845A CN115754845A (zh) 2023-03-07
CN115754845B true CN115754845B (zh) 2023-09-01

Family

ID=85344728

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211560528.5A Active CN115754845B (zh) 2022-12-07 2022-12-07 基于矢量光调制的原子磁强计空间磁场成像装置及方法

Country Status (1)

Country Link
CN (1) CN115754845B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001092949A1 (fr) * 2000-06-01 2001-12-06 Optware Corporation Modulateur spatial de lumiere
CN104777625A (zh) * 2015-04-21 2015-07-15 东南大学 一种生成椭圆偏振矢量光场的装置
CN105929458A (zh) * 2016-03-21 2016-09-07 吉林大学 航空磁场矢量检测装置及监测方法
JP2017040882A (ja) * 2015-08-21 2017-02-23 住友大阪セメント株式会社 光変調器
CN108061975A (zh) * 2017-12-07 2018-05-22 西北工业大学 一种高效产生任意矢量光场的方法及装置
CN111983526A (zh) * 2020-08-24 2020-11-24 清华大学 原子磁强计及磁场成像系统
CN113639883A (zh) * 2021-08-18 2021-11-12 之江实验室 一种碱金属原子磁强计自旋极化率空间分布原位测量系统及方法
CN113777106A (zh) * 2021-09-14 2021-12-10 之江实验室 一种原子磁强计的碱金属原子数密度空间分布均匀性的测试系统及方法
CN114488513A (zh) * 2022-02-12 2022-05-13 之江实验室 一种全矢量调制的单光纤高信噪比三维成像方法及装置
CN114527415A (zh) * 2022-02-22 2022-05-24 上海理工大学 一种多通道全光法极弱磁场检测装置
CN114996984A (zh) * 2022-03-04 2022-09-02 北京航空航天大学 一种基于三维极化率分布的极化梯度弛豫模拟计算方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001092949A1 (fr) * 2000-06-01 2001-12-06 Optware Corporation Modulateur spatial de lumiere
CN104777625A (zh) * 2015-04-21 2015-07-15 东南大学 一种生成椭圆偏振矢量光场的装置
JP2017040882A (ja) * 2015-08-21 2017-02-23 住友大阪セメント株式会社 光変調器
CN105929458A (zh) * 2016-03-21 2016-09-07 吉林大学 航空磁场矢量检测装置及监测方法
CN108061975A (zh) * 2017-12-07 2018-05-22 西北工业大学 一种高效产生任意矢量光场的方法及装置
CN111983526A (zh) * 2020-08-24 2020-11-24 清华大学 原子磁强计及磁场成像系统
CN113639883A (zh) * 2021-08-18 2021-11-12 之江实验室 一种碱金属原子磁强计自旋极化率空间分布原位测量系统及方法
CN113777106A (zh) * 2021-09-14 2021-12-10 之江实验室 一种原子磁强计的碱金属原子数密度空间分布均匀性的测试系统及方法
CN114488513A (zh) * 2022-02-12 2022-05-13 之江实验室 一种全矢量调制的单光纤高信噪比三维成像方法及装置
CN114527415A (zh) * 2022-02-22 2022-05-24 上海理工大学 一种多通道全光法极弱磁场检测装置
CN114996984A (zh) * 2022-03-04 2022-09-02 北京航空航天大学 一种基于三维极化率分布的极化梯度弛豫模拟计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
自由空间中偏振调制光场的传输及控制;刘圣,李鹏,章毅,韩磊,程华超,赵建林;《光学学报》;第36卷(第10期);全文 *

Also Published As

Publication number Publication date
CN115754845A (zh) 2023-03-07

Similar Documents

Publication Publication Date Title
US11360135B2 (en) Atom-based electromagnetic field sensing element and measurement system
CN105158709B (zh) 一种基于内嵌nv‑色心金刚石的磁场测量装置
CN105652223B (zh) 磁场测量方法以及磁场测量装置
JP6222974B2 (ja) 光ポンピング磁力計および磁気センシング方法
JP5707021B2 (ja) 磁場計測装置
CN104215553B (zh) 一种碱金属蒸汽的原子密度与极化率一体化测量装置
CN108287322B (zh) 一种无响应盲区的原子磁力仪及其测量外磁场的方法
CN106886000B (zh) 一种利用核磁共振实现磁场幅度稳定的装置及方法
US20130082700A1 (en) Nuclear magnetic resonance imaging apparatus and nuclear magnetic resonance imaging method
CN108519566B (zh) 一种基于光频移调制的serf原子磁强计装置及方法
CN105651649A (zh) 一种适用于原子磁强计的原子密度实时在线测量方法
CN113777106B (zh) 一种原子磁强计的碱金属原子数密度空间分布均匀性的测试系统及方法
CN111025206A (zh) 一种基于原子磁共振的静磁场空间分布测量系统及方法
US11313925B2 (en) Quantum sensor based on rare-earth-ion doped optical crystal and use thereof
US10024931B2 (en) Magnetic field measurement method and magnetic field measurement apparatus
CN115754845B (zh) 基于矢量光调制的原子磁强计空间磁场成像装置及方法
US2913665A (en) Optical microwave power measuring system
Long et al. Speedy in-situ magnetic field compensation algorithm for multiple-channel single-beam SERF atomic magnetometers
RU2650713C1 (ru) Способ измерения малых коэффициентов оптического поглощения нелинейно-оптических кристаллов
JP6660842B2 (ja) 緩和時間測定方法及び磁気共鳴測定装置
CN112229801B (zh) 一种测量磁场作用下镱原子团双折射效应的装置及方法
CN209784511U (zh) 基于43Ca+离子测量微弱高频交变磁场的装置
CN113341235A (zh) 一种测量装置
CN113341353B (zh) 基于数字微镜阵列的小尺度空间磁场梯度测量系统及方法
CN115754844B (zh) 基于方波磁场幅值的原子磁强计极化率操控方法及装置

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