CN115754844B - 基于方波磁场幅值的原子磁强计极化率操控方法及装置 - Google Patents
基于方波磁场幅值的原子磁强计极化率操控方法及装置 Download PDFInfo
- Publication number
- CN115754844B CN115754844B CN202211560465.3A CN202211560465A CN115754844B CN 115754844 B CN115754844 B CN 115754844B CN 202211560465 A CN202211560465 A CN 202211560465A CN 115754844 B CN115754844 B CN 115754844B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- square wave
- frequency
- alkali metal
- light path
- 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
Links
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种基于方波磁场幅值的原子磁强计极化率操控方法及装置,该方法在碱金属原子气室外部使用两套三轴亥姆霍兹线圈,一套用于施加高频方波调制磁场,控制碱金属原子的极化率,一套用于施加低频标定方波磁场,用减慢效应法实时测量极化率,上位机软件采集磁强计的输出信号并控制高频方波磁场的幅值,实现极化率的闭环控制。系统主要包括抽运光激光器、检测光激光器及光路调整元件、磁屏蔽筒、一套高频方波调制磁场线圈、一套低频方波标定磁场线圈、碱金属原子气室、隔热保温材料腔、无磁电加热装置、上位机数据处理软件,本发明可以在不改变原子磁强计光路的基础上将极化率稳定在设定值,操作简单。
Description
技术领域
本发明属于磁场测量领域,尤其涉及一种基于方波磁场幅值的原子磁强计极化率操控方法及装置。
背景技术
碱金属原子的自旋极化率对研究原子自旋交换驰豫状态有着非常重要的作用,且极化率沿着抽运光方向的投影含有外界带测磁场的信息,这使得原子磁强计光旋角的输出公式中包含了极化率的值,经试验证明,极化率的值在0.5附近,磁强计达到最佳工作状态,所以精密的测量并控制极化率的值,对提高磁强计灵敏度起着非常重要的作用。
目前,测量极化率的方法主要包括法拉第光旋角法、顺磁共振法、抽运光透过率法等,但前两种方法需要较大的外部磁场,第三种方法需要改变检测光路,使线偏振检测光成为圆偏振检测光。所以这几种方法都无法在不改变原子磁强计光路及极弱磁环境中直接进行极化率测量,也就无法实现极化率的控制。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了提供了一种基于方波磁场幅值的原子磁强计极化率操控方法及装置,可以在不改变原子磁强计光路的基础上将极化率稳定在设定值,操作简单。
本发明提供了一种基于方波磁场幅值的原子磁强计极化率操控方法,包括两套三轴亥姆霍兹线圈,一套用于施加三方向的高频方波调制磁场,通过改变调制磁场的幅值,改变极化率的大小;一套用于施加低频方波标定磁场,通过测量在方波标定磁场施加和关断过程中原子的瞬态响应,通过拟合布洛赫方程的瞬态解,计算出减慢因子,进而获得极化率的测量值;在上位机采集软件输入极化率的设定值,通过和测量值做对比,实时控制高频方波调制磁场的幅值改变极化率测量值,实现闭环控制,直到测量值达到设定值。
作为优选,所述高频方波调制磁场的频率大于屏蔽桶内碱金属原子在剩磁环境下的拉莫尔进动频率;低频方波标定磁场的频率避开原子自旋对磁场响应的敏感频率区域,同时保证磁场施加过程原子瞬态响应的有效性。
作为优选,本发明的操控方法包括以下步骤:
S1:开启无磁电加热系统对碱金属原子气室进行加热,使碱金属气体达到SERF态;
S2:调整抽运光光路和检测光光路;开启抽运激光器和检测激光器,使原子磁强计处于正常工作状态;
S3:在上位机数据采集软件上输入极化率的设定值;
S4:用低频标定磁场线圈施加方波磁场,用采集软件原子磁强计的瞬态响应信号,并用减慢效应法计算出极化率的测量值;
S5:上位机数据采集软件对比设定值和测量值,根据高频方波调制磁场幅值对应极化率的关系,改变磁场幅值;
S6:重复步骤S4和S5,直到极化率的测量值达到设定值。
作为优选,所述步骤S4中施加的方波磁场频率小于等于0.1Hz。
作为优选,两套三轴亥姆霍兹线圈的三轴磁场相互解耦采用柔性薄膜的方式粘贴在隔热保温材料腔的腔体上,或搭建两套三轴正交的线圈骨架。
本发明还提供了一种基于方波磁场幅值的原子磁强计极化率操控装置,包括碱金属原子气室、抽运光路、检测光路以及数据处理上位机,所述的碱金属原子气室外部由内向外依次设有无磁电加热装置、隔热保温材料腔、高频调制磁场线圈、低频标定磁场线圈、磁屏蔽系统,所述高频调制磁场线圈与数据处理上位机电连接,通过数据处理上位机控制高频调制磁场线圈的电流大小,所述的抽运光路从外部射入碱金属原子气室,所述的检测光路穿过碱金属原子气室并与数据处理上位机连接,所述的抽运光路和所述检测光路垂直。
作为优选,所述的抽运光路包括按光路方向依次设置的抽运光激光器、扩束镜、偏振片、抽运四分之一波片,所述抽运光激光器发出抽运光,抽运光经过扩束镜、偏振片、抽运四分之一波片后射入碱金属原子气室。
作为优选,所述的检测光路包括按光路方向依次设置的检测光激光器、起偏器、检测四分之一波片、PEM调制器、检偏器、光电探测器,所述检测光激光器发出检测光,检测光经过起偏器调整后射入碱金属原子气室并从另一端射出,射出的检测光再依次经过检测四分之一波片、PEM调制器、检偏器到达光电探测器,光电探测器与数据处理上位机电连接。
本发明通过施加低频方波标定磁场使原子磁强计的瞬态响应可以被采集到,然后用减慢效应法计算出极化率的值,这个值是极化率的测量值;通过上位采集软件输入极化率的设定值,将设定值和测量值做对比,通过改变高频方波调制磁场的幅值,达到闭环实时控制检测方向极化率值的目的,对稳定原子磁强计灵敏度有非常重要的作用。
附图说明
图1为本发明装置系统示意图;
图2为本发明的实施流程图;
图3为低频方波磁场时序和原子磁强计的瞬态响应图;
图4为高频磁场幅值和极化率的关系曲线;
图中:1-检测光激光器,2-起偏器,3-检测四分之一波片 4-PEM调制器,5-检偏器,6-光电探测器,7-磁屏蔽系统,8-低频标定磁场线圈,9-隔热保温材料腔,10-无磁电加热装置,11-碱金属原子气室,12-高频调制磁场线圈,13-抽运光激光器,14-扩束镜,15-偏振片,16-抽运四分之一波片,17-数据处理上位机。
具体实施方式
下面通过附图以及具体实施方式进一步说明本发明。
一种基于方波磁场幅值的原子磁强计极化率操控方法,包括两套三轴亥姆霍兹线圈,一套用于施加三方向的高频方波调制磁场,通过改变调制磁场的幅值,改变极化率的大小;一套用于施加低频方波标定磁场,通过测量在方波标定磁场施加和关断过程中原子的瞬态响应,通过拟合布洛赫方程的瞬态解,计算出减慢因子,进而获得极化率的测量值;在上位机采集软件输入极化率的设定值,通过和测量值做对比,实时控制高频方波调制磁场的幅值改变极化率测量值,实现闭环控制,直到测量值达到设定值。所述高频方波调制磁场的频率大于磁屏蔽系统内碱金属原子在剩磁环境下的拉莫尔进动频率;低频方波标定磁场的频率避开原子自旋对磁场响应的敏感频率区域。
如图2所示,本发明的具体实施步骤包括:
S1:开启无磁电加热系统对碱金属原子气室进行加热,使碱金属气体达到SERF状态;
S2:调整抽运光光路和检测光光路;开启抽运激光器和检测激光器,使原子磁强计处于正常工作状态;
S3:在上位机数据采集软件上输入极化率的设定值;
S4:用低频标定磁场线圈施加方波磁场,用采集软件原子磁强计的瞬态响应信号,并用减慢效应法计算出极化率的测量值;
S5:上位机数据采集软件对比设定值和测量值,根据高频方波调制磁场幅值对应极化率的关系,改变磁场幅值;
S6:重复步骤S4和S5,直到极化率的测量值达到设定值。
其中,所述步骤S4中施加的方波磁场频率小于等于0.1Hz。
两套三轴亥姆霍兹线圈的三轴磁场相互解耦采用柔性薄膜的方式粘贴在隔热保温材料腔的腔体上,或搭建两套三轴正交的线圈骨架。
实施例1
如图1所示,一种基于方波磁场幅值的原子磁强计极化率操控装置,包括碱金属原子气室11、抽运光路、检测光路以及数据处理上位机17,所述的碱金属原子气室11外部由内向外依次设有无磁电加热装置10、隔热保温材料腔9、高频调制磁场线圈12、低频标定磁场线圈8、磁屏蔽系统7,所述高频调制磁场线圈12与数据处理上位机17电连接,通过数据处理上位机17控制高频调制磁场线圈12的电流大小,所述的抽运光路从外部射入碱金属原子气室11,所述的检测光路穿过碱金属原子气室11并与数据处理上位机17连接,所述的抽运光路和所述检测光路垂直。
其中,所述的抽运光路包括按光路方向依次设置的抽运光激光器13、扩束镜14、偏振片15、抽运四分之一波片16,所述抽运光激光器13发出抽运光,抽运光经过扩束镜14、偏振片15、抽运四分之一波片16后射入碱金属原子气室11。
所述的检测光路包括按光路方向依次设置的检测光激光器1、起偏器2、检测四分之一波片3、PEM调制器4、检偏器5、光电探测器6,所述检测光激光器1发出检测光,检测光经过起偏器2调整后射入碱金属原子气室11并从另一端射出,射出的检测光再依次经过检测四分之一波片3、PEM调制器4、检偏器5到达光电探测器6,光电探测器6与数据处理上位机17电连接。
(1):开启无磁电加热系统,将碱金属原子气室加热,加热温度根据碱金属种类决定;
(2):开启抽运激光器和检测激光器,调整光路,使原子磁强计处于工作状态;抽运光激光器13输出激光的波长对应碱金属的fD1线,经过扩束镜14将激光的光斑扩大,然后经过偏振片15和抽运四分之一波片16后,保证以圆偏振形式进入碱金属原子气室11,抽运碱金属原子;检测光激光器1输出激光的波长对应碱金属的fD2线,经过起偏器2后,以线偏振的偏振态进入碱金属原子气室11,通过碱金属原子气室11后,经过检测四分之一波片3,进入PEM调制器4进行调制,经过检偏器5后,被光电探测器6接收;
(3):用低频标定磁场线圈8施加0.1Hz的方波磁场。如图3的(a)所示,通过低频标定磁场线圈8施加频率为0.1Hz的方波磁场信号,磁强计的瞬态响应如图3的(b)所示,每个磁场上升沿对应着瞬态响应信号的开启阶段,放大图如图3的(b)中点线所示。根据解算自旋极化率的瞬态响应,曲线的形式如下: 是稳定值,P0是初始值,Rtot是总驰豫率, />是减慢因子;
(4)在数据处理上位机17中输入极化率的设定值。如果设定值和步骤(3)测量出的极化率数值不同,用高频调制磁场线圈12调解极化率的值;
(5)在Y轴施加调制磁场下的布洛赫方程为:
(6)数据处理上位机17实时采集步骤(3)的极化率测量值并根据(4)的设定值按照步骤(5)进行调节,直到测量值和设定值相同。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围内。
Claims (6)
1.一种基于方波磁场幅值的原子磁强计极化率操控方法,其特征在于:包括两套三轴亥姆霍兹线圈,一套用于施加三方向的高频方波调制磁场,通过改变调制磁场的幅值,改变极化率的大小;一套用于施加低频方波标定磁场,通过测量在方波标定磁场施加和关断过程中原子的瞬态响应,通过拟合布洛赫方程的瞬态解,计算出减慢因子,进而获得极化率的测量值;在上位机采集软件输入极化率的设定值,通过和测量值做对比,实时控制高频方波调制磁场的幅值改变极化率测量值,实现闭环控制,直到测量值达到设定值,具体包括以下步骤:
S1:开启无磁电加热系统对碱金属原子气室进行加热,使碱金属气体达到SERF态;
S2:调整抽运光光路和检测光光路;开启抽运激光器和检测激光器,使原子磁强计处于正常工作状态;
S3:在上位机数据采集软件上输入极化率的设定值;
S4:用低频标定磁场线圈施加方波磁场,用采集软件原子磁强计的瞬态响应信号,并用减慢效应法计算出极化率的测量值,所述方波磁场频率小于等于0.1Hz;
S5:上位机数据采集软件对比设定值和测量值,根据高频方波调制磁场幅值对应极化率的关系,改变磁场幅值;
S6:重复步骤S4和S5,直到极化率的测量值达到设定值。
2.根据权利要求1所述的一种基于方波磁场幅值的原子磁强计极化率操控方法,其特征在于:所述高频方波调制磁场的频率大于屏蔽桶内碱金属原子在剩磁环境下的拉莫尔进动频率;低频方波调制磁场的频率避开原子自旋对磁场响应的敏感频率区域。
3.根据权利要求1所述的一种基于方波磁场幅值的原子磁强计极化率操控方法,其特征在于:两套三轴亥姆霍兹线圈的三轴磁场相互解耦采用柔性薄膜的方式粘贴在隔热保温材料腔的腔体上,或搭建两套三轴正交的线圈骨架。
4.根据权利要求1所述的一种基于方波磁场幅值的原子磁强计极化率操控方法的操控装置,包括碱金属原子气室、抽运光路、检测光路以及数据处理上位机,其特征在于:所述的碱金属原子气室外部由内向外依次设有无磁电加热装置、隔热保温材料腔、高频调制磁场线圈、低频标定磁场线圈、磁屏蔽系统,所述高频调制磁场线圈与数据处理上位机电连接,所述的抽运光路从外部射入碱金属原子气室,所述的检测光路穿过碱金属原子气室并与数据处理上位机连接,所述的抽运光路和所述检测光路垂直。
5.根据权利要求4所述的一种基于方波磁场幅值的原子磁强计极化率操控装置,其特征在于:所述的抽运光路包括按光路方向依次设置的抽运光激光器、扩束镜、偏振片、抽运四分之一波片,所述抽运光激光器发出抽运光,抽运光经过扩束镜、偏振片、抽运四分之一波片后射入碱金属原子气室。
6.根据权利要求4所述的一种基于方波磁场幅值的原子磁强计极化率操控装置,其特征在于:所述的检测光路包括按光路方向依次设置的检测光激光器、起偏器、检测四分之一波片、PEM调制器、检偏器、光电探测器,所述检测光激光器发出检测光,检测光经过起偏器调整后射入碱金属原子气室并从另一端射出,射出的检测光再依次经过检测四分之一波片、PEM调制器、检偏器到达光电探测器,光电探测器与数据处理上位机电连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211560465.3A CN115754844B (zh) | 2022-12-07 | 2022-12-07 | 基于方波磁场幅值的原子磁强计极化率操控方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211560465.3A CN115754844B (zh) | 2022-12-07 | 2022-12-07 | 基于方波磁场幅值的原子磁强计极化率操控方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115754844A CN115754844A (zh) | 2023-03-07 |
CN115754844B true CN115754844B (zh) | 2023-05-30 |
Family
ID=85344686
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211560465.3A Active CN115754844B (zh) | 2022-12-07 | 2022-12-07 | 基于方波磁场幅值的原子磁强计极化率操控方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115754844B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102901939A (zh) * | 2012-10-16 | 2013-01-30 | 北京航空航天大学 | 一种用于原子自旋器件稳定的原子自旋serf态的精密操控方法 |
RU2016136138A3 (zh) * | 2016-09-07 | 2018-03-15 | ||
CN108445428A (zh) * | 2018-04-11 | 2018-08-24 | 北京航空航天大学 | 一种serf原子磁强计电子极化率测量方法 |
US10274549B1 (en) * | 2016-07-13 | 2019-04-30 | AOSense, Inc. | Scalar atomic magnetometer with heading error suppression |
CN110244242A (zh) * | 2019-07-23 | 2019-09-17 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于相位延迟的碱金属原子自旋极化调控装置及方法 |
CN110426651A (zh) * | 2019-06-17 | 2019-11-08 | 北京航空航天大学 | 基于serf磁强计的三维原位磁线圈标定实验装置及方法 |
CN110879374A (zh) * | 2019-11-26 | 2020-03-13 | 北京航空航天大学 | 一种单光束自旋极化和检测方法 |
CN111060855A (zh) * | 2018-10-17 | 2020-04-24 | 北京自动化控制设备研究所 | 一种基于电子自旋磁共振线宽压窄的高灵敏磁场测量方法 |
CN112729269A (zh) * | 2020-12-30 | 2021-04-30 | 中国人民解放军国防科技大学 | 一种抑制碱金属与稀有气体原子耦合效应的工作方法 |
CN112924910A (zh) * | 2021-01-29 | 2021-06-08 | 北京航空航天大学 | 一种基于原位磁强计的屏蔽桶内剩磁测量方法 |
CN113639883A (zh) * | 2021-08-18 | 2021-11-12 | 之江实验室 | 一种碱金属原子磁强计自旋极化率空间分布原位测量系统及方法 |
CN114527414A (zh) * | 2022-02-16 | 2022-05-24 | 西安交通大学 | 自动补偿外干扰磁场的生物磁测量原子磁强计系统及方法 |
CN115389988A (zh) * | 2022-10-31 | 2022-11-25 | 之江实验室 | 一种原子磁强计交流磁场幅值测量装置及方法 |
-
2022
- 2022-12-07 CN CN202211560465.3A patent/CN115754844B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102901939A (zh) * | 2012-10-16 | 2013-01-30 | 北京航空航天大学 | 一种用于原子自旋器件稳定的原子自旋serf态的精密操控方法 |
US10274549B1 (en) * | 2016-07-13 | 2019-04-30 | AOSense, Inc. | Scalar atomic magnetometer with heading error suppression |
RU2016136138A3 (zh) * | 2016-09-07 | 2018-03-15 | ||
CN108445428A (zh) * | 2018-04-11 | 2018-08-24 | 北京航空航天大学 | 一种serf原子磁强计电子极化率测量方法 |
CN111060855A (zh) * | 2018-10-17 | 2020-04-24 | 北京自动化控制设备研究所 | 一种基于电子自旋磁共振线宽压窄的高灵敏磁场测量方法 |
CN110426651A (zh) * | 2019-06-17 | 2019-11-08 | 北京航空航天大学 | 基于serf磁强计的三维原位磁线圈标定实验装置及方法 |
CN110244242A (zh) * | 2019-07-23 | 2019-09-17 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于相位延迟的碱金属原子自旋极化调控装置及方法 |
CN110879374A (zh) * | 2019-11-26 | 2020-03-13 | 北京航空航天大学 | 一种单光束自旋极化和检测方法 |
CN112729269A (zh) * | 2020-12-30 | 2021-04-30 | 中国人民解放军国防科技大学 | 一种抑制碱金属与稀有气体原子耦合效应的工作方法 |
CN112924910A (zh) * | 2021-01-29 | 2021-06-08 | 北京航空航天大学 | 一种基于原位磁强计的屏蔽桶内剩磁测量方法 |
CN113639883A (zh) * | 2021-08-18 | 2021-11-12 | 之江实验室 | 一种碱金属原子磁强计自旋极化率空间分布原位测量系统及方法 |
CN114527414A (zh) * | 2022-02-16 | 2022-05-24 | 西安交通大学 | 自动补偿外干扰磁场的生物磁测量原子磁强计系统及方法 |
CN115389988A (zh) * | 2022-10-31 | 2022-11-25 | 之江实验室 | 一种原子磁强计交流磁场幅值测量装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115754844A (zh) | 2023-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108519566B (zh) | 一种基于光频移调制的serf原子磁强计装置及方法 | |
CN109188316A (zh) | 液晶相位补偿的自激式原子磁传感器及磁场测量方法 | |
CN106385283A (zh) | 一种应用于原子自旋进动检测的抽运光调制解调系统及方法 | |
CN103969604A (zh) | 射频原子磁力仪及其测量核磁共振信号的方法 | |
CN105301541A (zh) | 原子磁强计的磁线圈x、y轴非正交角的测量装置与方法 | |
CN110261796A (zh) | 一种基于零场共振的serf原子磁强计及磁补偿方法 | |
CN109358302B (zh) | 一种无被动磁屏蔽原子磁力计装置及测磁方法 | |
CN102859384A (zh) | 利用光诱发的超极化的核磁共振磁力计 | |
Lu et al. | Effects of AC magnetic field on spin-exchange relaxation of atomic magnetometer | |
CN112505595A (zh) | 一种高带宽高灵敏闭环serf原子磁力计装置 | |
CN110514193B (zh) | 核磁共振陀螺及探测方法 | |
US7315168B2 (en) | Shimming with MRI gradient | |
CN110045309A (zh) | 一种新型的光泵磁力仪精度自检测方法 | |
CN115754844B (zh) | 基于方波磁场幅值的原子磁强计极化率操控方法及装置 | |
US4847561A (en) | Domain and domain wall transition spectroscopy | |
Xing et al. | Optical rotation detection method based on acousto-optic modulation in an atomic spin co-magnetometer | |
RU2704391C1 (ru) | Способ управления атомарным магнитометрическим датчиком при работе в составе многоканальной диагностической системы | |
Long et al. | Speedy in-situ magnetic field compensation algorithm for multiple-channel single-beam SERF atomic magnetometers | |
JP5212972B2 (ja) | 計測装置および計測方法 | |
Li et al. | In situ simultaneous measurement of magnetic coil constants and nonorthogonal angles using atomic magnetometers | |
CN110411434A (zh) | 一种原子陀螺仪三维原位快速磁补偿方法 | |
Zhang et al. | Spin dynamics of magnetic resonance with parametric modulation in a potassium vapor cell | |
Cao et al. | In situ calibration of coils constant in a spin-exchange relaxation-free (SERF) co-magnetometer | |
CN113189528B (zh) | 应用于小型化原子磁力计的读出方法、装置、设备及介质 | |
JP4392941B2 (ja) | 磁気共鳴イメージング装置 |
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 |